Changeset 1353 for palm/trunk/SOURCE/advec_s_pw.f90
- Timestamp:
- Apr 8, 2014 3:21:23 PM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/advec_s_pw.f90
r1321 r1353 20 20 ! Current revisions: 21 21 ! ----------------- 22 ! 22 ! REAL constants provided with KIND-attribute 23 23 ! 24 24 ! Former revisions: … … 102 102 DO k = nzb_s_inner(j,i)+1, nzt 103 103 tend(k,j,i) = tend(k,j,i) & 104 -0.5 * ( ( u(k,j,i+1) - u_gtrans ) * ( sk(k,j,i+1) - sk(k,j,i) ) &105 - ( u(k,j,i) - u_gtrans ) * ( sk(k,j,i-1) - sk(k,j,i) ) &106 ) * ddx &107 -0.5 * ( ( v(k,j+1,i) - v_gtrans ) * ( sk(k,j+1,i) - sk(k,j,i) ) &108 - ( v(k,j,i) - v_gtrans ) * ( sk(k,j-1,i) - sk(k,j,i) ) &109 ) * ddy &110 - ( w(k,j,i) * ( sk(k+1,j,i) - sk(k,j,i) ) &111 - w(k-1,j,i) * ( sk(k-1,j,i) - sk(k,j,i) ) &112 ) * dd2zu(k)104 -0.5_wp * ( ( u(k,j,i+1) - u_gtrans ) * ( sk(k,j,i+1) - sk(k,j,i) ) & 105 - ( u(k,j,i) - u_gtrans ) * ( sk(k,j,i-1) - sk(k,j,i) ) & 106 ) * ddx & 107 -0.5_wp * ( ( v(k,j+1,i) - v_gtrans ) * ( sk(k,j+1,i) - sk(k,j,i) ) & 108 - ( v(k,j,i) - v_gtrans ) * ( sk(k,j-1,i) - sk(k,j,i) ) & 109 ) * ddy & 110 - ( w(k,j,i) * ( sk(k+1,j,i) - sk(k,j,i) ) & 111 - w(k-1,j,i) * ( sk(k-1,j,i) - sk(k,j,i) ) & 112 ) * dd2zu(k) 113 113 ENDDO 114 114 ENDDO … … 153 153 DO k = nzb_s_inner(j,i)+1, nzt 154 154 tend(k,j,i) = tend(k,j,i) & 155 -0.5 * ( ( u(k,j,i+1) - u_gtrans ) * ( sk(k,j,i+1) - sk(k,j,i) ) &156 - ( u(k,j,i) - u_gtrans ) * ( sk(k,j,i-1) - sk(k,j,i) ) &157 ) * ddx &158 -0.5 * ( ( v(k,j+1,i) - v_gtrans ) * ( sk(k,j+1,i) - sk(k,j,i) ) &159 - ( v(k,j,i) - v_gtrans ) * ( sk(k,j-1,i) - sk(k,j,i) ) &160 ) * ddy &161 - ( w(k,j,i) * ( sk(k+1,j,i) - sk(k,j,i) ) &162 - w(k-1,j,i) * ( sk(k-1,j,i) - sk(k,j,i) ) &163 ) * dd2zu(k)155 -0.5_wp * ( ( u(k,j,i+1) - u_gtrans ) * ( sk(k,j,i+1) - sk(k,j,i) ) & 156 - ( u(k,j,i) - u_gtrans ) * ( sk(k,j,i-1) - sk(k,j,i) ) & 157 ) * ddx & 158 -0.5_wp * ( ( v(k,j+1,i) - v_gtrans ) * ( sk(k,j+1,i) - sk(k,j,i) ) & 159 - ( v(k,j,i) - v_gtrans ) * ( sk(k,j-1,i) - sk(k,j,i) ) & 160 ) * ddy & 161 - ( w(k,j,i) * ( sk(k+1,j,i) - sk(k,j,i) ) & 162 - w(k-1,j,i) * ( sk(k-1,j,i) - sk(k,j,i) ) & 163 ) * dd2zu(k) 164 164 ENDDO 165 165
Note: See TracChangeset
for help on using the changeset viewer.