Changeset 1671


Ignore:
Timestamp:
Sep 25, 2015 3:29:37 AM (6 years ago)
Author:
raasch
Message:

bugfix: ghostpoint exchange for array diss in case that sgs velocities are used
for particles

File:
1 edited

Legend:

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

    r1586 r1671  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! bugfix: ghostpoint exchange for array diss in case that sgs velocities are used
     23! for particles
    2324!
    2425! Former revisions:
     
    214215
    215216    USE particle_attributes,                                                   &
    216         ONLY:  particle_advection, particle_advection_start, wang_kernel
     217        ONLY:  particle_advection, particle_advection_start,                   &
     218               use_sgs_for_particles, wang_kernel
    217219
    218220    USE pegrid
     
    446448                   CALL exchange_horiz( ql_vp, nbgp )
    447449                ENDIF
    448                 IF ( wang_kernel  .OR.  turbulence )  CALL exchange_horiz( diss, nbgp )
     450                IF ( wang_kernel  .OR.  turbulence  .OR.  use_sgs_for_particles ) &
     451                THEN
     452                   CALL exchange_horiz( diss, nbgp )
     453                ENDIF
    449454
    450455                IF ( numprocs == 1 )  THEN    ! workaround for single-core GPU runs
     
    510515                   CALL exchange_horiz( ql_vp, nbgp )
    511516                ENDIF
    512                 IF ( wang_kernel  .OR.  turbulence )  CALL exchange_horiz( diss, nbgp )
     517                IF ( wang_kernel  .OR.  turbulence  .OR.  use_sgs_for_particles ) &
     518                THEN
     519                   CALL exchange_horiz( diss, nbgp )
     520                ENDIF
    513521
    514522                IF ( numprocs == 1 )  THEN    ! workaround for single-core GPU runs
     
    597605                CALL exchange_horiz( ql_vp, nbgp )
    598606             ENDIF
    599              IF ( wang_kernel  .OR.  turbulence )  CALL exchange_horiz( diss, nbgp )
     607             IF ( wang_kernel  .OR.  turbulence  .OR.  use_sgs_for_particles ) &
     608             THEN
     609                CALL exchange_horiz( diss, nbgp )
     610             ENDIF
    600611
    601612             IF ( numprocs == 1 )  THEN    ! workaround for single-core GPU runs
Note: See TracChangeset for help on using the changeset viewer.