Changeset 2259 for palm


Ignore:
Timestamp:
Jun 8, 2017 9:09:11 AM (7 years ago)
Author:
gronemeier
Message:

Implemented synthetic turbulence generator

Location:
palm/trunk/SOURCE
Files:
1 added
10 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/Makefile

    r2256 r2259  
    2525# -----------------
    2626# $Id$
     27# Implemented synthetic turbulence generator
     28#
     29# 2256 2017-06-07 13:58:08Z suehring
    2730# Remove ring dependency in init_pegrid
    2831#
     
    366369        set_slicer_attributes_dvrp.f90 singleton_mod.f90 sor.f90 spectra_mod.f90 \
    367370        subsidence_mod.f90 sum_up_3d_data.f90 surface_mod.f90 \
    368         surface_coupler.f90 surface_layer_fluxes_mod.f90 swap_timelevel.f90 temperton_fft_mod.f90 \
     371        surface_coupler.f90 surface_layer_fluxes_mod.f90 swap_timelevel.f90 \
     372        synthetic_turbulence_generator_mod.f90 temperton_fft_mod.f90 \
    369373        time_integration.f90 time_to_string.f90 timestep.f90 \
    370374        timestep_scheme_steering.f90 transpose.f90 tridia_solver_mod.f90 \
     
    438442        netcdf_interface_mod.o posix_calls_from_fortran.o
    439443check_parameters.o: modules.o mod_kinds.o land_surface_model_mod.o \
    440         netcdf_interface_mod.o plant_canopy_model_mod.o pmc_interface_mod.o radiation_model_mod.o \
    441         spectra_mod.o subsidence_mod.o microphysics_mod.o wind_turbine_model_mod.o \
    442         urban_surface_mod.o
     444        microphysics_mod.o netcdf_interface_mod.o plant_canopy_model_mod.o pmc_interface_mod.o \
     445        radiation_model_mod.o spectra_mod.o subsidence_mod.o synthetic_turbulence_generator_mod.o \
     446        wind_turbine_model_mod.o urban_surface_mod.o
    443447close_file.o: modules.o mod_kinds.o netcdf_interface_mod.o
    444448compute_vpt.o: modules.o mod_kinds.o
     
    496500init_masks.o: modules.o mod_kinds.o netcdf_interface_mod.o
    497501init_ocean.o: modules.o eqn_state_seawater.o mod_kinds.o
    498 init_pegrid.o: modules.o mod_kinds.o spectra_mod.o
     502init_pegrid.o: modules.o mod_kinds.o spectra_mod.o synthetic_turbulence_generator_mod.o
    499503init_pt_anomaly.o: modules.o mod_kinds.o
    500504init_rankine.o: modules.o mod_kinds.o
     
    547551parin.o: modules.o cpulog_mod.o land_surface_model_mod.o mod_kinds.o netcdf_interface_mod.o \
    548552   plant_canopy_model_mod.o pmc_interface_mod.o progress_bar_mod.o spectra_mod.o \
    549    radiation_model_mod.o microphysics_mod.o virtual_flight_mod.o wind_turbine_model_mod.o
     553   synthetic_turbulence_generator_mod.o radiation_model_mod.o microphysics_mod.o \
     554   virtual_flight_mod.o wind_turbine_model_mod.o
    550555plant_canopy_model_mod.o: modules.o mod_kinds.o surface_mod.o
    551556pmc_interface_mod.o: modules.o mod_kinds.o pmc_child_mod.o pmc_general_mod.o \
     
    596601swap_timelevel.o: modules.o cpulog_mod.o mod_kinds.o land_surface_model_mod.o \
    597602   pmc_interface_mod.o urban_surface_mod.o
     603synthetic_turbulence_generator_mod.o: modules.o mod_kinds.o
    598604temperton_fft_mod.o: modules.o mod_kinds.o
    599605time_integration.o: modules.o advec_ws.o buoyancy.o calc_mean_profile.o \
     
    601607        ls_forcing_mod.o mod_kinds.o nudging_mod.o pmc_interface_mod.o production_e.o \
    602608        prognostic_equations.o progress_bar_mod.o radiation_model_mod.o \
    603         spectra_mod.o user_actions.o microphysics_mod.o surface_layer_fluxes_mod.o surface_mod.o \
     609        spectra_mod.o user_actions.o microphysics_mod.o synthetic_turbulence_generator_mod.o \
     610        surface_layer_fluxes_mod.o surface_mod.o \
    604611        urban_surface_mod.o virtual_flight_mod.o wind_turbine_model_mod.o
    605612time_to_string.o: mod_kinds.o
  • palm/trunk/SOURCE/check_parameters.f90

    r2251 r2259  
    2020! Current revisions:
    2121! -----------------
    22 !  
     22!
    2323!
    2424! Former revisions:
    2525! -----------------
    2626! $Id$
     27! Implemented synthetic turbulence generator
     28!
     29! 2251 2017-06-06 15:10:46Z Giersch
    2730!
    2831! 2250 2017-06-06 15:07:12Z Giersch
     
    499502    USE subsidence_mod
    500503    USE statistics
     504    USE synthetic_turbulence_generator_mod,                                    &
     505        ONLY:  stg_check_parameters
    501506    USE transpose_indices
    502507    USE urban_surface_mod,                                                     &
     
    12071212    ENDIF
    12081213
     1214!
     1215!-- When synthetic turbulence generator is used, peform addtional checks
     1216    IF ( synthetic_turbulence_generator )  CALL stg_check_parameters
    12091217!
    12101218!-- When plant canopy model is used, peform addtional checks
  • palm/trunk/SOURCE/header.f90

    r2258 r2259  
    2525! -----------------
    2626! $Id$
     27! Implemented synthetic turbulence generator
     28!
     29! 2258 2017-06-08 07:55:13Z suehring
    2730! Bugfix, add pre-preprocessor directives to enable non-parrallel mode
    2831!
     
    356359    USE plant_canopy_model_mod,                                                &
    357360        ONLY:  pcm_header, plant_canopy
    358        
    359 #if defined( __parallel )
     361
    360362    USE pmc_handle_communicator,                                               &
    361363        ONLY:  pmc_get_model_info
    362 #endif
    363364
    364365    USE pmc_interface,                                                         &
     
    373374    USE surface_mod,                                                           &
    374375        ONLY:  surf_def_h
     376
     377    USE synthetic_turbulence_generator_mod,                                    &
     378        ONLY:  stg_header
    375379
    376380    IMPLICIT NONE
     
    552556       WRITE ( io, 600 )  TRIM( nesting_mode ),                                &
    553557                          TRIM( nesting_datatransfer_mode )
    554 #if defined( __parallel )
    555558       CALL pmc_get_model_info( ncpl = ncpl, cpl_id = my_cpl_id )
    556559
     
    570573                             TRIM( cpl_name )
    571574       ENDDO
    572 #endif
    573575    ENDIF
    574576    WRITE ( io, 99 )
     
    988990       ENDIF
    989991    ENDIF
     992
     993    IF ( synthetic_turbulence_generator )  CALL stg_header ( io )
    990994
    991995    IF ( plant_canopy )  CALL pcm_header ( io )
  • palm/trunk/SOURCE/init_3d_model.f90

    r2252 r2259  
    2525! -----------------
    2626! $Id$
     27! Implemented synthetic turbulence generator
     28!
     29! 2252 2017-06-07 09:35:37Z knoop
    2730! rho_air now depending on surface_pressure even in Boussinesq mode
    2831!
     
    390393               sums_l_l, sums_up_fraction_l, sums_wsts_bc_l, ts_value,         &
    391394               weight_pres, weight_substep
    392  
     395
     396    USE synthetic_turbulence_generator_mod,                                    &
     397        ONLY:  stg_init, use_synthetic_turbulence_generator
     398
    393399    USE surface_layer_fluxes_mod,                                              &
    394400        ONLY:  init_surface_layer_fluxes
     
    10571063
    10581064!
     1065!--       Overwrite initial profiles in case of synthetic turbulence generator
     1066          IF( use_synthetic_turbulence_generator ) THEN
     1067             CALL stg_init
     1068          ENDIF
     1069
     1070!
    10591071!--       Use constructed initial profiles (velocity constant with height,
    10601072!--       temperature profile with constant gradient)
     
    15381550       IF ( passive_scalar )  ts_m  = 0.0_wp
    15391551       IF ( ocean          )  tsa_m = 0.0_wp
     1552!
     1553!--    Initialize synthetic turbulence generator in case of restart.
     1554       IF ( TRIM( initializing_actions ) == 'read_restart_data'  .AND.         &
     1555            use_synthetic_turbulence_generator )  CALL stg_init
    15401556
    15411557       CALL location_message( 'finished', .TRUE. )
  • palm/trunk/SOURCE/init_pegrid.f90

    r2238 r2259  
    2525! -----------------
    2626! $Id$
     27! Implemented synthetic turbulence generator
     28!
     29! 2238 2017-05-31 16:49:16Z suehring
    2730! Remove unnecessary module load of pmc_interface
    2831!
     
    209212    USE spectra_mod,                                                           &
    210213        ONLY:  calculate_spectra, dt_dosp
     214
     215    USE synthetic_turbulence_generator_mod,                                    &
     216        ONLY:  use_synthetic_turbulence_generator
    211217
    212218    USE transpose_indices,                                                     &
     
    423429!-- 1. transposition  z --> x
    424430!-- This transposition is not neccessary in case of a 1d-decomposition along x
    425     IF ( psolver == 'poisfft'  .OR.  calculate_spectra )  THEN
     431    IF ( psolver == 'poisfft'  .OR.  calculate_spectra  .OR.                   &
     432         use_synthetic_turbulence_generator )  THEN
    426433
    427434       IF ( pdims(2) /= 1 )  THEN
  • palm/trunk/SOURCE/modules.f90

    r2256 r2259  
    2525! -----------------
    2626! $Id$
     27! Implemented synthetic turbulence generator
     28!
     29! 2256 2017-06-07 13:58:08Z suehring
    2730! Change default value of zeta_min to -20
    2831! Increase dimension for wall_heatflux, etc.
     
    11651168    LOGICAL ::  stop_dt = .FALSE.                            !<
    11661169    LOGICAL ::  synchronous_exchange = .FALSE.               !<
     1170    LOGICAL ::  synthetic_turbulence_generator = .FALSE.     !< flag for synthetic turbulence generator module
    11671171    LOGICAL ::  terminate_run = .FALSE.                      !<
    11681172    LOGICAL ::  transpose_compute_overlap = .FALSE.          !<
  • palm/trunk/SOURCE/parin.f90

    r2233 r2259  
    2525! -----------------
    2626! $Id$
     27! Implemented synthetic turbulence generator
     28!
     29! 2233 2017-05-30 18:08:54Z suehring
    2730!
    2831! 2232 2017-05-30 17:47:52Z suehring
     
    315318    USE statistics,                                                            &
    316319        ONLY:  hom, hom_sum, pr_palm, region, statistic_regions
     320
     321    USE synthetic_turbulence_generator_mod,                                    &
     322        ONLY:  stg_parin
    317323
    318324    USE urban_surface_mod,                                                     &
     
    605611!--       if required
    606612          CALL flight_parin
    607 
     613!
     614!--       Check if synthetic turbulence generator is used and read stg_par if
     615!--       required
     616          CALL stg_parin
    608617!
    609618!--       Read user-defined variables
  • palm/trunk/SOURCE/read_var_list.f90

    r2233 r2259  
    2525! -----------------
    2626! $Id$
     27! Implemented synthetic turbulence generator
     28!
     29! 2233 2017-05-30 18:08:54Z suehring
    2730!
    2831! 2232 2017-05-30 17:47:52Z suehring
    2932! Replace wall_qflux, wall_sflux by wall_humidityflux and wall_scalarflux; add
    3033! wall_salinityflux
    31 ! +tunnel_height, tunnel_length, tunnel_width_x, tunnel_width_y, 
     34! +tunnel_height, tunnel_length, tunnel_width_x, tunnel_width_y,
    3235!  tunnel_wall_depth
    3336!
    3437! 2042 2016-11-02 13:47:31Z suehring
    3538! Bugfix, read restart data for wall_heatflux, wall_qflux and wall_sflux
    36 ! 
     39!
    3740! 2000 2016-08-20 18:09:15Z knoop
    3841! Forced header and separation lines into 80 columns
    39 ! 
     42!
    4043! 1992 2016-08-12 15:14:59Z suehring
    4144! top_scalarflux added
    42 ! 
     45!
    4346! 1960 2016-07-12 16:34:24Z suehring
    4447! Separate humidity and passive scalar
     
    6669! 1699 2015-10-29 08:02:35Z maronga
    6770! Bugfix: update of binary version from 3.9b to 4.0 was missing
    68 ! 
     71!
    6972! 1691 2015-10-26 16:17:44Z maronga
    70 ! Added output of most_method, constant_flux_layer, zeta_min, zeta_max. Removed 
     73! Added output of most_method, constant_flux_layer, zeta_min, zeta_max. Removed
    7174! output of prandtl_layer and rif_min, rif_max.
    72 ! 
     75!
    7376! 1682 2015-10-07 23:56:08Z knoop
    7477! Code annotations made doxygen readable
    75 ! 
     78!
    7679! 1615 2015-07-08 18:49:19Z suehring
    7780! Enable turbulent inflow for passive_scalar and humidity
     
    7982! 1585 2015-04-30 07:05:52Z maronga
    8083! Adapted for RRTMG
    81 ! 
     84!
    8285! 1560 2015-03-06 10:48:54Z keck
    8386! +recycling_yshift
    84 ! 
     87!
    8588! 1522 2015-01-14 10:53:12Z keck
    8689! added query for checking if the advection scheme in the restart run is the
    8790! same as the advection scheme in the corresponding initial run
    88 ! 
     91!
    8992! 1502 2014-12-03 18:22:31Z kanani
    90 ! Canopy module and parameters removed (parameters are always read from 
     93! Canopy module and parameters removed (parameters are always read from
    9194! canopy_par NAMELIST for initial and restart runs)
    92 ! 
     95!
    9396! 1496 2014-12-02 17:25:50Z maronga
    9497! Renamed "radiation" -> "cloud_top_radiation"
    95 ! 
     98!
    9699! 1484 2014-10-21 10:53:05Z kanani
    97100! Changes in the course of the canopy-model modularization:
     
    118121!
    119122! 1253 2013-11-07 10:48:12Z fricke
    120 ! Bugfix: add ref_state to read_parts_of_var_list, otherwise ref_state 
     123! Bugfix: add ref_state to read_parts_of_var_list, otherwise ref_state
    121124! is zero for initializing_actions = 'cyclic_fill'
    122125!
     
    142145!
    143146! 1053 2012-11-13 17:11:03Z hoffmann
    144 ! necessary expansions according to the two new prognostic equations (nr, qr) 
     147! necessary expansions according to the two new prognostic equations (nr, qr)
    145148! of the two-moment cloud physics scheme:
    146149! +bc_*_b, +bc_*_t, +bc_*_t_val, *_init, *_surface, *_surface_initial_change,
     
    148151! +surface_waterflux_*
    149152!
    150 ! in addition, steering parameters parameters of the two-moment cloud physics 
    151 ! scheme:   
     153! in addition, steering parameters parameters of the two-moment cloud physics
     154! scheme:
    152155! +cloud_scheme, +drizzle, +mu_constant, +mu_constant_value, +ventilation_effect
    153156!
     
    191194!------------------------------------------------------------------------------!
    192195 SUBROUTINE read_var_list
    193  
     196
    194197
    195198    USE arrays_3d,                                                             &
     
    201204    USE flight_mod,                                                            &
    202205        ONLY:  flight_read_restart_data
    203    
     206
    204207    USE grid_variables,                                                        &
    205208        ONLY:  dx, dy
     
    234237               v_max, v_max_ijk, w_max, w_max_ijk
    235238
     239    USE synthetic_turbulence_generator_mod,                                    &
     240        ONLY:  stg_read_restart_data
    236241
    237242    IMPLICIT NONE
     
    628633             READ ( 13 )  s_vertical_gradient_level
    629634          CASE ( 's_vertical_gradient_level_ind' )
    630              READ ( 13 )  s_vertical_gradient_level_ind 
     635             READ ( 13 )  s_vertical_gradient_level_ind
    631636          CASE ( 'sa_init' )
    632637             READ ( 13 )  sa_init
     
    646651             READ ( 13 )  surface_pressure
    647652          CASE ( 'surface_scalarflux' )
    648              READ ( 13 )  surface_scalarflux             
     653             READ ( 13 )  surface_scalarflux
    649654          CASE ( 'surface_waterflux' )
    650              READ ( 13 )  surface_waterflux     
     655             READ ( 13 )  surface_waterflux
    651656          CASE ( 'time_coupling' )
    652657             READ ( 13 )  time_coupling
     
    770775             READ ( 13 )  vg_vertical_gradient_level_ind
    771776          CASE ( 'virtual_flight' )
    772              READ ( 13 ) virtual_flight 
     777             READ ( 13 ) virtual_flight
    773778          CASE ( 'volume_flow_area' )
    774779             READ ( 13 )  volume_flow_area
     
    805810             WRITE( message_string, * ) 'unknown variable named "',         &
    806811                                        TRIM( variable_chr ), '" found in', &
    807                                         ' data from prior run on PE ', myid 
     812                                        ' data from prior run on PE ', myid
    808813             CALL message( 'read_var_list', 'PA0302', 1, 2, 0, 6, 0 )
    809814        END SELECT
     
    813818
    814819    ENDDO
    815    
     820
    816821    IF ( virtual_flight )  CALL flight_read_restart_data
    817822
     823    IF ( synthetic_turbulence_generator )  CALL stg_read_restart_data
    818824
    819825 END SUBROUTINE read_var_list
     
    828834!> run which used a smaller total domain or/and a different domain decomposition.
    829835!------------------------------------------------------------------------------!
    830  
     836
    831837 SUBROUTINE read_parts_of_var_list
    832838
     
    938944    READ ( 13 )  statistic_regions_on_file
    939945    IF ( statistic_regions_on_file /= statistic_regions )  THEN
    940        WRITE( message_string, * ) 'statistic regions on restart data file ',& 
     946       WRITE( message_string, * ) 'statistic regions on restart data file ',&
    941947                                  'differ from the current run',            &
    942948                                  '&statistic regions on file    = "',      &
     
    993999                CALL message( 'read_parts_of_var_list', 'PA0100', &
    9941000                                                                 1, 2, 0, 6, 0 )
    995              END IF                           
    996              
     1001             END IF
     1002
    9971003          CASE ( 'nx' )
    9981004             READ ( 13 )  nx_on_file
     
    10121018                CALL message( 'read_parts_of_var_list', 'PA0101', &
    10131019                                                                 1, 2, 0, 6, 0 )
    1014              END IF             
    1015              
     1020             END IF
     1021
    10161022          CASE DEFAULT
    10171023
     
    10401046!> Skipping the global control variables from restart-file (binary format)
    10411047!------------------------------------------------------------------------------!
    1042  
     1048
    10431049 SUBROUTINE skip_var_list
    10441050
    10451051    USE control_parameters,                                                    &
    10461052        ONLY:  virtual_flight
    1047  
     1053
    10481054    USE flight_mod,                                                            &
    10491055        ONLY:  flight_skip_var_list
    1050  
     1056
    10511057    IMPLICIT NONE
    10521058
     
    10681074    ENDDO
    10691075!
    1070 !-- In case of virtual flights, skip also variables related to 
     1076!-- In case of virtual flights, skip also variables related to
    10711077!-- this module.
    10721078    IF ( virtual_flight )  CALL flight_skip_var_list
    1073    
     1079
    10741080
    10751081 END SUBROUTINE skip_var_list
  • palm/trunk/SOURCE/time_integration.f90

    r2233 r2259  
    2525! -----------------
    2626! $Id$
     27! Implemented synthetic turbulence generator
     28!
     29! 2233 2017-05-30 18:08:54Z suehring
    2730!
    2831! 2232 2017-05-30 17:47:52Z suehring
     
    373376    USE urban_surface_mod,                                                     &
    374377        ONLY:  usm_material_heat_model, usm_material_model,                    &
    375                usm_radiation, usm_surface_energy_balance             
     378               usm_radiation, usm_surface_energy_balance
     379
     380    USE synthetic_turbulence_generator_mod,                                    &
     381        ONLY:  stg_main, use_synthetic_turbulence_generator
    376382
    377383    USE user_actions_mod,                                                      &
     
    665671!--       Impose a turbulent inflow using the recycling method
    666672          IF ( turbulent_inflow )  CALL  inflow_turbulence
     673
     674!
     675!--       Impose a turbulent inflow using synthetic generated turbulence
     676         IF ( use_synthetic_turbulence_generator ) THEN
     677            CALL  stg_main
     678         ENDIF
    667679
    668680!
  • palm/trunk/SOURCE/write_var_list.f90

    r2233 r2259  
    2525! -----------------
    2626! $Id$
     27! Implemented synthetic turbulence generator
     28!
     29! 2233 2017-05-30 18:08:54Z suehring
    2730!
    2831! 2232 2017-05-30 17:47:52Z suehring
    2932! Replace wall_qflux, wall_sflux by wall_humidityflux and wall_scalarflux; add
    3033! wall_salinityflux
    31 ! +tunnel_height, tunnel_lenth, tunnel_width_x, tunnel_width_y, 
     34! +tunnel_height, tunnel_lenth, tunnel_width_x, tunnel_width_y,
    3235! tunnel_wall_depth
    33 ! 
     36!
    3437! 2042 2016-11-02 13:47:31Z suehring
    3538! Bugfix, write restart data for wall_heatflux, wall_qflux and wall_sflux
    36 ! 
     39!
    3740! 2000 2016-08-20 18:09:15Z knoop
    3841! Forced header and separation lines into 80 columns
    39 ! 
     42!
    4043! 1992 2016-08-12 15:14:59Z suehring
    41 ! top scalarflux added 
     44! top scalarflux added
    4245!
    4346! 1957 2016-07-07 10:43:48Z suehring
    4447! flight module added
    45 ! 
     48!
    4649! 1849 2016-04-08 11:33:18Z hoffmann
    4750! Adapted for modularization of microphysics
     
    5053! spectra_mod added
    5154!
    52 ! 1831 2016-04-07 13:15:51Z hoffmann 
    53 ! turbulence renamed collision_turbulence, drizzle renamed 
     55! 1831 2016-04-07 13:15:51Z hoffmann
     56! turbulence renamed collision_turbulence, drizzle renamed
    5457! cloud_water_sedimentation
    5558!
     
    5962! 1705 2015-11-02 14:28:56Z maronga
    6063! Bugfix: two lines required swapping
    61 ! 
     64!
    6265! 1691 2015-10-26 16:17:44Z maronga
    63 ! Added output of most_method, constant_flux_layer, zeta_min, zeta_max. Removed 
     66! Added output of most_method, constant_flux_layer, zeta_min, zeta_max. Removed
    6467! output of prandtl_layer and rif_min, rif_max.
    65 ! 
     68!
    6669! 1682 2015-10-07 23:56:08Z knoop
    67 ! Code annotations made doxygen readable 
    68 ! 
     70! Code annotations made doxygen readable
     71!
    6972! 1585 2015-04-30 07:05:52Z maronga
    7073! Adapted for RRTMG
    71 ! 
     74!
    7275! 1551 2015-03-03 14:18:16Z maronga
    7376! Typo removed
    74 ! 
     77!
    7578! 1502 2014-12-03 18:22:31Z kanani
    76 ! Canopy module and parameters removed (parameters are always read from 
     79! Canopy module and parameters removed (parameters are always read from
    7780! canopy_par NAMELIST for initial and restart runs),
    78 ! Bugfix: added blanks in "cloud_top_radiation"-string to a total of 30 
     81! Bugfix: added blanks in "cloud_top_radiation"-string to a total of 30
    7982! characters
    80 ! 
     83!
    8184! 1496 2014-12-02 17:25:50Z maronga
    8285! Renamed "radiation" -> "cloud_top_radiation"
    83 ! 
     86!
    8487! 1484 2014-10-21 10:53:05Z kanani
    8588! Changes in the course of the canopy-model modularization:
     
    9093!
    9194! 1324 2014-03-21 09:13:16Z suehring
    92 ! Bugfix: ONLY statement for module netcdf_control removed 
     95! Bugfix: ONLY statement for module netcdf_control removed
    9396!
    9497! 1320 2014-03-20 08:40:49Z raasch
    9598! revision history before 2012 removed,
    9699! comment fields (!:) to be used for variable explanations added to
    97 ! all variable declaration statements 
     100! all variable declaration statements
    98101!
    99102! 1308 2014-03-13 14:58:42Z fricke
     
    116119!
    117120! 1053 2012-11-13 17:11:03Z hoffmann
    118 ! necessary expansions according to the two new prognostic equations (nr, qr) 
     121! necessary expansions according to the two new prognostic equations (nr, qr)
    119122! of the two-moment cloud physics scheme:
    120123! +bc_*_b, bc_*_t, bc_*_t_val, *_init, *_surface, *_surface_initial_change,
     
    122125! +surface_waterflux_*
    123126!
    124 ! in addition, steering parameters parameters of the two-moment cloud physics 
    125 ! scheme:   
     127! in addition, steering parameters parameters of the two-moment cloud physics
     128! scheme:
    126129! +cloud_scheme, +drizzle, +mu_constant, +mu_constant_value, +ventilation_effect
    127130!
     
    166169!------------------------------------------------------------------------------!
    167170 SUBROUTINE write_var_list
    168  
     171
    169172
    170173    USE arrays_3d,                                                             &
     
    173176
    174177    USE control_parameters
    175    
     178
    176179    USE flight_mod,                                                            &
    177180        ONLY:  flight_write_restart_data
    178    
     181
    179182    USE grid_variables,                                                        &
    180183        ONLY:  dx, dy
    181    
     184
    182185    USE indices,                                                               &
    183186        ONLY:  nz, nx, ny
     
    190193    USE model_1d,                                                              &
    191194        ONLY:  damp_level_1d, dt_pr_1d, dt_run_control_1d, end_time_1d
    192    
     195
    193196    USE netcdf_interface,                                                      &
    194197        ONLY:  netcdf_precision, output_for_t0
    195    
     198
    196199    USE particle_attributes,                                                   &
    197200        ONLY:  curvature_solution_effects, time_sort_particles
    198    
     201
    199202    USE pegrid
    200203
     
    204207    USE spectra_mod,                                                           &
    205208        ONLY:  average_count_sp
     209
     210    USE synthetic_turbulence_generator_mod,                                    &
     211        ONLY:  stg_write_restart_data
    206212
    207213    USE statistics,                                                            &
     
    209215               v_max_ijk, w_max, w_max_ijk
    210216
    211    
     217
    212218    IMPLICIT NONE
    213219
    214     CHARACTER (LEN=10) ::  binary_version   !< 
     220    CHARACTER (LEN=10) ::  binary_version   !<
    215221
    216222
     
    531537    WRITE ( 14 )  s_vertical_gradient_level
    532538    WRITE ( 14 )  's_vertical_gradient_level_ind '
    533     WRITE ( 14 )  s_vertical_gradient_level_ind 
     539    WRITE ( 14 )  s_vertical_gradient_level_ind
    534540    WRITE ( 14 )  'sa_init                       '
    535541    WRITE ( 14 )  sa_init
     
    549555    WRITE ( 14 )  surface_pressure
    550556    WRITE ( 14 )  'surface_scalarflux            '
    551     WRITE ( 14 )  surface_scalarflux   
     557    WRITE ( 14 )  surface_scalarflux
    552558    WRITE ( 14 )  'surface_waterflux             '
    553     WRITE ( 14 )  surface_waterflux   
     559    WRITE ( 14 )  surface_waterflux
    554560    WRITE ( 14 )  's_surface                     '
    555561    WRITE ( 14 )  s_surface
     
    711717    WRITE ( 14 )  zeta_min
    712718    WRITE ( 14 )  'z0h_factor                    '
    713     WRITE ( 14 )  z0h_factor 
     719    WRITE ( 14 )  z0h_factor
    714720
    715721!
     
    718724!
    719725!-- If required, write restart data for virtual measurements.
    720     IF ( virtual_flight )  CALL flight_write_restart_data 
    721    
    722    
     726    IF ( virtual_flight )  CALL flight_write_restart_data
     727
     728!
     729!-- If required, write restart data for virtual measurements.
     730    IF ( synthetic_turbulence_generator )  CALL stg_write_restart_data
     731
     732
    723733 END SUBROUTINE write_var_list
Note: See TracChangeset for help on using the changeset viewer.