Changeset 2298 for palm/trunk/SOURCE/pres.f90
- Timestamp:
- Jun 29, 2017 9:28:18 AM (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/pres.f90
r2233 r2298 25 25 ! ----------------- 26 26 ! $Id$ 27 ! comment changed + some formatting 28 ! 29 ! 2233 2017-05-30 18:08:54Z suehring 27 30 ! 28 31 ! 2232 2017-05-30 17:47:52Z suehring … … 537 540 !$OMP PARALLEL DO PRIVATE( i, j, k ) 538 541 DO m = 1, bc_h(0)%ns 539 i = bc_h(0)%i(m)540 j 541 k 542 i = bc_h(0)%i(m) 543 j = bc_h(0)%j(m) 544 k = bc_h(0)%k(m) 542 545 tend(k-1,j,i) = tend(k,j,i) 543 546 ENDDO … … 546 549 !$OMP PARALLEL DO PRIVATE( i, j, k ) 547 550 DO m = 1, bc_h(1)%ns 548 i = bc_h(1)%i(m)549 j 550 k 551 i = bc_h(1)%i(m) 552 j = bc_h(1)%j(m) 553 k = bc_h(1)%k(m) 551 554 tend(k+1,j,i) = tend(k,j,i) 552 555 ENDDO … … 559 562 !$OMP PARALLEL DO PRIVATE( i, j, k ) 560 563 DO m = 1, bc_h(0)%ns 561 i = bc_h(0)%i(m)562 j 563 k 564 i = bc_h(0)%i(m) 565 j = bc_h(0)%j(m) 566 k = bc_h(0)%k(m) 564 567 tend(k-1,j,i) = 0.0_wp 565 568 ENDDO … … 568 571 !$OMP PARALLEL DO PRIVATE( i, j, k ) 569 572 DO m = 1, bc_h(1)%ns 570 i = bc_h(1)%i(m)571 j 572 k 573 i = bc_h(1)%i(m) 574 j = bc_h(1)%j(m) 575 k = bc_h(1)%k(m) 573 576 tend(k+1,j,i) = 0.0_wp 574 577 ENDDO … … 616 619 ! 617 620 !-- Solve Poisson equation for perturbation pressure using Multigrid scheme, 618 !-- array tend is used to store the residuals, logical exchange_mg is used 619 !-- to discern data exchange in multigrid ( 1 ghostpoint ) and normal grid 620 !-- ( nbgp ghost points ). 621 !-- array tend is used to store the residuals. 621 622 622 623 !-- If the number of grid points of the gathered grid, which is collected
Note: See TracChangeset
for help on using the changeset viewer.