Changeset 9488 for branches/NERC

Timestamp:

2018-04-19T18:19:59+02:00 (6 years ago)

Author:

jpalmier

Message:

start with hard coded under-ice relax value

File:

• branches/NERC/dev_r5518_GO6_under_ice_relax/NEMOGCM/NEMO/OPA_SRC/SBC/sbcssr.F90 (modified) (2 diffs)

branches/NERC/dev_r5518_GO6_under_ice_relax/NEMOGCM/NEMO/OPA_SRC/SBC/sbcssr.F90

@@ -77 +77 @@
       REAL(wp) ::   zsrp     ! local scalar for unit conversion of rn_deds factor
       REAL(wp) ::   zerp_bnd ! local scalar for unit conversion of rn_epr_max factor
+      !! AXY (16-04-2018) scaling under-sea-ice relaxation
+      REAL(wp) ::   zice     ! local scalar to permit differential under sea-ice relaxation
+      REAL(wp) ::   zsrp_ice     !  local scalar for unit conversion of rn_deds factor - under ice
+      REAL(wp), PARAMETER ::   rn_deds_ice = -333.333     ! under ice relaxation term
       INTEGER  ::   ierror   ! return error code
       !!
@@ -107 +111 @@
             IF( nn_sssr == 1 ) THEN                                   !* Salinity damping term (salt flux only (sfx))
                zsrp = rn_deds / rday                                  ! from [mm/day] to [kg/m2/s]
+               zsrp_ice = rn_deds_ice / rday                                  ! from [mm/day] to [kg/m2/s]
 !CDIR COLLAPSE
                DO jj = 1, jpj
                   DO ji = 1, jpi
-                     zerp = zsrp * ( 1. - 2.*rnfmsk(ji,jj) )   &      ! No damping in vicinity of river mouths
+                     !! AXY (16-04-2018) Calculate ice-scaled zsrp term - from
+                     !!      rn_deds (ice-free) to rn_deds_ice (100% ice)
+                     zice = ((1. - fr_i(ji,jj)) * zsrp) +        &
+                        &   (fr_i(ji,jj) * zsrp_ice) 
+                     zerp = zice * ( 1. - 2.*rnfmsk(ji,jj) )   &      ! No damping in vicinity of river mouths
                         &        * ( sss_m(ji,jj) - sf_sss(1)%fnow(ji,jj,1) ) 
                      sfx(ji,jj) = sfx(ji,jj) + zerp                 ! salt flux