diff --git a/aeolis/constants.py b/aeolis/constants.py index 3502edb2..e7a4a601 100644 --- a/aeolis/constants.py +++ b/aeolis/constants.py @@ -109,6 +109,7 @@ 'germinate', # vegetation germination 'lateral', # vegetation lateral expansion 'vegfac', # Vegetation factor to modify shear stress by according to Raupach 1993 + 'fence_height', # Fence height 'R', # [m] wave runup 'eta', # [m] wave setup 'sigma_s', # [m] swash diff --git a/aeolis/wind.py b/aeolis/wind.py index cfdbdc4f..c4f51d9c 100644 --- a/aeolis/wind.py +++ b/aeolis/wind.py @@ -121,10 +121,7 @@ def interpolate(s, p, t): s['uws'] = - s['uw'] * np.sin((-p['alfa'] + s['udir']) / 180. * np.pi) # alfa [deg] is real world grid cell orientation (clockwise) s['uwn'] = - s['uw'] * np.cos((-p['alfa'] + s['udir']) / 180. * np.pi) - - if p['ny'] == 0: - s['uwn'][:,:] = 0. - + s['uw'] = np.abs(s['uw']) # Compute wind shear velocity