Ignore:
Timestamp:
Jun 9, 2016 4:25:25 PM (5 years ago)
Author:
suehring
Message:

several bugfixes in particle model and serial mode

File:
1 edited

Legend:

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

    r1823 r1929  
    1919! Current revisions:
    2020! ------------------
    21 !
     21! Call wall boundary conditions only if particles are in the vertical range of
     22! topography.
    2223!
    2324! Former revisions:
     
    111112
    112113    USE indices,                                                               &
    113         ONLY: nxl, nxr, nys, nyn, nzb, nzt
     114        ONLY: nxl, nxr, nys, nyn, nzb, nzb_max, nzt, nzb_w_inner
    114115
    115116    USE kinds
     
    300301                CALL lpm_advec(i,j,k)
    301302!
    302 !--             Particle reflection from walls
    303                 IF ( topography /= 'flat' )  THEN
     303!--             Particle reflection from walls. Only applied if the particles
     304!--             are in the vertical range of the topography. (Here, some
     305!--             optimization is still possible.)
     306                IF ( topography /= 'flat' .AND. k < nzb_max + 2 )  THEN
    304307                   CALL lpm_boundary_conds( 'walls' )
    305308                ENDIF
     
    334337       steps = steps + 1
    335338       dt_3d_reached_l = ALL(grid_particles(:,:,:)%time_loop_done)
    336 
    337339!
    338340!--    Find out, if all particles on every PE have completed the LES timestep
     
    359361!--    Horizontal boundary conditions including exchange between subdmains
    360362       CALL lpm_exchange_horiz
    361 
    362363!
    363364!--    Pack particles (eliminate those marked for deletion),
Note: See TracChangeset for help on using the changeset viewer.