- Timestamp:
- 2017-11-20T17:29:54+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r8126_ROBUST08_no_ghost/NEMOGCM/NEMO/OPA_SRC/nemogcm.F90
r8314 r8758 239 239 INTEGER :: ji ! dummy loop indices 240 240 INTEGER :: ios, ilocal_comm ! local integer 241 INTEGER :: iiarea, ijarea ! local integers 242 INTEGER :: iirest, ijrest ! local integers 241 243 CHARACTER(len=120), DIMENSION(30) :: cltxt, cltxt2, clnam 242 244 ! … … 282 284 ENDIF 283 285 ! 284 jpk = jpkglo285 !286 #if defined key_agrif287 IF( .NOT. Agrif_Root() ) THEN ! AGRIF children: specific setting (cf. agrif_user.F90)288 jpiglo = nbcellsx + 2 + 2*nbghostcells289 jpjglo = nbcellsy + 2 + 2*nbghostcells290 jpi = ( jpiglo-2*jpreci + (jpni-1+0) ) / jpni + 2*jpreci291 jpj = ( jpjglo-2*jprecj + (jpnj-1+0) ) / jpnj + 2*jprecj292 nperio = 0293 jperio = 0294 ln_use_jattr = .false.295 ENDIF296 #endif297 286 ! 298 287 ! !--------------------------------------------! … … 353 342 #endif 354 343 ENDIF 344 ! 345 #if defined key_agrif 346 IF( .NOT. Agrif_Root() ) THEN ! AGRIF children: specific setting (cf. agrif_user.F90) 347 jpiglo = nbcellsx + 2 + 2*nbghostcells 348 jpjglo = nbcellsy + 2 + 2*nbghostcells 349 jpi = ( jpiglo-2*nn_hls + (jpni-1+0) ) / jpni + 2*nn_hls 350 jpj = ( jpjglo-2*nn_hls + (jpnj-1+0) ) / jpnj + 2*nn_hls 351 jpimax = jpi 352 jpjmax = jpj 353 nperio = 0 354 jperio = 0 355 ln_use_jattr = .false. 356 ENDIF 357 #endif 355 358 356 359 IF( Agrif_Root() ) THEN ! AGRIF mother: specific setting from jpni and jpnj 360 iiarea = 1 + MOD( narea - 1 , jpni ) 361 ijarea = 1 + ( narea - 1 ) / jpni 362 iirest = 1 + MOD( jpiglo - 2*nn_hls - 1 , jpni ) 363 ijrest = 1 + MOD( jpjglo - 2*nn_hls - 1 , jpnj ) 357 364 #if defined key_nemocice_decomp 358 jpi = ( nx_global+2-2*jpreci + (jpni-1) ) / jpni + 2*jpreci ! first dim. 359 jpj = ( ny_global+2-2*jprecj + (jpnj-1) ) / jpnj + 2*jprecj ! second dim. 365 jpi = ( nx_global+2-2*nn_hls + (jpni-1) ) / jpni + 2*nn_hls ! first dim. 366 jpj = ( ny_global+2-2*nn_hls + (jpnj-1) ) / jpnj + 2*nn_hls ! second dim. 367 jpimax = jpi 368 jpjmax = jpj 369 IF( iiarea == jpni ) jpi = jpiglo - (jpni - 1) * (jpi - 2*nn_hls) 370 IF( ijarea == jpnj ) jpj = jpjglo - (jpnj - 1) * (jpj - 2*nn_hls) 360 371 #else 361 jpi = ( jpiglo -2*jpreci + (jpni-1) ) / jpni + 2*jpreci ! first dim. 362 jpj = ( jpjglo -2*jprecj + (jpnj-1) ) / jpnj + 2*jprecj ! second dim. 363 #endif 364 ENDIF 365 366 !!gm ??? why here it has already been done in line 301 ! 372 jpi = ( jpiglo -2*nn_hls + (jpni-1) ) / jpni + 2*nn_hls ! first dim. 373 jpj = ( jpjglo -2*nn_hls + (jpnj-1) ) / jpnj + 2*nn_hls ! second dim. 374 jpimax = jpi 375 jpjmax = jpj 376 IF( iiarea > iirest ) jpi = jpi - 1 377 IF( ijarea > ijrest ) jpj = jpj - 1 378 #endif 379 ENDIF 380 367 381 jpk = jpkglo ! third dim 368 !!gm end369 382 370 383 #if defined key_agrif … … 414 427 ! ! Domain decomposition 415 428 CALL mpp_init 429 IF( ln_nnogather ) CALL nemo_northcomms! northfold neighbour lists (must be done after the masks are defined) 416 430 ! 417 431 IF( nn_timing == 1 ) CALL timing_init … … 424 438 CALL dom_init ! Domain 425 439 IF( ln_crs ) CALL crs_init ! coarsened grid: domain initialization 426 IF( ln_nnogather ) CALL nemo_northcomms! northfold neighbour lists (must be done after the masks are defined)440 !IF( ln_nnogather ) CALL nemo_northcomms! northfold neighbour lists (must be done after the masks are defined) 427 441 IF( ln_ctl ) CALL prt_ctl_init ! Print control 428 442
Note: See TracChangeset
for help on using the changeset viewer.