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/sum_up_3d_data.f90

    r3943 r3994  
    2525! -----------------
    2626! $Id$
     27! output of turbulence intensity added
     28!
     29! 3943 2019-05-02 09:50:41Z maronga
    2730! Added output of qsws_av for green roofs.
    2831!
     
    283286    USE cpulog,                                                                &
    284287        ONLY:  cpu_log, log_point
     288
     289    USE diagnostic_output_quantities_mod,                                      &
     290        ONLY:  ti, ti_av
    285291
    286292    USE indices,                                                               &
     
    472478                ENDIF
    473479                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
    474486
    475487             CASE ( 'u' )
     
    10191031             ENDIF
    10201032
     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
    10211044          CASE ( 'u' )
    10221045             IF ( ALLOCATED( u_av ) ) THEN
Note: See TracChangeset for help on using the changeset viewer.