- Timestamp:
- 2020-11-27T17:26:33+01:00 (4 years ago)
- Location:
- NEMO/branches/2020/tickets_icb_1900
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/tickets_icb_1900
- Property svn:externals
-
NEMO/branches/2020/tickets_icb_1900/src/OCE/DYN/wet_dry.F90
r13237 r13899 57 57 REAL(wp), PUBLIC :: ssh_ref !: height of z=0 with respect to the geoid; 58 58 59 LOGICAL, PUBLIC :: ll_wd !: Wetting/drying activation switch if either ln_wd_il or ln_wd_dl59 LOGICAL, PUBLIC :: ll_wd = .FALSE. !: Wetting/drying activation switch (ln_wd_il or ln_wd_dl) <- default def if wad_init not called 60 60 61 61 PUBLIC wad_init ! initialisation routine called by step.F90 … … 111 111 112 112 r_rn_wdmin1 = 1 / rn_wdmin1 113 ll_wd = .FALSE.114 113 IF( ln_wd_il .OR. ln_wd_dl ) THEN 115 114 ll_wd = .TRUE. … … 174 173 ! 175 174 wdmask(:,:) = 1._wp 176 DO_2D _01_01175 DO_2D( 0, 1, 0, 1 ) 177 176 ! 178 177 IF( tmask(ji,jj,1) < 0.5_wp ) CYCLE ! we don't care about land cells … … 198 197 wdramp(:,:) = min((ht_0(:,:) + psshb1(:,:) - rn_wdmin1)/(rn_wdmin0 - rn_wdmin1),1.0_wp) 199 198 !jth assume don't need a lbc_lnk here 200 DO_2D _10_10199 DO_2D( 1, 0, 1, 0 ) 201 200 wdrampu(ji,jj) = MIN( wdramp(ji,jj) , wdramp(ji+1,jj) ) 202 201 wdrampv(ji,jj) = MIN( wdramp(ji,jj) , wdramp(ji,jj+1) ) … … 211 210 jflag = 0 ! flag indicating if any further iterations are needed 212 211 ! 213 DO_2D _01_01212 DO_2D( 0, 1, 0, 1 ) 214 213 IF( tmask(ji, jj, 1) < 0.5_wp ) CYCLE 215 214 IF( ht_0(ji,jj) > zdepwd ) CYCLE … … 307 306 zwdlmtv(:,:) = 1._wp 308 307 ! 309 DO_2D _01_01308 DO_2D( 0, 1, 0, 1 ) ! Horizontal Flux in u and v direction 310 309 ! 311 310 IF( tmask(ji, jj, 1 ) < 0.5_wp) CYCLE ! we don't care about land cells … … 333 332 jflag = 0 ! flag indicating if any further iterations are needed 334 333 ! 335 DO_2D _01_01334 DO_2D( 0, 1, 0, 1 ) 336 335 ! 337 336 IF( tmask(ji, jj, 1 ) < 0.5_wp ) CYCLE
Note: See TracChangeset
for help on using the changeset viewer.