m2cce/internship
m2cce/internship
1
Fork 0
Code Issues Activity

Updated projection without breakwater

Browse source
This commit is contained in:
Edgar P. Burkhart 2022-03-30 12:42:09 +02:00
parent dbced643ca
commit 0a122f76f3
Signed by: edpibu
GPG key ID: 9833D3C5A25BD227
2 changed files with 5 additions and 3 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

3
data/config-nb.ini
View file

@ -8,12 +8,13 @@ psize=Psize.dat
hires_step=0.5
[out]
margin=0.005
no_breakwater=True
root=out_nb
sub=bathy_sub.npy
out=bathy.npy
step=1
left=0
left=-300
right=150
[artha]

5
data/processing/projection.py
View file

@ -91,13 +91,14 @@ x_min_hires = x[z_crossing] + (bathy_hires[:, 0].min() - bathy_hires[hires_cross
x_max_hires = x[z_crossing] + (bathy_hires[:, 0].max() - bathy_hires[hires_crossing, 0])
log.debug(f"Replacing range: [{x_min_hires},{x_max_hires}]")
flt_x = (x > x_min_hires) & (x < x_max_hires)
hstru = np.zeros(z.shape)
poro = np.zeros(z.shape)
psize = np.zeros(z.shape)
if config.getboolean("out", "no_breakwater", fallback=False):
z[flt_x] = z[flt_x][-1]
flt_x = np.abs(x) < 250
z[flt_x] = np.linspace(z[flt_x][0], z[flt_x][-1], flt_x.sum())
else:
flt_x = (x > x_min_hires) & (x < x_max_hires)
z[flt_x] = interpolate.griddata(
(bathy_hires[:, 0],),
bathy_hires[:, 1],