Ignore:
Timestamp:
Mar 20, 2014 4:38:49 PM (7 years ago)
Author:
raasch
Message:

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

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/timestep_scheme_steering.f90

    r1321 r1322  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! REAL constants defined as wp-kind
    2323!
    2424! Former revisions:
     
    4848        ONLY:  intermediate_timestep_count, timestep_scheme, tsc
    4949
     50    USE kinds
     51
    5052    IMPLICIT NONE
    5153
     
    5759       IF ( timestep_scheme == 'runge-kutta-2' )  THEN
    5860          IF ( intermediate_timestep_count == 1 )  THEN
    59              tsc(1:5) = (/ 1.0, 1.0,  0.0, 0.0, 0.0 /)
     61             tsc(1:5) = (/ 1.0_wp, 1.0_wp,  0.0_wp, 0.0_wp, 0.0_wp /)
    6062          ELSE
    61              tsc(1:5) = (/ 1.0, 0.5, -0.5, 0.0, 1.0 /)
     63             tsc(1:5) = (/ 1.0_wp, 0.5_wp, -0.5_wp, 0.0_wp, 1.0_wp /)
    6264          ENDIF
    6365       ELSE
    6466          IF ( intermediate_timestep_count == 1 )  THEN
    65              tsc(1:5) = (/ 1.0,   1.0/3.0,        0.0, 0.0, 0.0 /)
     67             tsc(1:5) = (/ 1.0_wp,  1.0_wp /  3.0_wp,           0.0_wp, 0.0_wp, 0.0_wp /)
    6668          ELSEIF ( intermediate_timestep_count == 2 )  THEN
    67              tsc(1:5) = (/ 1.0, 15.0/16.0, -25.0/48.0, 0.0, 0.0 /)
     69             tsc(1:5) = (/ 1.0_wp, 15.0_wp / 16.0_wp, -25.0_wp/48.0_wp, 0.0_wp, 0.0_wp /)
    6870          ELSE
    69              tsc(1:5) = (/ 1.0,  8.0/15.0,   1.0/15.0, 0.0, 1.0 /)
     71             tsc(1:5) = (/ 1.0_wp,  8.0_wp / 15.0_wp,   1.0_wp/15.0_wp, 0.0_wp, 1.0_wp /)
    7072          ENDIF         
    7173       ENDIF
     
    7476!
    7577!--    Euler scheme
    76        tsc(1:5) = (/ 1.0, 1.0, 0.0, 0.0, 1.0 /)
     78       tsc(1:5) = (/ 1.0_wp, 1.0_wp, 0.0_wp, 0.0_wp, 1.0_wp /)
    7779
    7880    ENDIF
Note: See TracChangeset for help on using the changeset viewer.