 &initialization_parameters
 
         nx = 191, ny = 191, nz = 64,
         dx = 50.0, dy = 50.0, dz = 50.0,

         fft_method = 'fftw',

         humidity = .T., 

         large_scale_forcing = .T.,
         large_scale_subsidence = .T.,
         use_subsidence_tendencies = .F.,

         nudging = .T.,

         momentum_advec = 'ws-scheme',
         scalar_advec   = 'ws-scheme',

         initializing_actions = 'set_constant_profiles',

         bc_uv_b = 'dirichlet',
         bc_pt_b = 'dirichlet',
         bc_q_b  = 'dirichlet',

         latitude = 50.92,

         pt_reference = 290.0,
         reference_state = 'single_value', 

         constant_flux_layer = .T.,
         use_surface_fluxes = .T., 
         roughness_length = 0.01,

         rayleigh_damping_factor = 0.05,
         rayleigh_damping_height = 2700.0,
         scalar_rayleigh_damping = .T., /


 &runtime_parameters
 
         end_time = 21600.0, 

         termination_time_needed = 300.0,

         npex = 64,
         npey = 1, 

         create_disturbances = .T.,

         netcdf_data_format = 2,

         dt_run_control = 30.0,
         dt_dots        = 30.0,

         dt_dopr = 1800.0,
         averaging_interval_pr = 1800.0, 
         dt_averaging_input_pr = 30.0,

         data_output_pr = '#u','#v','w','#theta','ug','vg',
                          '#thetal','#q','#qv','ql',
                          'e','e*','p','#km','#kh','#l',
                          'w"u"','w*u*','wu', 'w"v"','w*v*','wv',
                          'w"theta"', 'w*theta*', 'wtheta', 
                          'w"thetav"', 'w*thetav*', 'wthetav', 
                          'w"q"', 'w*q*', 'wq', 
                          'w"qv"', 'w*qv*', 'wqv', 
                          'w*e*', 'u*2', 'v*2', 'w*2', 'theta*2', 'q*2',
                          'w*3', 'Sw', 'w*2theta*', 'w*theta*2',
                          'w*u*u*:dz','w*p*:dz', 'w"e:dz', 'hyp',
                          'td_lsa_thetal', 'td_lsa_q', 'td_sub_thetal', 'td_sub_q', 
                          'td_nud_thetal', 'td_nud_q', 'td_nud_u', 'td_nud_v',

         data_output = 'u', 'v', 'w', 'theta', 'thetal', 'p', 'q', 'ql'
                       'u_xy', 'v_xy', 'w_xy', 'theta_xy', 'thetal_xy',
                       'p_xy', 'q_xy', 'ql_xy', 
                       'u_xz', 'v_xz', 'w_xz', 'theta_xz', 'thetal_xz',
                       'p_xz', 'q_xz', 'ql_xz', 
                       'shf*_xy', 'qsws*_xy', 'lwp*_xy',

         nz_do3d = 32, 

         dt_do3d        = 3600.0,
         dt_do2d_xy     = 3600.0,
         dt_do2d_xz     = 3600.0,
         dt_data_output = 3600.0,

         do2d_at_begin = .T., !->check perturbations

         section_xy = 2, 4, 10, 15, 20, 25,  
         section_xz = 96,

         data_output_2d_on_each_pe = .T.,  /
         
 &bulk_cloud_parameters
 
          bulk_cloud_model = .T., /
