Changeset 879 for palm


Ignore:
Timestamp:
Apr 12, 2012 3:38:37 PM (13 years ago)
Author:
raasch
Message:

Bugfix: preprocessor statements for parallel execution added

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/poismg.f90

    r779 r879  
    77! Current revisions:
    88! -----------------
    9 !
     9! Bugfix: preprocessor statements for parallel execution added
    1010!
    1111! Former revisions:
     
    8888!-- Initialize arrays and variables used in this subroutine
    8989
    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 )
    9394    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(                   &
    9798                    mg_switch_to_pe0_level)+1) )
    9899    ELSE
    99100       ALLOCATE ( p3(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
    100101    ENDIF
     102#else
     103   ALLOCATE ( p3(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
     104#endif
     105 
    101106!
    102107!-- Ghost boundaries have to be added to divergence array.
     
    12721277!--    from the coarse to the next finer grid.
    12731278       IF ( grid_level == mg_switch_to_pe0_level+1 )  THEN
     1279
     1280#if defined( __parallel )
    12741281!
    12751282!--       At this level, the new residual first has to be scattered from
     
    13521359
    13531360          DEALLOCATE( p2_sub )
     1361#endif
    13541362
    13551363       ELSE
     1364
    13561365          CALL prolong( p2, p3 )
    13571366
Note: See TracChangeset for help on using the changeset viewer.