make sure that melt_f limits are not exceeded in dynamic calibration

OGGM · Feb 2, 2024 · 0b82ca3 · 0b82ca3
1 parent 5f8a8c5
commit 0b82ca3
Showing 1 changed file with 4 additions and 2 deletions.
diff --git a/oggm/core/dynamic_spinup.py b/oggm/core/dynamic_spinup.py
@@ -1978,8 +1978,10 @@ def c_fun(melt_f):
     def minimise_with_spline_fit(fct_to_minimise, melt_f_guess, mismatch):
         # defines limits of melt_f in accordance to maximal allowed change
         # between iterations
-        melt_f_limits = [melt_f_initial - melt_f_max_step_length,
-                         melt_f_initial + melt_f_max_step_length]
+        melt_f_limits = [max(melt_f_initial - melt_f_max_step_length,
+                             melt_f_min),
+                         min(melt_f_initial + melt_f_max_step_length,
+                             melt_f_max)]
 
         # this two variables indicate that the limits were already adapted to
         # avoid an error