&initialization_parameters !-- grid setup nx = 39, ny = 39, nz = 40, dx = 2.0, dy = 2.0, dz = 2.0, !-- 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 !-- tracer setup passive_scalar = .T., s_surface = 0.0008, !equivalent to 420ppm CO2 as background concentration surface_scalarflux = 0.02, !kg emission per square meter and second bc_s_b = 'neumann', !-- 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', topography_grid_convention = 'cell_center', ! default / &runtime_parameters !-- simulation time end_time = 1800.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 = 6.0, dt_dots = 6.0, !-- profiles skip_time_dopr = 0.0, dt_dopr = 1800.0, averaging_interval_pr = 1800.0, dt_averaging_input_pr = 6.0, do3d_at_begin = .TRUE., do2d_at_begin = .TRUE., dt_do2d_xy = 1800.0, dt_do2d_xz = 1800.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', '#s', 'ws', 'w*s*', 'w"s"', !-- output variables data_output = 'u', 'u_av', 'u_xy', 'u_xz', 'u_xz_av', 'v', 'v_av', 'w', 'w_av', 's', 's_av', 'ssws*_xy', 'pc', 'pc_av', 'pc_xy', 'pc_xz', 'pc_xz_av', section_xy = -1,21,22,23, section_xz = -1,17,18,19, / &particle_parameters !-- Boundary conditions for particles bc_par_lr = 'absorb', bc_par_b = 'reflect', bc_par_t = 'reflect', !-- Define point source psb = 45.0, pst = 45.0, pss = 40.0, psn = 40.0, psl = 25.0, psr = 25.0, !-- Time after new particle set is released dt_prel = 1.0, !-- Time steps after array are deallocated step_dealloc = 100, !-- Particles per point (default = 1) particles_per_point = 100, !-- Maximum time for particles particle_maximum_age = 200.0, !-- Add random fluctuations for initial particle position random_start_position = .FALSE., !-- SGS-velocities for particles use_sgs_for_particles = .TRUE., !-- Particle time series interval dt_dopts = 10.0, / &user_parameters emission_stripe_width_y = 6.0, !in meters emission_stripe_orig_y = 12.0, 62.0, !in meters s_ts_pos_x = 17.0, 61.0, 75.0, s_ts_pos_y = 25.0, 15.0, 41.0, /