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

palm/trunk/SOURCE/timestep_scheme_steering.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
+! REAL constants defined as wp-kind
 ! 
 ! Former revisions:
@@ -48 +48 @@
         ONLY:  intermediate_timestep_count, timestep_scheme, tsc
 
+    USE kinds
+
     IMPLICIT NONE
 
@@ -57 +59 @@
        IF ( timestep_scheme == 'runge-kutta-2' )  THEN
           IF ( intermediate_timestep_count == 1 )  THEN
-             tsc(1:5) = (/ 1.0, 1.0,  0.0, 0.0, 0.0 /)
+             tsc(1:5) = (/ 1.0_wp, 1.0_wp,  0.0_wp, 0.0_wp, 0.0_wp /)
           ELSE
-             tsc(1:5) = (/ 1.0, 0.5, -0.5, 0.0, 1.0 /)
+             tsc(1:5) = (/ 1.0_wp, 0.5_wp, -0.5_wp, 0.0_wp, 1.0_wp /)
           ENDIF
        ELSE
           IF ( intermediate_timestep_count == 1 )  THEN
-             tsc(1:5) = (/ 1.0,   1.0/3.0,        0.0, 0.0, 0.0 /)
+             tsc(1:5) = (/ 1.0_wp,  1.0_wp /  3.0_wp,           0.0_wp, 0.0_wp, 0.0_wp /)
           ELSEIF ( intermediate_timestep_count == 2 )  THEN
-             tsc(1:5) = (/ 1.0, 15.0/16.0, -25.0/48.0, 0.0, 0.0 /)
+             tsc(1:5) = (/ 1.0_wp, 15.0_wp / 16.0_wp, -25.0_wp/48.0_wp, 0.0_wp, 0.0_wp /)
           ELSE
-             tsc(1:5) = (/ 1.0,  8.0/15.0,   1.0/15.0, 0.0, 1.0 /)
+             tsc(1:5) = (/ 1.0_wp,  8.0_wp / 15.0_wp,   1.0_wp/15.0_wp, 0.0_wp, 1.0_wp /)
           ENDIF          
        ENDIF
@@ -74 +76 @@
 !
 !--    Euler scheme
-       tsc(1:5) = (/ 1.0, 1.0, 0.0, 0.0, 1.0 /)
+       tsc(1:5) = (/ 1.0_wp, 1.0_wp, 0.0_wp, 0.0_wp, 1.0_wp /)
 
     ENDIF