Ignore:
Timestamp:
May 22, 2019 6:08:09 PM (2 years ago)
Author:
suehring
Message:

new module for diagnostic output quantities added + output of turbulence intensity

File:
1 edited

Legend:

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

    r3988 r3994  
    2525! -----------------
    2626! $Id$
     27! output of turbulence intensity added
     28!
     29! 3988 2019-05-22 11:32:37Z kanani
    2730! Implement steerable output interval for virtual measurements
    2831!
     
    561564    USE date_and_time_mod,                                                                         &
    562565        ONLY:  calc_date_and_time, hour_call_emis, hour_of_year
     566
     567    USE diagnostic_output_quantities_mod,                                                          &
     568        ONLY:  diagnostic_output_quantities_calculate,                                             &
     569               timestep_number_at_prev_calc
    563570
    564571    USE flight_mod,                                                                                &
     
    16551662          time_do_sla = time_do_sla + dt_3d
    16561663          IF ( time_do_sla >= dt_averaging_input )  THEN
     1664             IF ( current_timestep_number > timestep_number_at_prev_calc )                         &
     1665                CALL diagnostic_output_quantities_calculate
     1666
    16571667             CALL sum_up_3d_data
    16581668             average_count_3d = average_count_3d + 1
     
    17331743!--    2d-data output (cross-sections)
    17341744       IF ( time_do2d_xy >= dt_do2d_xy  .AND.  time_since_reference_point >= skip_time_do2d_xy )  THEN
     1745          IF ( current_timestep_number > timestep_number_at_prev_calc )                            &
     1746             CALL diagnostic_output_quantities_calculate
     1747
    17351748          CALL data_output_2d( 'xy', 0 )
    17361749          time_do2d_xy = MOD( time_do2d_xy, MAX( dt_do2d_xy, dt_3d ) )
    17371750       ENDIF
    17381751       IF ( time_do2d_xz >= dt_do2d_xz  .AND.  time_since_reference_point >= skip_time_do2d_xz )  THEN
     1752          IF ( current_timestep_number > timestep_number_at_prev_calc )                            &
     1753
     1754             CALL diagnostic_output_quantities_calculate
    17391755          CALL data_output_2d( 'xz', 0 )
    17401756          time_do2d_xz = MOD( time_do2d_xz, MAX( dt_do2d_xz, dt_3d ) )
    17411757       ENDIF
    17421758       IF ( time_do2d_yz >= dt_do2d_yz  .AND.  time_since_reference_point >= skip_time_do2d_yz )  THEN
     1759          IF ( current_timestep_number > timestep_number_at_prev_calc )                            &
     1760             CALL diagnostic_output_quantities_calculate
     1761
    17431762          CALL data_output_2d( 'yz', 0 )
    17441763          time_do2d_yz = MOD( time_do2d_yz, MAX( dt_do2d_yz, dt_3d ) )
     
    17481767!--    3d-data output (volume data)
    17491768       IF ( time_do3d >= dt_do3d  .AND.  time_since_reference_point >= skip_time_do3d )  THEN
     1769          IF ( current_timestep_number > timestep_number_at_prev_calc )                            &
     1770             CALL diagnostic_output_quantities_calculate
     1771
    17501772          CALL data_output_3d( 0 )
    17511773          time_do3d = MOD( time_do3d, MAX( dt_do3d, dt_3d ) )
     
    17571779          IF ( time_domask(mid) >= dt_domask(mid)                                                  &
    17581780               .AND.  time_since_reference_point >= skip_time_domask(mid) )  THEN
     1781             IF ( current_timestep_number > timestep_number_at_prev_calc )                         &
     1782                CALL diagnostic_output_quantities_calculate
     1783
    17591784             CALL data_output_mask( 0 )
    17601785             time_domask(mid) = MOD( time_domask(mid), MAX( dt_domask(mid), dt_3d ) )
Note: See TracChangeset for help on using the changeset viewer.