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

palm/trunk/SOURCE/advec_s_up.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-!
+! REAL constants provided with KIND-attribute
 !
 ! Former revisions:
@@ -106 +106 @@
 !
 !--             x-direction
-                ukomp = 0.5 * ( u(k,j,i) + u(k,j,i+1) ) - u_gtrans
-                IF ( ukomp > 0.0 )  THEN
+                ukomp = 0.5_wp * ( u(k,j,i) + u(k,j,i+1) ) - u_gtrans
+                IF ( ukomp > 0.0_wp )  THEN
                    tend(k,j,i) = tend(k,j,i) - ukomp *                         &
                                          ( sk(k,j,i) - sk(k,j,i-1) ) * ddx
@@ -116 +116 @@
 !
 !--             y-direction
-                vkomp = 0.5 * ( v(k,j,i) + v(k,j+1,i) ) - v_gtrans
-                IF ( vkomp > 0.0 )  THEN
+                vkomp = 0.5_wp * ( v(k,j,i) + v(k,j+1,i) ) - v_gtrans
+                IF ( vkomp > 0.0_wp )  THEN
                    tend(k,j,i) = tend(k,j,i) - vkomp *                         &
                                          ( sk(k,j,i) - sk(k,j-1,i) ) * ddy
@@ -126 +126 @@
 !
 !--             z-direction
-                wkomp = 0.5 * ( w(k,j,i) + w(k-1,j,i) )
-                IF ( wkomp > 0.0 )  THEN
+                wkomp = 0.5_wp * ( w(k,j,i) + w(k-1,j,i) )
+                IF ( wkomp > 0.0_wp )  THEN
                    tend(k,j,i) = tend(k,j,i) - wkomp *                         &
                                          ( sk(k,j,i) - sk(k-1,j,i) ) * ddzu(k)
@@ -182 +182 @@
 !
 !--       x-direction
-          ukomp = 0.5 * ( u(k,j,i) + u(k,j,i+1) ) - u_gtrans
-          IF ( ukomp > 0.0 )  THEN
+          ukomp = 0.5_wp * ( u(k,j,i) + u(k,j,i+1) ) - u_gtrans
+          IF ( ukomp > 0.0_wp )  THEN
              tend(k,j,i) = tend(k,j,i) - ukomp *                               &
                                          ( sk(k,j,i) - sk(k,j,i-1) ) * ddx
@@ -192 +192 @@
 !
 !--       y-direction
-          vkomp = 0.5 * ( v(k,j,i) + v(k,j+1,i) ) - v_gtrans
-          IF ( vkomp > 0.0 )  THEN
+          vkomp = 0.5_wp * ( v(k,j,i) + v(k,j+1,i) ) - v_gtrans
+          IF ( vkomp > 0.0_wp )  THEN
              tend(k,j,i) = tend(k,j,i) - vkomp *                               &
                                          ( sk(k,j,i) - sk(k,j-1,i) ) * ddy
@@ -202 +202 @@
 !
 !--       z-direction
-          wkomp = 0.5 * ( w(k,j,i) + w(k-1,j,i) )
-          IF ( wkomp > 0.0 )  THEN
+          wkomp = 0.5_wp * ( w(k,j,i) + w(k-1,j,i) )
+          IF ( wkomp > 0.0_wp )  THEN
              tend(k,j,i) = tend(k,j,i) - wkomp *                               &
                                          ( sk(k,j,i) - sk(k-1,j,i) ) * ddzu(k)