fix(fortran): take into account case 'dG = 0' in 'frun_CemaNeige'

Refs: #178

fix(fortran): take into account case 'dG = 0' in 'frun_CemaNeige'
Refs: #178
15de9190 · Delaigue Olivier · a8bcc2ae · 15de9190
Commit 15de9190 authored 1 year ago by Delaigue Olivier
Hide whitespace changes
Inline Side-by-side

Showing

with 1 addition and 1 deletion
+1 -1
--- 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