Changeset 12114 for NEMO/branches/2019
- Timestamp:
- 2019-12-07T15:06:58+01:00 (4 years ago)
- Location:
- NEMO/branches/2019/dev_r11643_ENHANCE-11_CEthe_Shaconemo_diags
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11643_ENHANCE-11_CEthe_Shaconemo_diags/cfgs/SHARED/namelist_ice_ref
r11998 r12114 60 60 rn_icebfr = 15. ! maximum bottom stress per unit volume [N/m3] 61 61 rn_lfrelax = 1.e-5 ! relaxation time scale to reach static friction [s-1] 62 rn_tensile = 0. 05! isotropic tensile strength [0-0.5??]62 rn_tensile = 0.2 ! isotropic tensile strength [0-0.5??] 63 63 / 64 64 !------------------------------------------------------------------------------ … … 91 91 !------------------------------------------------------------------------------ 92 92 ln_rhg_EVP = .true. ! EVP rheology 93 ln_aEVP = . true.! adaptive rheology (Kimmritz et al. 2016 & 2017)93 ln_aEVP = .false. ! adaptive rheology (Kimmritz et al. 2016 & 2017) 94 94 rn_creepl = 2.0e-9 ! creep limit [1/s] 95 95 rn_ecc = 2.0 ! eccentricity of the elliptical yield curve -
NEMO/branches/2019/dev_r11643_ENHANCE-11_CEthe_Shaconemo_diags/cfgs/SHARED/namelist_ref
r11998 r12114 392 392 ln_sssr_bnd = .true. ! flag to bound erp term (associated with nn_sssr=2) 393 393 rn_sssr_bnd = 4.e0 ! ABS(Max/Min) value of the damping erp term [mm/day] 394 nn_sssr_ice = 1 ! control of sea surface restoring under sea-ice 395 ! 0 = no restoration under ice : * (1-icefrac) 396 ! 1 = restoration everywhere 397 ! >1 = enhanced restoration under ice : 1+(nn_icedmp-1)*icefrac 394 398 395 399 cn_dir = './' ! root directory for the SST/SSS data location -
NEMO/branches/2019/dev_r11643_ENHANCE-11_CEthe_Shaconemo_diags/src/OCE/SBC/sbcssr.F90
r11995 r12114 43 43 LOGICAL :: ln_sssr_bnd ! flag to bound erp term 44 44 REAL(wp) :: rn_sssr_bnd ! ABS(Max./Min.) value of erp term [mm/day] 45 INTEGER :: nn_ icedmp! Control of restoring under ice45 INTEGER :: nn_sssr_ice ! Control of restoring under ice 46 46 47 47 REAL(wp) , ALLOCATABLE, DIMENSION(:) :: buffer ! Temporary buffer for exchange … … 102 102 ENDIF 103 103 ! 104 IF( nn_sssr /= 0 .AND. nn_ icedmp/= 1 ) THEN105 ! use fraction of ice ( fr_i ) to adjust relaxation under ice if nn_ icedmp.ne. 1106 ! n.b. coefice is initialised and fixed to 1._wp if nn_ icedmp= 1104 IF( nn_sssr /= 0 .AND. nn_sssr_ice /= 1 ) THEN 105 ! use fraction of ice ( fr_i ) to adjust relaxation under ice if nn_sssr_ice .ne. 1 106 ! n.b. coefice is initialised and fixed to 1._wp if nn_sssr_ice = 1 107 107 DO jj = 1, jpj 108 108 DO ji = 1, jpi 109 SELECT CASE ( nn_ icedmp)109 SELECT CASE ( nn_sssr_ice ) 110 110 CASE ( 0 ) ; coefice(ji,jj) = 1._wp - fr_i(ji,jj) ! no/reduced damping under ice 111 CASE DEFAULT ; coefice(ji,jj) = 1._wp + (nn_icedmp-1)*fr_i(ji,jj) ! reinforced damping (x nn_icedmp) under ice )111 CASE DEFAULT ; coefice(ji,jj) = 1._wp + ( nn_sssr_ice - 1 ) * fr_i(ji,jj) ! reinforced damping (x nn_sssr_ice) under ice ) 112 112 END SELECT 113 113 END DO … … 170 170 TYPE(FLD_N) :: sn_sst, sn_sss ! informations about the fields to be read 171 171 NAMELIST/namsbc_ssr/ cn_dir, nn_sstr, nn_sssr, rn_dqdt, rn_deds, sn_sst, & 172 & sn_sss, ln_sssr_bnd, rn_sssr_bnd, nn_ icedmp172 & sn_sss, ln_sssr_bnd, rn_sssr_bnd, nn_sssr_ice 173 173 INTEGER :: ios 174 174 !!---------------------------------------------------------------------- … … 198 198 WRITE(numout,*) ' flag to bound erp term ln_sssr_bnd = ', ln_sssr_bnd 199 199 WRITE(numout,*) ' ABS(Max./Min.) erp threshold rn_sssr_bnd = ', rn_sssr_bnd, ' mm/day' 200 WRITE(numout,*) ' Cntrl of surface restoration under ice nn_ icedmp = ', nn_icedmp200 WRITE(numout,*) ' Cntrl of surface restoration under ice nn_sssr_ice = ', nn_sssr_ice 201 201 WRITE(numout,*) ' ( 0 = no restoration under ice)' 202 202 WRITE(numout,*) ' ( 1 = restoration everywhere )' … … 232 232 ENDIF 233 233 ! 234 coefice(:,:) = 1._wp ! Initialise coefice to 1._wp ; will not need to be changed if nn_ icedmp=1234 coefice(:,:) = 1._wp ! Initialise coefice to 1._wp ; will not need to be changed if nn_sssr_ice=1 235 235 ! !* Initialize qrp and erp if no restoring 236 236 IF( nn_sstr /= 1 ) qrp(:,:) = 0._wp
Note: See TracChangeset
for help on using the changeset viewer.