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/sum_up_3d_data.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!
     
    286289    USE cpulog,                                                                &
    287290        ONLY:  cpu_log, log_point
    288 
    289     USE diagnostic_output_quantities_mod,                                      &
    290         ONLY:  ti, ti_av
    291291
    292292    USE indices,                                                               &
     
    478478                ENDIF
    479479                tsurf_av = 0.0_wp
    480 
    481              CASE ( 'ti' )
    482                 IF ( .NOT. ALLOCATED( ti_av ) )  THEN
    483                    ALLOCATE( ti_av(nzb:nzt+1,nys:nyn,nxl:nxr) )
    484                 ENDIF
    485                 ti_av = 0.0_wp
    486480
    487481             CASE ( 'u' )
     
    10311025             ENDIF
    10321026
    1033           CASE ( 'ti' )
    1034              IF ( ALLOCATED( ti_av ) ) THEN
    1035                 DO  i = nxl, nxr
    1036                    DO  j = nys, nyn
    1037                       DO  k = nzb, nzt+1
    1038                          ti_av(k,j,i) = ti_av(k,j,i) + ti(k,j,i)
    1039                       ENDDO
    1040                    ENDDO
    1041                 ENDDO
    1042              ENDIF
    1043 
    10441027          CASE ( 'u' )
    10451028             IF ( ALLOCATED( u_av ) ) THEN
Note: See TracChangeset for help on using the changeset viewer.