Add Jan's calving to sandbox

Co-authored-by: Jan Malles <[email protected]>

oggm/core/flowline.py

@@ -99,6 +99,14 @@ def __init__(self, line=None, dx=1, map_dx=None,
     def has_ice(self):
         return np.any(self.thick > 0)
 
+    @utils.lazy_property
+    def water_depth(self):
+        return utils.clip_min(self.water_level - self.bed_h, 0)
+
+    @utils.lazy_property
+    def bed_below_wl(self):
+        return self.bed_h < self.water_level
+
     @Centerline.widths.getter
     def widths(self):
         """Compute the widths out of H and shape"""

oggm/params.cfg

@@ -362,6 +362,10 @@ dynamic_spinup_min_ice_thick = 10.
 # 4: Marine-terminating, Lake-terminating, Shelf-terminating
 tidewater_type = 2
 
+# ocean water density in kg m-3; should be >= ice density
+# for lake-terminating glaciers this could be changed to 1000 kg m-3
+ocean_density = 1028
+
 # Should we switch on the k-calving parameterisation for tidewater glaciers?
 use_kcalving_for_inversion = False
 # Its possible to use kcalving_for_run but not during the inversion.
@@ -391,6 +395,10 @@ free_board_marine_terminating = 10, 50
 # For lake terminating glaciers, we have no way to know the water level,
 # so we set an arbitrary free board value
 free_board_lake_terminating = 10
+## Stretch distance in hydrostatic pressure balance calculations for terminal
+# water-terminating cliffs (in meters) - See Malles et al.
+# (the current value of 8000m might be too high)
+max_calving_stretch_distance = 8000
 # We extend the calving glaciers by an arbitrary number of grid points,
 # and following an arbitrary slope.
 # How many grid points should we extend the calving front with?