Changeset 1053 for palm/trunk/SOURCE/init_cloud_physics.f90

Timestamp:

Nov 13, 2012 5:11:03 PM (11 years ago)

Author:

hoffmann

Message:

two-moment cloud physics implemented

Location:

palm/trunk

Files:

• . (modified) (1 prop)
• SOURCE (modified) (1 prop)
• SOURCE/init_cloud_physics.f90 (modified) (4 diffs)

palm/trunk/SOURCE/init_cloud_physics.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! calculation of the maximum timestep according to the terminal velocity of rain 
+! drops in the two moment cloud scheme
+!
+! calculation of frequently used constants (pirho_l, dpirho_l, schmidt_p_1d3, 
+! hyrho) 
 !
 ! Former revisions:
@@ -66 +70 @@
     REAL    ::  t_surface
 
-    ALLOCATE( hyp(nzb:nzt+1), pt_d_t(nzb:nzt+1), t_d_pt(nzb:nzt+1) )
+    ALLOCATE( hyp(nzb:nzt+1), pt_d_t(nzb:nzt+1), t_d_pt(nzb:nzt+1),  &
+              hyrho(nzb:nzt+1) )
 
 !
@@ -73 +78 @@
     l_d_r  = l_v / r_d
     l_d_rv = l_v / r_v
+
+    schmidt_p_1d3 = schmidt**( 1.0 / 3.0 )
+
+    pirho_l  = pi * rho_l / 6.0
+    dpirho_l = 1.0 / pirho_l 
+!
+!-- Calculate timestep according to precipitation
+    IF ( icloud_scheme == 0  .AND.  precipitation )  THEN
+       dt_precipitation = MINVAL( dzu(nzb+2:nzt) ) / w_precipitation
+    ENDIF
 !
 !-- Calculate factor used in equation for droplet growth by condensation
@@ -94 +109 @@
                    ( (t_surface - g/cp * zu(k)) / t_surface )**(1.0/0.286)
        pt_d_t(k) = ( 100000.0 / hyp(k) )**0.286
-       t_d_pt(k) = 1.0 / pt_d_t(k)       
+       t_d_pt(k) = 1.0 / pt_d_t(k)
+       hyrho(k)  = hyp(k) / ( r_d * t_d_pt(k) * pt_init(k) )       
     ENDDO