Ignore:
Timestamp:
Mar 27, 2015 9:56:27 AM (6 years ago)
Author:
raasch
Message:

optimized multigrid method installed, new parameter seed_follows_topography for particle release, small adjustment in subjob for HLRN

File:
1 edited

Legend:

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

    r1566 r1575  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! adjustments for psolver-queries, calculation of ngp_xz added
    2323!
    2424! Former revisions:
     
    761761!-- Calculate number of grid levels necessary for the multigrid poisson solver
    762762!-- as well as the gridpoint indices on each level
    763     IF ( psolver == 'multigrid' )  THEN
     763    IF ( psolver(1:9) == 'multigrid' )  THEN
    764764
    765765!
     
    977977!-- Do these calculations for the model grid and (if necessary) also
    978978!-- for the coarser grid levels used in the multigrid method
    979     ALLOCATE ( ngp_yz(0:maximum_grid_level), type_xz(0:maximum_grid_level),&
    980                type_yz(0:maximum_grid_level) )
     979    ALLOCATE ( ngp_xz(0:maximum_grid_level), ngp_yz(0:maximum_grid_level),     &
     980               type_xz(0:maximum_grid_level), type_yz(0:maximum_grid_level) )
    981981
    982982    nxl_l = nxl; nxr_l = nxr; nys_l = nys; nyn_l = nyn; nzb_l = nzb; nzt_l = nzt
     
    10071007!
    10081008!-- Definition of MPI-datatypes for multigrid method (coarser level grids)
    1009     IF ( psolver == 'multigrid' )  THEN
     1009    IF ( psolver(1:9) == 'multigrid' )  THEN
    10101010!   
    10111011!--    Definition of MPI-datatyoe as above, but only 1 ghost level is used
    10121012       DO  i = maximum_grid_level, 1 , -1
    10131013
     1014          ngp_xz(i) = (nzt_l - nzb_l + 2) * (nxr_l - nxl_l + 3)
    10141015          ngp_yz(i) = (nzt_l - nzb_l + 2) * (nyn_l - nys_l + 3)
    10151016
     
    11291130!
    11301131!-- Allocate wall flag arrays used in the multigrid solver
    1131     IF ( psolver == 'multigrid' )  THEN
     1132    IF ( psolver(1:9) == 'multigrid' )  THEN
    11321133
    11331134       DO  i = maximum_grid_level, 1, -1
Note: See TracChangeset for help on using the changeset viewer.