Changeset 10170 for NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/WAD/MY_SRC/usrdef_zgr.F90
- Timestamp:
- 2018-10-03T16:49:50+02:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/WAD/MY_SRC/usrdef_zgr.F90
r10074 r10170 234 234 zhu(ji,:) = 0.5_wp * ( zht(ji,:) + zht(ji+1,:) ) 235 235 END DO 236 CALL lbc_lnk( zhu, 'U', 1. ) ! boundary condition: this mask the surrounding grid-points236 CALL lbc_lnk( 'usrdef_zgr', zhu, 'U', 1. ) ! boundary condition: this mask the surrounding grid-points 237 237 ! ! ==>>> set by hand non-zero value on first/last columns & rows 238 238 DO ji = mi0(1), mi1(1) ! first row of global domain only … … 247 247 zhv(:,jj) = 0.5_wp * ( zht(:,jj) + zht(:,jj+1) ) 248 248 END DO 249 CALL lbc_lnk( zhv, 'V', 1. ) ! boundary condition: this mask the surrounding grid-points249 CALL lbc_lnk( 'usrdef_zgr', zhv, 'V', 1. ) ! boundary condition: this mask the surrounding grid-points 250 250 DO jj = mj0(1), mj1(1) ! first row of global domain only 251 251 zhv(:,jj) = zht(:,jj) … … 272 272 273 273 274 CALL lbc_lnk( z2d, 'T', 1. ) ! closed basin since jperio = 0 (see userdef_nam.F90)274 CALL lbc_lnk( 'usrdef_zgr', z2d, 'T', 1. ) ! closed basin since jperio = 0 (see userdef_nam.F90) 275 275 k_top(:,:) = NINT( z2d(:,:) ) 276 276 ! … … 311 311 END DO 312 312 END DO 313 CALL lbc_lnk( pdept, 'T', 1. )314 CALL lbc_lnk( pdepw, 'T', 1. )315 CALL lbc_lnk( pe3t , 'T', 1. )316 CALL lbc_lnk( pe3w , 'T', 1. )317 CALL lbc_lnk( pe3u , 'U', 1. )318 CALL lbc_lnk( pe3uw, 'U', 1. )319 CALL lbc_lnk( pe3f , 'F', 1. )320 CALL lbc_lnk( pe3v , 'V', 1. )321 CALL lbc_lnk( pe3vw, 'V', 1. )313 CALL lbc_lnk( 'usrdef_zgr', pdept, 'T', 1. ) 314 CALL lbc_lnk( 'usrdef_zgr', pdepw, 'T', 1. ) 315 CALL lbc_lnk( 'usrdef_zgr', pe3t , 'T', 1. ) 316 CALL lbc_lnk( 'usrdef_zgr', pe3w , 'T', 1. ) 317 CALL lbc_lnk( 'usrdef_zgr', pe3u , 'U', 1. ) 318 CALL lbc_lnk( 'usrdef_zgr', pe3uw, 'U', 1. ) 319 CALL lbc_lnk( 'usrdef_zgr', pe3f , 'F', 1. ) 320 CALL lbc_lnk( 'usrdef_zgr', pe3v , 'V', 1. ) 321 CALL lbc_lnk( 'usrdef_zgr', pe3vw, 'V', 1. ) 322 322 WHERE( pe3t (:,:,:) == 0._wp ) pe3t (:,:,:) = 1._wp 323 323 WHERE( pe3u (:,:,:) == 0._wp ) pe3u (:,:,:) = 1._wp
Note: See TracChangeset
for help on using the changeset viewer.