From 15de9190642b4b966b6bd3baa3d2a85ecfb55919 Mon Sep 17 00:00:00 2001
From: Delaigue Olivier <olivier.delaigue@irstea.fr>
Date: Tue, 5 Sep 2023 09:51:43 +0200
Subject: [PATCH] fix(fortran): take into account case 'dG = 0' in
'frun_CemaNeige' Refs: #178
---
src/frun_CEMANEIGE.f90 | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/src/frun_CEMANEIGE.f90 b/src/frun_CEMANEIGE.f90
index 047f2c32..959dff22 100644
--- a/src/frun_CEMANEIGE.f90
+++ b/src/frun_CEMANEIGE.f90
@@ -184,7 +184,7 @@
Gratio = MIN(Gratio+(Psol-Melt)/Gacc,1.d0) !Psol - Melt = dG
IF (Gratio.EQ.1.) Glocalmax = Gthreshold
ENDIF
- IF (dG.LT.0.) THEN
+ IF (dG.LE.0.) THEN
Gratio=MIN(G/Glocalmax,1.d0)
ENDIF
ELSE
--
GitLab