diff --git a/aeolis/threshold.py b/aeolis/threshold.py index c4bc2e84..db8f688d 100644 --- a/aeolis/threshold.py +++ b/aeolis/threshold.py @@ -238,22 +238,6 @@ def compute_moisture(s, p): return s -def dry_layer(s, p): - - Tdry_top = 12. * 3600. - zdry_max = 0.05 - - s['dzdry'] = (zdry_max - s['zdry']) * (p['dt_opt'] * p['accfac']) / Tdry_top + s['dzb'] - s['zdry'] += s['dzdry'] - s['zdry'] = np.minimum(np.maximum(s['zdry'], 0), zdry_max) - - ix = (s['zdry'] == 0.) - - s['uth'][ix] = np.inf - - - return s - #REMOVE?? CH # def compute_humidity(s, p):