Changeset 14251 for NEMO

Timestamp:

2020-12-24T08:42:55+01:00 (3 years ago)

Author:

vancop

Message:

Fix pond lids in topo ponds for christmas

File:

• NEMO/branches/2020/SI3-05_MeltPonds_topo/src/ICE/icethd_pnd.F90 (modified) (6 diffs)

NEMO/branches/2020/SI3-05_MeltPonds_topo/src/ICE/icethd_pnd.F90

@@ -44 +44 @@
    REAL(wp), PARAMETER :: &   ! shared parameters for topographic melt ponds
       zhi_min = 0.1_wp        , & ! minimum ice thickness with ponds (m) 
-      zTd     = 0.15_wp       , & ! temperature difference for freeze-up (C)
+      zTd     = 273.0_wp      , & ! temperature difference for freeze-up (C)
       zvp_min = 1.e-4_wp          ! minimum pond volume (m)
 
@@ -531 +531 @@
          zv_mlt          ! total amount of meltwater produced
 
-      REAL(wp), DIMENSION(jpi,jpj) ::   zvolp, &     !! total melt pond water available before redistribution and 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
                                         
@@ -624 +625 @@
          END DO ! jj
       END DO ! ji
+
+      zvolp_ini(:,:) = zvolp(:,:)
                   
       !--------------------------------------------------------------
@@ -639 +642 @@
             DO ji = 1, jpi
 
-            IF ( at_i(ji,jj) > 0.01 .AND. hm_i(ji,jj) > zhi_min .AND. vt_ip(ji,jj) > zvp_min * at_i(ji,jj) ) THEN
+            IF ( at_i(ji,jj) > 0.01 .AND. hm_i(ji,jj) > zhi_min .AND. zvolp_ini(ji,jj) > zvp_min * at_i(ji,jj) ) THEN
                   
                !--------------------------
@@ -690 +693 @@
 
                         ! 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  ! line valid for temp in C 
+                        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 ) &
@@ -726 +730 @@
                      !!! we have the category-dependent flux, let us use it ?
                      zfsurf = qns_ice(ji,jj,jl) + qsr_ice(ji,jj,jl)                     
-                     zdHui  = MAX ( -zfsurf * rdt_ice/zrhoi_L , 0._wp )
+                     zdHui  = MAX ( -zfsurf * rdt_ice / zrhoi_L , 0._wp )
                      zdvice = MIN ( zdHui * a_ip(ji,jj,jl) , v_ip(ji,jj,jl) )
+
                      IF ( zdvice > epsi10 ) THEN
+
                         v_il (ji,jj,jl)  = v_il(ji,jj,jl)   + zdvice
                         v_ip(ji,jj,jl)   = v_ip(ji,jj,jl)   - zdvice
-                        
                         diag_dvpn_lid(ji,jj) = diag_dvpn_lid(ji,jj) - zdvice      ! diag
+
 !                       zvolp(ji,jj)     = zvolp(ji,jj)     - zdvice
 !                       zfpond(ji,jj)    = zfpond(ji,jj)    - zdvice
+
                         h_ip(ji,jj,jl)   = v_ip(ji,jj,jl) / a_ip(ji,jj,jl) ! MV - in principle, this is useless as h_ip is computed in icevar
                      ENDIF