Changeset 15270 for NEMO/trunk/src/OCE/DOM/domain.F90
- Timestamp:
- 2021-09-17T16:27:55+02:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/OCE/DOM/domain.F90
r15267 r15270 20 20 !!---------------------------------------------------------------------- 21 21 !! dom_init : initialize the space and time domain 22 !! dom_glo : initialize global domain <--> local domain indices23 22 !! dom_nam : read and contral domain namelists 24 23 !! dom_ctl : control print for the ocean domain … … 124 123 ! !== Reference coordinate system ==! 125 124 ! 126 ! CALL dom_glo ! global domain versus local domain127 125 CALL dom_nam ! read namelist ( namrun, namdom ) 128 126 CALL dom_tile_init ! Tile domain … … 244 242 ! 245 243 END SUBROUTINE dom_init 246 247 248 SUBROUTINE dom_glo249 !!----------------------------------------------------------------------250 !! *** ROUTINE dom_glo ***251 !!252 !! ** Purpose : initialization of global domain <--> local domain indices253 !!254 !! ** Method :255 !!256 !! ** Action : - mig , mjg : local domain indices ==> global domain, including halos, indices257 !! - mig0, mjg0: local domain indices ==> global domain, excluding halos, indices258 !! - mi0 , mi1 : global domain indices ==> local domain indices259 !! - mj0 , mj1 (if global point not in the local domain ==> mi0>mi1 and/or mj0>mj1)260 !!----------------------------------------------------------------------261 INTEGER :: ji, jj ! dummy loop argument262 !!----------------------------------------------------------------------263 !264 DO ji = 1, jpi ! local domain indices ==> global domain indices, including halos265 mig(ji) = ji + nimpp - 1266 END DO267 DO jj = 1, jpj268 mjg(jj) = jj + njmpp - 1269 END DO270 ! ! local domain indices ==> global domain indices, excluding halos271 !272 mig0(:) = mig(:) - nn_hls273 mjg0(:) = mjg(:) - nn_hls274 ! ! global domain, including halos, indices ==> local domain indices275 ! ! (return (m.0,m.1)=(1,0) if data domain gridpoint is to the west/south of the276 ! ! local domain, or (m.0,m.1)=(jp.+1,jp.) to the east/north of local domain.277 DO ji = 1, jpiglo278 mi0(ji) = MAX( 1 , MIN( ji - nimpp + 1, jpi+1 ) )279 mi1(ji) = MAX( 0 , MIN( ji - nimpp + 1, jpi ) )280 END DO281 DO jj = 1, jpjglo282 mj0(jj) = MAX( 1 , MIN( jj - njmpp + 1, jpj+1 ) )283 mj1(jj) = MAX( 0 , MIN( jj - njmpp + 1, jpj ) )284 END DO285 IF(lwp) THEN ! control print286 WRITE(numout,*)287 WRITE(numout,*) 'dom_glo : domain: global <<==>> local '288 WRITE(numout,*) '~~~~~~~ '289 WRITE(numout,*) ' global domain: jpiglo = ', jpiglo, ' jpjglo = ', jpjglo, ' jpkglo = ', jpkglo290 WRITE(numout,*) ' local domain: jpi = ', jpi , ' jpj = ', jpj , ' jpk = ', jpk291 WRITE(numout,*)292 ENDIF293 !294 END SUBROUTINE dom_glo295 244 296 245
Note: See TracChangeset
for help on using the changeset viewer.