Changeset 1003 for palm/trunk/SOURCE/pres.f90
- Timestamp:
- Sep 14, 2012 2:35:53 PM (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/pres.f90
r779 r1003 4 4 ! Current revisions: 5 5 ! ----------------- 6 ! 6 ! adjustment of array tend for cases with unequal subdomain sizes removed 7 7 ! 8 8 ! Former revisions: … … 279 279 ELSE 280 280 !$OMP PARALLEL DO SCHEDULE( STATIC ) 281 DO i = nxl, nxr a282 DO j = nys, nyn a283 DO k = nzb+1, nzt a281 DO i = nxl, nxr 282 DO j = nys, nyn 283 DO k = nzb+1, nzt 284 284 d(k,j,i) = 0.0 285 285 ENDDO … … 409 409 410 410 ! 411 !-- Enlarge the size of tend, used as a working array for the transpositions412 IF ( nxra > nxr .OR. nyna > nyn .OR. nza > nz ) THEN413 DEALLOCATE( tend )414 ALLOCATE( tend(1:nza,nys:nyna,nxl:nxra) )415 ENDIF416 417 !418 411 !-- Solve Poisson equation via FFT and solution of tridiagonal matrices 419 412 IF ( psolver == 'poisfft' ) THEN … … 427 420 !-- are some optimization problems in poisfft 428 421 CALL poisfft_hybrid( d ) 429 ENDIF430 431 !432 !-- Resize tend to its normal size433 IF ( nxra > nxr .OR. nyna > nyn .OR. nza > nz ) THEN434 DEALLOCATE( tend )435 ALLOCATE( tend(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )436 422 ENDIF 437 423
Note: See TracChangeset
for help on using the changeset viewer.