Changeset 4741 for palm/trunk/SOURCE

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

Location:

palm/trunk/SOURCE

Files:

• calc_mean_profile.f90 (modified) (4 diffs)
• radiation_model_mod.f90 (modified) (3 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
 
 !

palm/trunk/SOURCE/radiation_model_mod.f90

@@ -28 +28 @@
 ! -----------------
 ! $Id$
+! Add option to force calculation of horizontal mean profiles independent on data output
+!
+! 4717 2020-09-30 22:27:40Z pavelkrc
 ! Fixes and optimizations of OpenMP parallelization, formatting of OpenMP
 ! directives (J. Resler)
@@ -4240 +4243 @@
 !
 !--       Calculate mean pt profile.
-          CALL calc_mean_profile( pt, 4 )
+          CALL calc_mean_profile( pt, 4, .TRUE. )
           pt_av = hom(:, 1, 4, 0)
 
           IF ( humidity )  THEN
-             CALL calc_mean_profile( q, 41 )
+             CALL calc_mean_profile( q, 41, .TRUE. )
              q_av  = hom(:, 1, 41, 0)
           ENDIF
@@ -4253 +4256 @@
           IF ( bulk_cloud_model )  THEN
 
-             CALL calc_mean_profile( ql, 54 )
+             CALL calc_mean_profile( ql, 54, .TRUE. )
              ! average ql is now in hom(:, 1, 54, 0)
              ql_av = hom(:, 1, 54, 0)