- Timestamp:
- Sep 10, 2019 12:49:41 PM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/TUTORIALS/cases/dispersion_eulerian_and_lpm_extended/INPUT/dispersion_eulerian_and_lpm_extended_p3d
r4098 r4224 1 !------------------------------------------------------------------------------- 2 !-- INITIALIZATION PARAMETER NAMELIST 3 ! Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/inipar 4 !------------------------------------------------------------------------------- 1 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) 2 12 3 !-- grid setup 4 nx = 39, 5 ny = 39, 6 nz = 40, 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 7 20 8 dx = 2.0, 9 dy = 2.0, 10 dz = 2.0, 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 11 48 12 !-- init-action 13 initializing_actions = 'set_constant_profiles',49 passive_scalar = .TRUE., ! switch on prognostic equation for pas. 50 ! scalar 14 51 15 !-- advection schemes 16 momentum_advec = 'ws-scheme', ! default advection scheme 17 scalar_advec = 'ws-scheme', 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 18 68 19 !-- wind speed setup 20 ug_surface = 1.0, 21 vg_surface = 0.0, 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 22 78 23 !-- temperature setup 24 neutral = .T., ! strictly neutral flow 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 25 85 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', 86 averaging_interval = 3600.0, ! averaging interval general data 87 dt_averaging_input = 6.0, ! general data sampling rate for averaging 31 88 32 !-- boundary conditions 33 bc_uv_t = 'neumann', ! free-slip boundary condition 34 roughness_length = 0.1, ! applies to all surfaces locally 89 dt_run_control = 6.0, ! output interval for run control 90 dt_dots = 6.0, ! output interval for time series 35 91 36 !-- special methods 37 dp_external = .T., ! use horizontal pressure gradient 38 dpdxy = -0.0002, 0.0, ! for forcing 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 39 97 40 !-- special setups 41 omega = 0.0, ! no Coriolis force98 do3d_at_begin = .TRUE., ! output of 3d volume data at the beginning 99 do2d_at_begin = .TRUE., ! output of 2d data at the beginning 42 100 43 topography = 'read_from_file', 44 topography_grid_convention = 'cell_center', ! default 45 / 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 46 192 47 193 48 &runtime_parameters49 50 !-- simulation time51 end_time = 10800.0,52 termination_time_needed = 900.0,53 54 !-- others55 create_disturbances = .T.,56 57 !-- global output settings58 skip_time_data_output = 0.0,59 dt_data_output = 3600.0,60 61 averaging_interval = 3600.0,62 dt_averaging_input = 6.0,63 64 !-- time series / run control65 dt_run_control = 6.0,66 dt_dots = 6.0,67 68 !-- profiles69 skip_time_dopr = 0.0,70 dt_dopr = 3600.0,71 averaging_interval_pr = 3600.0,72 dt_averaging_input_pr = 6.0,73 74 do3d_at_begin = .TRUE.,75 do2d_at_begin = .TRUE.,76 77 dt_do2d_xy = 3600.0,78 dt_do2d_xz = 3600.0,79 80 data_output_pr = '#u', 'u*2', 'wu', 'w*u*', 'w"u"',81 '#v', 'v*2', 'wv', 'w*v*', 'w"v"',82 'w', 'w*2',83 'e', 'e*',84 '#km',85 '#l',86 '#s', 'ws', 'w*s*', 'w"s"',87 88 !-- output variables89 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_parameters104 105 !-- Boundary conditions for particles106 bc_par_lr = 'absorb',107 bc_par_b = 'reflect',108 bc_par_t = 'reflect',109 110 !-- Define point source111 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 released119 dt_prel = 1.0,120 121 !-- Time steps after array are deallocated122 step_dealloc = 50,123 124 !-- Particles per point (default = 1)125 particles_per_point = 100,126 127 !-- Maximum time for particles128 particle_maximum_age = 200.0,129 130 !-- Add random fluctuations for initial particle position131 random_start_position = .FALSE.,132 133 !-- SGS-velocities for particles134 use_sgs_for_particles = .TRUE.,135 136 !-- Particle time series interval137 dt_dopts = 10.0,138 139 /140 141 &user_parameters142 emission_stripe_width_y = 6.0, !in meters143 emission_stripe_orig_y = 12.0, 62.0, !in meters144 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
Note: See TracChangeset
for help on using the changeset viewer.
