Changeset 3994 for palm/trunk/SOURCE/data_output_2d.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/data_output_2d.f90 (modified) (3 diffs)

palm/trunk/SOURCE/data_output_2d.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! output of turbulence intensity added
+! 
+! 3987 2019-05-22 09:52:13Z kanani
 ! Introduce alternative switch for debug output during timestepping
 ! 
@@ -321 +324 @@
     USE cpulog,                                                                &
         ONLY:  cpu_log, log_point
+
+    USE diagnostic_output_quantities_mod,                                      &
+        ONLY:  ti, ti_av
 
     USE indices,                                                               &
@@ -1261 +1267 @@
                 level_z(nzb+1) = zu(nzb+1)
                 
-                
-                
+             CASE ( 'ti_xy', 'ti_xz', 'ti_yz' )
+                IF ( av == 0 )  THEN
+                   to_be_resorted => ti
+                ELSE
+                   IF ( .NOT. ALLOCATED( ti_av ) ) THEN
+                      ALLOCATE( ti_av(nzb:nzt+1,nys:nyn,nxl:nxr) )
+                      ti_av = REAL( fill_value, KIND = wp )
+                   ENDIF
+                   to_be_resorted => ti_av
+                ENDIF
+                IF ( mode == 'xy' )  level_z = zu
+
              CASE ( 'w_xy', 'w_xz', 'w_yz' )
                 flag_nr = 3