&initialization_parameters !-- grid setup nx = 39, ny = 39, nz = 40, dx = 2.0, dy = 2.0, dz = 2.0, !-- pressure solver fft_method = 'temperton-algorithm', ! fast, but restrictions apply !-- init-action initializing_actions = 'set_constant_profiles', !-- 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 roughness_length = 0.1, ! applies to all surfaces locally !-- special methods dp_external = .T., ! use horizontal pressure gradient dpdxy = -0.0002, 0.0, ! for forcing !-- special setups omega = 0.0, ! no Coriolis force topography = 'read_from_file', / &runtime_parameters !-- simulation time end_time = 2700.0, termination_time_needed = 900.0, !-- others create_disturbances = .T., !-- global output settings skip_time_data_output = 0.0, dt_data_output = 1800.0, averaging_interval = 1800.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 = 1800.0, averaging_interval_pr = 1800.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', /