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

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:

• palm/trunk/SOURCE/average_3d_data.f90 (modified) (3 diffs)

palm/trunk/SOURCE/average_3d_data.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! output of turbulence intensity added
+! 
+! 3933 2019-04-25 12:33:20Z kanani
 ! Bugfix in CASE theta_2m*, removal of redundant code
 ! 
@@ -188 +191 @@
         ONLY:  cpu_log, log_point
 
+    USE diagnostic_output_quantities_mod,                                      &
+        ONLY:  ti_av
+
     USE indices,                                                               &
         ONLY:  nbgp, nxl, nxlg, nxr, nxrg, nyn, nyng, nys, nysg, nzb, nzt
@@ -515 +521 @@
              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) / REAL( average_count_3d, KIND=wp )
+                      ENDDO
+                   ENDDO
+                ENDDO
+             ENDIF
+
           CASE ( 'w' )
              IF ( ALLOCATED( w_av ) ) THEN