Changeset 1322 for palm/trunk/SOURCE/microphysics.f90

Timestamp:

Mar 20, 2014 4:38:49 PM (11 years ago)

Author:

raasch

Message:

REAL functions and a lot of REAL constants provided with KIND-attribute,
some small bugfixes

File:

• palm/trunk/SOURCE/microphysics.f90 (modified) (11 diffs)

palm/trunk/SOURCE/microphysics.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-!
+! REAL constants defined as wp-kind
 !
 ! Former revisions:
@@ -522 +522 @@
 !
 !--             Weight averaged radius of cloud droplets:
-                rc = 0.5 * ( xc * dpirho_l )**( 1.0 / 3.0 )
+                rc = 0.5 * ( xc * dpirho_l )**( 1.0 / 3.0_wp )
 
                 alpha_cc = ( a_1 + a_2 * nu_c ) / ( 1.0 + a_3 * nu_c )
@@ -529 +529 @@
 !
 !--             Mixing length (neglecting distance to ground and stratification)
-                l_mix = ( dx * dy * dzu(k) )**( 1.0 / 3.0 )
+                l_mix = ( dx * dy * dzu(k) )**( 1.0 / 3.0_wp )
 !
 !--             Limit dissipation rate according to Seifert, Nuijens and 
@@ -536 +536 @@
 !
 !--             Compute Taylor-microscale Reynolds number:
-                re_lambda = 6.0 / 11.0 * ( l_mix / c_const )**( 2.0 / 3.0 ) *  &
-                            SQRT( 15.0 / kin_vis_air ) * epsilon**( 1.0 / 6.0 )
+                re_lambda = 6.0 / 11.0 * ( l_mix / c_const )**( 2.0 / 3.0_wp ) *  &
+                            SQRT( 15.0 / kin_vis_air ) * epsilon**( 1.0 / 6.0_wp )
 !
 !--             The factor of 1.0E4 is needed to convert the dissipation rate
@@ -664 +664 @@
 !
 !--          Weight averaged diameter of rain drops:
-             dr = ( hyrho(k) * qr_1d(k) / nr_1d(k) * dpirho_l )**( 1.0 / 3.0 )
+             dr = ( hyrho(k) * qr_1d(k) / nr_1d(k) * dpirho_l )**( 1.0 / 3.0_wp )
 !
 !--          Collisional breakup rate (Seifert, 2008):
@@ -760 +760 @@
 !
 !--          Weight averaged diameter of rain drops:
-             dr = ( xr * dpirho_l )**( 1.0 / 3.0 )
+             dr = ( xr * dpirho_l )**( 1.0 / 3.0_wp )
 !
 !--          Compute ventilation factor and intercept parameter 
@@ -772 +772 @@
 !--             Slope parameter of gamma distribution (Seifert, 2008):
                 lambda_r = ( ( mu_r + 3.0 ) * ( mu_r + 2.0 ) *       &
-                             ( mu_r + 1.0 ) )**( 1.0 / 3.0 ) / dr
+                             ( mu_r + 1.0 ) )**( 1.0 / 3.0_wp ) / dr
 
                 mu_r_2   = mu_r + 2.0
@@ -849 +849 @@
 !
 !--    Sedimentation of cloud droplets (Heus et al., 2010):
-       sed_qc_const = k_st * ( 3.0 / ( 4.0 * pi * rho_l ))**( 2.0 / 3.0 ) *   &
+       sed_qc_const = k_st * ( 3.0 / ( 4.0 * pi * rho_l ))**( 2.0 / 3.0_wp ) *   &
                      EXP( 5.0 * LOG( sigma_gc )**2 )
 
@@ -856 +856 @@
        DO  k = nzt, nzb_s_inner(j,i)+1, -1
           IF ( qc_1d(k) > eps_sb )  THEN
-             sed_qc(k) = sed_qc_const * nc_1d(k)**( -2.0 / 3.0 ) * &
-                        ( qc_1d(k) * hyrho(k) )**( 5.0 / 3.0 )
+             sed_qc(k) = sed_qc_const * nc_1d(k)**( -2.0 / 3.0_wp ) * &
+                        ( qc_1d(k) * hyrho(k) )**( 5.0 / 3.0_wp )
           ELSE
              sed_qc(k) = 0.0
@@ -940 +940 @@
 !
 !--          Weight averaged diameter of rain drops:
-             dr = ( hyrho(k) * qr_1d(k) / nr_1d(k) * dpirho_l )**( 1.0 / 3.0 )
+             dr = ( hyrho(k) * qr_1d(k) / nr_1d(k) * dpirho_l )**( 1.0 / 3.0_wp )
 !
 !--          Shape parameter of gamma distribution (Milbrandt and Yau, 2005;
@@ -948 +948 @@
 !--          Slope parameter of gamma distribution (Seifert, 2008):
              lambda_r = ( ( mu_r + 3.0 ) * ( mu_r + 2.0 ) *          &
-                        ( mu_r + 1.0 ) )**( 1.0 / 3.0 ) / dr
+                        ( mu_r + 1.0 ) )**( 1.0 / 3.0_wp ) / dr
 
              w_nr(k) = MAX( 0.1, MIN( 20.0, a_term - b_term * ( 1.0 +          &