Changeset 550

Timestamp:

Aug 24, 2010 9:17:48 AM (14 years ago)

Author:

maronga

Message:

bugfix: calculation of subgrid-scale flux at nzt for use_top_fluxes=.T.

File:

• palm/trunk/SOURCE/flow_statistics.f90 (modified) (1 diff)

palm/trunk/SOURCE/flow_statistics.f90

@@ -519 +519 @@
 !--          Subgridscale fluxes at the top surface
              IF ( use_top_fluxes )  THEN
-                sums_l(nzt,12,tn) = sums_l(nzt,12,tn) + &
+                sums_l(nzt:nzt+1,12,tn) = sums_l(nzt:nzt+1,12,tn) + &
                                     uswst(j,i) * rmask(j,i,sr)    ! w"u"
-                sums_l(nzt,14,tn) = sums_l(nzt,14,tn) + &
+                sums_l(nzt:nzt+1,14,tn) = sums_l(nzt:nzt+1,14,tn) + &
                                     vswst(j,i) * rmask(j,i,sr)    ! w"v"
-                sums_l(nzt,16,tn) = sums_l(nzt,16,tn) + &
+                sums_l(nzt:nzt+1,16,tn) = sums_l(nzt:nzt+1,16,tn) + &
                                     tswst(j,i)  * rmask(j,i,sr)   ! w"pt"
-                sums_l(nzt,58,tn) = sums_l(nzt,58,tn) + &
+                sums_l(nzt:nzt+1,58,tn) = sums_l(nzt:nzt+1,58,tn) + &
                                     0.0 * rmask(j,i,sr)           ! u"pt"
-                sums_l(nzt,61,tn) = sums_l(nzt,61,tn) + &
-                                    0.0 * rmask(j,i,sr)           ! v"pt"
+                sums_l(nzt:nzt+1,61,tn) = sums_l(nzt:nzt+1,61,tn) + &
+                                    0.0 * rmask(j,i,sr)           ! v"pt" 
+
                 IF ( ocean )  THEN
                    sums_l(nzt,65,tn) = sums_l(nzt,65,tn) + &