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_w_up.f90

    r1321 r1353  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! REAL constants provided with KIND-attribute
    2323!
    2424! Former revisions:
     
    9494!
    9595!--             x-direction
    96                 ukomp = 0.25 * ( u(k,j,i)   + u(k,j,i+1) +                     &
     96                ukomp = 0.25_wp * ( u(k,j,i)   + u(k,j,i+1) +                  &
    9797                                 u(k+1,j,i) + u(k+1,j,i+1) ) - u_gtrans
    98                 IF ( ukomp > 0.0 )  THEN
     98                IF ( ukomp > 0.0_wp )  THEN
    9999                   tend(k,j,i) = tend(k,j,i) - ukomp *                         &
    100100                                         ( w(k,j,i) - w(k,j,i-1) ) * ddx
     
    105105!
    106106!--             y-direction
    107                 vkomp = 0.25 * ( v(k,j,i) + v(k,j+1,i) +                       &
     107                vkomp = 0.25_wp * ( v(k,j,i) + v(k,j+1,i) +                    &
    108108                                 v(k+1,j,i) + v(k+1,j+1,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                                         ( w(k,j,i) - w(k,j-1,i) ) * ddy
     
    116116!
    117117!--             z-direction
    118                 IF ( w(k,j,i) > 0.0 )  THEN
     118                IF ( w(k,j,i) > 0.0_wp )  THEN
    119119                   tend(k,j,i) = tend(k,j,i) - w(k,j,i) *                      &
    120120                                         ( w(k,j,i) - w(k-1,j,i) ) * ddzw(k)
     
    164164!
    165165!--       x-direction
    166           ukomp = 0.25 * ( u(k,j,i) + u(k,j,i+1) + u(k+1,j,i) + u(k+1,j,i+1)  &
     166          ukomp = 0.25_wp * ( u(k,j,i) + u(k,j,i+1) + u(k+1,j,i) + u(k+1,j,i+1) &
    167167                         ) - u_gtrans
    168           IF ( ukomp > 0.0 )  THEN
     168          IF ( ukomp > 0.0_wp )  THEN
    169169             tend(k,j,i) = tend(k,j,i) - ukomp *                               &
    170170                                         ( w(k,j,i) - w(k,j,i-1) ) * ddx
     
    175175!
    176176!--       y-direction
    177           vkomp = 0.25 * ( v(k,j,i) + v(k,j+1,i) + v(k+1,j,i) + v(k+1,j+1,i)  &
     177          vkomp = 0.25_wp * ( v(k,j,i) + v(k,j+1,i) + v(k+1,j,i) + v(k+1,j+1,i) &
    178178                         ) - v_gtrans
    179           IF ( vkomp > 0.0 )  THEN
     179          IF ( vkomp > 0.0_wp )  THEN
    180180             tend(k,j,i) = tend(k,j,i) - vkomp *                               &
    181181                                         ( w(k,j,i) - w(k,j-1,i) ) * ddy
     
    186186!
    187187!--       z-direction
    188           IF ( w(k,j,i) > 0.0 )  THEN
     188          IF ( w(k,j,i) > 0.0_wp )  THEN
    189189             tend(k,j,i) = tend(k,j,i) - w(k,j,i) *                            &
    190190                                         ( w(k,j,i) - w(k-1,j,i) ) * ddzw(k)
Note: See TracChangeset for help on using the changeset viewer.