Changeset 1498 for palm/trunk/SOURCE/flow_statistics.f90

Timestamp:

Dec 3, 2014 2:09:51 PM (9 years ago)

Author:

suehring

Message:

Comments added

File:

• palm/trunk/SOURCE/flow_statistics.f90 (modified) (3 diffs)

palm/trunk/SOURCE/flow_statistics.f90

@@ -21 +21 @@
 ! Current revisions:
 ! -----------------
-! 
+! Comments added
 !
 ! Former revisions:
@@ -219 +219 @@
 
 !
+!--    When calcuating horizontally-averaged total (resolved- plus subgrid-
+!--    scale) vertical fluxes and velocity variances by using commonly-
+!--    applied Reynolds-based methods ( e.g. <w'pt'> = (w-<w>)*(pt-<pt>) )
+!--    in combination with the 5th order advection scheme, pronounced 
+!--    artificial kinks could be observed in the vertical profiles near the 
+!--    surface. Please note: these kinks were not related to the model truth, 
+!--    i.e. these kinks are just related to an evaluation problem.   
+!--    In order avoid these kinks, vertical fluxes and horizontal as well 
+!--    vertical velocity variances are calculated directly within the advection
+!--    routines, according to the numerical discretization, to evaluate the 
+!--    statistical quantities as they will appear within the prognostic 
+!--    equations.
 !--    Copy the turbulent quantities, evaluated in the advection routines to 
-!--    the local array sums_l() for further computations
+!--    the local array sums_l() for further computations.
        IF ( ws_scheme_mom .AND. sr == 0 )  THEN
 
@@ -1517 +1529 @@
 
 !
+!--    When calcuating horizontally-averaged total (resolved- plus subgrid-
+!--    scale) vertical fluxes and velocity variances by using commonly-
+!--    applied Reynolds-based methods ( e.g. <w'pt'> = (w-<w>)*(pt-<pt>) )
+!--    in combination with the 5th order advection scheme, pronounced 
+!--    artificial kinks could be observed in the vertical profiles near the 
+!--    surface. Please note: these kinks were not related to the model truth, 
+!--    i.e. these kinks are just related to an evaluation problem.   
+!--    In order avoid these kinks, vertical fluxes and horizontal as well 
+!--    vertical velocity variances are calculated directly within the advection
+!--    routines, according to the numerical discretization, to evaluate the 
+!--    statistical quantities as they will appear within the prognostic 
+!--    equations. 
 !--    Copy the turbulent quantities, evaluated in the advection routines to
-!--    the local array sums_l() for further computations
+!--    the local array sums_l() for further computations.
        IF ( ws_scheme_mom .AND. sr == 0 )  THEN