- Timestamp:
- Apr 12, 2012 3:38:37 PM (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/poismg.f90
r779 r879 7 7 ! Current revisions: 8 8 ! ----------------- 9 ! 9 ! Bugfix: preprocessor statements for parallel execution added 10 10 ! 11 11 ! Former revisions: … … 88 88 !-- Initialize arrays and variables used in this subroutine 89 89 90 !-- If the number of grid points of the gathered grid, which is collected 91 !-- on PE0, is larger than the number of grid points of an PE, than array 92 !-- p3 will be enlarged. 90 !-- If the number of grid points of the gathered grid, which is collected 91 !-- on PE0, is larger than the number of grid points of an PE, than array 92 !-- p3 will be enlarged. 93 #if defined( __parallel ) 93 94 IF ( gathered_size > subdomain_size ) THEN 94 ALLOCATE( p3(nzb:nzt_mg(mg_switch_to_pe0_level)+1,nys_mg( &95 mg_switch_to_pe0_level)-1:nyn_mg(mg_switch_to_pe0_level)+1, 96 nxl_mg(mg_switch_to_pe0_level)-1:nxr_mg( 95 ALLOCATE( p3(nzb:nzt_mg(mg_switch_to_pe0_level)+1,nys_mg( & 96 mg_switch_to_pe0_level)-1:nyn_mg(mg_switch_to_pe0_level)+1,& 97 nxl_mg(mg_switch_to_pe0_level)-1:nxr_mg( & 97 98 mg_switch_to_pe0_level)+1) ) 98 99 ELSE 99 100 ALLOCATE ( p3(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) ) 100 101 ENDIF 102 #else 103 ALLOCATE ( p3(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) ) 104 #endif 105 101 106 ! 102 107 !-- Ghost boundaries have to be added to divergence array. … … 1272 1277 !-- from the coarse to the next finer grid. 1273 1278 IF ( grid_level == mg_switch_to_pe0_level+1 ) THEN 1279 1280 #if defined( __parallel ) 1274 1281 ! 1275 1282 !-- At this level, the new residual first has to be scattered from … … 1352 1359 1353 1360 DEALLOCATE( p2_sub ) 1361 #endif 1354 1362 1355 1363 ELSE 1364 1356 1365 CALL prolong( p2, p3 ) 1357 1366
Note: See TracChangeset
for help on using the changeset viewer.