Ignore:
Timestamp:
Apr 8, 2014 3:21:23 PM (7 years ago)
Author:
heinze
Message:

REAL constants provided with KIND-attribute

File:
1 edited

Legend:

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

    r1347 r1353  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! REAL constants provided with KIND-attribute
    2323!
    2424! Former revisions:
     
    9393
    9494       delt  = 1.0_wp / REAL( intervals, KIND=wp )
    95        sterm = delt * 0.5
     95       sterm = delt * 0.5_wp
    9696
    9797       DO  i = 1, intervals
    9898
    99           IF ( sterm > 0.5 )  THEN
    100              dn = -5.0
     99          IF ( sterm > 0.5_wp )  THEN
     100             dn = -5.0_wp
    101101          ELSE
    102              dn = 5.0
     102             dn = 5.0_wp
    103103          ENDIF
    104104
    105105          DO  j = 1, 15
    106              ex1 = dn * EXP( -dn ) - EXP( 0.5 * dn ) + EXP( -0.5 * dn )
     106             ex1 = dn * EXP( -dn ) - EXP( 0.5_wp * dn ) + EXP( -0.5_wp * dn )
    107107             ex2 = EXP( dn ) - EXP( -dn )
    108              ex3 = EXP( -dn ) * ( 1.0 - dn ) - 0.5 * EXP(  0.5 * dn ) &
    109                                              - 0.5 * EXP( -0.5 * dn )
     108             ex3 = EXP( -dn ) * ( 1.0_wp - dn ) - 0.5_wp * EXP(  0.5_wp * dn ) &
     109                                                - 0.5_wp * EXP( -0.5_wp * dn )
    110110             ex4 = EXP( dn ) + EXP( -dn )
    111111             ex5 = dn * sterm + ex1 / ex2
     
    115115          ENDDO
    116116
    117           IF ( sterm < 0.5 )  dn = MAX(  2.95E-2_wp, dn )
    118           IF ( sterm > 0.5 )  dn = MIN( -2.95E-2_wp, dn )
     117          IF ( sterm < 0.5_wp )  dn = MAX(  2.95E-2_wp, dn )
     118          IF ( sterm > 0.5_wp )  dn = MIN( -2.95E-2_wp, dn )
    119119          ex1 = EXP( -dn )
    120120          ex2 = EXP( dn ) - ex1
    121121          aex(i) = -ex1 / ex2
    122           bex(i) = 1.0 / ex2
     122          bex(i) = 1.0_wp / ex2
    123123          dex(i) = dn
    124           eex(i) = EXP( dex(i) * 0.5 )
     124          eex(i) = EXP( dex(i) * 0.5_wp )
    125125          sterm = sterm + delt
    126126
Note: See TracChangeset for help on using the changeset viewer.