&initialization_parameters 

     !-- Numerical grid   
         nx = 49, 
         ny = 49, 
         nz = 100,

         dx = 20.0, 
         dy = 20.0, 
         dz = 20.0,

     !-- Pressure solver  
         fft_method = 'temperton-algorithm',

     !-- Initialization
!         initializing_actions = 'set_constant_profiles',

         initializing_actions = 'set_1d-model_profiles',
         end_time_1d = 864000.0,
         damp_level_1d = 1250.0,
         mixing_length_1d = 'blackadar',
         dissipation_1d = 'detering',

         dt_run_control_1d = 3600.0,
         dt_pr_1d = 86400.0,

     !-- Rayleigh damping
         rayleigh_damping_height = 1250.0,
         rayleigh_damping_factor = 0.01,
         galilei_transformation = .T.,

     !-- Roughness length         
         roughness_length = 0.1, 

     !-- Geostrophic wind (large-scale pressure gradient)
         ug_surface = 5.0, 
         vg_surface = 0.0,

     !-- Neutral stratification
         pt_surface = 300.0,
         pt_vertical_gradient = 0.0,
         pt_vertical_gradient_level = 0.0,
         neutral = .T.,

     !-- Boundary conditions
         bc_uv_b =  'dirichlet', 
         bc_uv_t = 'neumann',

         bc_p_b = 'neumann', 
         bc_p_t = 'dirichlet',

         constant_flux_layer = .T.,
         /

 &runtime_parameters  

     !-- Simulation time 1.5h
         end_time = 5400.0,

     !-- Disturbances
         create_disturbances = .T.,
         dt_disturb = 150.0, 
         disturbance_energy_limit = 0.01,
 
     !-- Run control
         dt_run_control = 60.0,

     !-- Time series
         dt_dopr = 5400.0,
         
     !-- Vertical profiles
         averaging_interval_pr = 5400.0, 
         dt_averaging_input_pr = 50.0,
         data_output_pr = '#u', '#v', '#theta', 
                          'w"u"', 'w*u*', 'wu', 'w"v"', 'w*v*', 'wv','w*2', 'theta*2',
                          
         /
         
 &spectra_parameters
 
     !-- Output settings for spectra
         data_output_sp =  'u', 'v', 'w',
         dt_dosp = 5400.0,
         comp_spectra_level =  2, 5, 10, 15, 20,
         spectra_direction =  'x' , 'y' , 'x',
         averaging_interval_sp = 5400.0,
         
         /