Changeset 1353 for palm/trunk/SOURCE/advec_v_pw.f90
- Timestamp:
- Apr 8, 2014 3:21:23 PM (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/advec_v_pw.f90
r1321 r1353 20 20 ! Current revisions: 21 21 ! ----------------- 22 ! 22 ! REAL constants provided with KIND-attribute 23 23 ! 24 24 ! Former revisions: … … 90 90 91 91 92 gu = 2.0 * u_gtrans93 gv = 2.0 * v_gtrans92 gu = 2.0_wp * u_gtrans 93 gv = 2.0_wp * v_gtrans 94 94 DO i = nxl, nxr 95 95 DO j = nysv, nyn 96 96 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 * ( & 98 98 ( v(k,j,i+1) * ( u(k,j-1,i+1) + u(k,j,i+1) - gu ) & 99 99 - v(k,j,i-1) * ( u(k,j-1,i) + u(k,j,i) - gu ) ) * ddx & … … 103 103 - v(k-1,j,i) * ( w(k-1,j-1,i) + w(k-1,j,i) ) ) & 104 104 * ddzw(k) & 105 )105 ) 106 106 ENDDO 107 107 ENDDO … … 141 141 142 142 143 gu = 2.0 * u_gtrans144 gv = 2.0 * v_gtrans143 gu = 2.0_wp * u_gtrans 144 gv = 2.0_wp * v_gtrans 145 145 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 * ( & 147 147 ( v(k,j,i+1) * ( u(k,j-1,i+1) + u(k,j,i+1) - gu ) & 148 148 - v(k,j,i-1) * ( u(k,j-1,i) + u(k,j,i) - gu ) ) * ddx & … … 152 152 - v(k-1,j,i) * ( w(k-1,j-1,i) + w(k-1,j,i) ) ) & 153 153 * ddzw(k) & 154 )154 ) 155 155 ENDDO 156 156
Note: See TracChangeset
for help on using the changeset viewer.