Changeset 3039


Ignore:
Timestamp:
May 24, 2018 1:13:11 PM (3 years ago)
Author:
schwenkel
Message:

bugfix for lcm with grid stretching

Location:
palm/trunk/SOURCE
Files:
4 edited

Legend:

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

    r2718 r3039  
    2525! -----------------
    2626! $Id$
     27! bugfix for lcm with grid stretching
     28!
     29! 2718 2018-01-02 08:49:38Z maronga
    2730! Corrected "Former revisions" section
    2831!
     
    6770
    6871    USE arrays_3d,                                                             &
    69         ONLY:  ql, ql_v, ql_vp
     72        ONLY:  dzw, ql, ql_v, ql_vp
    7073
    7174    USE cloud_parameters,                                                      &
     
    127130             IF ( ql_v(k,j,i) /= 0.0_wp )  THEN
    128131                ql(k,j,i) = ql(k,j,i) + rho_l * 1.33333333_wp * pi *           &
    129                                         ql_v(k,j,i) /                       &
    130                                         ( rho_surface * dx * dy * dz )
     132                                        ql_v(k,j,i) /                          &
     133                                        ( rho_surface * dx * dy * dzw(k) )
    131134
    132135                IF ( ql(k,j,i) < 0.0_wp )  THEN
  • palm/trunk/SOURCE/lpm_droplet_collision.f90

    r2718 r3039  
    2525! -----------------
    2626! $Id$
     27! bugfix for lcm with grid stretching
     28!
     29! 2718 2018-01-02 08:49:38Z maronga
    2730! Corrected "Former revisions" section
    2831!
     
    114117
    115118    USE arrays_3d,                                                             &
    116         ONLY:  diss, ql_v, ql_vp
     119        ONLY:  diss, dzw, ql_v, ql_vp
    117120
    118121    USE cloud_parameters,                                                      &
     
    174177    number_of_particles   = prt_count(k,j,i)
    175178    factor_volume_to_mass = 4.0_wp / 3.0_wp * pi * rho_l
    176     ddV                   = 1.0_wp / ( dx * dy * dz )
     179    ddV                   = 1.0_wp / ( dx * dy * dzw(k) )
    177180!
    178181!-- Collision requires at least one super droplet inside the box
  • palm/trunk/SOURCE/lpm_droplet_condensation.f90

    r2718 r3039  
    2525! -----------------
    2626! $Id$
     27! bugfix for lcm with grid stretching
     28!
     29! 2718 2018-01-02 08:49:38Z maronga
    2730! Corrected "Former revisions" section
    2831!
     
    124127
    125128    USE arrays_3d,                                                             &
    126         ONLY:  hyp, pt, q, ql_c, ql_v
     129        ONLY:  dzw, hyp, pt, q, ql_c, ql_v
    127130
    128131    USE cloud_parameters,                                                      &
     
    372375                                   rho_l * 1.33333333_wp * pi *                &
    373376                                   ( new_r(n)**3 - particles(n)%radius**3 ) /  &
    374                                    ( rho_surface * dx * dy * dz )
     377                                   ( rho_surface * dx * dy * dzw(kp) )
    375378!
    376379!--    Check if the increase in liqid water is not too big. If this is the case,
  • palm/trunk/SOURCE/lpm_init.f90

    r2967 r3039  
    2525! -----------------
    2626! $Id$
     27! bugfix for lcm with grid stretching
     28!
     29! 2967 2018-04-13 11:22:08Z raasch
    2730! nesting routine is only called if nesting is switched on
    2831!
     
    209212
    210213    USE arrays_3d,                                                             &
    211         ONLY:  de_dx, de_dy, de_dz, zu, zw
     214        ONLY:  de_dx, de_dy, de_dz, dzw, zu, zw
    212215
    213216    USE control_parameters,                                                    &
     
    11611164                     na(3) / ( SQRT( 2.0 * pi ) * log_sigma(3) ) *                     &
    11621165                     EXP( - LOG10( r_mid / rm(3) )**2 / ( 2.0 * log_sigma(3)**2 ) )    &
    1163                    ) * ( LOG10(r_r) - LOG10(r_l) ) * ( dx * dy * dz )
     1166                   ) * ( LOG10(r_r) - LOG10(r_l) ) * ( dx * dy * dzw(kp) )
    11641167
    11651168!
Note: See TracChangeset for help on using the changeset viewer.