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

REAL constants provided with KIND-attribute

File:
1 edited

Legend:

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

    r1321 r1353  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! REAL constants provided with KIND-attribute
    2323!
    2424! Former revisions:
     
    8989       REAL(wp)    ::  gv !:
    9090 
    91        gu = 2.0 * u_gtrans
    92        gv = 2.0 * v_gtrans
     91       gu = 2.0_wp * u_gtrans
     92       gv = 2.0_wp * v_gtrans
    9393       DO  i = nxlu, nxr
    9494          DO  j = nys, nyn
    9595             DO  k = nzb_u_inner(j,i)+1, nzt
    96                 tend(k,j,i) = tend(k,j,i) - 0.25 * (                           &
     96                tend(k,j,i) = tend(k,j,i) - 0.25_wp * (                        &
    9797                         ( u(k,j,i+1) * ( u(k,j,i+1) + u(k,j,i) - gu )         &
    9898                         - u(k,j,i-1) * ( u(k,j,i) + u(k,j,i-1) - gu ) ) * ddx &
     
    102102                         - u(k-1,j,i) * ( w(k-1,j,i) + w(k-1,j,i-1) ) )        &
    103103                                                                  * ddzw(k)    &
    104                                                    )
     104                                                      )
    105105             ENDDO
    106106          ENDDO
     
    139139       REAL(wp)    ::  gv !:
    140140
    141        gu = 2.0 * u_gtrans
    142        gv = 2.0 * v_gtrans
     141       gu = 2.0_wp * u_gtrans
     142       gv = 2.0_wp * v_gtrans
    143143       DO  k = nzb_u_inner(j,i)+1, nzt
    144           tend(k,j,i) = tend(k,j,i) - 0.25 * (                                 &
     144          tend(k,j,i) = tend(k,j,i) - 0.25_wp * (                              &
    145145                         ( u(k,j,i+1) * ( u(k,j,i+1) + u(k,j,i) - gu )         &
    146146                         - u(k,j,i-1) * ( u(k,j,i) + u(k,j,i-1) - gu ) ) * ddx &
     
    150150                         - u(k-1,j,i) * ( w(k-1,j,i) + w(k-1,j,i-1) ) )        &
    151151                                                                  * ddzw(k)    &
    152                                              )
     152                                                )
    153153       ENDDO
    154154
Note: See TracChangeset for help on using the changeset viewer.