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

Timestamp:

Jan 18, 2011 4:19:48 PM (14 years ago)

Author:

suehring

Message:

Right computation of the pressure using Runge-Kutta weighting coefficients. Consideration of the pressure gradient during the time integration removed. Removed bugfix concerning velocity variances.

File:

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

palm/trunk/SOURCE/flow_statistics.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
-! When advection is computed with ws-scheme, turbulent fluxes are already
-! computed in the respective advection routines and buffered in arrays
-! sums_xx_ws_l(). This is due to a consistent treatment of statistics with the
-! numerics and to avoid unphysical kinks near the surface. 
-! So some if requests has to be done to dicern between fluxes from ws-scheme
-! other advection schemes. 
-! Furthermore the computation of z_i is only done if the heat flux exceeds a
-! minimum value. This affects only simulations of a neutral boundary layer and
-! is due to reasons of computations in the advection scheme.
-!
+! Top bc for the horizontal velocity variances added for ocean runs.
+! Setting the corresponding bottom bc moved to advec_ws.
 !
 ! Former revisions:
@@ -157 +149 @@
        IF ( ws_scheme_mom )  THEN
 !       
-!--       Boundary condition for u'u' and v'v', because below the surface no
-!--       computation for these quantities is done.
-          DO  i = nxl, nxr
-             DO  j =  nys, nyn
-                sums_us2_ws_l(nzb_u_inner(j,i),sr) =                          & 
-                    sums_us2_ws_l(nzb_u_inner(j,i)+1,sr)
-                sums_vs2_ws_l(nzb_v_inner(j,i),sr) =                          &  
-                    sums_vs2_ws_l(nzb_v_inner(j,i)+1,sr)
-             ENDDO
-          ENDDO
+!--       According to the Neumann bc for the horizontal velocity components,
+!--       the corresponding fluxes has to satisfiy the same bc.
+          IF ( ocean )  THEN
+             sums_us2_ws_l(nzt+1,sr) = sums_us2_ws_l(nzt,sr)
+             sums_vs2_ws_l(nzt+1,sr) = sums_vs2_ws_l(nzt,sr)    
+          ENDIF
 !          
 !--       Swap the turbulent quantities evaluated in advec_ws.
@@ -184 +172 @@
        ENDIF
        IF ( ws_scheme_sca )  THEN
-          sums_l(:,17,0) = sums_wspts_ws_l(:,sr)      ! w*pts* from advec_s_ws
+          sums_l(:,17,0) = sums_wspts_ws_l(:,sr)      ! w*pt* from advec_s_ws
           IF ( ocean ) sums_l(:,66,0) = sums_wssas_ws_l(:,sr) ! w*sa*
           IF ( humidity  .OR.  passive_scalar ) sums_l(:,49,0) =              &