Changeset 14685

Timestamp:

2021-04-08T17:35:36+02:00 (3 years ago)

Author:

clem

Message:

4.0-HEAD: solve negative sublimation problems (ticket #2649)

File:

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

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

@@ -224 +224 @@
       zevap_rema(1:npti) = 0._wp
       DO ji = 1, npti
-         IF( evap_ice_1d(ji) > 0._wp ) THEN
-            zdeltah   (ji) = MAX( - evap_ice_1d(ji) * r1_rhos * rdt_ice, - h_s_1d(ji) )   ! amount of snw that sublimates, < 0            
-            zevap_rema(ji) = MAX( 0._wp, evap_ice_1d(ji) * rdt_ice + zdeltah(ji) * rhos ) ! remaining evap in kg.m-2 (used for ice sublimation later on)
-         ENDIF
+         zdeltah   (ji) = MAX( - evap_ice_1d(ji) * r1_rhos * rdt_ice, - h_s_1d(ji) )   ! amount of snw that sublimates, < 0           
+         zevap_rema(ji) = evap_ice_1d(ji) * rdt_ice + zdeltah(ji) * rhos               ! remaining evap in kg.m-2 (used for ice sublimation later on)
       END DO