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

palm/trunk/SOURCE/subsidence.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! REAL constants provided with KIND-attribute
 ! 
 ! Former revisions:
@@ -89 +89 @@
        IF ( .NOT. ALLOCATED( w_subs )) THEN
           ALLOCATE( w_subs(nzb:nzt+1) )
-          w_subs = 0.0
+          w_subs = 0.0_wp
        ENDIF 
 
@@ -102 +102 @@
 !--   using the given gradients
       i = 1
-      gradient = 0.0
-      ws_surface = 0.0
+      gradient = 0.0_wp
+      ws_surface = 0.0_wp
       
 
@@ -110 +110 @@
          IF ( i < 11 ) THEN
             IF ( subs_vertical_gradient_level(i) < zu(k)  .AND. &
-                 subs_vertical_gradient_level(i) >= 0.0 )  THEN
-               gradient = subs_vertical_gradient(i) / 100.0
+                 subs_vertical_gradient_level(i) >= 0.0_wp )  THEN
+               gradient = subs_vertical_gradient(i) / 100.0_wp
                subs_vertical_gradient_level_i(i) = k - 1
                i = i + 1
             ENDIF
          ENDIF
-         IF ( gradient /= 0.0 )  THEN
+         IF ( gradient /= 0.0_wp )  THEN
             IF ( k /= 1 )  THEN
                w_subs(k) = w_subs(k-1) + dzu(k) * gradient
             ELSE
-               w_subs(k) = ws_surface   + 0.5 * dzu(k) * gradient
+               w_subs(k) = ws_surface   + 0.5_wp * dzu(k) * gradient
             ENDIF
          ELSE
@@ -130 +130 @@
 !--   In case of no given gradients for the subsidence/ascent velocity,
 !--   choose zero gradient
-      IF ( subs_vertical_gradient_level(1) == -9999999.9 )  THEN
-         subs_vertical_gradient_level(1) = 0.0
+      IF ( subs_vertical_gradient_level(1) == -9999999.9_wp )  THEN
+         subs_vertical_gradient_level(1) = 0.0_wp
       ENDIF
 
@@ -171 +171 @@
           DO  j = nys, nyn
              DO  k = nzb_s_inner(j,i)+1, nzt 
-                IF ( w_subs(k) < 0.0 ) THEN    ! large-scale subsidence
+                IF ( w_subs(k) < 0.0_wp ) THEN    ! large-scale subsidence
                    tendency(k,j,i) = tendency(k,j,i) - w_subs(k) *  &
                                      ( var(k+1,j,i) - var(k,j,i) ) * ddzu(k+1)
@@ -187 +187 @@
 
        DO  k = nzb, nzt
-          IF ( w_subs(k) < 0.0 ) THEN      ! large-scale subsidence
+          IF ( w_subs(k) < 0.0_wp ) THEN      ! large-scale subsidence
              var_mod(k) = var_init(k) - dt_3d * w_subs(k) *  &
                                ( var_init(k+1) - var_init(k) ) * ddzu(k+1)
@@ -195 +195 @@
 !--   At the upper boundary, the initial profile is shifted with aid of 
 !--   the gradient tmp_grad. (This is ok if the gradients are linear.)
-      IF ( w_subs(nzt) < 0.0 ) THEN
+      IF ( w_subs(nzt) < 0.0_wp ) THEN
          tmp_grad = ( var_init(nzt+1) - var_init(nzt) ) * ddzu(nzt+1)
          var_mod(nzt+1) = var_init(nzt+1) -  &
@@ -203 +203 @@
 
       DO  k = nzt+1, nzb+1, -1
-         IF ( w_subs(k) >= 0.0 ) THEN  ! large-scale ascent
+         IF ( w_subs(k) >= 0.0_wp ) THEN  ! large-scale ascent
             var_mod(k) = var_init(k) - dt_3d * w_subs(k) *  &
                                ( var_init(k) - var_init(k-1) ) * ddzu(k) 
@@ -213 +213 @@
 
 
-      IF ( w_subs(nzb+1) >= 0.0 ) THEN
+      IF ( w_subs(nzb+1) >= 0.0_wp ) THEN
          var_mod(nzb) = var_init(nzb)
       ENDIF
@@ -253 +253 @@
 !--    Influence of w_subsidence on the current tendency term
        DO  k = nzb_s_inner(j,i)+1, nzt 
-          IF ( w_subs(k) < 0.0 ) THEN      ! large-scale subsidence
+          IF ( w_subs(k) < 0.0_wp ) THEN      ! large-scale subsidence
              tendency(k,j,i) = tendency(k,j,i) - w_subs(k) *  &
                                ( var(k+1,j,i) - var(k,j,i) ) * ddzu(k+1)
@@ -269 +269 @@
 
           DO  k = nzb, nzt
-             IF ( w_subs(k) < 0.0 ) THEN      ! large-scale subsidence
+             IF ( w_subs(k) < 0.0_wp ) THEN      ! large-scale subsidence
                 var_mod(k) = var_init(k) - dt_3d * w_subs(k) *  &
                                   ( var_init(k+1) - var_init(k) ) * ddzu(k+1)
@@ -277 +277 @@
 !--       At the upper boundary, the initial profile is shifted with aid of 
 !--       the gradient tmp_grad. (This is ok if the gradients are linear.)
-          IF ( w_subs(nzt) < 0.0 ) THEN
+          IF ( w_subs(nzt) < 0.0_wp ) THEN
              tmp_grad = ( var_init(nzt+1) - var_init(nzt) ) * ddzu(nzt+1)
              var_mod(nzt+1) = var_init(nzt+1) -  &
@@ -285 +285 @@
 
           DO  k = nzt+1, nzb+1, -1
-             IF ( w_subs(k) >= 0.0 ) THEN  ! large-scale ascent
+             IF ( w_subs(k) >= 0.0_wp ) THEN  ! large-scale ascent
                 var_mod(k) = var_init(k) - dt_3d * w_subs(k) *  &
                                    ( var_init(k) - var_init(k-1) ) * ddzu(k)
@@ -295 +295 @@
 
 
-          IF ( w_subs(nzb+1) >= 0.0 ) THEN
+          IF ( w_subs(nzb+1) >= 0.0_wp ) THEN
              var_mod(nzb) = var_init(nzb)
           ENDIF