Changeset 3207 for palm


Ignore:
Timestamp:
Aug 27, 2018 12:55:33 PM (6 years ago)
Author:
schwenkel
Message:

Minor bugfix for sgs-velocities in case of cloud droplets

File:
1 edited

Legend:

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

    r3189 r3207  
    2525! -----------------
    2626! $Id$
     27! Minor bugfix for sgs-velocities in case of cloud droplets
     28!
     29! 3189 2018-08-06 13:18:55Z raasch
    2730! Bugfix: Index of the array dzw has to be k+1 during the interpolation.
    2831! Otherwise k=0 causes an abortion because dzw is allocated from 1 to nzt+1
     
    977980                   IF ( use_sgs_for_particles )  THEN
    978981                      lagr_timescale(n) = km(kp,jp,ip) / MAX( e(kp,jp,ip), 1.0E-20_wp )
    979                       RL             = EXP( -1.0_wp * dt_3d / lagr_timescale(n) )
     982                      RL             = EXP( -1.0_wp * dt_3d / MAX( lagr_timescale(n), &
     983                                             1.0E-20_wp ) )
    980984                      sigma          = SQRT( e(kp,jp,ip) )
    981985
     
    10491053                IF ( use_sgs_for_particles )  THEN
    10501054                    lagr_timescale(n) = km(kp,jp,ip) / MAX( e(kp,jp,ip), 1.0E-20_wp )
    1051                     RL             = EXP( -1.0_wp * dt_3d / lagr_timescale(n) )
     1055                     RL             = EXP( -1.0_wp * dt_3d / MAX( lagr_timescale(n), &
     1056                                             1.0E-20_wp ) )
    10521057                    sigma          = SQRT( e(kp,jp,ip) )
    10531058
Note: See TracChangeset for help on using the changeset viewer.