Changeset 3919 for branches/2012
- Timestamp:
- 2013-06-13T12:52:53+02:00 (11 years ago)
- Location:
- branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO
- Files:
-
- 4 edited
-
NST_SRC/agrif_opa_sponge.F90 (modified) (6 diffs)
-
NST_SRC/agrif_opa_update.F90 (modified) (5 diffs)
-
NST_SRC/agrif_top_sponge.F90 (modified) (3 diffs)
-
OPA_SRC/LBC/lib_mpp.F90 (modified) (23 diffs)
Legend:
- Unmodified
- Added
- Removed
-
branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/NST_SRC/agrif_opa_sponge.F90
r3917 r3919 35 35 REAL(wp), POINTER, DIMENSION(:,: ) :: localviscsponge 36 36 REAL(wp), POINTER, DIMENSION(:,: ) :: ztu, ztv 37 REAL(wp), POINTER, DIMENSION(:,:,:,:) :: ztab 4d37 REAL(wp), POINTER, DIMENSION(:,:,:,:) :: ztab 38 38 REAL(wp), POINTER, DIMENSION(:,:,:,:) :: tsbdiff 39 39 40 40 #if defined SPONGE 41 41 CALL wrk_alloc( jpi, jpj, localviscsponge, ztu, ztv ) 42 CALL wrk_alloc( jpi, jpj, jpk, jpts, ztab 4d, tsbdiff )42 CALL wrk_alloc( jpi, jpj, jpk, jpts, ztab, tsbdiff ) 43 43 44 44 timecoeff = REAL(Agrif_NbStepint(),wp)/Agrif_rhot() … … 46 46 Agrif_SpecialValue=0. 47 47 Agrif_UseSpecialValue = .TRUE. 48 ztab 4d= 0.e049 CALL Agrif_Bc_Variable(ztab 4d, tsa_id,calledweight=timecoeff,procname=interptsn)48 ztab = 0.e0 49 CALL Agrif_Bc_Variable(ztab, tsa_id,calledweight=timecoeff,procname=interptsn) 50 50 Agrif_UseSpecialValue = .FALSE. 51 51 52 tsbdiff(:,:,:,:) = tsb(:,:,:,:) - ztab 4d(:,:,:,:)52 tsbdiff(:,:,:,:) = tsb(:,:,:,:) - ztab(:,:,:,:) 53 53 54 54 spongearea = 2 + 2 * Agrif_irhox() … … 148 148 149 149 CALL wrk_dealloc( jpi, jpj, localviscsponge, ztu, ztv ) 150 CALL wrk_dealloc( jpi, jpj, jpk, jpts, ztab 4d, tsbdiff )150 CALL wrk_dealloc( jpi, jpj, jpk, jpts, ztab, tsbdiff ) 151 151 #endif 152 152 … … 166 166 REAL(wp), POINTER, DIMENSION(:,:,:) :: ubdiff, vbdiff 167 167 REAL(wp), POINTER, DIMENSION(:,:,:) :: rotdiff, hdivdiff 168 REAL(wp), POINTER, DIMENSION(:,:,:) :: ztab 3d168 REAL(wp), POINTER, DIMENSION(:,:,:) :: ztab 169 169 170 170 #if defined SPONGE 171 171 CALL wrk_alloc( jpi, jpj, localviscsponge ) 172 CALL wrk_alloc( jpi, jpj, jpk, ztab 3d, ubdiff, vbdiff, rotdiff, hdivdiff )172 CALL wrk_alloc( jpi, jpj, jpk, ztab, ubdiff, vbdiff, rotdiff, hdivdiff ) 173 173 174 174 timecoeff = REAL(Agrif_NbStepint(),wp)/Agrif_rhot() … … 176 176 Agrif_SpecialValue=0. 177 177 Agrif_UseSpecialValue = ln_spc_dyn 178 ztab 3d= 0.e0179 CALL Agrif_Bc_Variable(ztab 3d, ua_id,calledweight=timecoeff,procname=interpun)178 ztab = 0.e0 179 CALL Agrif_Bc_Variable(ztab, ua_id,calledweight=timecoeff,procname=interpun) 180 180 Agrif_UseSpecialValue = .FALSE. 181 181 182 ubdiff(:,:,:) = (ub(:,:,:) - ztab 3d(:,:,:))*umask(:,:,:)183 184 ztab 3d= 0.e0182 ubdiff(:,:,:) = (ub(:,:,:) - ztab(:,:,:))*umask(:,:,:) 183 184 ztab = 0.e0 185 185 Agrif_SpecialValue=0. 186 186 Agrif_UseSpecialValue = ln_spc_dyn 187 CALL Agrif_Bc_Variable(ztab 3d, va_id,calledweight=timecoeff,procname=interpvn)187 CALL Agrif_Bc_Variable(ztab, va_id,calledweight=timecoeff,procname=interpvn) 188 188 Agrif_UseSpecialValue = .FALSE. 189 189 190 vbdiff(:,:,:) = (vb(:,:,:) - ztab 3d(:,:,:))*vmask(:,:,:)190 vbdiff(:,:,:) = (vb(:,:,:) - ztab(:,:,:))*vmask(:,:,:) 191 191 192 192 spongearea = 2 + 2 * Agrif_irhox() … … 318 318 ! ! =============== 319 319 CALL wrk_dealloc( jpi, jpj, localviscsponge ) 320 CALL wrk_dealloc( jpi, jpj, jpk, ztab 3d, ubdiff, vbdiff, rotdiff, hdivdiff )320 CALL wrk_dealloc( jpi, jpj, jpk, ztab, ubdiff, vbdiff, rotdiff, hdivdiff ) 321 321 322 322 #endif -
branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/NST_SRC/agrif_opa_update.F90
r3917 r3919 32 32 !! 33 33 INTEGER, INTENT(in) :: kt 34 REAL(wp), POINTER, DIMENSION(:,:,:,:) :: ztab 4d34 REAL(wp), POINTER, DIMENSION(:,:,:,:) :: ztab 35 35 36 36 37 37 IF((Agrif_NbStepint() .NE. (Agrif_irhot()-1)).AND.(kt /= 0)) RETURN 38 38 #if defined TWO_WAY 39 CALL wrk_alloc( jpi, jpj, jpk, jpts, ztab 4d)39 CALL wrk_alloc( jpi, jpj, jpk, jpts, ztab ) 40 40 41 41 Agrif_UseSpecialValueInUpdate = .TRUE. … … 43 43 44 44 IF (MOD(nbcline,nbclineupdate) == 0) THEN 45 CALL Agrif_Update_Variable(ztab 4d,tsn_id, procname=updateTS)46 ELSE 47 CALL Agrif_Update_Variable(ztab 4d,tsn_id,locupdate=(/0,2/), procname=updateTS)45 CALL Agrif_Update_Variable(ztab,tsn_id, procname=updateTS) 46 ELSE 47 CALL Agrif_Update_Variable(ztab,tsn_id,locupdate=(/0,2/), procname=updateTS) 48 48 ENDIF 49 49 50 50 Agrif_UseSpecialValueInUpdate = .FALSE. 51 51 52 CALL wrk_dealloc( jpi, jpj, jpk, jpts, ztab 4d)52 CALL wrk_dealloc( jpi, jpj, jpk, jpts, ztab ) 53 53 #endif 54 54 … … 62 62 INTEGER, INTENT(in) :: kt 63 63 REAL(wp), POINTER, DIMENSION(:,:) :: ztab2d 64 REAL(wp), POINTER, DIMENSION(:,:,:) :: ztab 3d64 REAL(wp), POINTER, DIMENSION(:,:,:) :: ztab 65 65 66 66 … … 68 68 #if defined TWO_WAY 69 69 CALL wrk_alloc( jpi, jpj, ztab2d ) 70 CALL wrk_alloc( jpi, jpj, jpk, ztab 3d)70 CALL wrk_alloc( jpi, jpj, jpk, ztab ) 71 71 72 72 IF (mod(nbcline,nbclineupdate) == 0) THEN 73 CALL Agrif_Update_Variable(ztab 3d,un_id,procname = updateU)74 CALL Agrif_Update_Variable(ztab 3d,vn_id,procname = updateV)75 ELSE 76 CALL Agrif_Update_Variable(ztab 3d,un_id,locupdate=(/0,1/),procname = updateU)77 CALL Agrif_Update_Variable(ztab 3d,vn_id,locupdate=(/0,1/),procname = updateV)73 CALL Agrif_Update_Variable(ztab,un_id,procname = updateU) 74 CALL Agrif_Update_Variable(ztab,vn_id,procname = updateV) 75 ELSE 76 CALL Agrif_Update_Variable(ztab,un_id,locupdate=(/0,1/),procname = updateU) 77 CALL Agrif_Update_Variable(ztab,vn_id,locupdate=(/0,1/),procname = updateV) 78 78 ENDIF 79 79 … … 89 89 90 90 CALL wrk_dealloc( jpi, jpj, ztab2d ) 91 CALL wrk_dealloc( jpi, jpj, jpk, ztab 3d)91 CALL wrk_dealloc( jpi, jpj, jpk, ztab ) 92 92 93 93 !Done in step -
branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/NST_SRC/agrif_top_sponge.F90
r3917 r3919 36 36 REAL(wp) :: ztra, zabe1, zabe2, zbtr 37 37 REAL(wp), POINTER, DIMENSION(:,:) :: localviscsponge 38 REAL(wp), POINTER, DIMENSION(:,:,:,:) :: trbdiff, ztru, ztrv, ztab 4d38 REAL(wp), POINTER, DIMENSION(:,:,:,:) :: trbdiff, ztru, ztrv, ztab 39 39 40 40 #if defined SPONGE_TOP 41 41 CALL wrk_alloc( jpi, jpj, localviscsponge ) 42 CALL wrk_alloc( jpi, jpj, jpk, jptra, trbdiff, ztru, ztrv, ztab 4d)42 CALL wrk_alloc( jpi, jpj, jpk, jptra, trbdiff, ztru, ztrv, ztab ) 43 43 44 44 timecoeff = REAL(Agrif_NbStepint(),wp)/Agrif_rhot() … … 46 46 Agrif_SpecialValue=0. 47 47 Agrif_UseSpecialValue = .TRUE. 48 ztab 4d= 0.e049 CALL Agrif_Bc_Variable(ztab 4d, tra_id,calledweight=timecoeff,procname=interptrn)48 ztab = 0.e0 49 CALL Agrif_Bc_Variable(ztab, tra_id,calledweight=timecoeff,procname=interptrn) 50 50 Agrif_UseSpecialValue = .FALSE. 51 51 52 trbdiff(:,:,:,:) = trb(:,:,:,:) - ztab 4d(:,:,:,:)52 trbdiff(:,:,:,:) = trb(:,:,:,:) - ztab(:,:,:,:) 53 53 54 54 spongearea = 2 + 2 * Agrif_irhox() … … 147 147 148 148 CALL wrk_dealloc( jpi, jpj, localviscsponge ) 149 CALL wrk_dealloc( jpi, jpj, jpk, jptra, trbdiff, ztru, ztrv, ztab 4d)149 CALL wrk_dealloc( jpi, jpj, jpk, jptra, trbdiff, ztru, ztrv, ztab ) 150 150 151 151 #endif -
branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/OPA_SRC/LBC/lib_mpp.F90
r3800 r3919 159 159 160 160 ! Arrays used in mpp_lbc_north_3d() 161 REAL(wp), DIMENSION(:,:,:) , ALLOCATABLE, SAVE :: ztab, znorthloc162 REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE, SAVE :: znorthgloio163 REAL(wp), DIMENSION(:,:,:) , ALLOCATABLE, SAVE :: zfoldwk ! Workspace for message transfers avoiding mpi_allgather161 REAL(wp), DIMENSION(:,:,:) , ALLOCATABLE, SAVE :: tab_3d, xnorthloc 162 REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE, SAVE :: xnorthgloio 163 REAL(wp), DIMENSION(:,:,:) , ALLOCATABLE, SAVE :: foldwk ! Workspace for message transfers avoiding mpi_allgather 164 164 165 165 ! Arrays used in mpp_lbc_north_2d() 166 REAL(wp), DIMENSION(:,:) , ALLOCATABLE, SAVE :: ztab_2d, znorthloc_2d167 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, SAVE :: znorthgloio_2d168 REAL(wp), DIMENSION(:,:) , ALLOCATABLE, SAVE :: zfoldwk_2d ! Workspace for message transfers avoiding mpi_allgather166 REAL(wp), DIMENSION(:,:) , ALLOCATABLE, SAVE :: tab_2d, xnorthloc_2d 167 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, SAVE :: xnorthgloio_2d 168 REAL(wp), DIMENSION(:,:) , ALLOCATABLE, SAVE :: foldwk_2d ! Workspace for message transfers avoiding mpi_allgather 169 169 170 170 ! Arrays used in mpp_lbc_north_e() 171 REAL(wp), DIMENSION(:,:) , ALLOCATABLE, SAVE :: ztab_e, znorthloc_e172 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, SAVE :: znorthgloio_e171 REAL(wp), DIMENSION(:,:) , ALLOCATABLE, SAVE :: tab_e, xnorthloc_e 172 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, SAVE :: xnorthgloio_e 173 173 174 174 ! North fold arrays used to minimise the use of allgather operations. Set in nemo_northcomms (nemogcm) so need to be public … … 209 209 & tr2we(1-jpr2dj:jpj+jpr2dj,jpreci+jpr2di,2) , & 210 210 ! 211 & ztab(jpiglo,4,jpk) , znorthloc(jpi,4,jpk) , znorthgloio(jpi,4,jpk,jpni) , &212 & zfoldwk(jpi,4,jpk) , &213 ! 214 & ztab_2d(jpiglo,4) , znorthloc_2d(jpi,4) , znorthgloio_2d(jpi,4,jpni) , &215 & zfoldwk_2d(jpi,4) , &216 ! 217 & ztab_e(jpiglo,4+2*jpr2dj) , znorthloc_e(jpi,4+2*jpr2dj) , znorthgloio_e(jpi,4+2*jpr2dj,jpni) , &211 & tab_3d(jpiglo,4,jpk) , xnorthloc(jpi,4,jpk) , xnorthgloio(jpi,4,jpk,jpni) , & 212 & foldwk(jpi,4,jpk) , & 213 ! 214 & tab_2d(jpiglo,4) , xnorthloc_2d(jpi,4) , xnorthgloio_2d(jpi,4,jpni) , & 215 & foldwk_2d(jpi,4) , & 216 ! 217 & tab_e(jpiglo,4+2*jpr2dj) , xnorthloc_e(jpi,4+2*jpr2dj) , xnorthgloio_e(jpi,4+2*jpr2dj,jpni) , & 218 218 ! 219 219 & STAT=lib_mpp_alloc ) … … 2229 2229 ityp = -1 2230 2230 ijpjm1 = 3 2231 ztab(:,:,:) = 0.e02232 ! 2233 DO jj = nlcj - ijpj +1, nlcj ! put in znorthloc the last 4 jlines of pt3d2231 tab_3d(:,:,:) = 0.e0 2232 ! 2233 DO jj = nlcj - ijpj +1, nlcj ! put in xnorthloc the last 4 jlines of pt3d 2234 2234 ij = jj - nlcj + ijpj 2235 znorthloc(:,ij,:) = pt3d(:,jj,:)2235 xnorthloc(:,ij,:) = pt3d(:,jj,:) 2236 2236 END DO 2237 2237 ! 2238 ! ! Build in procs of ncomm_north the znorthgloio2238 ! ! Build in procs of ncomm_north the xnorthgloio 2239 2239 itaille = jpi * jpk * ijpj 2240 2240 IF ( l_north_nogather ) THEN … … 2246 2246 ij = jj - nlcj + ijpj 2247 2247 DO ji = 1, nlci 2248 ztab(ji+nimpp-1,ij,:) = pt3d(ji,jj,:)2248 tab_3d(ji+nimpp-1,ij,:) = pt3d(ji,jj,:) 2249 2249 END DO 2250 2250 END DO … … 2271 2271 2272 2272 DO jr = 1,nsndto(ityp) 2273 CALL mppsend(5, znorthloc, itaille, isendto(jr,ityp), ml_req_nf(jr) )2273 CALL mppsend(5, xnorthloc, itaille, isendto(jr,ityp), ml_req_nf(jr) ) 2274 2274 END DO 2275 2275 DO jr = 1,nsndto(ityp) 2276 CALL mpprecv(5, zfoldwk, itaille, isendto(jr,ityp))2276 CALL mpprecv(5, foldwk, itaille, isendto(jr,ityp)) 2277 2277 iproc = isendto(jr,ityp) + 1 2278 2278 ildi = nldit (iproc) … … 2281 2281 DO jj = 1, ijpj 2282 2282 DO ji = ildi, ilei 2283 ztab(ji+iilb-1,jj,:) = zfoldwk(ji,jj,:)2283 tab_3d(ji+iilb-1,jj,:) = foldwk(ji,jj,:) 2284 2284 END DO 2285 2285 END DO … … 2296 2296 2297 2297 IF ( ityp .lt. 0 ) THEN 2298 CALL MPI_ALLGATHER( znorthloc , itaille, MPI_DOUBLE_PRECISION, &2299 & znorthgloio, itaille, MPI_DOUBLE_PRECISION, ncomm_north, ierr )2298 CALL MPI_ALLGATHER( xnorthloc , itaille, MPI_DOUBLE_PRECISION, & 2299 & xnorthgloio, itaille, MPI_DOUBLE_PRECISION, ncomm_north, ierr ) 2300 2300 ! 2301 2301 DO jr = 1, ndim_rank_north ! recover the global north array … … 2306 2306 DO jj = 1, ijpj 2307 2307 DO ji = ildi, ilei 2308 ztab(ji+iilb-1,jj,:) = znorthgloio(ji,jj,:,jr)2308 tab_3d(ji+iilb-1,jj,:) = xnorthgloio(ji,jj,:,jr) 2309 2309 END DO 2310 2310 END DO … … 2312 2312 ENDIF 2313 2313 ! 2314 ! The ztabarray has been either:2314 ! The tab_3d array has been either: 2315 2315 ! a. Fully populated by the mpi_allgather operation or 2316 2316 ! b. Had the active points for this domain and northern neighbours populated … … 2319 2319 ! this domain will be identical. 2320 2320 ! 2321 CALL lbc_nfd( ztab, cd_type, psgn ) ! North fold boundary condition2321 CALL lbc_nfd( tab_3d, cd_type, psgn ) ! North fold boundary condition 2322 2322 ! 2323 2323 DO jj = nlcj-ijpj+1, nlcj ! Scatter back to pt3d 2324 2324 ij = jj - nlcj + ijpj 2325 2325 DO ji= 1, nlci 2326 pt3d(ji,jj,:) = ztab(ji+nimpp-1,ij,:)2326 pt3d(ji,jj,:) = tab_3d(ji+nimpp-1,ij,:) 2327 2327 END DO 2328 2328 END DO … … 2361 2361 ityp = -1 2362 2362 ijpjm1 = 3 2363 ztab_2d(:,:) = 0.e02364 ! 2365 DO jj = nlcj-ijpj+1, nlcj ! put in znorthloc_2d the last 4 jlines of pt2d2363 tab_2d(:,:) = 0.e0 2364 ! 2365 DO jj = nlcj-ijpj+1, nlcj ! put in xnorthloc_2d the last 4 jlines of pt2d 2366 2366 ij = jj - nlcj + ijpj 2367 znorthloc_2d(:,ij) = pt2d(:,jj)2367 xnorthloc_2d(:,ij) = pt2d(:,jj) 2368 2368 END DO 2369 2369 2370 ! ! Build in procs of ncomm_north the znorthgloio_2d2370 ! ! Build in procs of ncomm_north the xnorthgloio_2d 2371 2371 itaille = jpi * ijpj 2372 2372 IF ( l_north_nogather ) THEN … … 2378 2378 ij = jj - nlcj + ijpj 2379 2379 DO ji = 1, nlci 2380 ztab_2d(ji+nimpp-1,ij) = pt2d(ji,jj)2380 tab_2d(ji+nimpp-1,ij) = pt2d(ji,jj) 2381 2381 END DO 2382 2382 END DO … … 2404 2404 2405 2405 DO jr = 1,nsndto(ityp) 2406 CALL mppsend(5, znorthloc_2d, itaille, isendto(jr,ityp), ml_req_nf(jr) )2406 CALL mppsend(5, xnorthloc_2d, itaille, isendto(jr,ityp), ml_req_nf(jr) ) 2407 2407 END DO 2408 2408 DO jr = 1,nsndto(ityp) 2409 CALL mpprecv(5, zfoldwk_2d, itaille, isendto(jr,ityp))2409 CALL mpprecv(5, foldwk_2d, itaille, isendto(jr,ityp)) 2410 2410 iproc = isendto(jr,ityp) + 1 2411 2411 ildi = nldit (iproc) … … 2414 2414 DO jj = 1, ijpj 2415 2415 DO ji = ildi, ilei 2416 ztab_2d(ji+iilb-1,jj) = zfoldwk_2d(ji,jj)2416 tab_2d(ji+iilb-1,jj) = foldwk_2d(ji,jj) 2417 2417 END DO 2418 2418 END DO … … 2429 2429 2430 2430 IF ( ityp .lt. 0 ) THEN 2431 CALL MPI_ALLGATHER( znorthloc_2d , itaille, MPI_DOUBLE_PRECISION, &2432 & znorthgloio_2d, itaille, MPI_DOUBLE_PRECISION, ncomm_north, ierr )2431 CALL MPI_ALLGATHER( xnorthloc_2d , itaille, MPI_DOUBLE_PRECISION, & 2432 & xnorthgloio_2d, itaille, MPI_DOUBLE_PRECISION, ncomm_north, ierr ) 2433 2433 ! 2434 2434 DO jr = 1, ndim_rank_north ! recover the global north array … … 2439 2439 DO jj = 1, ijpj 2440 2440 DO ji = ildi, ilei 2441 ztab_2d(ji+iilb-1,jj) = znorthgloio_2d(ji,jj,jr)2441 tab_2d(ji+iilb-1,jj) = xnorthgloio_2d(ji,jj,jr) 2442 2442 END DO 2443 2443 END DO … … 2445 2445 ENDIF 2446 2446 ! 2447 ! The ztab array has been either:2447 ! The tab array has been either: 2448 2448 ! a. Fully populated by the mpi_allgather operation or 2449 2449 ! b. Had the active points for this domain and northern neighbours populated … … 2452 2452 ! this domain will be identical. 2453 2453 ! 2454 CALL lbc_nfd( ztab_2d, cd_type, psgn ) ! North fold boundary condition2454 CALL lbc_nfd( tab_2d, cd_type, psgn ) ! North fold boundary condition 2455 2455 ! 2456 2456 ! … … 2458 2458 ij = jj - nlcj + ijpj 2459 2459 DO ji = 1, nlci 2460 pt2d(ji,jj) = ztab_2d(ji+nimpp-1,ij)2460 pt2d(ji,jj) = tab_2d(ji+nimpp-1,ij) 2461 2461 END DO 2462 2462 END DO … … 2491 2491 ! 2492 2492 ijpj=4 2493 ztab_e(:,:) = 0.e02493 tab_e(:,:) = 0.e0 2494 2494 2495 2495 ij=0 2496 ! put in znorthloc_e the last 4 jlines of pt2d2496 ! put in xnorthloc_e the last 4 jlines of pt2d 2497 2497 DO jj = nlcj - ijpj + 1 - jpr2dj, nlcj +jpr2dj 2498 2498 ij = ij + 1 2499 2499 DO ji = 1, jpi 2500 znorthloc_e(ji,ij)=pt2d(ji,jj)2500 xnorthloc_e(ji,ij)=pt2d(ji,jj) 2501 2501 END DO 2502 2502 END DO 2503 2503 ! 2504 2504 itaille = jpi * ( ijpj + 2 * jpr2dj ) 2505 CALL MPI_ALLGATHER( znorthloc_e(1,1) , itaille, MPI_DOUBLE_PRECISION, &2506 & znorthgloio_e(1,1,1), itaille, MPI_DOUBLE_PRECISION, ncomm_north, ierr )2505 CALL MPI_ALLGATHER( xnorthloc_e(1,1) , itaille, MPI_DOUBLE_PRECISION, & 2506 & xnorthgloio_e(1,1,1), itaille, MPI_DOUBLE_PRECISION, ncomm_north, ierr ) 2507 2507 ! 2508 2508 DO jr = 1, ndim_rank_north ! recover the global north array … … 2513 2513 DO jj = 1, ijpj+2*jpr2dj 2514 2514 DO ji = ildi, ilei 2515 ztab_e(ji+iilb-1,jj) = znorthgloio_e(ji,jj,jr)2515 tab_e(ji+iilb-1,jj) = xnorthgloio_e(ji,jj,jr) 2516 2516 END DO 2517 2517 END DO … … 2521 2521 ! 2. North-Fold boundary conditions 2522 2522 ! ---------------------------------- 2523 CALL lbc_nfd( ztab_e(:,:), cd_type, psgn, pr2dj = jpr2dj )2523 CALL lbc_nfd( tab_e(:,:), cd_type, psgn, pr2dj = jpr2dj ) 2524 2524 2525 2525 ij = jpr2dj … … 2528 2528 ij = ij +1 2529 2529 DO ji= 1, nlci 2530 pt2d(ji,jj) = ztab_e(ji+nimpp-1,ij)2530 pt2d(ji,jj) = tab_e(ji+nimpp-1,ij) 2531 2531 END DO 2532 2532 END DO
Note: See TracChangeset
for help on using the changeset viewer.
