- 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/tests/WAD/MY_SRC/usrdef_zgr.F90
r12740 r13899 15 15 !!--------------------------------------------------------------------- 16 16 USE oce ! ocean variables 17 USE dom_oce , ONLY: ht_0, mi0, mi1, nimpp, njmpp, & 18 & mj0, mj1, glamt, gphit ! ocean space and time domain 17 USE dom_oce , ONLY: ht_0, mi0, mi1, mj0, mj1, glamt, gphit ! ocean space and time domain 19 18 USE usrdef_nam ! User defined : namelist variables 20 19 USE wet_dry , ONLY: rn_wdmin1, rn_wdmin2, rn_wdld ! Wetting and drying … … 244 243 ! at v-point: averaging zht 245 244 zhv = 0._wp 246 DO_2D _00_00245 DO_2D( 0, 0, 0, 0 ) 247 246 zhv(ji,jj) = 0.5_wp * ( zht(ji,jj) + zht(ji,jj+1) ) 248 247 END_2D … … 281 280 ht_0 = zht 282 281 k_bot(:,:) = jpkm1 * k_top(:,:) !* bottom ocean = jpk-1 (here use k_top as a land mask) 283 DO_2D _11_11282 DO_2D( 1, 1, 1, 1 ) 284 283 IF( zht(ji,jj) <= -(rn_wdld - rn_wdmin2)) THEN 285 284 k_bot(ji,jj) = 0 … … 290 289 ! !* terrain-following coordinate with e3.(k)=cst) 291 290 ! ! OVERFLOW case : identical with j-index (T=V, U=F) 292 DO_2D _00_00291 DO_2D( 0, 0, 0, 0 ) 293 292 z1_jpkm1 = 1._wp / REAL( k_bot(ji,jj) - k_top(ji,jj) + 1 , wp) 294 293 DO jk = 1, jpk
Note: See TracChangeset
for help on using the changeset viewer.