| 1 | !------------------------------------------------------------------------------- |
|---|
| 2 | !-- INITIALIZATION PARAMETER NAMELIST |
|---|
| 3 | ! Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/inipar |
|---|
| 4 | !------------------------------------------------------------------------------- |
|---|
| 5 | &initialization_parameters |
|---|
| 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) |
|---|
| 12 | |
|---|
| 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 |
|---|
| 20 | |
|---|
| 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 | ! |
|---|
| 27 | !-- boundary conditions |
|---|
| 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 |
|---|
| 48 | |
|---|
| 49 | passive_scalar = .TRUE., ! switch on prognostic equation for pas. |
|---|
| 50 | ! scalar |
|---|
| 51 | |
|---|
| 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 |
|---|
| 68 | |
|---|
| 69 | !------------------------------------------------------------------------------- |
|---|
| 70 | !-- RUNTIME PARAMETER NAMELIST |
|---|
| 71 | ! Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/d3par |
|---|
| 72 | !------------------------------------------------------------------------------- |
|---|
| 73 | &runtime_parameters |
|---|
| 74 | ! |
|---|
| 75 | !-- run steering |
|---|
| 76 | !------------------------------------------------------------------------------- |
|---|
| 77 | end_time = 10800.0, ! simulation time of the 3D model |
|---|
| 78 | |
|---|
| 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 |
|---|
| 85 | |
|---|
| 86 | averaging_interval = 3600.0, ! averaging interval general data |
|---|
| 87 | dt_averaging_input = 6.0, ! general data sampling rate for averaging |
|---|
| 88 | |
|---|
| 89 | dt_run_control = 6.0, ! output interval for run control |
|---|
| 90 | dt_dots = 6.0, ! output interval for time series |
|---|
| 91 | |
|---|
| 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 |
|---|
| 97 | |
|---|
| 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 |
|---|
| 100 | |
|---|
| 101 | dt_do2d_xy = 3600.0, ! output interval for xy cross sections |
|---|
| 102 | dt_do2d_xz = 3600.0, ! output interval for xz cross sections |
|---|
| 103 | |
|---|
| 104 | data_output_pr = '#u', 'u*2', 'wu', 'w*u*', 'w"u"', |
|---|
| 105 | '#v', 'v*2', 'wv', 'w*v*', 'w"v"', |
|---|
| 106 | 'w', 'w*2', |
|---|
| 107 | 'e', 'e*', |
|---|
| 108 | '#km', |
|---|
| 109 | '#l', |
|---|
| 110 | '#s', 'ws', 'w*s*', 'w"s"', ! profile output |
|---|
| 111 | |
|---|
| 112 | data_output = 'u', 'u_av', |
|---|
| 113 | 'u_xy', 'u_xz', 'u_xz_av', |
|---|
| 114 | 'v', 'v_av', |
|---|
| 115 | 'w', 'w_av', |
|---|
| 116 | 's', 's_av', |
|---|
| 117 | 'ssws*_xy', |
|---|
| 118 | 'pc', 'pc_av', |
|---|
| 119 | 'pc_xy', 'pc_xz', 'pc_xz_av', ! 2d and/or 3d |
|---|
| 120 | ! output |
|---|
| 121 | |
|---|
| 122 | section_xy = -1, 21, 22, 23, ! grid index for 2D XY cross |
|---|
| 123 | ! sections |
|---|
| 124 | section_xz = -1, 17, 18, 19, ! grid index for 2D XZ cross |
|---|
| 125 | ! sections |
|---|
| 126 | |
|---|
| 127 | / ! end of runtime parameter namelist |
|---|
| 128 | |
|---|
| 129 | !------------------------------------------------------------------------------- |
|---|
| 130 | !-- PARTICLE PARAMETER NAMELIST |
|---|
| 131 | ! Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/parpar |
|---|
| 132 | !------------------------------------------------------------------------------- |
|---|
| 133 | &particle_parameters |
|---|
| 134 | ! |
|---|
| 135 | !-- boundary conditions for particle transport |
|---|
| 136 | !------------------------------------------------------------------------------- |
|---|
| 137 | bc_par_lr = 'absorb', ! lateral boundary condition |
|---|
| 138 | bc_par_b = 'reflect', ! bottom boundary condition |
|---|
| 139 | bc_par_t = 'reflect', ! top boundary condition |
|---|
| 140 | ! |
|---|
| 141 | !-- define point source |
|---|
| 142 | !------------------------------------------------------------------------------- |
|---|
| 143 | psb = 45.0, ! bottom edge of a particle source |
|---|
| 144 | pst = 45.0, ! top edge of a particle source |
|---|
| 145 | pss = 40.0, ! front ("south") edge of a particle source |
|---|
| 146 | psn = 40.0, ! rear ("north") edge of a particle source |
|---|
| 147 | psl = 25.0, ! left edge of a particle source |
|---|
| 148 | psr = 25.0, ! right edge of a particle source |
|---|
| 149 | |
|---|
| 150 | particles_per_point = 100, ! particles per point (default = 1) |
|---|
| 151 | |
|---|
| 152 | random_start_position = .FALSE., ! add random fluctuations for initial |
|---|
| 153 | ! particle position |
|---|
| 154 | ! |
|---|
| 155 | !-- define time variables of particle relaese |
|---|
| 156 | !------------------------------------------------------------------------------- |
|---|
| 157 | dt_prel = 1.0, ! time after new particle set is released |
|---|
| 158 | |
|---|
| 159 | step_dealloc = 50, ! time steps after arrays are deallocated |
|---|
| 160 | |
|---|
| 161 | particle_maximum_age = 200.0, ! maximum allowed age of particles in s |
|---|
| 162 | ! |
|---|
| 163 | !-- physics of particles |
|---|
| 164 | !------------------------------------------------------------------------------- |
|---|
| 165 | use_sgs_for_particles = .TRUE., ! SGS-velocities for particles |
|---|
| 166 | ! |
|---|
| 167 | !-- particle output |
|---|
| 168 | !------------------------------------------------------------------------------- |
|---|
| 169 | dt_dopts = 10.0, ! output interval for particle quantities |
|---|
| 170 | |
|---|
| 171 | / ! end of particle parameter namelist |
|---|
| 172 | |
|---|
| 173 | !------------------------------------------------------------------------------- |
|---|
| 174 | !-- USER PARAMETER NAMELIST |
|---|
| 175 | ! Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/userpar |
|---|
| 176 | !------------------------------------------------------------------------------- |
|---|
| 177 | &user_parameters |
|---|
| 178 | ! |
|---|
| 179 | !-- definition of emission stripe |
|---|
| 180 | !------------------------------------------------------------------------------- |
|---|
| 181 | emission_stripe_width_y = 6.0, ! width of emission stripe in meters |
|---|
| 182 | emission_stripe_orig_y = 12.0, 62.0, ! origin of emission stripe in meters |
|---|
| 183 | ! |
|---|
| 184 | !-- defition of time series output |
|---|
| 185 | !------------------------------------------------------------------------------- |
|---|
| 186 | s_ts_pos_x = 17.0, 61.0, 75.0, ! x position of scalar time |
|---|
| 187 | ! series output |
|---|
| 188 | s_ts_pos_y = 25.0, 15.0, 41.0, ! y-position of scalar time |
|---|
| 189 | ! series output |
|---|
| 190 | |
|---|
| 191 | / ! end of user parameter namelist |
|---|
| 192 | |
|---|
| 193 | |
|---|
