[4224] | 1 | !------------------------------------------------------------------------------- |
---|
| 2 | !-- INITIALIZATION PARAMETER NAMELIST |
---|
| 3 | ! Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/inipar |
---|
| 4 | !------------------------------------------------------------------------------- |
---|
[4002] | 5 | &initialization_parameters |
---|
[4224] | 6 | ! |
---|
| 7 | !-- grid parameters |
---|
| 8 | !------------------------------------------------------------------------------- |
---|
| 9 | nx = 39, ! Number of gridboxes in x-direction (nx+1) |
---|
| 10 | ny = 39, ! Number of gridboxes in y-direction (ny+1) |
---|
| 11 | nz = 40, ! Number of gridboxes in z-direction (nz) |
---|
[4002] | 12 | |
---|
[4224] | 13 | dx = 2.0, ! Size of single gridbox in x-direction |
---|
| 14 | dy = 2.0, ! Size of single gridbox in y-direction |
---|
| 15 | dz = 2.0, ! Size of single gridbox in z-direction |
---|
| 16 | ! |
---|
| 17 | !-- initialization |
---|
| 18 | !------------------------------------------------------------------------------- |
---|
| 19 | initializing_actions = 'set_constant_profiles', ! initial conditions |
---|
[4002] | 20 | |
---|
[4224] | 21 | ug_surface = 1.0, ! u-comp of geostrophic wind at surface |
---|
| 22 | vg_surface = 0.0, ! v-comp of geostrophic wind at surface |
---|
| 23 | |
---|
| 24 | s_surface = 0.0008, ! equivalent to 420ppm CO2 as background |
---|
| 25 | ! concentration |
---|
| 26 | ! |
---|
[4002] | 27 | !-- boundary conditions |
---|
[4224] | 28 | !------------------------------------------------------------------------------- |
---|
| 29 | bc_uv_t = 'neumann', ! free-slip boundary condition |
---|
| 30 | |
---|
| 31 | surface_scalarflux = 0.02, ! kg emission per square meter and second |
---|
| 32 | bc_s_b = 'neumann', ! bottom boundary condition for scalar |
---|
| 33 | ! concentration. Required with scalar |
---|
| 34 | !flux |
---|
| 35 | |
---|
| 36 | roughness_length = 0.1, ! roughness length (in m) for momentum |
---|
| 37 | ! |
---|
| 38 | !-- numerics |
---|
| 39 | !------------------------------------------------------------------------------- |
---|
| 40 | fft_method = 'temperton-algorithm', ! build-in fft method |
---|
| 41 | |
---|
| 42 | momentum_advec = 'ws-scheme', ! default advection scheme |
---|
| 43 | scalar_advec = 'ws-scheme', ! default advection scheme |
---|
| 44 | ! |
---|
| 45 | !-- mode |
---|
| 46 | !------------------------------------------------------------------------------- |
---|
| 47 | neutral = .TRUE., ! strictly neutral flow |
---|
[4002] | 48 | |
---|
[4224] | 49 | passive_scalar = .TRUE., ! switch on prognostic equation for pas. |
---|
| 50 | ! scalar |
---|
[4002] | 51 | |
---|
[4224] | 52 | dp_external = .TRUE., ! use ext. pressure gradient as driving force |
---|
| 53 | dpdxy = -0.0002, 0.0, ! external pressure gradient applied |
---|
| 54 | ! in x- and y-direction (Pa/m) |
---|
| 55 | ! |
---|
| 56 | !-- physics |
---|
| 57 | !------------------------------------------------------------------------------- |
---|
| 58 | omega = 0.0, ! angular velocity of the rotating system, |
---|
| 59 | ! no Coriolis force |
---|
| 60 | ! |
---|
| 61 | !-- topography |
---|
| 62 | !------------------------------------------------------------------------------- |
---|
| 63 | topography = 'read_from_file', ! ASCII Input file contains the |
---|
| 64 | ! arbitrary topography height |
---|
| 65 | ! information in m. |
---|
| 66 | |
---|
| 67 | / ! end of initialization parameter namelist |
---|
[4002] | 68 | |
---|
[4224] | 69 | !------------------------------------------------------------------------------- |
---|
| 70 | !-- RUNTIME PARAMETER NAMELIST |
---|
| 71 | ! Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/d3par |
---|
| 72 | !------------------------------------------------------------------------------- |
---|
[4002] | 73 | &runtime_parameters |
---|
[4224] | 74 | ! |
---|
| 75 | !-- run steering |
---|
| 76 | !------------------------------------------------------------------------------- |
---|
| 77 | end_time = 10800.0, ! simulation time of the 3D model |
---|
[4002] | 78 | |
---|
[4224] | 79 | create_disturbances = .TRUE., ! randomly perturbate horiz. velocity |
---|
| 80 | ! |
---|
| 81 | !-- data output |
---|
| 82 | !------------------------------------------------------------------------------- |
---|
| 83 | skip_time_data_output = 0.0, ! no data output before this time interval |
---|
| 84 | dt_data_output = 3600.0, ! output interval for general data |
---|
[4002] | 85 | |
---|
[4224] | 86 | averaging_interval = 3600.0, ! averaging interval general data |
---|
| 87 | dt_averaging_input = 6.0, ! general data sampling rate for averaging |
---|
[4002] | 88 | |
---|
[4224] | 89 | dt_run_control = 6.0, ! output interval for run control |
---|
| 90 | dt_dots = 6.0, ! output interval for time series |
---|
[4002] | 91 | |
---|
[4224] | 92 | skip_time_dopr = 0.0, ! no profile data output before this time |
---|
| 93 | ! interval |
---|
| 94 | dt_dopr = 3600.0, ! output interval for profile data |
---|
| 95 | averaging_interval_pr = 3600.0, ! averaging interval profile data |
---|
| 96 | dt_averaging_input_pr = 6.0, ! profile data sampling rate for averaging |
---|
[4002] | 97 | |
---|
[4224] | 98 | do3d_at_begin = .TRUE., ! output of 3d volume data at the beginning |
---|
| 99 | do2d_at_begin = .TRUE., ! output of 2d data at the beginning |
---|
[4002] | 100 | |
---|
[4224] | 101 | dt_do2d_xy = 3600.0, ! output interval for xy cross sections |
---|
| 102 | dt_do2d_xz = 3600.0, ! output interval for xz cross sections |
---|
[4002] | 103 | |
---|
[4409] | 104 | data_output_pr = '#u', 'u*2', 'wu', 'w*u*', 'w"u"', '#v', 'v*2', 'wv', 'w*v*', 'w"v"', 'w', 'w*2', 'e', 'e*', '#km', '#l', '#s', 'ws', 'w*s*', 'w"s"', |
---|
[4002] | 105 | |
---|
[4409] | 106 | 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', |
---|
[4002] | 107 | |
---|
[4224] | 108 | section_xy = -1, 21, 22, 23, ! grid index for 2D XY cross |
---|
| 109 | ! sections |
---|
| 110 | section_xz = -1, 17, 18, 19, ! grid index for 2D XZ cross |
---|
| 111 | ! sections |
---|
[4002] | 112 | |
---|
[4224] | 113 | / ! end of runtime parameter namelist |
---|
[4002] | 114 | |
---|
[4224] | 115 | !------------------------------------------------------------------------------- |
---|
| 116 | !-- PARTICLE PARAMETER NAMELIST |
---|
| 117 | ! Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/parpar |
---|
| 118 | !------------------------------------------------------------------------------- |
---|
[4002] | 119 | &particle_parameters |
---|
[4224] | 120 | ! |
---|
| 121 | !-- boundary conditions for particle transport |
---|
| 122 | !------------------------------------------------------------------------------- |
---|
| 123 | bc_par_lr = 'absorb', ! lateral boundary condition |
---|
| 124 | bc_par_b = 'reflect', ! bottom boundary condition |
---|
| 125 | bc_par_t = 'reflect', ! top boundary condition |
---|
| 126 | ! |
---|
| 127 | !-- define point source |
---|
| 128 | !------------------------------------------------------------------------------- |
---|
| 129 | psb = 45.0, ! bottom edge of a particle source |
---|
| 130 | pst = 45.0, ! top edge of a particle source |
---|
| 131 | pss = 40.0, ! front ("south") edge of a particle source |
---|
| 132 | psn = 40.0, ! rear ("north") edge of a particle source |
---|
| 133 | psl = 25.0, ! left edge of a particle source |
---|
| 134 | psr = 25.0, ! right edge of a particle source |
---|
| 135 | |
---|
| 136 | particles_per_point = 100, ! particles per point (default = 1) |
---|
| 137 | |
---|
| 138 | random_start_position = .FALSE., ! add random fluctuations for initial |
---|
| 139 | ! particle position |
---|
| 140 | ! |
---|
| 141 | !-- define time variables of particle relaese |
---|
| 142 | !------------------------------------------------------------------------------- |
---|
| 143 | dt_prel = 1.0, ! time after new particle set is released |
---|
| 144 | |
---|
| 145 | step_dealloc = 50, ! time steps after arrays are deallocated |
---|
| 146 | |
---|
| 147 | particle_maximum_age = 200.0, ! maximum allowed age of particles in s |
---|
| 148 | ! |
---|
| 149 | !-- physics of particles |
---|
| 150 | !------------------------------------------------------------------------------- |
---|
| 151 | use_sgs_for_particles = .TRUE., ! SGS-velocities for particles |
---|
| 152 | ! |
---|
| 153 | !-- particle output |
---|
| 154 | !------------------------------------------------------------------------------- |
---|
| 155 | dt_dopts = 10.0, ! output interval for particle quantities |
---|
| 156 | |
---|
| 157 | / ! end of particle parameter namelist |
---|
[4002] | 158 | |
---|
[4224] | 159 | !------------------------------------------------------------------------------- |
---|
| 160 | !-- USER PARAMETER NAMELIST |
---|
| 161 | ! Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/userpar |
---|
| 162 | !------------------------------------------------------------------------------- |
---|
[4002] | 163 | &user_parameters |
---|
[4224] | 164 | ! |
---|
| 165 | !-- definition of emission stripe |
---|
| 166 | !------------------------------------------------------------------------------- |
---|
| 167 | emission_stripe_width_y = 6.0, ! width of emission stripe in meters |
---|
| 168 | emission_stripe_orig_y = 12.0, 62.0, ! origin of emission stripe in meters |
---|
| 169 | ! |
---|
| 170 | !-- defition of time series output |
---|
| 171 | !------------------------------------------------------------------------------- |
---|
| 172 | s_ts_pos_x = 17.0, 61.0, 75.0, ! x position of scalar time |
---|
| 173 | ! series output |
---|
| 174 | s_ts_pos_y = 25.0, 15.0, 41.0, ! y-position of scalar time |
---|
| 175 | ! series output |
---|
[4002] | 176 | |
---|
[4224] | 177 | / ! end of user parameter namelist |
---|
[4002] | 178 | |
---|
| 179 | |
---|