Changeset 12371 for NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/ICE/icevar.F90

Timestamp:

2020-02-12T13:13:26+01:00 (4 years ago)

Author:

dancopsey

Message:

Add changes to match CICE settings.

File:

• NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/ICE/icevar.F90 (modified) (3 diffs)

NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/ICE/icevar.F90

@@ -317 +317 @@
       REAL(wp), PARAMETER :: zsi0 = 3.5_wp
       REAL(wp), PARAMETER :: zsi1 = 4.5_wp
+
+      REAL(wp), PARAMETER :: min_salin = 0.1_wp     ! Minimum salinity in sea ice
+      REAL(wp) ::            zn                     ! Fraction of ice depth
+      REAL(wp), PARAMETER :: msal    = 0.573_wp
+      REAL(wp), PARAMETER :: nsal    = 0.407_wp
+      REAL(wp), PARAMETER :: saltmax = 9.6_wp
+      ! REAL(wp), PARAMETER, DIMENSION(4) :: weights = (/0.06363419, 0.25545967, 0.33155539, 0.34935076/)
+      REAL(wp), PARAMETER, DIMENSION(4) :: weights = (/0.0703805, 0.25526203, 0.32860314, 0.34575433/)
+      
       !!-------------------------------------------------------------------
 
@@ -350 +359 @@
             DO jj = 1, jpj
                DO ji = 1, jpi
-                  zalpha(ji,jj,jl) = MAX(  0._wp , MIN( ( zsi1 - s_i(ji,jj,jl) ) * z1_dS , 1._wp )  )
+                  zalpha(ji,jj,jl) = 1._wp
                   !                             ! force a constant profile when SSS too low (Baltic Sea)
                   IF( 2._wp * s_i(ji,jj,jl) >= sss_m(ji,jj) )   zalpha(ji,jj,jl) = 0._wp  
@@ -363 +372 @@
                DO jj = 1, jpj
                   DO ji = 1, jpi
-                     !                          ! linear profile with 0 surface value
-                     zs0 = z_slope_s(ji,jj,jl) * ( REAL(jk,wp) - 0.5_wp ) * h_i(ji,jj,jl) * r1_nlay_i
+
+                     !  Copied over from CICE
+                     ! zn = (real(jk,kind=dbl_kind)-0.5_wp) / real(nlay_i,kind=dbl_kind)
+                     ! sz_i(ji,jj,jk,jl)=(saltmax/2.0_wp)*(1.0_wp-cos(pi*zn**(nsal/(msal+zn))))
+                     ! sz_i(ji,jj,jk,jl) = max(sz_i(ji,jj,jk,jl), min_salin)
+
+                     ! Weighting method to match CICE
+                     zs0 = (s_i(ji,jj,jl)*4.0_wp) * weights(jk)
                      zs  = zalpha(ji,jj,jl) * zs0 + ( 1._wp - zalpha(ji,jj,jl) ) * s_i(ji,jj,jl)     ! weighting the profile
                      sz_i(ji,jj,jk,jl) = MIN( rn_simax, MAX( zs, rn_simin ) )
+
+
+                     ! zs0 = z_slope_s(ji,jj,jl) * ( REAL(jk,wp) - 0.5_wp ) * h_i(ji,jj,jl) * r1_nlay_i
+                     ! zs  = zalpha(ji,jj,jl) * zs0 + ( 1._wp - zalpha(ji,jj,jl) ) * s_i(ji,jj,jl)     ! weighting the profile
+                     ! sz_i(ji,jj,jk,jl) = MIN( rn_simax, MAX( zs, rn_simin ) )
 
                      IF ( jl == 1 .AND. jj == 26 .AND. ji == 42 ) THEN
                        write(numout,*) 'icevar: jk, sz_i(ji,jj,jk,jl), s_i(ji,jj,jl), weights(jk) = ',jk, ' ',sz_i(ji,jj,jk,jl), ' ',s_i(ji,jj,jl), ' ',weights(jk)
                      ENDIF
+
                   END DO
                END DO