&initialization_parameters !-- grid setup nx = 79, ny = 39, nz = 40, dx = 2.0, dy = 2.0, dz = 2.0, !-- pressure solver psolver = 'multigrid', ! required for non-cyclic bc mg_cycles = 2, ! needs to be checked !-- init-action initializing_actions = 'cyclic_fill', !-- advection schemes momentum_advec = 'ws-scheme', ! default advection scheme scalar_advec = 'ws-scheme', !-- wind speed setup ug_surface = 1.0, vg_surface = 0.0, !-- temperature setup neutral = .T., ! strictly neutral flow !-- boundary conditions bc_uv_t = 'neumann', ! free-slip boundary condition bc_lr = 'dirichlet/radiation', ! non-cyclic left/right boundary roughness_length = 0.1, ! applies to all surfaces locally !-- special methods turbulent_inflow = .TRUE., recycling_width = 80.0, inflow_damping_height = 80.0, !-- special setups omega = 0.0, ! no Coriolis force topography = 'read_from_file', topography_grid_convention = 'cell_center', ! default / &runtime_parameters !-- simulation time end_time = 600.0, termination_time_needed = 600.0, !-- others !-- global output settings skip_time_data_output = 0.0, dt_data_output = 600.0, averaging_interval = 600.0, dt_averaging_input = 6.0, !-- time series / run control dt_run_control = 60.0, dt_dots = 60.0, !-- profiles skip_time_dopr = 0.0, dt_dopr = 600.0, averaging_interval_pr = 600.0, dt_averaging_input_pr = 6.0, data_output_pr = '#u', 'u*2', 'wu', 'w*u*', 'w"u"', '#v', 'v*2', 'wv', 'w*v*', 'w"v"', 'w*2', 'e', 'e*', '#km', '#l', !-- output variables data_output = 'u', 'u_av', 'v', 'v_av', 'w', 'w_av', /