Changeset 15045

Timestamp:

2021-06-23T12:43:59+02:00 (3 years ago)

Author:

clem

Message:

4.0-HEAD: bug found by ecmwf: do not do ice remapping if, by some mystery and in very rare occasion, hR=hL=hice.

File:

• NEMO/releases/r4.0/r4.0-HEAD/src/ICE/iceitd.F90 (modified) (1 diff)

NEMO/releases/r4.0/r4.0-HEAD/src/ICE/iceitd.F90

@@ -368 +368 @@
             !
             ! Compute coefficients of g(eta) = g0 + g1*eta
-            zdhr = 1._wp / (phR(ji) - phL(ji))
+            IF( phR(ji) > phL(ji) ) THEN   ;   zdhr = 1._wp / (phR(ji) - phL(ji))
+            ELSE                           ;   zdhr = 0._wp ! if hR=hL=hice => no remapping
+            ENDIF
+            !!zdhr = 1._wp / (phR(ji) - phL(ji))
             zwk1 = 6._wp * paice(ji) * zdhr
             zwk2 = ( phice(ji) - phL(ji) ) * zdhr