- 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/SAS_SRC/nemogcm.F90
r8314 r8758 74 74 ! 75 75 #if defined key_agrif 76 CALL Agrif_Init_Grids() ! AGRIF: set the meshes76 CALL Agrif_Init_Grids() ! AGRIF: set the meshes 77 77 #endif 78 78 ! … … 167 167 INTEGER :: ilocal_comm ! local integer 168 168 INTEGER :: ios, inum ! - - 169 INTEGER :: iiarea, ijarea ! local integers 170 INTEGER :: iirest, ijrest ! local integers 169 171 CHARACTER(len=120), DIMENSION(30) :: cltxt, cltxt2, clnam 170 172 CHARACTER(len=80) :: clname … … 220 222 ENDIF 221 223 ! 222 jpk = jpkglo223 !224 #if defined key_agrif225 IF( .NOT. Agrif_Root() ) THEN ! AGRIF children: specific setting (cf. agrif_user.F90)226 jpiglo = nbcellsx + 2 + 2*nbghostcells227 jpjglo = nbcellsy + 2 + 2*nbghostcells228 jpi = ( jpiglo-2*jpreci + (jpni-1+0) ) / jpni + 2*jpreci229 jpj = ( jpjglo-2*jprecj + (jpnj-1+0) ) / jpnj + 2*jprecj230 nperio = 0231 jperio = 0232 ln_use_jattr = .false.233 ENDIF234 #endif235 224 ! 236 225 ! !--------------------------------------------! … … 289 278 #endif 290 279 ENDIF 280 ! 281 #if defined key_agrif 282 IF( .NOT. Agrif_Root() ) THEN ! AGRIF children: specific setting (cf. agrif_user.F90) 283 jpiglo = nbcellsx + 2 + 2*nbghostcells 284 jpjglo = nbcellsy + 2 + 2*nbghostcells 285 jpi = ( jpiglo-2*nn_hls + (jpni-1+0) ) / jpni + 2*nn_hls 286 jpj = ( jpjglo-2*nn_hls + (jpnj-1+0) ) / jpnj + 2*nn_hls 287 jpimax = jpi 288 jpjmax = jpj 289 nperio = 0 290 jperio = 0 291 ln_use_jattr = .false. 292 ENDIF 293 #endif 291 294 292 295 IF( Agrif_Root() ) THEN ! AGRIF mother: specific setting from jpni and jpnj 296 iiarea = 1 + MOD( narea - 1 , jpni ) 297 ijarea = 1 + ( narea - 1 ) / jpni 298 iirest = 1 + MOD( jpiglo - 2*nn_hls - 1 , jpni ) 299 ijrest = 1 + MOD( jpjglo - 2*nn_hls - 1 , jpnj ) 293 300 #if defined key_nemocice_decomp 294 jpi = ( nx_global+2-2*jpreci + (jpni-1) ) / jpni + 2*jpreci ! first dim. 295 jpj = ( ny_global+2-2*jprecj + (jpnj-1) ) / jpnj + 2*jprecj ! second dim. 301 jpi = ( nx_global+2-2*nn_hls + (jpni-1) ) / jpni + 2*nn_hls ! first dim. 302 jpj = ( ny_global+2-2*nn_hls + (jpnj-1) ) / jpnj + 2*nn_hls ! second dim. 303 jpimax = jpi 304 jpjmax = jpj 305 IF( iiarea == jpni ) jpi = jpiglo - (jpni - 1) * (jpi - 2*nn_hls) 306 IF( ijarea == jpnj ) jpj = jpjglo - (jpnj - 1) * (jpj - 2*nn_hls) 296 307 #else 297 jpi = ( jpiglo -2*jpreci + (jpni-1) ) / jpni + 2*jpreci ! first dim. 298 jpj = ( jpjglo -2*jprecj + (jpnj-1) ) / jpnj + 2*jprecj ! second dim. 299 #endif 300 ENDIF 308 jpi = ( jpiglo -2*nn_hls + (jpni-1) ) / jpni + 2*nn_hls ! first dim. 309 jpj = ( jpjglo -2*nn_hls + (jpnj-1) ) / jpnj + 2*nn_hls ! second dim. 310 jpimax = jpi 311 jpjmax = jpj 312 IF( iiarea > iirest ) jpi = jpi - 1 313 IF( ijarea > ijrest ) jpj = jpj - 1 314 #endif 315 ENDIF 316 317 jpk = jpkglo ! third dim 301 318 302 319 #if defined key_agrif … … 676 693 IF ((sxM .gt. sxT) .AND. (sxM .lt. dxT)) THEN 677 694 nsndto = nsndto + 1 678 695 isendto(nsndto) = jn 679 696 ELSEIF ((sxM .le. sxT) .AND. (dxM .ge. dxT)) THEN 680 697 nsndto = nsndto + 1
Note: See TracChangeset
for help on using the changeset viewer.