Merge pull request #1608 from knutfrode/dev

[run-ex] OpenBerg: setting property moving to 0 for grounded icebergs…

Author: knutfrode

examples/example_openberg.py

@@ -32,6 +32,7 @@
 # Simulating drift for 48 hours
 o = OpenBerg()
 o.set_config('drift:vertical_profile', False)
+o.set_config('drift:horizontal_diffusivity', 100)
 o.add_readers_from_list(forcing)
 o.seed_elements(**icebergs)
 o.run(duration=timedelta(days=2))

opendrift/models/openberg.py

@@ -474,7 +474,8 @@ def dynamic(t,iceb_vel, water_vel, wind_vel, wave_height, wave_direction, Ao,
         hwall = draft - water_depth
         grounded = np.logical_and(hwall >= 0, grounding)
         if any(grounded) and grounding:
-            logger.info(f"Grounding condition : Icebergs grounded = {len(hwall[hwall>0])}, hwall={np.round(hwall[hwall>0],3)} meters")
+            logger.debug(f"Grounding condition : Icebergs grounded = {len(hwall[hwall>0])}, hwall={np.round(hwall[hwall>0],3)} meters")
+            self.elements.moving[grounded] = 0  # Grounded icebergs shall not move, also not with diffusivity
 
         sol = solve_ivp(dynamic, [0, self.time_step.total_seconds()], V0,
                         args=(water_vel, wind_vel, wave_height, wave_direction, Ao, Aa, rho_water,
@@ -537,7 +538,7 @@ def roll_over(self):
         W, L = np.min([L, W], axis=0), np.max([L, W], axis=0)
         mask = (W / H) < crit
         if any(mask):
-            logger.info(f"Rolling over : {np.sum(mask)} icebergs ...")
+            logger.debug(f"Rolling over : {np.sum(mask)} icebergs ...")
             nL, nW, nH = (np.max([L[mask], H[mask]], axis=0), np.min([L[mask], H[mask]], axis=0), W[mask])
             L[mask], W[mask], H[mask] = nL, nW, nH
         depthib = H * alpha