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

Timestamp:

Mar 27, 2014 1:18:20 PM (10 years ago)

Author:

heinze

Message:

Bugfix: REAL constants provided with KIND-attribute especially in call of intrinsic function like MAX, MIN, SIGN

File:

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

palm/trunk/SOURCE/microphysics.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
+! Bugfix: REAL constants provided with KIND-attribute especially in call of 
+! intrinsic function like MAX, MIN, SIGN
 ! 
 ! Former revisions:
@@ -988 +989 @@
              d_max  = MAX( qr_1d(k+1), qr_1d(k), qr_1d(k-1) ) - qr_1d(k)
 
-             qr_slope(k) = SIGN(1.0, d_mean) * MIN ( 2.0 * d_min, 2.0 * d_max, &
+             qr_slope(k) = SIGN(1.0_wp, d_mean) * MIN ( 2.0 * d_min, 2.0 * d_max, &
                                                      ABS( d_mean ) )
 
@@ -995 +996 @@
              d_max  = MAX( nr_1d(k+1), nr_1d(k), nr_1d(k-1) ) - nr_1d(k)
 
-             nr_slope(k) = SIGN(1.0, d_mean) * MIN ( 2.0 * d_min, 2.0 * d_max, &
+             nr_slope(k) = SIGN(1.0_wp, d_mean) * MIN ( 2.0 * d_min, 2.0 * d_max, &
                                                      ABS( d_mean ) )
           ENDDO
@@ -1017 +1018 @@
           z_run = 0.0 ! height above z(k)
           k_run = k
-          c_run = MIN( 1.0, c_nr(k) )
+          c_run = MIN( 1.0_wp, c_nr(k) )
           DO WHILE ( c_run > 0.0  .AND.  k_run <= nzt )
              flux  = flux + hyrho(k_run) *                                    &
@@ -1024 +1025 @@
              z_run = z_run + dzu(k_run)
              k_run = k_run + 1
-             c_run = MIN( 1.0, c_nr(k_run) - z_run * ddzu(k_run) )
+             c_run = MIN( 1.0_wp, c_nr(k_run) - z_run * ddzu(k_run) )
           ENDDO
 !
@@ -1041 +1042 @@
           z_run = 0.0 ! height above z(k)
           k_run = k
-          c_run = MIN( 1.0, c_qr(k) )
+          c_run = MIN( 1.0_wp, c_qr(k) )
 
           DO WHILE ( c_run > 0.0  .AND.  k_run <= nzt-1 )
@@ -1050 +1051 @@
              z_run = z_run + dzu(k_run)
              k_run = k_run + 1
-             c_run = MIN( 1.0, c_qr(k_run) - z_run * ddzu(k_run) )
+             c_run = MIN( 1.0_wp, c_qr(k_run) - z_run * ddzu(k_run) )
 
           ENDDO