Changeset 183 for palm/trunk/SOURCE/diffusion_s.f90
- Timestamp:
- Aug 4, 2008 3:39:12 PM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/diffusion_s.f90
r139 r183 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! 6 ! bugfix: calculation of fluxes at vertical surfaces 7 7 ! 8 8 ! Former revisions: … … 90 90 91 91 tend(k,j,i) = tend(k,j,i) & 92 + 0.5 * ( fwxp(j,i) *&92 + ( fwxp(j,i) * 0.5 * & 93 93 ( kh(k,j,i) + kh(k,j,i+1) ) * ( s(k,j,i+1)-s(k,j,i) ) & 94 94 + ( 1.0 - fwxp(j,i) ) * wall_s_flux(1) & 95 -fwxm(j,i) * 95 -fwxm(j,i) * 0.5 * & 96 96 ( kh(k,j,i) + kh(k,j,i-1) ) * ( s(k,j,i)-s(k,j,i-1) ) & 97 97 + ( 1.0 - fwxm(j,i) ) * wall_s_flux(2) & 98 98 ) * ddx2 & 99 + 0.5 * ( fwyp(j,i) *&99 + ( fwyp(j,i) * 0.5 * & 100 100 ( kh(k,j,i) + kh(k,j+1,i) ) * ( s(k,j+1,i)-s(k,j,i) ) & 101 101 + ( 1.0 - fwyp(j,i) ) * wall_s_flux(3) & 102 -fwym(j,i) * 102 -fwym(j,i) * 0.5 * & 103 103 ( kh(k,j,i) + kh(k,j-1,i) ) * ( s(k,j,i)-s(k,j-1,i) ) & 104 104 + ( 1.0 - fwym(j,i) ) * wall_s_flux(4) & … … 200 200 201 201 tend(k,j,i) = tend(k,j,i) & 202 + 0.5 * ( fwxp(j,i) *&202 + ( fwxp(j,i) * 0.5 * & 203 203 ( kh(k,j,i) + kh(k,j,i+1) ) * ( s(k,j,i+1)-s(k,j,i) ) & 204 204 + ( 1.0 - fwxp(j,i) ) * wall_s_flux(1) & 205 -fwxm(j,i) * 205 -fwxm(j,i) * 0.5 * & 206 206 ( kh(k,j,i) + kh(k,j,i-1) ) * ( s(k,j,i)-s(k,j,i-1) ) & 207 207 + ( 1.0 - fwxm(j,i) ) * wall_s_flux(2) & 208 208 ) * ddx2 & 209 + 0.5 * ( fwyp(j,i) *&209 + ( fwyp(j,i) * 0.5 * & 210 210 ( kh(k,j,i) + kh(k,j+1,i) ) * ( s(k,j+1,i)-s(k,j,i) ) & 211 211 + ( 1.0 - fwyp(j,i) ) * wall_s_flux(3) & 212 -fwym(j,i) * 212 -fwym(j,i) * 0.5 * & 213 213 ( kh(k,j,i) + kh(k,j-1,i) ) * ( s(k,j,i)-s(k,j-1,i) ) & 214 214 + ( 1.0 - fwym(j,i) ) * wall_s_flux(4) &
Note: See TracChangeset
for help on using the changeset viewer.