Changeset 12740 for NEMO/trunk/tests/WAD
- Timestamp:
- 2020-04-12T11:03:06+02:00 (4 years ago)
- Location:
- NEMO/trunk/tests/WAD/MY_SRC
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/tests/WAD/MY_SRC/usrdef_hgr.F90
r10074 r12740 26 26 PUBLIC usr_def_hgr ! called by domhgr.F90 27 27 28 !! * Substitutions 29 # include "do_loop_substitute.h90" 28 30 !!---------------------------------------------------------------------- 29 31 !! NEMO/OCE 4.0 , NEMO Consortium (2018) … … 72 74 ! !== grid point position ==! (in kilometers) 73 75 zfact = rn_dx * 1.e-3 ! conversion in km 74 DO jj = 1, jpj 75 DO ji = 1, jpi ! longitude 76 plamt(ji,jj) = zfact * ( - 0.5 + REAL( ji-1 + nimpp-1 , wp ) ) 77 plamu(ji,jj) = zfact * ( REAL( ji-1 + nimpp-1 , wp ) ) 78 plamv(ji,jj) = plamt(ji,jj) 79 plamf(ji,jj) = plamu(ji,jj) 80 ! ! latitude 81 pphit(ji,jj) = zfact * ( - 0.5 + REAL( jj-1 + njmpp-1 , wp ) ) 82 pphiu(ji,jj) = pphit(ji,jj) 83 pphiv(ji,jj) = zfact * ( REAL( jj-1 + njmpp-1 , wp ) ) 84 pphif(ji,jj) = pphiv(ji,jj) 85 END DO 86 END DO 76 DO_2D_11_11 77 ! ! longitude 78 plamt(ji,jj) = zfact * ( - 0.5 + REAL( ji-1 + nimpp-1 , wp ) ) 79 plamu(ji,jj) = zfact * ( REAL( ji-1 + nimpp-1 , wp ) ) 80 plamv(ji,jj) = plamt(ji,jj) 81 plamf(ji,jj) = plamu(ji,jj) 82 ! ! latitude 83 pphit(ji,jj) = zfact * ( - 0.5 + REAL( jj-1 + njmpp-1 , wp ) ) 84 pphiu(ji,jj) = pphit(ji,jj) 85 pphiv(ji,jj) = zfact * ( REAL( jj-1 + njmpp-1 , wp ) ) 86 pphif(ji,jj) = pphiv(ji,jj) 87 END_2D 87 88 ! 88 89 ! !== Horizontal scale factors ==! (in meters) -
NEMO/trunk/tests/WAD/MY_SRC/usrdef_istate.F90
r10074 r12740 26 26 PUBLIC usr_def_istate ! called in istate.F90 27 27 28 !! * Substitutions 29 # include "do_loop_substitute.h90" 28 30 !!---------------------------------------------------------------------- 29 31 !! NEMO/OCE 4.0 , NEMO Consortium (2018) … … 174 176 ! Apply minimum wetdepth criterion 175 177 ! 176 do jj = 1,jpj 177 do ji = 1,jpi 178 IF( ht_0(ji,jj) + pssh(ji,jj) < rn_wdmin1 ) THEN 179 pssh(ji,jj) = ptmask(ji,jj,1)*( rn_wdmin1 - ht_0(ji,jj) ) 180 ENDIF 181 end do 182 end do 178 DO_2D_11_11 179 IF( ht_0(ji,jj) + pssh(ji,jj) < rn_wdmin1 ) THEN 180 pssh(ji,jj) = ptmask(ji,jj,1)*( rn_wdmin1 - ht_0(ji,jj) ) 181 ENDIF 182 END_2D 183 183 ! 184 184 END SUBROUTINE usr_def_istate -
NEMO/trunk/tests/WAD/MY_SRC/usrdef_zgr.F90
r12377 r12740 29 29 PUBLIC usr_def_zgr ! called by domzgr.F90 30 30 31 !! * Substitutions 32 # include "do_loop_substitute.h90" 31 33 !!---------------------------------------------------------------------- 32 34 !! NEMO/OCE 4.0 , NEMO Consortium (2018) … … 242 244 ! at v-point: averaging zht 243 245 zhv = 0._wp 244 DO jj = 1, jpjm1245 zhv( :,jj) = 0.5_wp * ( zht(:,jj) + zht(:,jj+1) )246 END DO246 DO_2D_00_00 247 zhv(ji,jj) = 0.5_wp * ( zht(ji,jj) + zht(ji,jj+1) ) 248 END_2D 247 249 CALL lbc_lnk( 'usrdef_zgr', zhv, 'V', 1. ) ! boundary condition: this mask the surrounding grid-points 248 250 DO jj = mj0(1), mj1(1) ! first row of global domain only … … 279 281 ht_0 = zht 280 282 k_bot(:,:) = jpkm1 * k_top(:,:) !* bottom ocean = jpk-1 (here use k_top as a land mask) 281 DO jj = 1, jpj 282 DO ji = 1, jpi 283 IF( zht(ji,jj) <= -(rn_wdld - rn_wdmin2)) THEN 284 k_bot(ji,jj) = 0 285 k_top(ji,jj) = 0 286 ENDIF 287 END DO 288 END DO 283 DO_2D_11_11 284 IF( zht(ji,jj) <= -(rn_wdld - rn_wdmin2)) THEN 285 k_bot(ji,jj) = 0 286 k_top(ji,jj) = 0 287 ENDIF 288 END_2D 289 289 ! 290 290 ! !* terrain-following coordinate with e3.(k)=cst) 291 291 ! ! OVERFLOW case : identical with j-index (T=V, U=F) 292 DO jj = 1, jpjm1 293 DO ji = 1, jpim1 294 z1_jpkm1 = 1._wp / REAL( k_bot(ji,jj) - k_top(ji,jj) + 1 , wp) 295 DO jk = 1, jpk 296 zwet = MAX( zht(ji,jj), rn_wdmin1 ) 297 pdept(ji,jj,jk) = zwet * z1_jpkm1 * ( REAL( jk , wp ) - 0.5_wp ) 298 pdepw(ji,jj,jk) = zwet * z1_jpkm1 * ( REAL( jk-1 , wp ) ) 299 pe3t (ji,jj,jk) = zwet * z1_jpkm1 300 pe3w (ji,jj,jk) = zwet * z1_jpkm1 301 zwet = MAX( zhu(ji,jj), rn_wdmin1 ) 302 pe3u (ji,jj,jk) = zwet * z1_jpkm1 303 pe3uw(ji,jj,jk) = zwet * z1_jpkm1 304 pe3f (ji,jj,jk) = zwet * z1_jpkm1 305 zwet = MAX( zhv(ji,jj), rn_wdmin1 ) 306 pe3v (ji,jj,jk) = zwet * z1_jpkm1 307 pe3vw(ji,jj,jk) = zwet * z1_jpkm1 308 END DO 309 END DO 310 END DO 292 DO_2D_00_00 293 z1_jpkm1 = 1._wp / REAL( k_bot(ji,jj) - k_top(ji,jj) + 1 , wp) 294 DO jk = 1, jpk 295 zwet = MAX( zht(ji,jj), rn_wdmin1 ) 296 pdept(ji,jj,jk) = zwet * z1_jpkm1 * ( REAL( jk , wp ) - 0.5_wp ) 297 pdepw(ji,jj,jk) = zwet * z1_jpkm1 * ( REAL( jk-1 , wp ) ) 298 pe3t (ji,jj,jk) = zwet * z1_jpkm1 299 pe3w (ji,jj,jk) = zwet * z1_jpkm1 300 zwet = MAX( zhu(ji,jj), rn_wdmin1 ) 301 pe3u (ji,jj,jk) = zwet * z1_jpkm1 302 pe3uw(ji,jj,jk) = zwet * z1_jpkm1 303 pe3f (ji,jj,jk) = zwet * z1_jpkm1 304 zwet = MAX( zhv(ji,jj), rn_wdmin1 ) 305 pe3v (ji,jj,jk) = zwet * z1_jpkm1 306 pe3vw(ji,jj,jk) = zwet * z1_jpkm1 307 END DO 308 END_2D 311 309 CALL lbc_lnk( 'usrdef_zgr', pdept, 'T', 1. ) 312 310 CALL lbc_lnk( 'usrdef_zgr', pdepw, 'T', 1. )
Note: See TracChangeset
for help on using the changeset viewer.