Changeset 12636 for NEMO/branches/UKMO/NEMO_4.0.1_meltpond_improvements/src/OCE/BDY/bdyice.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/BDY/bdyice.F90 (modified) (5 diffs)

NEMO/branches/UKMO/NEMO_4.0.1_meltpond_improvements/src/OCE/BDY/bdyice.F90

@@ -98 +98 @@
                  &                      , v_i , 'T', 1., v_s , 'T', 1., sv_i, 'T', 1.                &
                  &                      , kfillmode=jpfillnothing ,lsend=llsend1, lrecv=llrecv1      )
+            CALL lbc_lnk_multi( 'bdyice', lh_ip,'T', 1., kfillmode=jpfillnothing ,lsend=llsend1, lrecv=llrecv1      )
             ! exchange 4d arrays :   third dimension = 1   and then   third dimension = jpk
             CALL lbc_lnk_multi( 'bdyice', t_s , 'T', 1., e_s , 'T', 1., kfillmode=jpfillnothing ,lsend=llsend1, lrecv=llrecv1 )
@@ -163 +164 @@
             a_ip(ji,jj,  jl) = ( a_ip(ji,jj,  jl) * zwgt1 + dta%aip(i_bdy,jl) * zwgt ) * tmask(ji,jj,1)  ! Ice  pond concentration
             h_ip(ji,jj,  jl) = ( h_ip(ji,jj,  jl) * zwgt1 + dta%hip(i_bdy,jl) * zwgt ) * tmask(ji,jj,1)  ! Ice  pond depth
+            lh_ip(ji,jj, jl) = ( lh_ip(ji,jj, jl) * zwgt1 + dta%hip(i_bdy,jl) * zwgt ) * tmask(ji,jj,1)  ! Ice  pond lid thickness
             !
             sz_i(ji,jj,:,jl) = s_i(ji,jj,jl)
@@ -170 +172 @@
                a_ip(ji,jj,jl) = 0._wp
                h_ip(ji,jj,jl) = 0._wp
+               lh_ip(ji,jj,jl) = 0._wp
             ENDIF
             !
@@ -231 +234 @@
                a_ip(ji,jj,  jl) = a_ip(ib,jb,  jl)
                h_ip(ji,jj,  jl) = h_ip(ib,jb,  jl)
+               lh_ip(ji,jj, jl) = lh_ip(ib,jb, jl)
                !
                sz_i(ji,jj,:,jl) = sz_i(ib,jb,:,jl)
@@ -280 +284 @@
                a_ip(ji,jj,  jl) = 0._wp
                v_ip(ji,jj,  jl) = 0._wp
+               lh_ip(ji,jj, jl) = 0._wp
                t_su(ji,jj,  jl) = rt0
                t_s (ji,jj,:,jl) = rt0