From ea46fb0272fadb103711e8bc5bf1f28c4a5b10c2 Mon Sep 17 00:00:00 2001
From: Quentin Rodier <quentin.rodier@meteo.fr>
Date: Wed, 13 Apr 2022 11:08:55 +0200
Subject: [PATCH] Quentin 13/04/2022: bugfix infinite loop sedimentation init
 inner values of ZREMAINT (bug in mesonh only)

---
 src/common/micro/mode_ice4_sedimentation_split.F90 | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/src/common/micro/mode_ice4_sedimentation_split.F90 b/src/common/micro/mode_ice4_sedimentation_split.F90
index 439022136..4864c1472 100644
--- a/src/common/micro/mode_ice4_sedimentation_split.F90
+++ b/src/common/micro/mode_ice4_sedimentation_split.F90
@@ -336,7 +336,8 @@ END IF
 PINPRX(:,:) = 0.
 ZINVTSTEP=1./PTSTEP
 ZRSMIN(:) = ICED%XRTMIN(:) * ZINVTSTEP
-ZREMAINT(:,:) = PTSTEP
+ZREMAINT(:,:) = 0.
+ZREMAINT(D%NIB:D%NIE,D%NJB:D%NJE) = PTSTEP
 !
 DO WHILE (ANY(ZREMAINT>0.))
   ISEDIM = 0
-- 
GitLab