Changeset 14777

Timestamp:

2021-05-03T10:57:28+02:00 (3 years ago)

Author:

clem

Message:

4.0-HEAD: avoid roundoff error (it does not change results, it is only important for online diags)

File:

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

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

@@ -121 +121 @@
          DO ji = 1, npti
             rswitch      = MAX( 0._wp , SIGN( 1._wp , zhnew(ji) - epsi20 ) ) 
-            qnew(ji,jk1) = rswitch * ( zeh_cum1(ji,jk1) - zeh_cum1(ji,jk1-1) ) / MAX( zhnew(ji), epsi20 )
+            qnew(ji,jk1) = rswitch * MAX( 0._wp, zeh_cum1(ji,jk1) - zeh_cum1(ji,jk1-1) ) / MAX( zhnew(ji), epsi20 ) ! max for roundoff error
          END DO
       END DO