Ignore:
Timestamp:
May 27, 2016 2:35:57 PM (5 years ago)
Author:
raasch
Message:

bugfixes for calculating run control quantities, bugfix for calculating pressure with fft-method in case of Neumann conditions both at bottom and top, steering of pres modified, ocean mode now uses initial density profile as reference in the buoyancy term

File:
1 edited

Legend:

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

    r1854 r1918  
    1919! Current revisions:
    2020! -----------------
    21 !
     21! in case of Wicker-Skamarock scheme, calculate disturbance kinetic energy here,
     22! if flow_statistics is called before the first initial time step
    2223!
    2324! Former revisions:
     
    630631!--    Computation of statistics when ws-scheme is not used. Else these
    631632!--    quantities are evaluated in the advection routines.
    632        IF ( .NOT. ws_scheme_mom .OR. sr /= 0 )  THEN
     633       IF ( .NOT. ws_scheme_mom .OR. sr /= 0 .OR. simulated_time == 0.0_wp )   &
     634       THEN
    633635          !$OMP DO
    634636          DO  i = nxl, nxr
     
    16231625
    16241626!
    1625 !-- If required, sum up horizontal averages for subsequent time averaging
    1626     IF ( do_sum )  THEN
     1627!-- If required, sum up horizontal averages for subsequent time averaging.
     1628!-- Do not sum, if flow statistics is called before the first initial time step.
     1629    IF ( do_sum  .AND.  simulated_time /= 0.0_wp )  THEN
    16271630       IF ( average_count_pr == 0 )  hom_sum = 0.0_wp
    16281631       hom_sum = hom_sum + hom(:,1,:,:)
     
    22302233!--    Computation of statistics when ws-scheme is not used. Else these
    22312234!--    quantities are evaluated in the advection routines.
    2232        IF ( .NOT. ws_scheme_mom  .OR.  sr /= 0 )  THEN
     2235       IF ( .NOT. ws_scheme_mom .OR. sr /= 0 .OR. simulated_time == 0.0_wp )   &
     2236       THEN
    22332237
    22342238          !$OMP DO
     
    35663570!
    35673571!-- If required, sum up horizontal averages for subsequent time averaging
    3568     IF ( do_sum )  THEN
     3572!-- Do not sum, if flow statistics is called before the first initial time step.
     3573    IF ( do_sum  .AND.  simulated_time /= 0.0_wp )  THEN
    35693574       IF ( average_count_pr == 0 )  hom_sum = 0.0_wp
    35703575       hom_sum = hom_sum + hom(:,1,:,:)
Note: See TracChangeset for help on using the changeset viewer.