1
Fork 0

Fixed shift of bathy

This commit is contained in:
Edgar P. Burkhart 2022-03-04 12:00:35 +01:00
parent 7686aeffe5
commit 78a3f7fb54
Signed by: edpibu
GPG key ID: 9833D3C5A25BD227

View file

@ -34,9 +34,9 @@ poro = np.loadtxt(root.joinpath(config.get("data", "poro")))
psize = np.loadtxt(root.joinpath(config.get("data", "psize")))
log.info("Generating grid")
x_hires = np.arange(-0.5 * bathy_hires.size, 0, 0.5)
x_lores = np.arange(-1 * bathy_lores.size, 0, 1)
x_hstru = np.arange(-0.5 * hstru.size, 0, 0.5)
x_hires = -np.arange(0, 0.5 * bathy_hires.size, 0.5)[::-1]
x_lores = -np.arange(0, 1 * bathy_lores.size, 1)[::-1]
x_hstru = -np.arange(0, 0.5 * hstru.size, 0.5)[::-1]
log.info("Generating output data")
bathy_hires_pd = pd.Series(bathy_hires.copy(), index=x_hires)
@ -84,9 +84,7 @@ if config.getboolean("proc", "plot", fallback=False):
ax.plot(x_hires, bathy_hires, label="High-res")
ax.plot(x_lores, bathy_lores, label="Low-res")
ax.plot(bathy.index, bathy.z, ls="-.", c="k", label="Combined")
ax.plot(bathy.index, bathy.hstru, label="Hstru")
ax.plot(bathy.index, bathy.poro, label="Poro")
ax.plot(bathy.index, bathy.psize, label="Psize")
ax.plot(bathy.index, bathy.z + bathy.hstru, label="Hstru")
ax.grid()
ax.legend()