Changeset 12402 for NEMO/branches/UKMO/NEMO_4.0.1_fix_cpl/src/OCE/BDY/bdyice.F90

Timestamp:

2020-02-18T16:53:58+01:00 (4 years ago)

Author:

dancopsey

Message:

Add melt pond lid code.

File:

• NEMO/branches/UKMO/NEMO_4.0.1_fix_cpl/src/OCE/BDY/bdyice.F90 (modified) (5 diffs)

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

@@ -95 +95 @@
             ! exchange 3d arrays
             CALL lbc_lnk_multi( 'bdyice', a_i , 'T', 1., h_i , 'T', 1., h_s , 'T', 1., oa_i, 'T', 1. &
-                 &                      , a_ip, 'T', 1., v_ip, 'T', 1., s_i , 'T', 1., t_su, 'T', 1. &
+                 &                      , a_ip, 'T', 1., v_ip, 'T', 1., lh_ip,'T', 1., s_i , 'T', 1., t_su, 'T', 1. &
                  &                      , v_i , 'T', 1., v_s , 'T', 1., sv_i, 'T', 1.                &
                  &                      , kfillmode=jpfillnothing ,lsend=llsend1, lrecv=llrecv1      )
@@ -163 +163 @@
             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 +171 @@
                a_ip(ji,jj,jl) = 0._wp
                h_ip(ji,jj,jl) = 0._wp
+               lh_ip(ji,jj,jl) = 0._wp
             ENDIF
             !
@@ -231 +233 @@
                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 +283 @@
                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