Changeset 14252

Timestamp:

2020-12-24T09:07:31+01:00 (3 years ago)

Author:

vancop

Message:

bug fixes in topo ponds for trunk

File:

• NEMO/trunk/src/ICE/icethd_pnd.F90 (modified) (5 diffs)

NEMO/trunk/src/ICE/icethd_pnd.F90

@@ -511 +511 @@
       !!-------------------------------------------------------------------
       REAL(wp), PARAMETER :: &   ! shared parameters for topographic melt ponds
-         zTd     = 0.15_wp       , & ! temperature difference for freeze-up (C)
+         zTd     = 273._wp       , & ! temperature difference for freeze-up (K)
          zvp_min = 1.e-4_wp          ! minimum pond volume (m)
 
@@ -530 +530 @@
          zv_mlt          ! total amount of meltwater produced
 
-      REAL(wp), DIMENSION(jpi,jpj) ::   zvolp, &     !! total melt pond water available before redistribution and drainage
-                                        zvolp_res    !! remaining melt pond water available after drainage
+      REAL(wp), DIMENSION(jpi,jpj) ::   zvolp_ini , &   !! total melt pond water available before redistribution and drainage
+                                        zvolp     , &   !! total melt pond water volume
+                                        zvolp_res       !! remaining melt pond water available after drainage
 
       REAL(wp), DIMENSION(jpi,jpj,jpl) ::   z1_a_i
@@ -624 +625 @@
       END DO ! ji
 
+      zvolp_ini(:,:) = zvolp(:,:)
+
       !--------------------------------------------------------------
       ! Redistribute and drain water from ponds
@@ -637 +640 @@
          DO_2D( 1, 1, 1, 1 )
 
-            IF ( at_i(ji,jj) > 0.01 .AND. hm_i(ji,jj) > rn_himin .AND. vt_ip(ji,jj) > zvp_min * at_i(ji,jj) ) THEN
+            IF ( at_i(ji,jj) > 0.01 .AND. hm_i(ji,jj) > rn_himin .AND. zvolp_ini(ji,jj) > zvp_min * at_i(ji,jj) ) THEN
 
                !--------------------------
@@ -688 +691 @@
 
                         ! differential growth of base of surface floating ice layer
-                        zdTice = MAX( - t_su(ji,jj,jl) - zTd , 0._wp ) ! > 0
+                        zdTice = MAX( - ( t_su(ji,jj,jl) - zTd ) , 0._wp ) ! > 0
                         zomega = rcnd_i * zdTice / zrhoi_L
                         zdHui  = SQRT( 2._wp * zomega * rDt_ice + ( v_il(ji,jj,jl) / a_i(ji,jj,jl) )**2 ) &