Changeset 1353 for palm/trunk/SOURCE/calc_liquid_water_content.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/calc_liquid_water_content.f90 (modified) (2 diffs)

palm/trunk/SOURCE/calc_liquid_water_content.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-!
+! REAL constants provided with KIND-attribute
 !
 ! Former revisions:
@@ -97 +97 @@
 !
 !--          Compute saturation water vapor pressure at t_l
-             e_s = 610.78 * EXP( 17.269 * ( t_l - 273.16 ) /                   &
-                                          ( t_l - 35.86 ) )
+             e_s = 610.78_wp * EXP( 17.269_wp * ( t_l - 273.16_wp ) /          &
+                                                ( t_l - 35.86_wp ) )
 
 !
 !--          Compute approximation of saturation humidity
-             q_s = 0.622 * e_s / ( hyp(k) - 0.378 * e_s )
+             q_s = 0.622_wp * e_s / ( hyp(k) - 0.378_wp * e_s )
 
 !
 !--          Correction factor
-             alpha = 0.622 * l_d_r * l_d_cp / ( t_l * t_l )
+             alpha = 0.622_wp * l_d_r * l_d_cp / ( t_l * t_l )
 
 !
 !--          Correction of the approximated value
 !--          (see: Cuijpers + Duynkerke, 1993, JAS, 23)
-             q_s = q_s * ( 1.0 + alpha * q(k,j,i) ) / ( 1.0 + alpha * q_s )
+             q_s = q_s * ( 1.0_wp + alpha * q(k,j,i) ) / ( 1.0_wp + alpha * q_s )
 
 !
 !--          Compute the liquid water content
              IF ( icloud_scheme == 0  .AND.  precipitation)  THEN
-                IF ( ( q(k,j,i) - q_s - qr(k,j,i) ) > 0.0 ) THEN
+                IF ( ( q(k,j,i) - q_s - qr(k,j,i) ) > 0.0_wp ) THEN
                    qc(k,j,i) = q(k,j,i) - q_s - qr(k,j,i)
                    ql(k,j,i) = qc(k,j,i) + qr(k,j,i)
                 ELSE
                    IF ( q(k,j,i) < qr(k,j,i) )  q(k,j,i) = qr(k,j,i)
-                   qc(k,j,i) = 0.0 
+                   qc(k,j,i) = 0.0_wp 
                    ql(k,j,i) = qr(k,j,i)
                 ENDIF
              ELSEIF ( icloud_scheme == 0  .AND.  .NOT. precipitation )  THEN
-                IF ( ( q(k,j,i) - q_s ) > 0.0 ) THEN
+                IF ( ( q(k,j,i) - q_s ) > 0.0_wp ) THEN
                    qc(k,j,i) = q(k,j,i) - q_s
                    ql(k,j,i) = qc(k,j,i)
                 ELSE
-                   qc(k,j,i) = 0.0 
-                   ql(k,j,i) = 0.0
+                   qc(k,j,i) = 0.0_wp 
+                   ql(k,j,i) = 0.0_wp
                 ENDIF
              ELSE
-                IF ( ( q(k,j,i) - q_s ) > 0.0 ) THEN
+                IF ( ( q(k,j,i) - q_s ) > 0.0_wp ) THEN
                    ql(k,j,i) = q(k,j,i) - q_s
                 ELSE
-                   ql(k,j,i) = 0.0 
+                   ql(k,j,i) = 0.0_wp 
                 ENDIF
              ENDIF