Changeset 3004 for palm/trunk/SOURCE/flow_statistics.f90

Timestamp:

Apr 27, 2018 12:33:25 PM (6 years ago)

Author:

Giersch

Message:

precipitation_rate removed, further allocation checks for data output of averaged quantities implemented, double CALL of flow_statistics at the beginning of time_integration removed, further minor bugfixes, comments added

File:

• palm/trunk/SOURCE/flow_statistics.f90 (modified) (2 diffs)

palm/trunk/SOURCE/flow_statistics.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Comments related to the calculation of the inversion height expanded
+! 
+! 3003 2018-04-23 10:22:58Z Giersch
 ! The inversion height will not be calcuated before the first timestep in 
 ! case of restarts.
@@ -2007 +2010 @@
 !--    The corresponding height is assumed as the boundary layer height, if it
 !--    is less than 1.5 times the height where the heat flux becomes negative
-!--    (positive) for the first time.
+!--    (positive) for the first time. Attention: the resolved vertical sensible 
+!--    heat flux (hom(:,1,17,sr) = w*pt*) is not known at the beginning because
+!--    the calculation happens in advec_s_ws which is called after 
+!--    flow_statistics. Therefore z_i is directly taken from restart data at 
+!--    the beginning of restart runs.  
        IF ( TRIM( initializing_actions ) /= 'read_restart_data' .OR.           &
             simulated_time_at_begin /= simulated_time ) THEN