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

palm/trunk/SOURCE/init_rankine.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-!
+! REAL constants provided with KIND-attribute
 !
 ! Former revisions:
@@ -89 +89 @@
 !-- Default: eddy radius rc, eddy strength z,
 !--          position of eddy centre: ic, jc, kc1, kc2
-    rc  =  4.0 * dx
+    rc  =  4.0_wp * dx
     ic  =  ( nx+1 ) / 2
     jc  =  ic
@@ -111 +111 @@
     DO  i = nxl, nxr
        DO  j = nys, nyn
-          x = ( i - ic - 0.5 ) * dx
-          y = ( j - jc ) * dy
+          x = ( i - ic - 0.5_wp ) * dx
+          y = ( j - jc          ) * dy
           radius = SQRT( x**2 + y**2 )
-          IF ( radius <= 2.0 * rc )  THEN
-             betrag = radius / ( 2.0 * rc ) * 0.08_wp
-          ELSEIF ( radius > 2.0 * rc  .AND.  radius < 8.0 * rc )  THEN
-             betrag = 0.08 * EXP( -( radius - 2.0_wp * rc ) / 2.0_wp )
-          ELSE
-             betrag = 0.0
-          ENDIF
-
-          IF ( x == 0.0 )  THEN
-             IF ( y > 0.0 )  THEN
+          IF ( radius <= 2.0_wp * rc )  THEN
+             betrag = radius / ( 2.0_wp * rc ) * 0.08_wp
+          ELSEIF ( radius > 2.0_wp * rc  .AND.  radius < 8.0_wp * rc )  THEN
+             betrag = 0.08_wp * EXP( -( radius - 2.0_wp * rc ) / 2.0_wp )
+          ELSE
+             betrag = 0.0_wp
+          ENDIF
+
+          IF ( x == 0.0_wp )  THEN
+             IF ( y > 0.0_wp )  THEN
                 alpha = pi / 2.0_wp
-             ELSEIF ( y < 0.0 )  THEN
-                alpha = 3.0 * pi / 2.0_wp
-             ENDIF
-          ELSE
-             IF ( x < 0.0 )  THEN
+             ELSEIF ( y < 0.0_wp )  THEN
+                alpha = 3.0_wp * pi / 2.0_wp
+             ENDIF
+          ELSE
+             IF ( x < 0.0_wp )  THEN
                 alpha = ATAN( y / x ) + pi
              ELSE
-                IF ( y < 0.0 )  THEN
+                IF ( y < 0.0_wp )  THEN
                    alpha = ATAN( y / x ) + 2.0_wp * pi
                 ELSE
@@ -152 +152 @@
     DO  i = nxl, nxr
        DO  j = nys, nyn
-          x = ( i - ic ) * dx
-          y = ( j - jc - 0.5) * dy
+          x = ( i - ic          ) * dx
+          y = ( j - jc - 0.5_wp ) * dy
           radius = SQRT( x**2 + y**2 )
-          IF ( radius <= 2.0 * rc )  THEN
+          IF ( radius <= 2.0_wp * rc )  THEN
              betrag = radius / ( 2.0_wp * rc ) * 0.08_wp
-          ELSEIF ( radius > 2.0 * rc  .AND.  radius < 8.0 * rc )  THEN
-             betrag = 0.08 * EXP( -( radius - 2.0_wp * rc ) / 2.0_wp )
-          ELSE
-             betrag = 0.0
-          ENDIF
-
-          IF ( x == 0.0 )  THEN
-             IF ( y > 0.0 )  THEN
+          ELSEIF ( radius > 2.0_wp * rc  .AND.  radius < 8.0_wp * rc )  THEN
+             betrag = 0.08_wp * EXP( -( radius - 2.0_wp * rc ) / 2.0_wp )
+          ELSE
+             betrag = 0.0_wp
+          ENDIF
+
+          IF ( x == 0.0_wp )  THEN
+             IF ( y > 0.0_wp )  THEN
                 alpha = pi / 2.0_wp
-             ELSEIF ( y < 0.0 )  THEN
-                alpha = 3.0 * pi / 2.0_wp
-             ENDIF
-          ELSE
-             IF ( x < 0.0 )  THEN
+             ELSEIF ( y < 0.0_wp )  THEN
+                alpha = 3.0_wp * pi / 2.0_wp
+             ENDIF
+          ELSE
+             IF ( x < 0.0_wp )  THEN
                 alpha = ATAN( y / x ) + pi
              ELSE
-                IF ( y < 0.0 )  THEN
+                IF ( y < 0.0_wp )  THEN
                    alpha = ATAN( y / x ) + 2.0_wp * pi
                 ELSE