Ignore:
Timestamp:
Mar 20, 2014 4:38:49 PM (8 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/init_rankine.f90

    r1321 r1322  
    2020! Current revisions:
    2121! -----------------
     22! REAL constants defined as wp_kind
    2223!
    2324! Former revisions:
     
    111112          radius = SQRT( x**2 + y**2 )
    112113          IF ( radius <= 2.0 * rc )  THEN
    113              betrag = radius / ( 2.0 * rc ) * 0.08
     114             betrag = radius / ( 2.0 * rc ) * 0.08_wp
    114115          ELSEIF ( radius > 2.0 * rc  .AND.  radius < 8.0 * rc )  THEN
    115              betrag = 0.08 * EXP( -( radius - 2.0 * rc ) / 2.0 )
     116             betrag = 0.08 * EXP( -( radius - 2.0_wp * rc ) / 2.0_wp )
    116117          ELSE
    117118             betrag = 0.0
     
    120121          IF ( x == 0.0 )  THEN
    121122             IF ( y > 0.0 )  THEN
    122                 alpha = pi / 2.0
     123                alpha = pi / 2.0_wp
    123124             ELSEIF ( y < 0.0 )  THEN
    124                 alpha = 3.0 * pi / 2.0
     125                alpha = 3.0 * pi / 2.0_wp
    125126             ENDIF
    126127          ELSE
     
    129130             ELSE
    130131                IF ( y < 0.0 )  THEN
    131                    alpha = ATAN( y / x ) + 2.0 * pi
     132                   alpha = ATAN( y / x ) + 2.0_wp * pi
    132133                ELSE
    133134                   alpha = ATAN( y / x )
     
    152153          radius = SQRT( x**2 + y**2 )
    153154          IF ( radius <= 2.0 * rc )  THEN
    154              betrag = radius / ( 2.0 * rc ) * 0.08
     155             betrag = radius / ( 2.0_wp * rc ) * 0.08_wp
    155156          ELSEIF ( radius > 2.0 * rc  .AND.  radius < 8.0 * rc )  THEN
    156              betrag = 0.08 * EXP( -( radius - 2.0 * rc ) / 2.0 )
     157             betrag = 0.08 * EXP( -( radius - 2.0_wp * rc ) / 2.0_wp )
    157158          ELSE
    158159             betrag = 0.0
     
    161162          IF ( x == 0.0 )  THEN
    162163             IF ( y > 0.0 )  THEN
    163                 alpha = pi / 2.0
     164                alpha = pi / 2.0_wp
    164165             ELSEIF ( y < 0.0 )  THEN
    165                 alpha = 3.0 * pi / 2.0
     166                alpha = 3.0 * pi / 2.0_wp
    166167             ENDIF
    167168          ELSE
     
    170171             ELSE
    171172                IF ( y < 0.0 )  THEN
    172                    alpha = ATAN( y / x ) + 2.0 * pi
     173                   alpha = ATAN( y / x ) + 2.0_wp * pi
    173174                ELSE
    174175                   alpha = ATAN( y / x )
Note: See TracChangeset for help on using the changeset viewer.