Changeset 678 for palm/trunk/SOURCE/flow_statistics.f90
 Timestamp:
 Feb 2, 2011 2:31:56 PM (12 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SOURCE/flow_statistics.f90
r674 r678 4 4 ! Current revisions: 5 5 !  6 ! Bugfix in calculation of divergence of vertical flux of resolved scale 7 ! energy, pressure fluctuations, and flux of pressure fluctuation itself 6 8 ! 7 9 ! Former revisions: … … 802 804 sums_ll(k,1) = sums_ll(k,1) + 0.5 * w(k,j,i) * ( & 803 805 ( 0.25 * ( u(k,j,i)+u(k+1,j,i)+u(k,j,i+1)+u(k+1,j,i+1) & 804  2.0* ( hom(k,1,1,sr) + hom(k+1,1,1,sr) ) &806  0.5 * ( hom(k,1,1,sr) + hom(k+1,1,1,sr) ) & 805 807 ) )**2 & 806 808 + ( 0.25 * ( v(k,j,i)+v(k+1,j,i)+v(k,j+1,i)+v(k+1,j+1,i) & 807  2.0* ( hom(k,1,2,sr) + hom(k+1,1,2,sr) ) &809  0.5 * ( hom(k,1,2,sr) + hom(k+1,1,2,sr) ) & 808 810 ) )**2 & 809 811 + w(k,j,i)**2 ) … … 820 822 sums_ll(nzt+1,2) = 0.0 821 823 822 DO k = nzb _s_inner(j,i)+1, nzt824 DO k = nzb+1, nzt 823 825 sums_l(k,55,tn) = ( sums_ll(k,1)  sums_ll(k1,1) ) * ddzw(k) 824 826 sums_l(k,56,tn) = ( sums_ll(k,2)  sums_ll(k1,2) ) * ddzw(k)
Note: See TracChangeset
for help on using the changeset viewer.