Changeset 4071 for palm

Timestamp:

Jul 3, 2019 8:02:00 PM (5 years ago)

Author:

suehring

Message:

Bugfix in synthetic turbulence generator in non-nested cases (offline or self nesting) when no cyclic fill is used: in order to get correct initial inflow profiles at the left boundary the array mean_inflow_profiles need to be initialized accordingly.

File:

• palm/trunk/SOURCE/synthetic_turbulence_generator_mod.f90 (modified) (5 diffs)

palm/trunk/SOURCE/synthetic_turbulence_generator_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfix, initialize mean_inflow_profiles in case turbulence and inflow
+! information is not read from file. 
+! 
+! 4022 2019-06-12 11:52:39Z suehring
 ! Several bugfixes and improvements
 ! - revise bias correction of the imposed perturbations (correction via volume
@@ -224 +228 @@
                mean_inflow_profiles,                                           &
                q,                                                              &
-               pt,                                                             & 
+               q_init,                                                         &
+               pt,                                                             &
+               pt_init,                                                        &
                u,                                                              &
                u_init,                                                         &
@@ -246 +252 @@
                dt_3d,                                                          &
                e_init,                                                         &
+               humidity,                                                       &
                initializing_actions,                                           &
                intermediate_timestep_count,                                    &
@@ -252 +259 @@
                message_string,                                                 &
                nesting_offline,                                                &
+               neutral,                                                        &
                num_mean_inflow_profiles,                                       &
                rans_mode,                                                      &
@@ -757 +765 @@
     CALL RANDOM_SEED(put=seed)
 !
-!-- Allocate required arrays
-!-- mean_inflow profiles must not be allocated in offline nesting
+!-- Allocate required arrays.
+!-- In case no offline nesting or self nesting is used, the arrary
+!-- mean_inflow profiles is required. Check if it is already allocated, else
+!-- allocate and initialize it appropriately. Note, in case turbulence and 
+!-- inflow information is read from file, mean_inflow_profiles is already
+!-- allocated and initialized appropriately. 
     IF ( .NOT. nesting_offline  .AND.  .NOT.  child_domain )  THEN 
-       IF ( .NOT. ALLOCATED( mean_inflow_profiles ) )                          &
+       IF ( .NOT. ALLOCATED( mean_inflow_profiles ) )  THEN
           ALLOCATE( mean_inflow_profiles(nzb:nzt+1,1:num_mean_inflow_profiles) )
+          mean_inflow_profiles = 0.0_wp
+          mean_inflow_profiles(:,1) = u_init
+          mean_inflow_profiles(:,2) = v_init
+!
+!--       Even though potential temperature and humidity are not perturbed,
+!--       they need to be initialized appropriately. 
+          IF ( .NOT. neutral )                                                 &
+             mean_inflow_profiles(:,4) = pt_init
+          IF ( humidity )                                                      &
+             mean_inflow_profiles(:,6) = q_init       
+       ENDIF   
     ENDIF