Changeset 9488


Ignore:
Timestamp:
2018-04-19T18:19:59+02:00 (3 years ago)
Author:
jpalmier
Message:

start with hard coded under-ice relax value

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/NERC/dev_r5518_GO6_under_ice_relax/NEMOGCM/NEMO/OPA_SRC/SBC/sbcssr.F90

    r6486 r9488  
    7777      REAL(wp) ::   zsrp     ! local scalar for unit conversion of rn_deds factor 
    7878      REAL(wp) ::   zerp_bnd ! local scalar for unit conversion of rn_epr_max factor 
     79      !! AXY (16-04-2018) scaling under-sea-ice relaxation 
     80      REAL(wp) ::   zice     ! local scalar to permit differential under sea-ice relaxation 
     81      REAL(wp) ::   zsrp_ice     !  local scalar for unit conversion of rn_deds factor - under ice 
     82      REAL(wp), PARAMETER ::   rn_deds_ice = -333.333     ! under ice relaxation term 
    7983      INTEGER  ::   ierror   ! return error code 
    8084      !! 
     
    107111            IF( nn_sssr == 1 ) THEN                                   !* Salinity damping term (salt flux only (sfx)) 
    108112               zsrp = rn_deds / rday                                  ! from [mm/day] to [kg/m2/s] 
     113               zsrp_ice = rn_deds_ice / rday                                  ! from [mm/day] to [kg/m2/s] 
    109114!CDIR COLLAPSE 
    110115               DO jj = 1, jpj 
    111116                  DO ji = 1, jpi 
    112                      zerp = zsrp * ( 1. - 2.*rnfmsk(ji,jj) )   &      ! No damping in vicinity of river mouths 
     117                     !! AXY (16-04-2018) Calculate ice-scaled zsrp term - from 
     118                     !!      rn_deds (ice-free) to rn_deds_ice (100% ice) 
     119                     zice = ((1. - fr_i(ji,jj)) * zsrp) +        & 
     120                        &   (fr_i(ji,jj) * zsrp_ice)  
     121                     zerp = zice * ( 1. - 2.*rnfmsk(ji,jj) )   &      ! No damping in vicinity of river mouths 
    113122                        &        * ( sss_m(ji,jj) - sf_sss(1)%fnow(ji,jj,1) )  
    114123                     sfx(ji,jj) = sfx(ji,jj) + zerp                 ! salt flux 
Note: See TracChangeset for help on using the changeset viewer.