Changeset 3994 for palm/trunk/SOURCE/sum_up_3d_data.f90

Timestamp:

May 22, 2019 6:08:09 PM (5 years ago)

Author:

suehring

Message:

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

File:

• palm/trunk/SOURCE/sum_up_3d_data.f90 (modified) (4 diffs)

palm/trunk/SOURCE/sum_up_3d_data.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! output of turbulence intensity added
+! 
+! 3943 2019-05-02 09:50:41Z maronga
 ! Added output of qsws_av for green roofs.
 ! 
@@ -283 +286 @@
     USE cpulog,                                                                &
         ONLY:  cpu_log, log_point
+
+    USE diagnostic_output_quantities_mod,                                      &
+        ONLY:  ti, ti_av
 
     USE indices,                                                               &
@@ -472 +478 @@
                 ENDIF
                 tsurf_av = 0.0_wp
+
+             CASE ( 'ti' )
+                IF ( .NOT. ALLOCATED( ti_av ) )  THEN
+                   ALLOCATE( ti_av(nzb:nzt+1,nys:nyn,nxl:nxr) )
+                ENDIF
+                ti_av = 0.0_wp
 
              CASE ( 'u' )
@@ -1019 +1031 @@
              ENDIF
 
+          CASE ( 'ti' )
+             IF ( ALLOCATED( ti_av ) ) THEN
+                DO  i = nxl, nxr
+                   DO  j = nys, nyn
+                      DO  k = nzb, nzt+1
+                         ti_av(k,j,i) = ti_av(k,j,i) + ti(k,j,i)
+                      ENDDO
+                   ENDDO
+                ENDDO
+             ENDIF
+
           CASE ( 'u' )
              IF ( ALLOCATED( u_av ) ) THEN