Changeset 3937 for palm/trunk/SOURCE/init_3d_model.f90

Timestamp:

Apr 29, 2019 3:09:07 PM (5 years ago)

Author:

suehring

Message:

Move initialization of STG behind initialization of offline nesting; Bugfix in STG in case of very early restart; calculation of scaling parameters used for parametrization of synthetic turbulence profiles improved; in offline nesting, set boundary value at upper-left and upper-south grid points for u- and v-component, respectively

File:

• palm/trunk/SOURCE/init_3d_model.f90 (modified) (7 diffs)

palm/trunk/SOURCE/init_3d_model.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Move initialization of synthetic turbulence generator behind initialization
+! of offline nesting. Remove call for stg_adjust, as this is now already done
+! in stg_init. 
+! 
+! 3900 2019-04-16 15:17:43Z suehring
 ! Fix problem with LOD = 2 initialization
 ! 
@@ -666 +671 @@
 
     USE synthetic_turbulence_generator_mod,                                    &
-        ONLY:  parametrize_inflow_turbulence, stg_adjust, stg_init,            &
-               use_syn_turb_gen
+        ONLY:  parametrize_inflow_turbulence, stg_init, use_syn_turb_gen
 
     USE surface_layer_fluxes_mod,                                              &
@@ -1305 +1309 @@
 !--       fluxes, etc. 
           CALL init_surfaces
-!
-!--       Initialize turbulence generator
-          IF( use_syn_turb_gen )  CALL stg_init
 
           CALL location_message( 'initializing with INIFOR', 'finished' )
@@ -1375 +1376 @@
 !--       fluxes, etc. 
           CALL init_surfaces
-!
-!--       Initialize synthetic turbulence generator if required 
-          IF( use_syn_turb_gen )  CALL stg_init
 
           CALL location_message( 'initializing with 1D model profiles', 'finished' )
@@ -1446 +1444 @@
 !--       fluxes, etc. 
           CALL init_surfaces
-!
-!--       Initialize synthetic turbulence generator if required 
-          IF( use_syn_turb_gen )  CALL stg_init
           
           CALL location_message( 'initializing with constant profiles', 'finished' )
@@ -1868 +1863 @@
        ENDIF
        IF ( passive_scalar )  ts_m  = 0.0_wp
-!
-!--    Initialize synthetic turbulence generator in case of restart.
-       IF ( TRIM( initializing_actions ) == 'read_restart_data'  .AND.         &
-            use_syn_turb_gen )  CALL stg_init
 
        CALL location_message( 'initializing in case of restart / cyclic_fill', 'finished' )
@@ -2315 +2306 @@
     IF ( agents_active )  CALL mas_init
 !
-!-- In case the synthetic turbulence generator does not have any information
-!-- about the inflow turbulence, these information will be parametrized 
-!-- depending on the initial atmospheric conditions and surface properties.
-!-- Please note, within pre-determined time intervals these turbulence 
-!-- information can be updated if desired. 
-    IF ( use_syn_turb_gen  .AND.  parametrize_inflow_turbulence )  THEN
-       CALL stg_adjust
-    ENDIF
+!-- Initialization of synthetic turbulence generator
+    IF ( use_syn_turb_gen )  CALL stg_init
 !
 !-- Initializing actions for all other modules