Merge pull request #82 from openearth/74-purpose-of-th_drylayer

Delete dry_layer function

Author: Sierd

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):