Changeset 2259 for palm/trunk/SOURCE/read_var_list.f90

Timestamp:

Jun 8, 2017 9:09:11 AM (7 years ago)

Author:

gronemeier

Message:

Implemented synthetic turbulence generator

File:

• palm/trunk/SOURCE/read_var_list.f90 (modified) (20 diffs)

palm/trunk/SOURCE/read_var_list.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Implemented synthetic turbulence generator
+!
+! 2233 2017-05-30 18:08:54Z suehring
 !
 ! 2232 2017-05-30 17:47:52Z suehring
 ! Replace wall_qflux, wall_sflux by wall_humidityflux and wall_scalarflux; add
 ! wall_salinityflux
-! +tunnel_height, tunnel_length, tunnel_width_x, tunnel_width_y, 
+! +tunnel_height, tunnel_length, tunnel_width_x, tunnel_width_y,
 !  tunnel_wall_depth
 !
 ! 2042 2016-11-02 13:47:31Z suehring
 ! Bugfix, read restart data for wall_heatflux, wall_qflux and wall_sflux
-! 
+!
 ! 2000 2016-08-20 18:09:15Z knoop
 ! Forced header and separation lines into 80 columns
-! 
+!
 ! 1992 2016-08-12 15:14:59Z suehring
 ! top_scalarflux added
-! 
+!
 ! 1960 2016-07-12 16:34:24Z suehring
 ! Separate humidity and passive scalar
@@ -66 +69 @@
 ! 1699 2015-10-29 08:02:35Z maronga
 ! Bugfix: update of binary version from 3.9b to 4.0 was missing
-! 
+!
 ! 1691 2015-10-26 16:17:44Z maronga
-! Added output of most_method, constant_flux_layer, zeta_min, zeta_max. Removed 
+! Added output of most_method, constant_flux_layer, zeta_min, zeta_max. Removed
 ! output of prandtl_layer and rif_min, rif_max.
-! 
+!
 ! 1682 2015-10-07 23:56:08Z knoop
 ! Code annotations made doxygen readable
-! 
+!
 ! 1615 2015-07-08 18:49:19Z suehring
 ! Enable turbulent inflow for passive_scalar and humidity
@@ -79 +82 @@
 ! 1585 2015-04-30 07:05:52Z maronga
 ! Adapted for RRTMG
-! 
+!
 ! 1560 2015-03-06 10:48:54Z keck
 ! +recycling_yshift
-! 
+!
 ! 1522 2015-01-14 10:53:12Z keck
 ! added query for checking if the advection scheme in the restart run is the
 ! same as the advection scheme in the corresponding initial run
-! 
+!
 ! 1502 2014-12-03 18:22:31Z kanani
-! Canopy module and parameters removed (parameters are always read from 
+! Canopy module and parameters removed (parameters are always read from
 ! canopy_par NAMELIST for initial and restart runs)
-! 
+!
 ! 1496 2014-12-02 17:25:50Z maronga
 ! Renamed "radiation" -> "cloud_top_radiation"
-! 
+!
 ! 1484 2014-10-21 10:53:05Z kanani
 ! Changes in the course of the canopy-model modularization:
@@ -118 +121 @@
 !
 ! 1253 2013-11-07 10:48:12Z fricke
-! Bugfix: add ref_state to read_parts_of_var_list, otherwise ref_state 
+! Bugfix: add ref_state to read_parts_of_var_list, otherwise ref_state
 ! is zero for initializing_actions = 'cyclic_fill'
 !
@@ -142 +145 @@
 !
 ! 1053 2012-11-13 17:11:03Z hoffmann
-! necessary expansions according to the two new prognostic equations (nr, qr) 
+! necessary expansions according to the two new prognostic equations (nr, qr)
 ! of the two-moment cloud physics scheme:
 ! +bc_*_b, +bc_*_t, +bc_*_t_val, *_init, *_surface, *_surface_initial_change,
@@ -148 +151 @@
 ! +surface_waterflux_*
 !
-! in addition, steering parameters parameters of the two-moment cloud physics 
-! scheme:   
+! in addition, steering parameters parameters of the two-moment cloud physics
+! scheme:
 ! +cloud_scheme, +drizzle, +mu_constant, +mu_constant_value, +ventilation_effect
 !
@@ -191 +194 @@
 !------------------------------------------------------------------------------!
  SUBROUTINE read_var_list
- 
+
 
     USE arrays_3d,                                                             &
@@ -201 +204 @@
     USE flight_mod,                                                            &
         ONLY:  flight_read_restart_data
-    
+
     USE grid_variables,                                                        &
         ONLY:  dx, dy
@@ -234 +237 @@
                v_max, v_max_ijk, w_max, w_max_ijk
 
+    USE synthetic_turbulence_generator_mod,                                    &
+        ONLY:  stg_read_restart_data
 
     IMPLICIT NONE
@@ -628 +633 @@
              READ ( 13 )  s_vertical_gradient_level
           CASE ( 's_vertical_gradient_level_ind' )
-             READ ( 13 )  s_vertical_gradient_level_ind 
+             READ ( 13 )  s_vertical_gradient_level_ind
           CASE ( 'sa_init' )
              READ ( 13 )  sa_init
@@ -646 +651 @@
              READ ( 13 )  surface_pressure
           CASE ( 'surface_scalarflux' )
-             READ ( 13 )  surface_scalarflux             
+             READ ( 13 )  surface_scalarflux
           CASE ( 'surface_waterflux' )
-             READ ( 13 )  surface_waterflux      
+             READ ( 13 )  surface_waterflux
           CASE ( 'time_coupling' )
              READ ( 13 )  time_coupling
@@ -770 +775 @@
              READ ( 13 )  vg_vertical_gradient_level_ind
           CASE ( 'virtual_flight' )
-             READ ( 13 ) virtual_flight  
+             READ ( 13 ) virtual_flight
           CASE ( 'volume_flow_area' )
              READ ( 13 )  volume_flow_area
@@ -805 +810 @@
              WRITE( message_string, * ) 'unknown variable named "',         &
                                         TRIM( variable_chr ), '" found in', &
-                                        ' data from prior run on PE ', myid 
+                                        ' data from prior run on PE ', myid
              CALL message( 'read_var_list', 'PA0302', 1, 2, 0, 6, 0 )
         END SELECT
@@ -813 +818 @@
 
     ENDDO
-    
+
     IF ( virtual_flight )  CALL flight_read_restart_data
 
+    IF ( synthetic_turbulence_generator )  CALL stg_read_restart_data
 
  END SUBROUTINE read_var_list
@@ -828 +834 @@
 !> run which used a smaller total domain or/and a different domain decomposition.
 !------------------------------------------------------------------------------!
- 
+
  SUBROUTINE read_parts_of_var_list
 
@@ -938 +944 @@
     READ ( 13 )  statistic_regions_on_file
     IF ( statistic_regions_on_file /= statistic_regions )  THEN
-       WRITE( message_string, * ) 'statistic regions on restart data file ',& 
+       WRITE( message_string, * ) 'statistic regions on restart data file ',&
                                   'differ from the current run',            &
                                   '&statistic regions on file    = "',      &
@@ -993 +999 @@
                 CALL message( 'read_parts_of_var_list', 'PA0100', &
                                                                  1, 2, 0, 6, 0 )
-             END IF                           
-             
+             END IF
+
           CASE ( 'nx' )
              READ ( 13 )  nx_on_file
@@ -1012 +1018 @@
                 CALL message( 'read_parts_of_var_list', 'PA0101', &
                                                                  1, 2, 0, 6, 0 )
-             END IF             
-             
+             END IF
+
           CASE DEFAULT
 
@@ -1040 +1046 @@
 !> Skipping the global control variables from restart-file (binary format)
 !------------------------------------------------------------------------------!
- 
+
  SUBROUTINE skip_var_list
 
     USE control_parameters,                                                    &
         ONLY:  virtual_flight
- 
+
     USE flight_mod,                                                            &
         ONLY:  flight_skip_var_list
- 
+
     IMPLICIT NONE
 
@@ -1068 +1074 @@
     ENDDO
 !
-!-- In case of virtual flights, skip also variables related to 
+!-- In case of virtual flights, skip also variables related to
 !-- this module.
     IF ( virtual_flight )  CALL flight_skip_var_list
-    
+
 
  END SUBROUTINE skip_var_list