Changeset 1353 for palm/trunk/SOURCE/advec_u_pw.f90
- Timestamp:
- Apr 8, 2014 3:21:23 PM (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/advec_u_pw.f90
r1321 r1353 20 20 ! Current revisions: 21 21 ! ----------------- 22 ! 22 ! REAL constants provided with KIND-attribute 23 23 ! 24 24 ! Former revisions: … … 89 89 REAL(wp) :: gv !: 90 90 91 gu = 2.0 * u_gtrans92 gv = 2.0 * v_gtrans91 gu = 2.0_wp * u_gtrans 92 gv = 2.0_wp * v_gtrans 93 93 DO i = nxlu, nxr 94 94 DO j = nys, nyn 95 95 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 * ( & 97 97 ( u(k,j,i+1) * ( u(k,j,i+1) + u(k,j,i) - gu ) & 98 98 - u(k,j,i-1) * ( u(k,j,i) + u(k,j,i-1) - gu ) ) * ddx & … … 102 102 - u(k-1,j,i) * ( w(k-1,j,i) + w(k-1,j,i-1) ) ) & 103 103 * ddzw(k) & 104 )104 ) 105 105 ENDDO 106 106 ENDDO … … 139 139 REAL(wp) :: gv !: 140 140 141 gu = 2.0 * u_gtrans142 gv = 2.0 * v_gtrans141 gu = 2.0_wp * u_gtrans 142 gv = 2.0_wp * v_gtrans 143 143 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 * ( & 145 145 ( u(k,j,i+1) * ( u(k,j,i+1) + u(k,j,i) - gu ) & 146 146 - u(k,j,i-1) * ( u(k,j,i) + u(k,j,i-1) - gu ) ) * ddx & … … 150 150 - u(k-1,j,i) * ( w(k-1,j,i) + w(k-1,j,i-1) ) ) & 151 151 * ddzw(k) & 152 )152 ) 153 153 ENDDO 154 154
Note: See TracChangeset
for help on using the changeset viewer.