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

palm/trunk/SOURCE/advec_u_up.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-!
+! REAL constants provided with KIND-attribute
 !
 ! Former revisions:
@@ -96 +96 @@
 !--             x-direction
                 ukomp = u(k,j,i) - u_gtrans
-                IF ( ukomp > 0.0 )  THEN
+                IF ( ukomp > 0.0_wp )  THEN
                    tend(k,j,i) = tend(k,j,i) - ukomp *                         &
                                          ( u(k,j,i) - u(k,j,i-1) ) * ddx
@@ -105 +105 @@
 !
 !--             y-direction
-                vkomp = 0.25 * ( v(k,j,i)   + v(k,j+1,i) +                     &
+                vkomp = 0.25_wp * ( v(k,j,i)   + v(k,j+1,i) +                  &
                                  v(k,j,i-1) + v(k,j+1,i-1) ) - v_gtrans
-                IF ( vkomp > 0.0 )  THEN
+                IF ( vkomp > 0.0_wp )  THEN
                    tend(k,j,i) = tend(k,j,i) - vkomp *                         &
                                          ( u(k,j,i) - u(k,j-1,i) ) * ddy
@@ -116 +116 @@
 !
 !--             z-direction
-                wkomp = 0.25 * ( w(k,j,i)   + w(k-1,j,i) +                     &
+                wkomp = 0.25_wp * ( w(k,j,i)   + w(k-1,j,i) +                  &
                                  w(k,j,i-1) + w(k-1,j,i-1) )
-                IF ( wkomp > 0.0 )  THEN
+                IF ( wkomp > 0.0_wp )  THEN
                    tend(k,j,i) = tend(k,j,i) - wkomp *                         &
                                          ( u(k,j,i) - u(k-1,j,i) ) * ddzu(k)
@@ -168 +168 @@
 !--       x-direction
           ukomp = u(k,j,i) - u_gtrans
-          IF ( ukomp > 0.0 )  THEN
+          IF ( ukomp > 0.0_wp )  THEN
              tend(k,j,i) = tend(k,j,i) - ukomp *                               &
                                          ( u(k,j,i) - u(k,j,i-1) ) * ddx
@@ -177 +177 @@
 !
 !--       y-direction
-          vkomp = 0.25 * ( v(k,j,i) + v(k,j+1,i) + v(k,j,i-1) + v(k,j+1,i-1)   &
+          vkomp = 0.25_wp * ( v(k,j,i) + v(k,j+1,i) + v(k,j,i-1) + v(k,j+1,i-1) &
                          ) - v_gtrans
-          IF ( vkomp > 0.0 )  THEN
+          IF ( vkomp > 0.0_wp )  THEN
              tend(k,j,i) = tend(k,j,i) - vkomp *                               &
                                          ( u(k,j,i) - u(k,j-1,i) ) * ddy
@@ -188 +188 @@
 !
 !--       z-direction
-          wkomp = 0.25 * ( w(k,j,i) + w(k-1,j,i) + w(k,j,i-1) + w(k-1,j,i-1) )
-          IF ( wkomp > 0.0 )  THEN
+          wkomp = 0.25_wp * ( w(k,j,i) + w(k-1,j,i) + w(k,j,i-1) + w(k-1,j,i-1) )
+          IF ( wkomp > 0.0_wp )  THEN
              tend(k,j,i) = tend(k,j,i) - wkomp *                               &
                                          ( u(k,j,i) - u(k-1,j,i) ) * ddzu(k)