Index: palm/trunk/SOURCE/advec_ws.f90
===================================================================
--- palm/trunk/SOURCE/advec_ws.f90	(revision 1572)
+++ palm/trunk/SOURCE/advec_ws.f90	(revision 1580)
@@ -20,5 +20,5 @@
 ! Current revisions:
 ! ------------------
-! 
+! Bugfix: statistical evaluation of scalar fluxes in case of monotonic limiter.
 ! 
 ! Former revisions:
@@ -1121,53 +1121,49 @@
 
 !
-!--    Evaluation of statistics. Please note, in case of using monotone limiter
-!--    vertical fluxes will be not calculated by the advective fluxes.
-       IF ( .NOT. monotonic_adjustment )  THEN
-
-          SELECT CASE ( sk_char )
-
-             CASE ( 'pt' )
-
-                DO  k = nzb, nzt
-                   sums_wspts_ws_l(k,tn) = sums_wspts_ws_l(k,tn) +            &
-                                          ( flux_t(k) + diss_t(k) )           &
-                                  * weight_substep(intermediate_timestep_count)
-                ENDDO
-             
-             CASE ( 'sa' )
-
-                DO  k = nzb, nzt
-                   sums_wssas_ws_l(k,tn) = sums_wssas_ws_l(k,tn) +            &
-                                          ( flux_t(k) + diss_t(k) )           &
-                                  * weight_substep(intermediate_timestep_count)
-                ENDDO
-             
-             CASE ( 'q' )
-
-                DO  k = nzb, nzt
-                   sums_wsqs_ws_l(k,tn)  = sums_wsqs_ws_l(k,tn) +             &
-                                         ( flux_t(k) + diss_t(k) )            &
-                                  * weight_substep(intermediate_timestep_count)
-                ENDDO
-
-             CASE ( 'qr' )
-
-                DO  k = nzb, nzt
-                   sums_wsqrs_ws_l(k,tn)  = sums_wsqrs_ws_l(k,tn) +           &
-                                         ( flux_t(k) + diss_t(k) )            &
-                                  * weight_substep(intermediate_timestep_count)
-                ENDDO
-
-             CASE ( 'nr' )
-
-                DO  k = nzb, nzt
-                   sums_wsnrs_ws_l(k,tn)  = sums_wsnrs_ws_l(k,tn) +           &
-                                         ( flux_t(k) + diss_t(k) )            &
-                                  * weight_substep(intermediate_timestep_count)
-                ENDDO
-
-            END SELECT
-         ENDIF
-
+!--    Evaluation of statistics.
+       SELECT CASE ( sk_char )
+
+          CASE ( 'pt' )
+
+             DO  k = nzb, nzt
+                sums_wspts_ws_l(k,tn) = sums_wspts_ws_l(k,tn) +               &
+                                       ( flux_t(k) + diss_t(k) )              &
+                               * weight_substep(intermediate_timestep_count)
+             ENDDO
+            
+          CASE ( 'sa' )
+
+             DO  k = nzb, nzt
+                sums_wssas_ws_l(k,tn) = sums_wssas_ws_l(k,tn) +               &
+                                       ( flux_t(k) + diss_t(k) )              &
+                               * weight_substep(intermediate_timestep_count)
+             ENDDO
+            
+          CASE ( 'q' )
+
+             DO  k = nzb, nzt
+                sums_wsqs_ws_l(k,tn)  = sums_wsqs_ws_l(k,tn) +                &
+                                      ( flux_t(k) + diss_t(k) )               &
+                               * weight_substep(intermediate_timestep_count)
+             ENDDO
+
+          CASE ( 'qr' )
+
+             DO  k = nzb, nzt
+                sums_wsqrs_ws_l(k,tn)  = sums_wsqrs_ws_l(k,tn) +              &
+                                      ( flux_t(k) + diss_t(k) )               &
+                               * weight_substep(intermediate_timestep_count)
+             ENDDO
+
+          CASE ( 'nr' )
+
+             DO  k = nzb, nzt
+                sums_wsnrs_ws_l(k,tn)  = sums_wsnrs_ws_l(k,tn) +              &
+                                      ( flux_t(k) + diss_t(k) )               &
+                               * weight_substep(intermediate_timestep_count)
+             ENDDO
+
+         END SELECT
+         
     END SUBROUTINE advec_s_ws_ij
 
@@ -3251,45 +3247,39 @@
              ENDDO
 !
-!--          Evaluation of statistics. Please note, in case of using monotone
-!--          limiter vertical fluxes will be not calculated by the advective 
-!--          fluxes.
-             IF ( .NOT. monotonic_adjustment )  THEN
-
-                SELECT CASE ( sk_char )
-
-                    CASE ( 'pt' )
-                       DO  k = nzb, nzt
-                          sums_wspts_ws_l(k,tn) = sums_wspts_ws_l(k,tn)       &
-                           + ( flux_t(k) + diss_t(k) )                        &
-                           *   weight_substep(intermediate_timestep_count)
-                       ENDDO
-                    CASE ( 'sa' )
-                       DO  k = nzb, nzt
-                          sums_wssas_ws_l(k,tn) = sums_wssas_ws_l(k,tn)       &
-                           + ( flux_t(k) + diss_t(k) )                        &
-                           *   weight_substep(intermediate_timestep_count)
-                       ENDDO
-                    CASE ( 'q' )
-                       DO  k = nzb, nzt
-                          sums_wsqs_ws_l(k,tn) = sums_wsqs_ws_l(k,tn)         &
-                           + ( flux_t(k) + diss_t(k) )                        &
-                           *   weight_substep(intermediate_timestep_count)
-                       ENDDO
-                    CASE ( 'qr' )
-                       DO  k = nzb, nzt
-                          sums_wsqrs_ws_l(k,tn) = sums_wsqrs_ws_l(k,tn)       &
-                           + ( flux_t(k) + diss_t(k) )                        &
-                           *   weight_substep(intermediate_timestep_count)
-                       ENDDO
-                    CASE ( 'nr' )
-                       DO  k = nzb, nzt
-                          sums_wsnrs_ws_l(k,tn) = sums_wsnrs_ws_l(k,tn)       &
-                           + ( flux_t(k) + diss_t(k) )                        &
-                           *   weight_substep(intermediate_timestep_count)
-                       ENDDO
-
-                 END SELECT
-
-              ENDIF
+!--          Evaluation of statistics. 
+             SELECT CASE ( sk_char )
+
+                 CASE ( 'pt' )
+                    DO  k = nzb, nzt
+                       sums_wspts_ws_l(k,tn) = sums_wspts_ws_l(k,tn)          &
+                        + ( flux_t(k) + diss_t(k) )                           &
+                        *   weight_substep(intermediate_timestep_count)
+                    ENDDO
+                 CASE ( 'sa' )
+                    DO  k = nzb, nzt
+                       sums_wssas_ws_l(k,tn) = sums_wssas_ws_l(k,tn)          &
+                        + ( flux_t(k) + diss_t(k) )                           &
+                        *   weight_substep(intermediate_timestep_count)
+                    ENDDO
+                 CASE ( 'q' )
+                    DO  k = nzb, nzt
+                       sums_wsqs_ws_l(k,tn) = sums_wsqs_ws_l(k,tn)            &
+                        + ( flux_t(k) + diss_t(k) )                           &
+                        *   weight_substep(intermediate_timestep_count)
+                    ENDDO
+                 CASE ( 'qr' )
+                    DO  k = nzb, nzt
+                       sums_wsqrs_ws_l(k,tn) = sums_wsqrs_ws_l(k,tn)          &
+                        + ( flux_t(k) + diss_t(k) )                           &
+                        *   weight_substep(intermediate_timestep_count)
+                    ENDDO
+                 CASE ( 'nr' )
+                    DO  k = nzb, nzt
+                       sums_wsnrs_ws_l(k,tn) = sums_wsnrs_ws_l(k,tn)          &
+                        + ( flux_t(k) + diss_t(k) )                           &
+                        *   weight_substep(intermediate_timestep_count)
+                    ENDDO
+
+              END SELECT
 
          ENDDO
