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

REAL constants provided with KIND-attribute

File:
1 edited

Legend:

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

    r1321 r1353  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! REAL constants provided with KIND-attribute
    2323!
    2424! Former revisions:
     
    9090 
    9191
    92        gu = 2.0 * u_gtrans
    93        gv = 2.0 * v_gtrans
     92       gu = 2.0_wp * u_gtrans
     93       gv = 2.0_wp * v_gtrans
    9494       DO  i = nxl, nxr
    9595          DO  j = nysv, nyn
    9696             DO  k = nzb_v_inner(j,i)+1, nzt
    97                 tend(k,j,i) = tend(k,j,i) - 0.25 * (                           &
     97                tend(k,j,i) = tend(k,j,i) - 0.25_wp * (                        &
    9898                         ( v(k,j,i+1) * ( u(k,j-1,i+1) + u(k,j,i+1) - gu )     &
    9999                         - v(k,j,i-1) * ( u(k,j-1,i) + u(k,j,i) - gu ) ) * ddx &
     
    103103                         - v(k-1,j,i) * ( w(k-1,j-1,i) + w(k-1,j,i) ) )        &
    104104                                                                  * ddzw(k)    &
    105                                                    )
     105                                                      )
    106106             ENDDO
    107107          ENDDO
     
    141141
    142142
    143        gu = 2.0 * u_gtrans
    144        gv = 2.0 * v_gtrans
     143       gu = 2.0_wp * u_gtrans
     144       gv = 2.0_wp * v_gtrans
    145145       DO  k = nzb_v_inner(j,i)+1, nzt
    146           tend(k,j,i) = tend(k,j,i) - 0.25 * (                                 &
     146          tend(k,j,i) = tend(k,j,i) - 0.25_wp * (                              &
    147147                         ( v(k,j,i+1) * ( u(k,j-1,i+1) + u(k,j,i+1) - gu )     &
    148148                         - v(k,j,i-1) * ( u(k,j-1,i) + u(k,j,i) - gu ) ) * ddx &
     
    152152                         - v(k-1,j,i) * ( w(k-1,j-1,i) + w(k-1,j,i) ) )        &
    153153                                                                  * ddzw(k)    &
    154                                              )
     154                                                )
    155155       ENDDO
    156156 
Note: See TracChangeset for help on using the changeset viewer.