Display layered results; fix read_swash for layered vectors
This commit is contained in:
parent
cb05a03a28
commit
12ad3b5c59
3 changed files with 86 additions and 10 deletions
|
@ -23,22 +23,29 @@ config.read("config.ini")
|
|||
inp = pathlib.Path(config.get("post", "inp"))
|
||||
root = pathlib.Path(config.get("swash", "out"))
|
||||
|
||||
data = np.load(inp.joinpath(config.get("post", "case")))
|
||||
bathy = pd.read_hdf(
|
||||
pathlib.Path(config.get("data", "out")).joinpath("bathy.h5"), "bathy"
|
||||
)
|
||||
|
||||
x = data["x"]
|
||||
t = data["t"]
|
||||
|
||||
wl = np.maximum(data["watl"], -data["botl"])
|
||||
print(x.size, -np.arange(0, 1 * bathy.hstru.size, 1)[::-1].size)
|
||||
def data(var):
|
||||
return np.load(inp.joinpath(f"{var}.npy"))
|
||||
|
||||
|
||||
x = data("xp")
|
||||
t = data("tsec")
|
||||
|
||||
watl = data("watl")
|
||||
botl = data("botl")
|
||||
|
||||
wl = np.maximum(watl, -botl)
|
||||
# print(x.size, -np.arange(0, 1 * bathy.hstru.size, 1)[::-1].size)
|
||||
|
||||
fig, ax = plt.subplots()
|
||||
ax.plot(x, -data["botl"], c="k")
|
||||
ax.fill_between(
|
||||
x, -data["botl"], -data["botl"] + bathy.hstru, color="k", alpha=0.2
|
||||
)
|
||||
ax.plot(x, -botl, c="k")
|
||||
# ax.fill_between(
|
||||
# x, -botl, -data["botl"] + bathy.hstru, color="k", alpha=0.2
|
||||
# )
|
||||
(line,) = ax.plot(x, wl[0])
|
||||
|
||||
|
||||
|
|
69
swash/processing/layers.py
Normal file
69
swash/processing/layers.py
Normal file
|
@ -0,0 +1,69 @@
|
|||
import argparse
|
||||
import configparser
|
||||
import logging
|
||||
import pathlib
|
||||
|
||||
import matplotlib.pyplot as plt
|
||||
import matplotlib.animation as animation
|
||||
import numpy as np
|
||||
import pandas as pd
|
||||
|
||||
|
||||
parser = argparse.ArgumentParser(description="Animate swash output")
|
||||
parser.add_argument("-v", "--verbose", action="count", default=0)
|
||||
args = parser.parse_args()
|
||||
|
||||
logging.basicConfig(level=max((10, 20 - 10 * args.verbose)))
|
||||
log = logging.getLogger("post")
|
||||
|
||||
log.info("Starting post-processing")
|
||||
config = configparser.ConfigParser()
|
||||
config.read("config.ini")
|
||||
|
||||
inp = pathlib.Path(config.get("post", "inp"))
|
||||
root = pathlib.Path(config.get("swash", "out"))
|
||||
|
||||
bathy = pd.read_hdf(
|
||||
pathlib.Path(config.get("data", "out")).joinpath("bathy.h5"), "bathy"
|
||||
)
|
||||
|
||||
|
||||
def data(var):
|
||||
return np.load(inp.joinpath(f"{var}.npy"))
|
||||
|
||||
|
||||
x = data("xp")
|
||||
t = data("tsec")
|
||||
|
||||
watl = data("watl")
|
||||
botl = data("botl")
|
||||
zk = data("zk")
|
||||
velk = data("velk")
|
||||
|
||||
wl = np.maximum(watl, -botl)
|
||||
# print(x.size, -np.arange(0, 1 * bathy.hstru.size, 1)[::-1].size)
|
||||
|
||||
fig, ax = plt.subplots()
|
||||
ax.plot(x, -botl, c="k")
|
||||
# ax.fill_between(
|
||||
# x, -botl, -data["botl"] + bathy.hstru, color="k", alpha=0.2
|
||||
# )
|
||||
lines = ax.plot(x, zk[0].T)
|
||||
|
||||
print(velk.shape)
|
||||
velk = velk.reshape((6001, 10, 2, 1251))
|
||||
vk = np.sqrt((velk[1000] ** 2).sum(axis=1))
|
||||
print(vk.shape)
|
||||
plt.imshow(vk)
|
||||
plt.colorbar()
|
||||
|
||||
# def animate(i):
|
||||
# for line, z in zip(lines, zk[i]):
|
||||
# line.set_ydata(z)
|
||||
# return lines
|
||||
#
|
||||
# ani = animation.FuncAnimation(
|
||||
# fig, animate, frames=wl[:, 0].size, interval=20, blit=True
|
||||
# )
|
||||
|
||||
plt.show(block=True)
|
|
@ -41,7 +41,7 @@ class ReadSwash:
|
|||
return self.read_nohead(path).reshape((self._n_t, -1, self._n_x))
|
||||
|
||||
def read_vector_lay(self, path):
|
||||
return self.read_nohead(path).reshape((self._n_t, 2, -1, self._n_x))
|
||||
return self.read_nohead(path).reshape((self._n_t, -1, 2, self._n_x))
|
||||
|
||||
@property
|
||||
def t(self):
|
||||
|
|
Loading…
Reference in a new issue