Ignore:
Timestamp:
Jun 18, 2019 10:32:41 AM (5 years ago)
Author:
suehring
Message:

diagnostic output: Modularize diagnostic output, rename subroutines; formatting adjustments; allocate arrays only when required; add output of uu, vv, ww to enable variance calculation via temporal EC method; radiation: bugfix in masked data output; flow_statistics: Correct conversion to kinematic vertical scalar fluxes in case of pw-scheme and statistic regions

File:
1 edited

Legend:

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

    r3994 r4039  
    2525! -----------------
    2626! $Id$
     27! modularize diagnostic output
     28!
     29! 3994 2019-05-22 18:08:09Z suehring
    2730! output of turbulence intensity added
    2831!
     
    324327    USE cpulog,                                                                &
    325328        ONLY:  cpu_log, log_point
    326 
    327     USE diagnostic_output_quantities_mod,                                      &
    328         ONLY:  ti, ti_av
    329329
    330330    USE indices,                                                               &
     
    11711171                   IF ( ibc_uv_b == 0 )  local_pf(:,:,nzb) = 0.0_wp
    11721172                ENDIF
    1173 
     1173               
    11741174             CASE ( 'us*_xy' )        ! 2d-array
    11751175                IF ( av == 0 )  THEN
     
    12661266                two_d = .TRUE.
    12671267                level_z(nzb+1) = zu(nzb+1)
    1268                
    1269              CASE ( 'ti_xy', 'ti_xz', 'ti_yz' )
    1270                 IF ( av == 0 )  THEN
    1271                    to_be_resorted => ti
    1272                 ELSE
    1273                    IF ( .NOT. ALLOCATED( ti_av ) ) THEN
    1274                       ALLOCATE( ti_av(nzb:nzt+1,nys:nyn,nxl:nxr) )
    1275                       ti_av = REAL( fill_value, KIND = wp )
    1276                    ENDIF
    1277                    to_be_resorted => ti_av
    1278                 ENDIF
    1279                 IF ( mode == 'xy' )  level_z = zu
    12801268
    12811269             CASE ( 'w_xy', 'w_xz', 'w_yz' )
Note: See TracChangeset for help on using the changeset viewer.