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

Timestamp:

Apr 8, 2014 3:21:23 PM (10 years ago)

Author:

heinze

Message:

REAL constants provided with KIND-attribute

File:

• palm/trunk/SOURCE/init_cloud_physics.f90 (modified) (5 diffs)

palm/trunk/SOURCE/init_cloud_physics.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! REAL constants provided with KIND-attribute
 ! 
 ! Former revisions:
@@ -93 +93 @@
     REAL(wp) ::  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) )
 
@@ -104 +104 @@
     schmidt_p_1d3 = schmidt**( 1.0_wp / 3.0_wp )
 
-    pirho_l  = pi * rho_l / 6.0
-    dpirho_l = 1.0 / pirho_l 
+    pirho_l  = pi * rho_l / 6.0_wp
+    dpirho_l = 1.0_wp / pirho_l 
 !
 !-- Calculate timestep according to precipitation
     IF ( icloud_scheme == 0  .AND.  precipitation )  THEN
-       dt_precipitation = c_sedimentation * MINVAL( dzu(nzb+2:nzt) ) /        &
+       dt_precipitation = c_sedimentation * MINVAL( dzu(nzb+2:nzt) ) /         &
                           w_precipitation
     ENDIF
 !
 !-- Calculate factor used in equation for droplet growth by condensation
-    bfactor = 3.0 * vanthoff * mass_of_solute * molecular_weight_of_water &
-              / ( 4.0 * pi * rho_l * molecular_weight_of_solute )
+    bfactor = 3.0_wp * vanthoff * mass_of_solute * molecular_weight_of_water   &
+              / ( 4.0_wp * pi * rho_l * molecular_weight_of_solute )
 !
 !-- Calculate:
@@ -121 +121 @@
 !-- t / pt : ratio of actual and potential temperature (t_d_pt)
 !-- p_0(z) : vertical profile of the hydrostatic pressure (hyp)
-    t_surface = pt_surface * ( surface_pressure / 1000.0 )**0.286_wp
+    t_surface = pt_surface * ( surface_pressure / 1000.0_wp )**0.286_wp
     DO  k = nzb, nzt+1
 !
 !--    Check temperature in case of too large domain height
-       IF ( ( t_surface - g/cp * zu(k) ) < 0.0 )  THEN
+       IF ( ( t_surface - g/cp * zu(k) ) < 0.0_wp )  THEN
           WRITE( message_string, * )  'absolute temperature < 0.0 at zu(', k, &
                                       ') = ', zu(k)
           CALL message( 'init_cloud_physics', 'PA0142', 1, 2, 0, 6, 0 )
        ENDIF
-       hyp(k)    = surface_pressure * 100.0 * &
+       hyp(k)    = surface_pressure * 100.0_wp * &
                    ( (t_surface - g/cp * zu(k)) / t_surface )**(1.0_wp/0.286_wp)
-       pt_d_t(k) = ( 100000.0 / hyp(k) )**0.286_wp
-       t_d_pt(k) = 1.0 / pt_d_t(k)
+       pt_d_t(k) = ( 100000.0_wp / hyp(k) )**0.286_wp
+       t_d_pt(k) = 1.0_wp / pt_d_t(k)
        hyrho(k)  = hyp(k) / ( r_d * t_d_pt(k) * pt_init(k) )       
     ENDDO
@@ -139 +139 @@
 !
 !-- Compute reference density
-    rho_surface = surface_pressure * 100.0 / ( r_d * t_surface )
+    rho_surface = surface_pressure * 100.0_wp / ( r_d * t_surface )