Changeset 3347 for palm/trunk/SOURCE/time_integration.f90

Timestamp:

Oct 15, 2018 2:21:08 PM (6 years ago)

Author:

suehring

Message:

Offline nesting revised and separated from large_scale_forcing_mod; Time-dependent synthetic turbulence generator; bugfixes in USM/LSM radiation model and init_pegrid

File:

• palm/trunk/SOURCE/time_integration.f90 (modified) (8 diffs)

palm/trunk/SOURCE/time_integration.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! - offline nesting separated from large-scale forcing module
+! - changes for synthetic turbulence generator
+! 
+! 3343 2018-10-15 10:38:52Z suehring
 ! - Formatting, clean-up, comments (kanani)
 ! - Added CALL to chem_emissions_setup (Russo)
@@ -459 +463 @@
 
     USE lsf_nudging_mod,                                                       &
-        ONLY:  calc_tnudge, ls_forcing_surf, ls_forcing_vert, nudge_ref,       &
-               lsf_nesting_offline, lsf_nesting_offline_mass_conservation
+        ONLY:  calc_tnudge, ls_forcing_surf, ls_forcing_vert, nudge_ref
 
     USE multi_agent_system_mod,                                                &
         ONLY:  agents_active, multi_agent_system
-
+        
+    USE nesting_offl_mod,                                                      &
+        ONLY:  nesting_offl_bc, nesting_offl_mass_conservation
+        
     USE netcdf_data_input_mod,                                                 &
-        ONLY:  chem_emis, chem_emis_att, nest_offl, netcdf_data_input_lsf
+        ONLY:  chem_emis, chem_emis_att, nest_offl,                            &
+               netcdf_data_input_offline_nesting
 
     USE ocean_mod,                                                             &
@@ -515 +522 @@
 
     USE synthetic_turbulence_generator_mod,                                    &
-        ONLY:  stg_main, use_syn_turb_gen
+        ONLY:  dt_stg_call, dt_stg_adjust, parametrize_inflow_turbulence,      &
+               stg_adjust, stg_main, time_stg_adjust, time_stg_call,           &
+               use_syn_turb_gen
 
     USE user_actions_mod,                                                      &
@@ -638 +647 @@
        IF ( nesting_offline )  THEN
           IF ( nest_offl%time(nest_offl%tind_p) <= time_since_reference_point )&
-             CALL netcdf_data_input_lsf
+             CALL netcdf_data_input_offline_nesting
        ENDIF
 
@@ -918 +927 @@
 !--       Impose a turbulent inflow using the recycling method
           IF ( turbulent_inflow )  CALL  inflow_turbulence
-
-!
-!--       Impose a turbulent inflow using synthetic generated turbulence
-          IF ( use_syn_turb_gen )  CALL  stg_main
 
 !
@@ -960 +965 @@
 !--       boundaries. 
           IF ( nesting_offline  .AND.  intermediate_timestep_count ==          &
-                                       intermediate_timestep_count_max  )  THEN
-             CALL lsf_nesting_offline
-!
-!--          Moreover, ensure mass conservation
-             CALL lsf_nesting_offline_mass_conservation
-          ENDIF
-
+                                       intermediate_timestep_count_max  )      &
+             CALL nesting_offl_bc
+!
+!--       Impose a turbulent inflow using synthetic generated turbulence, 
+!--       only once per time step.
+          IF ( use_syn_turb_gen  .AND.  time_stg_call >= dt_stg_call  .AND.    &
+             intermediate_timestep_count == intermediate_timestep_count_max )  THEN! &
+             CALL stg_main
+          ENDIF
+!
+!--       Ensure mass conservation. This need to be done after imposing
+!--       synthetic turbulence and top boundary condition for pressure is set to
+!--       Neumann conditions. 
+!--       Is this also required in case of Dirichlet?
+          IF ( nesting_offline )  CALL nesting_offl_mass_conservation
 !
 !--       Reduce the velocity divergence via the equation for perturbation
@@ -1213 +1226 @@
        time_dopr_listing          = time_dopr_listing        + dt_3d
        time_run_control   = time_run_control + dt_3d
+!
+!--    In case of synthetic turbulence generation and parametrized turbulence
+!--    information, update the time counter and if required, adjust the 
+!--    STG to new atmospheric conditions.
+       IF ( use_syn_turb_gen  )  THEN
+          IF ( parametrize_inflow_turbulence )  THEN
+             time_stg_adjust = time_stg_adjust + dt_3d
+             IF ( time_stg_adjust >= dt_stg_adjust )  CALL stg_adjust
+          ENDIF
+          time_stg_call = time_stg_call + dt_3d
+       ENDIF
 
 !
@@ -1472 +1496 @@
 
     ENDDO   ! time loop
-
+    
+!
 !-- Vertical nesting: Deallocate variables initialized for vertical nesting    
     IF ( vnest_init )  CALL vnest_deallocate