Changeset 1353 for palm/trunk/SOURCE/advec_w_pw.f90
- Timestamp:
- Apr 8, 2014 3:21:23 PM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/advec_w_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 = nys, nyn 96 96 DO k = nzb_w_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 ( w(k,j,i+1) * ( u(k+1,j,i+1) + u(k,j,i+1) - gu ) & 99 99 - w(k,j,i-1) * ( u(k+1,j,i) + u(k,j,i) - gu ) ) * ddx & … … 103 103 - w(k-1,j,i) * ( w(k,j,i) + w(k-1,j,i) ) ) & 104 104 * ddzu(k+1) & 105 )105 ) 106 106 ENDDO 107 107 ENDDO … … 140 140 REAL(wp) :: gv !: 141 141 142 gu = 2.0 * u_gtrans143 gv = 2.0 * v_gtrans142 gu = 2.0_wp * u_gtrans 143 gv = 2.0_wp * v_gtrans 144 144 DO k = nzb_w_inner(j,i)+1, nzt 145 tend(k,j,i) = tend(k,j,i) - 0.25 * (&145 tend(k,j,i) = tend(k,j,i) - 0.25_wp * ( & 146 146 ( w(k,j,i+1) * ( u(k+1,j,i+1) + u(k,j,i+1) - gu ) & 147 147 - w(k,j,i-1) * ( u(k+1,j,i) + u(k,j,i) - gu ) ) * ddx &
Note: See TracChangeset
for help on using the changeset viewer.