Changeset 3039

Timestamp:

May 24, 2018 1:13:11 PM (7 years ago)

Author:

schwenkel

Message:

bugfix for lcm with grid stretching

Location:

palm/trunk/SOURCE

Files:

• lpm_calc_liquid_water_content.f90 (modified) (3 diffs)
• lpm_droplet_collision.f90 (modified) (3 diffs)
• lpm_droplet_condensation.f90 (modified) (3 diffs)
• lpm_init.f90 (modified) (3 diffs)

palm/trunk/SOURCE/lpm_calc_liquid_water_content.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! bugfix for lcm with grid stretching
+! 
+! 2718 2018-01-02 08:49:38Z maronga
 ! Corrected "Former revisions" section
 ! 
@@ -67 +70 @@
 
     USE arrays_3d,                                                             &
-        ONLY:  ql, ql_v, ql_vp
+        ONLY:  dzw, ql, ql_v, ql_vp
 
     USE cloud_parameters,                                                      &
@@ -127 +130 @@
              IF ( ql_v(k,j,i) /= 0.0_wp )  THEN
                 ql(k,j,i) = ql(k,j,i) + rho_l * 1.33333333_wp * pi *           &
-                                        ql_v(k,j,i) /                       &
-                                        ( rho_surface * dx * dy * dz )
+                                        ql_v(k,j,i) /                          &
+                                        ( rho_surface * dx * dy * dzw(k) )
 
                 IF ( ql(k,j,i) < 0.0_wp )  THEN

palm/trunk/SOURCE/lpm_droplet_collision.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! bugfix for lcm with grid stretching
+! 
+! 2718 2018-01-02 08:49:38Z maronga
 ! Corrected "Former revisions" section
 ! 
@@ -114 +117 @@
 
     USE arrays_3d,                                                             &
-        ONLY:  diss, ql_v, ql_vp
+        ONLY:  diss, dzw, ql_v, ql_vp
 
     USE cloud_parameters,                                                      &
@@ -174 +177 @@
     number_of_particles   = prt_count(k,j,i)
     factor_volume_to_mass = 4.0_wp / 3.0_wp * pi * rho_l
-    ddV                   = 1.0_wp / ( dx * dy * dz )
+    ddV                   = 1.0_wp / ( dx * dy * dzw(k) )
 !
 !-- Collision requires at least one super droplet inside the box

palm/trunk/SOURCE/lpm_droplet_condensation.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! bugfix for lcm with grid stretching
+! 
+! 2718 2018-01-02 08:49:38Z maronga
 ! Corrected "Former revisions" section
 ! 
@@ -124 +127 @@
 
     USE arrays_3d,                                                             &
-        ONLY:  hyp, pt, q, ql_c, ql_v
+        ONLY:  dzw, hyp, pt, q, ql_c, ql_v
 
     USE cloud_parameters,                                                      &
@@ -372 +375 @@
                                    rho_l * 1.33333333_wp * pi *                &
                                    ( new_r(n)**3 - particles(n)%radius**3 ) /  &
-                                   ( rho_surface * dx * dy * dz )
+                                   ( rho_surface * dx * dy * dzw(kp) )
 !
 !--    Check if the increase in liqid water is not too big. If this is the case,

palm/trunk/SOURCE/lpm_init.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! bugfix for lcm with grid stretching
+! 
+! 2967 2018-04-13 11:22:08Z raasch
 ! nesting routine is only called if nesting is switched on
 ! 
@@ -209 +212 @@
 
     USE arrays_3d,                                                             &
-        ONLY:  de_dx, de_dy, de_dz, zu, zw
+        ONLY:  de_dx, de_dy, de_dz, dzw, zu, zw
 
     USE control_parameters,                                                    &
@@ -1161 +1164 @@
                      na(3) / ( SQRT( 2.0 * pi ) * log_sigma(3) ) *                     &
                      EXP( - LOG10( r_mid / rm(3) )**2 / ( 2.0 * log_sigma(3)**2 ) )    &
-                   ) * ( LOG10(r_r) - LOG10(r_l) ) * ( dx * dy * dz )
+                   ) * ( LOG10(r_r) - LOG10(r_l) ) * ( dx * dy * dzw(kp) )
 
 !