Changeset 1322 for palm/trunk/SOURCE/init_rankine.f90

Timestamp:

Mar 20, 2014 4:38:49 PM (10 years ago)

Author:

raasch

Message:

REAL functions and a lot of REAL constants provided with KIND-attribute,
some small bugfixes

File:

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

palm/trunk/SOURCE/init_rankine.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
+! REAL constants defined as wp_kind
 !
 ! Former revisions:
@@ -111 +112 @@
           radius = SQRT( x**2 + y**2 )
           IF ( radius <= 2.0 * rc )  THEN
-             betrag = radius / ( 2.0 * rc ) * 0.08
+             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 * rc ) / 2.0 )
+             betrag = 0.08 * EXP( -( radius - 2.0_wp * rc ) / 2.0_wp )
           ELSE
              betrag = 0.0
@@ -120 +121 @@
           IF ( x == 0.0 )  THEN
              IF ( y > 0.0 )  THEN
-                alpha = pi / 2.0
+                alpha = pi / 2.0_wp
              ELSEIF ( y < 0.0 )  THEN
-                alpha = 3.0 * pi / 2.0
+                alpha = 3.0 * pi / 2.0_wp
              ENDIF
           ELSE
@@ -129 +130 @@
              ELSE
                 IF ( y < 0.0 )  THEN
-                   alpha = ATAN( y / x ) + 2.0 * pi
+                   alpha = ATAN( y / x ) + 2.0_wp * pi
                 ELSE
                    alpha = ATAN( y / x )
@@ -152 +153 @@
           radius = SQRT( x**2 + y**2 )
           IF ( radius <= 2.0 * rc )  THEN
-             betrag = radius / ( 2.0 * rc ) * 0.08
+             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 * rc ) / 2.0 )
+             betrag = 0.08 * EXP( -( radius - 2.0_wp * rc ) / 2.0_wp )
           ELSE
              betrag = 0.0
@@ -161 +162 @@
           IF ( x == 0.0 )  THEN
              IF ( y > 0.0 )  THEN
-                alpha = pi / 2.0
+                alpha = pi / 2.0_wp
              ELSEIF ( y < 0.0 )  THEN
-                alpha = 3.0 * pi / 2.0
+                alpha = 3.0 * pi / 2.0_wp
              ENDIF
           ELSE
@@ -170 +171 @@
              ELSE
                 IF ( y < 0.0 )  THEN
-                   alpha = ATAN( y / x ) + 2.0 * pi
+                   alpha = ATAN( y / x ) + 2.0_wp * pi
                 ELSE
                    alpha = ATAN( y / x )