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/advec_v_up.f90

    r1321 r1353  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! REAL constants provided with KIND-attribute
    2323!
    2424! Former revisions:
     
    9595!
    9696!--             x-direction
    97                 ukomp = 0.25 * ( u(k,j,i)   + u(k,j-1,i) +                     &
     97                ukomp = 0.25_wp * ( u(k,j,i)   + u(k,j-1,i) +                  &
    9898                                 u(k,j,i+1) + u(k,j-1,i+1) ) - u_gtrans
    99                 IF ( ukomp > 0.0 )  THEN
     99                IF ( ukomp > 0.0_wp )  THEN
    100100                   tend(k,j,i) = tend(k,j,i) - ukomp *                         &
    101101                                         ( v(k,j,i) - v(k,j,i-1) ) * ddx
     
    107107!--             y-direction
    108108                vkomp = v(k,j,i) - v_gtrans
    109                 IF ( vkomp > 0.0 )  THEN
     109                IF ( vkomp > 0.0_wp )  THEN
    110110                   tend(k,j,i) = tend(k,j,i) - vkomp *                         &
    111111                                         ( v(k,j,i) - v(k,j-1,i) ) * ddy
     
    116116!
    117117!--             z-direction
    118                 wkomp = 0.25 * ( w(k,j,i)  + w(k-1,j,i) +                      &
     118                wkomp = 0.25_wp * ( w(k,j,i)  + w(k-1,j,i) +                   &
    119119                                 w(k,j-1,i) + w(k-1,j-1,i) )
    120                 IF ( wkomp > 0.0 )  THEN
     120                IF ( wkomp > 0.0_wp )  THEN
    121121                   tend(k,j,i) = tend(k,j,i) - wkomp *                         &
    122122                                         ( v(k,j,i) - v(k-1,j,i) ) * ddzu(k)
     
    167167!
    168168!--       x-direction
    169           ukomp = 0.25 * ( u(k,j,i) + u(k,j-1,i) + u(k,j,i+1) + u(k,j-1,i+1)  &
     169          ukomp = 0.25_wp * ( u(k,j,i) + u(k,j-1,i) + u(k,j,i+1) + u(k,j-1,i+1) &
    170170                         ) - u_gtrans
    171           IF ( ukomp > 0.0 )  THEN
     171          IF ( ukomp > 0.0_wp )  THEN
    172172             tend(k,j,i) = tend(k,j,i) - ukomp *                               &
    173173                                         ( v(k,j,i) - v(k,j,i-1) ) * ddx
     
    179179!--       y-direction
    180180          vkomp = v(k,j,i) - v_gtrans
    181           IF ( vkomp > 0.0 )  THEN
     181          IF ( vkomp > 0.0_wp )  THEN
    182182             tend(k,j,i) = tend(k,j,i) - vkomp *                               &
    183183                                         ( v(k,j,i) - v(k,j-1,i) ) * ddy
     
    188188!
    189189!--       z-direction
    190           wkomp = 0.25 * ( w(k,j,i) + w(k-1,j,i) + w(k,j-1,i) + w(k-1,j-1,i) )
    191           IF ( wkomp > 0.0 )  THEN
     190          wkomp = 0.25_wp * ( w(k,j,i) + w(k-1,j,i) + w(k,j-1,i) + w(k-1,j-1,i) )
     191          IF ( wkomp > 0.0_wp )  THEN
    192192             tend(k,j,i) = tend(k,j,i) - wkomp *                               &
    193193                                         ( v(k,j,i) - v(k-1,j,i) ) * ddzu(k)
Note: See TracChangeset for help on using the changeset viewer.