Changeset 12636 for NEMO/branches/UKMO/NEMO_4.0.1_meltpond_improvements/src/OCE/SBC/sbccpl.F90

Timestamp:

2020-04-01T09:44:12+02:00 (4 years ago)

Author:

dancopsey

Message:

Add:

• Melt pond lids
• Melt pond maximum area and thickness
• Melt pond vertical flushing
• Area contributing to melt ponds depends on total ice fraction

File:

• NEMO/branches/UKMO/NEMO_4.0.1_meltpond_improvements/src/OCE/SBC/sbccpl.F90 (modified) (2 diffs)

NEMO/branches/UKMO/NEMO_4.0.1_meltpond_improvements/src/OCE/SBC/sbccpl.F90

@@ -2047 +2047 @@
       INTEGER, INTENT(in) ::   kt
       !
+      !
       INTEGER ::   ji, jj, jl   ! dummy loop indices
       INTEGER ::   isec, info   ! local integer
       REAL(wp) ::   zumax, zvmax
       REAL(wp), DIMENSION(jpi,jpj)     ::   zfr_l, ztmp1, ztmp2, zotx1, zoty1, zotz1, zitx1, zity1, zitz1
-      REAL(wp), DIMENSION(jpi,jpj,jpl) ::   ztmp3, ztmp4   
+      REAL(wp), DIMENSION(jpi,jpj,jpl) ::   ztmp3, ztmp4
+
       !!----------------------------------------------------------------------
       !
@@ -2256 +2258 @@
             SELECT CASE( sn_snd_mpnd%clcat )  
             CASE( 'yes' )  
-               ztmp3(:,:,1:jpl) =  a_ip(:,:,1:jpl)
-               ztmp4(:,:,1:jpl) =  v_ip(:,:,1:jpl)  
+
+               ztmp3(:,:,1:jpl) =  a_ip_eff(:,:,1:jpl)
+               ztmp4(:,:,1:jpl) =  h_ip(:,:,1:jpl)
+
             CASE( 'no' )  
                ztmp3(:,:,:) = 0.0