Ignore:
Timestamp:
Nov 7, 2011 2:18:25 PM (10 years ago)
Author:
fricke
Message:

Modifications of the multigrid pressure solver

File:
1 edited

Legend:

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

    r720 r778  
    44! Current revisions:
    55! -----------------
    6 !
     6! New allocation of tend when multigrid is used and the collected field on PE0
     7! has more grid points than the subdomain of an PE.
    78!
    89! Former revisions:
     
    526527!--    to discern data exchange in multigrid ( 1 ghostpoint ) and normal grid
    527528!--    ( nbgp ghost points ).
     529
     530!--    If the number of grid points of the gathered grid, which is collected
     531!--    on PE0, is larger than the number of grid points of an PE, than array
     532!--    tend will be enlarged.
     533       IF ( gathered_size > subdomain_size )  THEN
     534          DEALLOCATE( tend )
     535          ALLOCATE( tend(nzb:nzt_mg(mg_switch_to_pe0_level)+1,nys_mg(          &
     536                    mg_switch_to_pe0_level)-1:nyn_mg(mg_switch_to_pe0_level)+1,&
     537                    nxl_mg(mg_switch_to_pe0_level)-1:nxr_mg(                   &
     538                    mg_switch_to_pe0_level)+1) )
     539       ENDIF
     540
    528541       CALL poismg( tend )
     542
     543       IF ( gathered_size > subdomain_size )  THEN
     544          DEALLOCATE( tend )
     545          ALLOCATE( tend(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     546       ENDIF
    529547
    530548!
Note: See TracChangeset for help on using the changeset viewer.