Changeset 12114 for NEMO/branches/2019/dev_r11643_ENHANCE-11_CEthe_Shaconemo_diags/src/OCE/SBC/sbcssr.F90
- Timestamp:
- 2019-12-07T15:06:58+01:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
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.