Ignore:
Timestamp:
Jan 18, 2011 4:19:48 PM (11 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:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/flow_statistics.f90

    r668 r673  
    44! Current revisions:
    55! -----------------
    6 ! When advection is computed with ws-scheme, turbulent fluxes are already
    7 ! computed in the respective advection routines and buffered in arrays
    8 ! sums_xx_ws_l(). This is due to a consistent treatment of statistics with the
    9 ! numerics and to avoid unphysical kinks near the surface.
    10 ! So some if requests has to be done to dicern between fluxes from ws-scheme
    11 ! other advection schemes.
    12 ! Furthermore the computation of z_i is only done if the heat flux exceeds a
    13 ! minimum value. This affects only simulations of a neutral boundary layer and
    14 ! is due to reasons of computations in the advection scheme.
    15 !
     6! Top bc for the horizontal velocity variances added for ocean runs.
     7! Setting the corresponding bottom bc moved to advec_ws.
    168!
    179! Former revisions:
     
    157149       IF ( ws_scheme_mom )  THEN
    158150!       
    159 !--       Boundary condition for u'u' and v'v', because below the surface no
    160 !--       computation for these quantities is done.
    161           DO  i = nxl, nxr
    162              DO  j =  nys, nyn
    163                 sums_us2_ws_l(nzb_u_inner(j,i),sr) =                          &
    164                     sums_us2_ws_l(nzb_u_inner(j,i)+1,sr)
    165                 sums_vs2_ws_l(nzb_v_inner(j,i),sr) =                          & 
    166                     sums_vs2_ws_l(nzb_v_inner(j,i)+1,sr)
    167              ENDDO
    168           ENDDO
     151!--       According to the Neumann bc for the horizontal velocity components,
     152!--       the corresponding fluxes has to satisfiy the same bc.
     153          IF ( ocean )  THEN
     154             sums_us2_ws_l(nzt+1,sr) = sums_us2_ws_l(nzt,sr)
     155             sums_vs2_ws_l(nzt+1,sr) = sums_vs2_ws_l(nzt,sr)   
     156          ENDIF
    169157!         
    170158!--       Swap the turbulent quantities evaluated in advec_ws.
     
    184172       ENDIF
    185173       IF ( ws_scheme_sca )  THEN
    186           sums_l(:,17,0) = sums_wspts_ws_l(:,sr)      ! w*pts* from advec_s_ws
     174          sums_l(:,17,0) = sums_wspts_ws_l(:,sr)      ! w*pt* from advec_s_ws
    187175          IF ( ocean ) sums_l(:,66,0) = sums_wssas_ws_l(:,sr) ! w*sa*
    188176          IF ( humidity  .OR.  passive_scalar ) sums_l(:,49,0) =              &
Note: See TracChangeset for help on using the changeset viewer.