[4002] | 1 | &initialization_parameters |
---|
| 2 | |
---|
| 3 | !-- grid setup |
---|
| 4 | nx = 39, |
---|
| 5 | ny = 39, |
---|
| 6 | nz = 40, |
---|
| 7 | |
---|
| 8 | dx = 2.0, |
---|
| 9 | dy = 2.0, |
---|
| 10 | dz = 2.0, |
---|
| 11 | |
---|
| 12 | !-- init-action |
---|
| 13 | initializing_actions = 'set_constant_profiles', |
---|
| 14 | |
---|
| 15 | !-- advection schemes |
---|
| 16 | momentum_advec = 'ws-scheme', ! default advection scheme |
---|
| 17 | scalar_advec = 'ws-scheme', |
---|
| 18 | |
---|
| 19 | !-- wind speed setup |
---|
| 20 | ug_surface = 1.0, |
---|
| 21 | vg_surface = 0.0, |
---|
| 22 | |
---|
| 23 | !-- temperature setup |
---|
| 24 | neutral = .T., ! strictly neutral flow |
---|
| 25 | |
---|
| 26 | !-- tracer setup |
---|
| 27 | passive_scalar = .T., |
---|
| 28 | s_surface = 0.0008, !equivalent to 420ppm CO2 as background concentration |
---|
| 29 | surface_scalarflux = 0.02, !kg emission per square meter and second |
---|
| 30 | bc_s_b = 'neumann', |
---|
| 31 | |
---|
| 32 | !-- boundary conditions |
---|
| 33 | bc_uv_t = 'neumann', ! free-slip boundary condition |
---|
| 34 | roughness_length = 0.1, ! applies to all surfaces locally |
---|
| 35 | |
---|
| 36 | !-- special methods |
---|
| 37 | dp_external = .T., ! use horizontal pressure gradient |
---|
| 38 | dpdxy = -0.0002, 0.0, ! for forcing |
---|
| 39 | |
---|
| 40 | !-- special setups |
---|
| 41 | omega = 0.0, ! no Coriolis force |
---|
| 42 | |
---|
| 43 | topography = 'read_from_file', |
---|
| 44 | topography_grid_convention = 'cell_center', ! default |
---|
| 45 | / |
---|
| 46 | |
---|
| 47 | |
---|
| 48 | &runtime_parameters |
---|
| 49 | |
---|
| 50 | !-- simulation time |
---|
[4098] | 51 | end_time = 10800.0, |
---|
| 52 | termination_time_needed = 900.0, |
---|
[4002] | 53 | |
---|
| 54 | !-- others |
---|
| 55 | create_disturbances = .T., |
---|
| 56 | |
---|
| 57 | !-- global output settings |
---|
| 58 | skip_time_data_output = 0.0, |
---|
[4098] | 59 | dt_data_output = 3600.0, |
---|
[4002] | 60 | |
---|
[4098] | 61 | averaging_interval = 3600.0, |
---|
[4002] | 62 | dt_averaging_input = 6.0, |
---|
| 63 | |
---|
| 64 | !-- time series / run control |
---|
| 65 | dt_run_control = 6.0, |
---|
| 66 | dt_dots = 6.0, |
---|
| 67 | |
---|
| 68 | !-- profiles |
---|
| 69 | skip_time_dopr = 0.0, |
---|
[4098] | 70 | dt_dopr = 3600.0, |
---|
| 71 | averaging_interval_pr = 3600.0, |
---|
[4002] | 72 | dt_averaging_input_pr = 6.0, |
---|
| 73 | |
---|
| 74 | do3d_at_begin = .TRUE., |
---|
| 75 | do2d_at_begin = .TRUE., |
---|
| 76 | |
---|
[4098] | 77 | dt_do2d_xy = 3600.0, |
---|
| 78 | dt_do2d_xz = 3600.0, |
---|
[4002] | 79 | |
---|
| 80 | data_output_pr = '#u', 'u*2', 'wu', 'w*u*', 'w"u"', |
---|
| 81 | '#v', 'v*2', 'wv', 'w*v*', 'w"v"', |
---|
[4098] | 82 | 'w', 'w*2', |
---|
[4002] | 83 | 'e', 'e*', |
---|
| 84 | '#km', |
---|
| 85 | '#l', |
---|
| 86 | '#s', 'ws', 'w*s*', 'w"s"', |
---|
| 87 | |
---|
| 88 | !-- output variables |
---|
| 89 | data_output = 'u', 'u_av', |
---|
| 90 | 'u_xy', 'u_xz', 'u_xz_av', |
---|
| 91 | 'v', 'v_av', |
---|
| 92 | 'w', 'w_av', |
---|
| 93 | 's', 's_av', |
---|
| 94 | 'ssws*_xy', |
---|
| 95 | 'pc', 'pc_av', 'pc_xy', 'pc_xz', 'pc_xz_av', |
---|
| 96 | |
---|
| 97 | |
---|
| 98 | section_xy = -1,21,22,23, |
---|
| 99 | section_xz = -1,17,18,19, |
---|
| 100 | |
---|
| 101 | / |
---|
| 102 | |
---|
| 103 | &particle_parameters |
---|
| 104 | |
---|
| 105 | !-- Boundary conditions for particles |
---|
| 106 | bc_par_lr = 'absorb', |
---|
| 107 | bc_par_b = 'reflect', |
---|
| 108 | bc_par_t = 'reflect', |
---|
| 109 | |
---|
| 110 | !-- Define point source |
---|
| 111 | psb = 45.0, |
---|
| 112 | pst = 45.0, |
---|
| 113 | pss = 40.0, |
---|
| 114 | psn = 40.0, |
---|
| 115 | psl = 25.0, |
---|
| 116 | psr = 25.0, |
---|
| 117 | |
---|
| 118 | !-- Time after new particle set is released |
---|
| 119 | dt_prel = 1.0, |
---|
| 120 | |
---|
| 121 | !-- Time steps after array are deallocated |
---|
[4098] | 122 | step_dealloc = 50, |
---|
[4002] | 123 | |
---|
| 124 | !-- Particles per point (default = 1) |
---|
| 125 | particles_per_point = 100, |
---|
| 126 | |
---|
| 127 | !-- Maximum time for particles |
---|
| 128 | particle_maximum_age = 200.0, |
---|
| 129 | |
---|
| 130 | !-- Add random fluctuations for initial particle position |
---|
| 131 | random_start_position = .FALSE., |
---|
| 132 | |
---|
| 133 | !-- SGS-velocities for particles |
---|
| 134 | use_sgs_for_particles = .TRUE., |
---|
| 135 | |
---|
| 136 | !-- Particle time series interval |
---|
| 137 | dt_dopts = 10.0, |
---|
| 138 | |
---|
| 139 | / |
---|
| 140 | |
---|
| 141 | &user_parameters |
---|
| 142 | emission_stripe_width_y = 6.0, !in meters |
---|
| 143 | emission_stripe_orig_y = 12.0, 62.0, !in meters |
---|
| 144 | s_ts_pos_x = 17.0, 61.0, 75.0, |
---|
| 145 | s_ts_pos_y = 25.0, 15.0, 41.0, |
---|
| 146 | |
---|
| 147 | / |
---|
| 148 | |
---|
| 149 | |
---|