Changeset 4741 for palm/trunk/SOURCE/calc_mean_profile.f90

Timestamp:

Oct 14, 2020 2:32:50 PM (4 years ago)

Author:

suehring

Message:

radiation: Add option to force calculation of horizontal mean profiles independent on data output

File:

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

palm/trunk/SOURCE/calc_mean_profile.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
-! 
+!
+!
 ! Former revisions:
 ! -----------------
 ! $Id$
+! Add option to force calculation of horizontal mean profiles independent on data output
+! 
+! 4542 2020-05-19 15:45:12Z raasch
 ! file re-formatted to follow the PALM coding standard
 !
@@ -66 +69 @@
 !> @todo Missing subroutine description.
 !--------------------------------------------------------------------------------------------------!
-    SUBROUTINE calc_mean_profile( var, pr )
+    SUBROUTINE calc_mean_profile( var, pr, force_calc )
 
        USE control_parameters,                                                                     &
@@ -84 +87 @@
        IMPLICIT NONE
 
+       LOGICAL, OPTIONAL ::  force_calc             !< passed flag used to force calculation of mean profiles independend on data output
+       LOGICAL           ::  force_calc_l = .FALSE. !< control flag
+
        INTEGER(iwp) ::  i                  !<
        INTEGER(iwp) ::  j                  !<
@@ -99 +105 @@
 !--    results with jobs which are calling flow_statistics at different time intervals. At
 !--    initialization, intermediate_timestep_count = 0 is considered as well.
+!--    Note, calc_mean_profile is also called from the radiation. Especially during the spinup when
+!--    no data output is called, the values in hom do not necessarily reflect horizontal mean.
+!--    The same is also for nested simulations during the spinup, where hom is initialized on a
+!--    ealier point in time than the initialization via the parent domain takes place, meaning
+!--    that hom does not necessarily reflect the horizontal mean during the spinup, too.
+!--    In order to force the computation of horizontal mean profiles independent on data output,
+!--    i.e. independent on flow_statistics, add a special control flag.
+       IF ( PRESENT( force_calc ) )  THEN
+          IF ( force_calc )  force_calc_l = .TRUE.
+       ENDIF
 
-       IF ( .NOT. flow_statistics_called  .AND.  intermediate_timestep_count <= 1 )  THEN
+       IF ( ( .NOT. flow_statistics_called  .AND.  intermediate_timestep_count <= 1 ) .OR.         &
+            force_calc_l )  THEN
 
 !