From f65fc83331f4b18857002bb503a4c28212a48c62 Mon Sep 17 00:00:00 2001
From: Christa van IJzendoorn
 <52737913+christavanijzendoorn@users.noreply.github.com>
Date: Fri, 26 Jul 2024 16:44:33 -0700
Subject: [PATCH 1/2] max_moist implementation

Convert hard coded moisture threshold at which no transport is forced by increasing uth to inf. Max_moist is now used to have a user-defined value.
---
 aeolis/constants.py | 1 +
 aeolis/threshold.py | 6 ++++--
 2 files changed, 5 insertions(+), 2 deletions(-)

diff --git a/aeolis/constants.py b/aeolis/constants.py
index ef57da34..7b2e8a23 100644
--- a/aeolis/constants.py
+++ b/aeolis/constants.py
@@ -291,6 +291,7 @@
     'Cl_gw'                         : 0.7,                # [m] Groundwater overheight due to runup
     'in_gw'                         : 0,                  # [m] Initial groundwater level
     'GW_stat'                       : 1,                  # [m] Landward static groundwater boundary (if static boundary is defined)
+    'max_moist'                     : 10.,           # NEWCH      # [%] Moisture content (volumetric in percent) above which the threshold shear velocity is set to infinity (no transport, default value Delgado-Fernandez, 2010)
     'theta_dyn'                     : 33.,                # [degrees] Initial Dynamic angle of repose, critical dynamic slope for avalanching
     'theta_stat'                    : 34.,                # [degrees] Initial Static angle of repose, critical static slope for avalanching
     'avg_time'                      : 86400.,             # [s] Indication of the time period over which the bed level change is averaged for vegetation growth
diff --git a/aeolis/threshold.py b/aeolis/threshold.py
index 126cf21e..f2fc0f87 100644
--- a/aeolis/threshold.py
+++ b/aeolis/threshold.py
@@ -234,9 +234,11 @@ def compute_moisture(s, p):
     else:
         logger.log_and_raise('Unknown moisture formulation [%s]' % p['method_moist'], exc=ValueError)
 
+    th_mg = (p['max_moist']/100 * p['rhow'] / (p['rhog'] * (1. - p['porosity'])))
     # should be .04 according to Pye and Tsoar
-    # should be .64 according to Delgado-Fernandez (10% vol.)
-    ix = mg > 0.064
+    # should be .064 according to Delgado-Fernandez (10% vol.)
+    # .097 corresponds to 15%
+    ix = mg > th_mg
     s['uth'][ix] = np.inf
     
     return s

From 8f60d47002fb09e4e42c2bffc391eaf4c1a3db48 Mon Sep 17 00:00:00 2001
From: Christa van IJzendoorn
 <52737913+christavanijzendoorn@users.noreply.github.com>
Date: Fri, 26 Jul 2024 16:54:02 -0700
Subject: [PATCH 2/2] Update constants.py

remove NEWCH
---
 aeolis/constants.py | 1 +
 1 file changed, 1 insertion(+)

diff --git a/aeolis/constants.py b/aeolis/constants.py
index 7b2e8a23..5c50e094 100644
--- a/aeolis/constants.py
+++ b/aeolis/constants.py
@@ -292,6 +292,7 @@
     'in_gw'                         : 0,                  # [m] Initial groundwater level
     'GW_stat'                       : 1,                  # [m] Landward static groundwater boundary (if static boundary is defined)
     'max_moist'                     : 10.,           # NEWCH      # [%] Moisture content (volumetric in percent) above which the threshold shear velocity is set to infinity (no transport, default value Delgado-Fernandez, 2010)
+    'max_moist'                     : 10.,                # [%] Moisture content (volumetric in percent) above which the threshold shear velocity is set to infinity (no transport, default value Delgado-Fernandez, 2010)
     'theta_dyn'                     : 33.,                # [degrees] Initial Dynamic angle of repose, critical dynamic slope for avalanching
     'theta_stat'                    : 34.,                # [degrees] Initial Static angle of repose, critical static slope for avalanching
     'avg_time'                      : 86400.,             # [s] Indication of the time period over which the bed level change is averaged for vegetation growth