1 | <!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN"> |
---|
2 | <html><head> |
---|
3 | <meta content="text/html; charset=ISO-8859-1" http-equiv="content-type"><title>PALM chapter 4.2</title></head> |
---|
4 | <body><h3 style="line-height: 100%;"><a name="Kapitel4.2"></a>4.2 <a href="#Laufparameter">Runtime |
---|
5 | parameters</a> and <a href="#Paketparameter">package |
---|
6 | parameters</a></h3> |
---|
7 | <h3 style="margin-bottom: 0cm; line-height: 100%;"><a name="Laufparameter"></a>Runtime parameters:</h3> |
---|
8 | <br><br><table style="text-align: left; width: 100%;" border="1" cellpadding="2" cellspacing="2"> <tbody> |
---|
9 | <tr> |
---|
10 | <td style="vertical-align: top;"><font size="4"><b>Parameter |
---|
11 | name</b></font></td> <td style="vertical-align: top;"><font size="4"><b>Type</b></font></td> |
---|
12 | <td style="vertical-align: top;"> <p><b><font size="4">Default</font></b> <br> <b><font size="4">value</font></b></p> </td> |
---|
13 | <td style="vertical-align: top;"><font size="4"><b>Explanation</b></font></td> |
---|
14 | </tr> <tr> <td style="vertical-align: top;"><a name="averaging_interval"></a><span style="font-weight: bold;">averaging_interval</span><br> |
---|
15 | </td> <td style="vertical-align: top;">R<br> </td> |
---|
16 | <td style="vertical-align: top;"><span style="font-style: italic;">0.0</span><br> </td> |
---|
17 | <td style="vertical-align: top;">Averaging interval for |
---|
18 | all output of temporally averaged data (in s).<br><br>This |
---|
19 | parameter defines the time interval length for temporally averaged data |
---|
20 | (vertical profiles, spectra, 2d cross-sections, 3d volume data). By |
---|
21 | default, data are not subject to temporal averaging. The |
---|
22 | interval |
---|
23 | length is limited by the parameter <a href="#dt_data_output_av">dt_data_output_av</a>. |
---|
24 | In any case, <span style="font-weight: bold;">averaging_interval</span> |
---|
25 | <= <span style="font-weight: bold;">dt_data_output_av</span> |
---|
26 | must hold.<br><br>If |
---|
27 | an interval is defined, then by default the average is calculated from |
---|
28 | the data values of all timesteps lying within this interval. The number |
---|
29 | of time levels entering into the average can be reduced with the |
---|
30 | parameter <a href="#dt_averaging_input">dt_averaging_input</a>.<br><br>If |
---|
31 | an averaging interval can not be completed at the end of a run, it |
---|
32 | will be finished at the beginning of the next restart run. Thus for |
---|
33 | restart runs, averaging intervals do not |
---|
34 | necessarily begin at the beginning of the run.<br><br>Parameters |
---|
35 | <a href="#averaging_interval_pr">averaging_interval_pr</a> |
---|
36 | and <a href="#averaging_interval_sp">averaging_interval_sp</a> |
---|
37 | can be used to define different averaging intervals for vertical |
---|
38 | profile data and spectra, respectively.<br> </td> </tr> |
---|
39 | <tr> <td style="vertical-align: top;"> <p><a name="averaging_interval_pr"></a><b>averaging_interval_pr</b></p> |
---|
40 | </td> <td style="vertical-align: top;">R<br> </td> |
---|
41 | <td style="vertical-align: top;"><span style="font-style: italic;">value of <a href="#averaging_interval">averaging_<br> |
---|
42 | interval</a><br> |
---|
43 | </span> </td> <td style="vertical-align: top;"><p>Averaging |
---|
44 | interval for output of vertical profiles to |
---|
45 | local |
---|
46 | file <font color="#000000"><font color="#000000"><a href="chapter_3.4.html#DATA_1D_PR_NETCDF">DATA_1D_PR_NETCDF</a> |
---|
47 | </font></font>and/or <a href="chapter_3.4.html#PLOT1D_DATA">PLOT1D_DATA</a> |
---|
48 | (in s). </p> <p>If |
---|
49 | this parameter is given a non-zero value, temporally |
---|
50 | averaged vertical profile data are output. By default, profile data |
---|
51 | data are not subject to temporal averaging. The interval length is |
---|
52 | limited by the parameter <a href="#dt_dopr">dt_dopr</a>. |
---|
53 | In any case <b>averaging_interval_pr</b> <= <b>dt_dopr |
---|
54 | </b>must |
---|
55 | hold.</p>If an interval is defined, then by default the average |
---|
56 | is calculated |
---|
57 | from the data values of all timesteps lying within this interval. The |
---|
58 | number of time levels entering into the average can be reduced with the |
---|
59 | parameter <a href="#dt_averaging_input_pr">dt_averaging_input_pr</a>. |
---|
60 | <p>If |
---|
61 | an averaging interval can not be completed at the end of a run, it will |
---|
62 | be finished at the beginning of the next restart run. Thus for restart |
---|
63 | runs, averaging intervals do not |
---|
64 | necessarily begin at the beginning of the run.</p> </td> </tr> |
---|
65 | <tr> <td style="vertical-align: top;"><a name="call_psolver_at_all_substeps"></a><span style="font-weight: bold;">call_psolver_at_all_<br> |
---|
66 | substeps</span></td> <td style="vertical-align: top;">L<br> |
---|
67 | </td> <td style="vertical-align: top;"><span style="font-style: italic;">.T.</span><br> </td> |
---|
68 | <td style="vertical-align: top;">Switch |
---|
69 | to steer the call of the pressure solver.<br> <br> |
---|
70 | In order to speed-up performance, the Poisson equation for perturbation |
---|
71 | pressure (see <a href="#psolver">psolver</a>) can |
---|
72 | be called only at the last substep of multistep Runge-Kutta |
---|
73 | timestep schemes (see <a href="chapter_4.1.html#timestep_scheme">timestep_scheme</a>) |
---|
74 | by setting <span style="font-weight: bold;">call_psolver_at_all_substeps</span> |
---|
75 | = <span style="font-style: italic;">.F.</span>. |
---|
76 | In many cases, this sufficiently reduces the divergence of the velocity |
---|
77 | field. Nevertheless, small-scale ripples (2-delta-x) may occur. In this |
---|
78 | case and in case |
---|
79 | of non-cyclic lateral boundary conditions, <span style="font-weight: bold;">call_psolver_at_all_timesteps</span> |
---|
80 | = <span style="font-style: italic;">.T.</span> |
---|
81 | should be used. <span style="font-weight: bold;"></span></td> |
---|
82 | </tr> <tr> <td style="vertical-align: top;"><p><a name="fcl_factor"></a><b>cfl_factor</b></p> |
---|
83 | </td> <td style="vertical-align: top;">R<br> </td> |
---|
84 | <td style="vertical-align: top;"> <p><i>0.1, |
---|
85 | 0.8 or 0.9</i> <br> <i>(see right)</i></p> |
---|
86 | </td> <td style="vertical-align: top;"> <p lang="en-GB">Time step limiting factor. </p> |
---|
87 | <p><span lang="en-GB">In the model, the <span lang="en-GB">maximum |
---|
88 | allowed </span>time step according to CFL and |
---|
89 | diffusion-criterion |
---|
90 | dt_max is reduced by </span><a href="chapter_4.1.html#dt"><span lang="en-GB">dt</span></a> <span lang="en-GB">= |
---|
91 | dt_max * <b>cfl_factor</b> |
---|
92 | in order to avoid stability problems which may arise in the vicinity of |
---|
93 | the maximum allowed timestep. The condition <i>0.0</i> |
---|
94 | < <b>cfl_factor</b> |
---|
95 | < <i>1.0 </i>applies.<br> </span></p> |
---|
96 | <p><span lang="en-GB">The default value of |
---|
97 | cfl_factor depends on |
---|
98 | the </span><a href="chapter_4.1.html#timestep_scheme"><span lang="en-GB">timestep_scheme</span></a><span lang="en-GB"> used:<br> </span></p> <p><span lang="en-GB">For the third order Runge-Kutta scheme it |
---|
99 | is <b>cfl_factor</b> = </span><span style="font-style: italic;">0.9</span><span lang="en-GB">.<br> </span></p> <p><span lang="en-GB">In case of the leapfrog scheme a quite |
---|
100 | restrictive value of <span style="font-weight: bold;">cfl_factor</span> |
---|
101 | = <span style="font-style: italic;">0.1 </span></span><span lang="en-GB">is used because for larger values the velocity |
---|
102 | divergence |
---|
103 | significantly effects the accuracy of the model results.</span><a href="chapter_4.1.html#scalar_advec"><span lang="en-GB"></span></a><span lang="en-GB"> Possibly larger values may |
---|
104 | be used with the leapfrog scheme but these are to be determined by |
---|
105 | appropriate test runs.<span style="font-family: times new roman;"><br> |
---|
106 | </span></span></p> <span lang="en-GB"><span style="font-family: times new roman;"></span><font face="Times New Roman">The default value for the Euler |
---|
107 | scheme is <span style="font-weight: bold;">cfl_factor</span> |
---|
108 | = <span style="font-style: italic;">0.8</span> .</font></span></td> |
---|
109 | </tr><tr> <td style="vertical-align: top;"> <p><a name="create_disturbances"></a><b>create_disturbances</b></p> |
---|
110 | </td> <td style="vertical-align: top;">L<br> </td> |
---|
111 | <td style="vertical-align: top;"><span style="font-style: italic;">.T.</span><br> </td> |
---|
112 | <td style="vertical-align: top;"> <p>Switch to |
---|
113 | impose random perturbations to the horizontal |
---|
114 | velocity field. </p> <p>With <b>create_disturbances</b> |
---|
115 | = <i>.T.,</i> random |
---|
116 | perturbations can be imposed to the horizontal velocity field at |
---|
117 | certain times e.g. in order to trigger off the onset of convection, |
---|
118 | etc..<br> </p> <p>The temporal interval between |
---|
119 | these times can be steered with <a href="#dt_disturb">dt_disturb</a>, |
---|
120 | the vertical range of the perturbations with <a href="#disturbance_level_b">disturbance_level_b</a> |
---|
121 | and <a href="#disturbance_level_t">disturbance_level_t</a>, |
---|
122 | and the perturbation amplitude with <a href="#disturbance_amplitude">disturbance_amplitude</a>. |
---|
123 | In case of non-cyclic lateral boundary conditions (see <a href="chapter_4.1.html#bc_lr">bc_lr</a> |
---|
124 | and <a href="chapter_4.1.html#bc_ns">bc_ns</a>), |
---|
125 | the horizontal range of the perturbations is determined by <a href="chapter_4.1.html#inflow_disturbance_begin">inflow_disturbance_begin</a> |
---|
126 | and <a href="chapter_4.1.html#inflow_disturbance_end">inflow_disturbance_end</a>. |
---|
127 | A perturbation is added to each grid point with its individual value |
---|
128 | determined by multiplying the disturbance amplitude with a uniformly |
---|
129 | distributed random number. |
---|
130 | After this, the arrays of u and v are smoothed by applying a |
---|
131 | Shuman-filter twice and made divergence-free by applying the pressure |
---|
132 | solver.<br> </p> <p>The random number generator to |
---|
133 | be used can be chosen with <a href="chapter_4.1.html#random_generator">random_generator</a>.<br> |
---|
134 | </p> <p>As soon as the desired flow features have |
---|
135 | developed |
---|
136 | (e.g. convection has started), further imposing of |
---|
137 | perturbations |
---|
138 | is not necessary and can be omitted (this does not hold for non-cyclic |
---|
139 | lateral boundaries!). This can be steered by assigning |
---|
140 | an upper limit value for the perturbation energy (the perturbation |
---|
141 | energy is defined by the deviation of the velocity from the mean flow) |
---|
142 | using the parameter <a href="#disturbance_energy_limit">disturbance_energy_limit</a>. |
---|
143 | As soon as the perturbation energy has exceeded this energy limit, no |
---|
144 | more random perturbations are assigned<br> |
---|
145 | . <br> |
---|
146 | Timesteps where a random perturbation has been imposed are marked in |
---|
147 | the local file <a href="chapter_3.4.html#RUN_CONTROL">RUN_CONTROL</a> |
---|
148 | by the character "D" appended to the values of the maximum horizontal |
---|
149 | velocities. </p> </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="cross_normalized_x"></a><b>cross_normalized_x</b></p> |
---|
150 | </td> <td style="vertical-align: top;">C*10 |
---|
151 | <br> (100)</td> <td style="vertical-align: top;"><i>100 * ' '</i></td> |
---|
152 | <td style="vertical-align: top;"> <p>Type of |
---|
153 | normalization applied to the x-coordinate of vertical |
---|
154 | profiles to be plotted with <span style="font-weight: bold;">profil</span>.</p> |
---|
155 | <p>This parameter only applies for <a href="chapter_4.2.html#data_output_format">data_output_format</a> |
---|
156 | = <span style="font-style: italic;">'profil'</span>.</p><p>If |
---|
157 | vertical profiles are plotted with the plot software <span style="font-weight: bold;">profil</span> (data on |
---|
158 | local file <a href="chapter_3.4.html#PLOT1D_DATA">PLOT1D_DATA</a>, |
---|
159 | parameters on local file <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/PLOT1D_PAR">PLOT1D_PAR</a>) |
---|
160 | the x-values of the data points can be normalized with respect to |
---|
161 | certain quantities (e.g. the near-surface heat flux) in order to ensure |
---|
162 | a better comparability. This type of normalization then applies to all |
---|
163 | profiles of one coordinate system (panel). The normalization quantities |
---|
164 | are re-calculated for each output time of each individual profile. If |
---|
165 | temporally averaged profiles are output (see <a href="#averaging_interval_pr">averaging_interval_pr</a>), |
---|
166 | then the normalization quantities are also temporally averaged |
---|
167 | accordingly. If the value of a normalization quantity becomes zero, |
---|
168 | then normalization for the total respective coordinate system (panel) |
---|
169 | is switched off automatically (also for the y-axis).<br> </p> |
---|
170 | <p>By default, the normalization quantities are calculated as the |
---|
171 | horizontal mean of the total model domain and and these values are also |
---|
172 | used for the normalization of profiles from subdomains (see <a href="chapter_4.1.html#statistic_regions">statistic_regions</a>). |
---|
173 | Instead of this, they can be calculated from the data of a certain |
---|
174 | subdomain by using the parameter <a href="#normalizing_region">normalizing_region</a> |
---|
175 | (however, they are used again for all subdomains and even for the total |
---|
176 | domain). </p> <p>The user can choose between |
---|
177 | the following normalization |
---|
178 | quantities: <br> </p> <table style="text-align: left; width: 100%;" cellpadding="2" cellspacing="2"> <tbody> <tr> <td style="vertical-align: top;"><i>'wpt0'</i></td> |
---|
179 | <td style="vertical-align: top;">Normalization with |
---|
180 | respect |
---|
181 | to the total surface sensible heat |
---|
182 | flux (k=0 ).</td> </tr> <tr> <td style="vertical-align: middle;"><i>'ws2'</i></td> |
---|
183 | <td style="vertical-align: top;">Normalization with |
---|
184 | respect |
---|
185 | to w<sub>*</sub> <sup>2</sup> |
---|
186 | (square of the characteristic vertical wind speed of the CBL)</td> |
---|
187 | </tr> <tr> <td style="vertical-align: top;"><i>'tsw2'</i></td> |
---|
188 | <td style="vertical-align: top;">Normalization with |
---|
189 | respect |
---|
190 | to the square of the characteristic |
---|
191 | temperature of the CBL theta<sub>*</sub> (this is defined |
---|
192 | as ratio of |
---|
193 | the surface heat flux and w<sub>*</sub>).</td> </tr> |
---|
194 | <tr> <td style="vertical-align: middle;"><i>'ws3'</i></td> |
---|
195 | <td style="vertical-align: top;">Normalization with |
---|
196 | respect |
---|
197 | to w<sub>*</sub> <sup>3</sup>.</td> </tr> |
---|
198 | <tr> <td style="vertical-align: middle;"><i>'ws2tsw'</i></td> |
---|
199 | <td style="vertical-align: top;">Normalization with |
---|
200 | respect |
---|
201 | to w<sub>*</sub><sup>2</sup>theta<sub>*</sub> |
---|
202 | (for definition of theta<sub>*</sub> see <span style="font-style: italic;">'tsw2'</span>).</td> |
---|
203 | </tr> <tr> <td style="vertical-align: middle;"><i>'wstsw2'</i></td> |
---|
204 | <td style="vertical-align: top;">Normalization with |
---|
205 | respect |
---|
206 | to w<sub>*</sub><sup>2 </sup>theta<sub>*</sub> |
---|
207 | (for definition of theta<sub>*</sub> see <span style="font-style: italic;">'tsw2'</span>).</td> |
---|
208 | </tr> </tbody> </table> <p>For each |
---|
209 | coordinate system (panel) to be drawn (see <a href="#cross_profiles">cross_profiles</a>) |
---|
210 | an individual normalization quantity can be assigned. For example: if <span style="font-weight: bold;">cross_normalized_x</span> = |
---|
211 | <span style="font-style: italic;">'ws2'</span><i>,'ws3'</i>, |
---|
212 | then the |
---|
213 | x-values in the 1st coordinate system are normalized with respect to w<sub>*</sub><sup>2</sup> |
---|
214 | and in the 2nd system with respect to w<sub>*</sub><sup>3</sup>. |
---|
215 | Data |
---|
216 | of the further coordinate systems (if any are to be drawn) are not |
---|
217 | normalized. </p> <p>Using a normalization |
---|
218 | leaves all vertical profile data on |
---|
219 | local file <a href="chapter_3.4.html#PLOT1D_DATA">PLOT1D_DATA</a> |
---|
220 | unaffected, it only affects the visualization. Within <span style="font-weight: bold;">profil</span>, the |
---|
221 | normalization is steered |
---|
222 | by parameter <a href="http://www.muk.uni-hannover.de/institut/software/profil_beschreibung.html#NORMX">normx</a> |
---|
223 | which may be changed subsequently by the user in the parameter file |
---|
224 | (local file <a href="chapter_3.4.html#PLOT1D_PAR">PLOT1D_PAR</a>).<br> |
---|
225 | <br> |
---|
226 | The assigned normalization quantity is noted in the axes labels of the |
---|
227 | respective coordinate systems (see <a href="#cross_xtext">cross_xtext</a>).</p> |
---|
228 | </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="cross_normalized_y"></a><b>cross_normalized_y</b></p> |
---|
229 | </td> <td style="vertical-align: top;">C*10 |
---|
230 | <br> (100)</td> <td style="vertical-align: top;"><i>100 * ' '</i></td> |
---|
231 | <td style="vertical-align: top;"> <p>Type of |
---|
232 | normalization applied to the y-coordinate of vertical |
---|
233 | profiles to be plotted with <span style="font-weight: bold;">profil</span>. </p> |
---|
234 | <p>This parameter only applies for <a href="chapter_4.2.html#data_output_format">data_output_format</a> |
---|
235 | = <span style="font-style: italic;">'profil'</span>.</p><p>If |
---|
236 | vertical profiles are plotted with the plot software <span style="font-weight: bold;">profil</span> (data on |
---|
237 | local file <a href="chapter_3.4.html#PLOT1D_DATA">PLOT1D_DATA</a>, |
---|
238 | parameter on local file <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/PLOT1D_PAR">PLOT1D_PAR</a>) |
---|
239 | the y-values of the data points can be normalized with respect to |
---|
240 | certain quantities (at present only the normalization with respect to |
---|
241 | the boundary layer height is possible) in order to ensure a better |
---|
242 | comparability. </p> <p>The user can choose between the |
---|
243 | following normalization |
---|
244 | quantities: <br> </p> <table style="text-align: left; width: 100%;" cellpadding="2" cellspacing="2"> <tbody> <tr> <td style="vertical-align: top;"><i>'z_i'</i></td> |
---|
245 | <td style="vertical-align: top;">Normalization with |
---|
246 | respect |
---|
247 | to the boundary layer height |
---|
248 | (determined from the height where the heat flux achieves its minimum |
---|
249 | value).</td> </tr> </tbody> </table> <p>For |
---|
250 | further explanations see <a href="#cross_normalized_x">cross_normalized_x.</a></p> |
---|
251 | </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="cross_profiles"></a><b>cross_profiles</b></p> |
---|
252 | </td> <td style="vertical-align: top;">C*100 |
---|
253 | <br> (100)</td> <td style="vertical-align: top;">see right<br> </td> |
---|
254 | <td style="vertical-align: top;"> <p>Determines |
---|
255 | which vertical profiles are to be presented in |
---|
256 | which coordinate system if the plot software <span style="font-weight: bold;">profil</span> is used. |
---|
257 | </p> <p>This parameter only applies for |
---|
258 | <a href="chapter_4.2.html#data_output_format">data_output_format</a> |
---|
259 | = <span style="font-style: italic;">'profil'</span>.</p><p>The |
---|
260 | default assignment is: </p> <p><b>cross_profiles</b> |
---|
261 | = </p> <ul> <p><span style="font-family: monospace; font-style: italic;">' |
---|
262 | u v ',</span><br> <span style="font-family: monospace; font-style: italic;">' pt |
---|
263 | ', </span><br style="font-family: monospace; font-style: italic;"> <span style="font-family: monospace; font-style: italic;">' |
---|
264 | w"pt" w*pt* w*pt*BC wpt wptBC ', </span><br style="font-family: monospace; font-style: italic;"> <span style="font-family: monospace; font-style: italic;">' |
---|
265 | w"u" w*u* wu w"v"w*v* wv ', </span><br style="font-family: monospace; font-style: italic;"> <span style="font-family: monospace; font-style: italic;">' km |
---|
266 | kh ',</span><br style="font-family: monospace; font-style: italic;"> <span style="font-family: monospace; font-style: italic;">' l ' |
---|
267 | ,</span><br> |
---|
268 | 14 * <span style="font-family: monospace; font-style: italic;">' |
---|
269 | '</span></p> </ul> <p>If plot output of |
---|
270 | vertical profiles is produced (see <a href="#data_output_pr">data_output_pr</a>), |
---|
271 | the appropriate data are written to the local file <a href="chapter_3.4.html#PLOT1D_DATA">PLOT1D_DATA</a>. |
---|
272 | Simultaneously, the model produces a parameter file (local name <a href="chapter_3.4.html#PLOT1D_PAR">PLOT1D_PAR</a>) |
---|
273 | which describes the layout for a plot to be generated with the plot |
---|
274 | program <span style="font-weight: bold;">profil</span>. |
---|
275 | The parameter <b>cross_profiles</b> |
---|
276 | determines how many coordinate systems (panels) the plot contains and |
---|
277 | which profiles are supposed to be drawn into which panel. <b>cross_profiles</b> |
---|
278 | expects a character string (up to 100 characters long) for each |
---|
279 | coordinate system, which consists of the names of the profiles to be |
---|
280 | drawn into this system (all available profiles and their respective |
---|
281 | names are described at parameter <a href="#data_output_pr">data_output_pr</a>). |
---|
282 | The single names have to be separated by one blank (' ') and a blank |
---|
283 | must be spent also at the beginning and at the end of the |
---|
284 | string. </p> <p>Example: </p> <ul> |
---|
285 | <p><b>cross_profiles</b> = <span style="font-family: monospace; font-style: italic;">' u v ', |
---|
286 | ' pt '</span></p> </ul> <p>In this case the |
---|
287 | plot consists of two coordinate systems |
---|
288 | (panels) with the first panel containing the profiles of the horizontal |
---|
289 | velocity components (<span style="font-style: italic;">'u'</span> |
---|
290 | and <span style="font-style: italic;">'v'</span>) |
---|
291 | of all output times (see <a href="#dt_dopr">dt_dopr</a>) |
---|
292 | and the second one containing the profiles of the potential temperature |
---|
293 | (<span style="font-style: italic;">'pt'</span>).<br> |
---|
294 | </p> <p>Whether the coordinate systems are actually drawn, |
---|
295 | depends on |
---|
296 | whether data of the appropriate profiles were output during the run |
---|
297 | (profiles to be output have to be selected with the parameter <a href="#data_output_pr">data_output_pr</a>). |
---|
298 | For example if <b>data_output_pr</b> = <span style="font-style: italic;">'u'</span>, <span style="font-style: italic;">'v'</span> was assigned, |
---|
299 | then |
---|
300 | the plot only consists of one panel, since no profiles of the potential |
---|
301 | temperature were output. On the other hand, if profiles were assigned |
---|
302 | to <b>data_output_pr </b>whose names do not appear in <b>cross_profiles</b>, |
---|
303 | then the respective profile data are output (<a href="chapter_3.4.html#PLOT1D_DATA">PLOT1D_DATA</a>) |
---|
304 | but they are not drawn in the plot. <br> </p> |
---|
305 | The arrangement of the panels in the plot can be controlled |
---|
306 | with the parameters <a href="#profile_columns">profile_columns</a> |
---|
307 | and <a href="#profile_rows">profile_rows</a>. |
---|
308 | Up to 100 panels systems are allowed in a plot (however, they may be |
---|
309 | distributed on several pages).</td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="cross_xtext"></a><b>cross_xtext</b></p> |
---|
310 | </td> <td style="vertical-align: top;">C*40 |
---|
311 | <br> (100)</td> <td style="vertical-align: top;">see right<br> </td> |
---|
312 | <td style="vertical-align: top;"> <p>x-axis labels |
---|
313 | of vertical profile coordinate systems to be |
---|
314 | plotted with <span style="font-weight: bold;">profil</span>. |
---|
315 | </p> <p>This parameter only applies for <a href="chapter_4.2.html#data_output_format">data_output_format</a> |
---|
316 | = <span style="font-style: italic;">'profil'</span>.</p><p>The |
---|
317 | default assignment is: </p> <p><b>cross_xtext</b> |
---|
318 | = </p> <ul> <p><span style="font-style: italic;">'wind speed in |
---|
319 | ms>->1', </span><br style="font-style: italic;"> <span style="font-style: italic;">'pot. temperature in |
---|
320 | K', </span><br style="font-style: italic;"> |
---|
321 | <span style="font-style: italic;">'heat flux in K |
---|
322 | ms>->1', </span><br style="font-style: italic;"> <span style="font-style: italic;">'momentum flux in |
---|
323 | m>2s>2', </span><br style="font-style: italic;"> <span style="font-style: italic;">'eddy diffusivity in |
---|
324 | m>2s>->1', </span><br style="font-style: italic;"> <span style="font-style: italic;">'mixing length in m',</span> |
---|
325 | <br>14 * <span style="font-style: italic;">' '</span></p> |
---|
326 | </ul> <p>This parameter can be used to assign x-axis |
---|
327 | labels to vertical |
---|
328 | profiles to be plotted with the plot software <span style="font-weight: bold;">profil </span>(for output |
---|
329 | of vertical |
---|
330 | profile data see <a href="#data_output_pr">data_output_pr</a>).<br> |
---|
331 | The labels are assigned to those coordinate systems (panels) defined by |
---|
332 | <a href="#cross_profiles">cross_profiles</a> |
---|
333 | according to their respective order (compare the default values of <b>cross_xtext</b> |
---|
334 | and <b>cross_profiles</b>). </p> <p>Umlauts |
---|
335 | are possible (write “ in front of, similar to TeX), as |
---|
336 | well as super- and subscripts (use ">" or "<" in front of |
---|
337 | each |
---|
338 | character), special characters etc. (see UNIRAS manuals) when using the |
---|
339 | plot software <a href="http://www.muk.uni-hannover.de/institut/software/profil_beschreibung.html#chapter3.2.6">profil</a>.</p> |
---|
340 | </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="cycle_mg"></a><b>cycle_mg</b></p> |
---|
341 | </td> <td style="vertical-align: top;">C*1</td> |
---|
342 | <td style="vertical-align: top;"><i>'w'</i></td> |
---|
343 | <td style="vertical-align: top;"> <p>Type of cycle |
---|
344 | to be used with the multi-grid method. </p> <p>This |
---|
345 | parameter determines which type of cycle is applied in |
---|
346 | the multi-grid method used for solving the Poisson equation for |
---|
347 | perturbation pressure (see <a href="#psolver">psolver</a>). |
---|
348 | It defines in which way it is switched between the fine and coarse |
---|
349 | grids. So-called v- and w-cycles are realized (i.e. <b>cycle_mg</b> |
---|
350 | may be assigned the values <i>'v'</i> or <i>'w'</i>). |
---|
351 | The |
---|
352 | computational cost of w-cycles is much higher than that of v-cycles, |
---|
353 | however, w-cycles give a much better convergence. </p> </td> |
---|
354 | </tr> <tr> <td style="vertical-align: top;"><p><a name="data_output"></a><b>data_output</b></p> |
---|
355 | </td> <td style="vertical-align: top;">C * 10 (100)<br> |
---|
356 | </td> <td style="vertical-align: top;"><span style="font-style: italic;">100 * ' '</span><br> |
---|
357 | </td> <td style="vertical-align: top;">Quantities |
---|
358 | for which 2d cross section and/or 3d volume data are to be output.<br><br>PALM |
---|
359 | allows the output of instantaneous data as well as of temporally |
---|
360 | averaged data which is steered by the strings assigned to this |
---|
361 | parameter (see below).<br><br>By default, cross section |
---|
362 | data are output (depending on the selected cross sections(s), see |
---|
363 | below) to local files <a href="chapter_3.4.html#DATA_2D_XY_NETCDF">DATA_2D_XY_NETCDF</a>, |
---|
364 | <a href="chapter_3.4.html#DATA_2D_XZ_NETCDF">DATA_2D_XZ_NETCDF</a> |
---|
365 | and/or <a href="chapter_3.4.html#DATA_2D_YZ_NETCDF">DATA_2D_YZ_NETCDF</a>. |
---|
366 | Volume data are output to file <a href="chapter_3.4.html#DATA_3D_NETCDF">DATA_3D_NETCDF</a>. |
---|
367 | If the user has switched on the output of temporally averaged data, |
---|
368 | these are written seperately to local files <a href="chapter_3.4.html#DATA_2D_XY_AV_NETCDF">DATA_2D_XY_AV_NETCDF</a>, |
---|
369 | <a href="chapter_3.4.html#DATA_2D_XZ_AV_NETCDF">DATA_2D_XZ_AV_NETCDF</a>, |
---|
370 | <a href="chapter_4.3.html#DATA_2D_YZ_AV_NETCDF">DATA_2D_YZ_AV_NETCDF</a>, |
---|
371 | and <a href="chapter_3.4.html#DATA_3D_AV_NETCDF">DATA_3D_AV_NETCDF</a>, |
---|
372 | respectively.<br><br>The |
---|
373 | filenames already suggest that all files have NetCDF format. |
---|
374 | Informations about the file content (kind of quantities, array |
---|
375 | dimensions and grid coordinates) are part of the self describing NetCDF |
---|
376 | format and can be extracted from the NetCDF files using the command |
---|
377 | "ncdump -c <filename>". See chapter <a href="chapter_4.5.1.html">4.5.1</a> about processing |
---|
378 | the PALM NetCDF data.<br><br>The following quantities are |
---|
379 | available for output by default:<br><br><table style="text-align: left; width: 576px; height: 481px;" border="1" cellpadding="2" cellspacing="2"><tbody><tr><td style="width: 106px;"><span style="font-weight: bold;">quantity |
---|
380 | name</span></td><td style="width: 196px;"><span style="font-weight: bold;">meaning</span></td><td><span style="font-weight: bold;">unit</span></td><td><span style="font-weight: bold;">remarks</span></td></tr><tr><td style="width: 106px;"><span style="font-style: italic;">e</span></td><td style="width: 196px;">SGS TKE</td><td>m<sup>2</sup>/s<sup>2</sup></td><td></td></tr><tr><td style="width: 106px; vertical-align: top;"><span style="font-style: italic;">lwp*</span></td><td style="width: 196px; vertical-align: top;">liquid water path</td><td style="vertical-align: top;">m</td><td style="vertical-align: top;">only horizontal cross section |
---|
381 | is allowed, requires <a href="chapter_4.1.html#cloud_physics">cloud_physics</a> |
---|
382 | = <span style="font-style: italic;">.TRUE.</span></td></tr><tr><td style="width: 106px; vertical-align: top;"><span style="font-style: italic;">p</span></td><td style="width: 196px; vertical-align: top;">perturpation |
---|
383 | pressure</td><td style="vertical-align: top;">N/m<sup>2</sup>, |
---|
384 | Pa</td><td style="vertical-align: top;"></td></tr><tr><td style="width: 106px; vertical-align: top;"><span style="font-style: italic;">pc</span></td><td style="width: 196px; vertical-align: top;">particle/droplet |
---|
385 | concentration</td><td style="vertical-align: top;">#/gridbox</td><td style="vertical-align: top;">requires that particle |
---|
386 | advection is switched on by <span style="font-weight: bold;">mrun</span>-option |
---|
387 | "-p particles"</td></tr><tr><td style="width: 106px; vertical-align: top;"><span style="font-style: italic;">pr</span></td><td style="width: 196px; vertical-align: top;">mean |
---|
388 | particle/droplet radius </td><td style="vertical-align: top;">m</td><td style="vertical-align: top;">requires that particle |
---|
389 | advection is switched on by <span style="font-weight: bold;">mrun</span>-option |
---|
390 | "-p particles"</td></tr><tr><td style="width: 106px; vertical-align: top;"><span style="font-style: italic;">pt</span></td><td style="width: 196px; vertical-align: top;">potential |
---|
391 | temperature<br></td><td style="vertical-align: top;">K</td><td style="vertical-align: top;"></td></tr><tr><td style="width: 106px; vertical-align: top;"><span style="font-style: italic;">q</span></td><td style="width: 196px; vertical-align: top;">specific humidity |
---|
392 | (or total water content, if cloud physics is switched on)</td><td style="vertical-align: top;">kg/kg</td><td style="vertical-align: top;">requires <a href="chapter_4.1.html#moisture">moisture</a> = <span style="font-style: italic;">.TRUE.</span></td></tr><tr><td style="width: 106px; vertical-align: top;"><span style="font-style: italic;">ql</span></td><td style="width: 196px; vertical-align: top;">liquid water |
---|
393 | content</td><td style="vertical-align: top;">kg/kg</td><td style="vertical-align: top;">requires <a href="chapter_4.1.html#cloud_physics">cloud_physics</a> |
---|
394 | = <span style="font-style: italic;">.TRUE.</span> |
---|
395 | or <a href="chapter_4.1.html#cloud_droplets">cloud_droplets</a> |
---|
396 | = <span style="font-style: italic;">.TRUE.</span></td></tr><tr><td style="width: 106px; vertical-align: top;"><span style="font-style: italic;">ql_c</span></td><td style="width: 196px; vertical-align: top;">change in liquid |
---|
397 | water content due to condensation/evaporation during last timestep</td><td style="vertical-align: top;">kg/kg</td><td style="vertical-align: top;">requires <a href="chapter_4.1.html#cloud_droplets">cloud_droplets</a> |
---|
398 | = <span style="font-style: italic;">.TRUE.</span></td></tr><tr><td style="width: 106px; vertical-align: top;"><span style="font-style: italic;">ql_v</span></td><td style="width: 196px; vertical-align: top;">volume of liquid |
---|
399 | water</td><td style="vertical-align: top;">m<sup>3</sup>/gridbox</td><td style="vertical-align: top;">requires <a href="chapter_4.1.html#cloud_droplets">cloud_droplets</a> |
---|
400 | = <span style="font-style: italic;">.TRUE.</span></td></tr><tr><td style="width: 106px; vertical-align: top;"><span style="font-style: italic;">ql_vp</span></td><td style="width: 196px; vertical-align: top;">weighting factor</td><td style="vertical-align: top;"></td><td style="vertical-align: top;">requires <a href="chapter_4.1.html#cloud_droplets">cloud_droplets</a> |
---|
401 | = <span style="font-style: italic;">.TRUE.</span></td></tr><tr><td style="width: 106px; vertical-align: top;"><span style="font-style: italic;">qv</span></td><td style="width: 196px; vertical-align: top;">water vapor |
---|
402 | content (specific humidity)</td><td style="vertical-align: top;">kg/kg</td><td style="vertical-align: top;">requires <a href="chapter_4.1.html#cloud_physics">cloud_physics</a> |
---|
403 | = <span style="font-style: italic;">.TRUE.</span></td></tr><tr><td style="width: 106px; vertical-align: top;"><span style="font-style: italic;">s</span></td><td style="width: 196px; vertical-align: top;">concentration of |
---|
404 | the scalar</td><td style="vertical-align: top;">1/m<sup>3</sup></td><td style="vertical-align: top;">requires <a href="chapter_4.1.html#passive_scalar">passive_scalar</a> |
---|
405 | = <span style="font-style: italic;">.TRUE.</span></td></tr><tr><td style="width: 106px; vertical-align: top;"><span style="font-style: italic;">t*</span></td><td style="width: 196px; vertical-align: top;">(near surface) |
---|
406 | characteristic temperature</td><td style="vertical-align: top;">K</td><td style="vertical-align: top;">only horizontal cross section |
---|
407 | is allowed</td></tr><tr><td style="width: 106px; vertical-align: top;"><span style="font-style: italic;">u</span></td><td style="width: 196px; vertical-align: top;">u-component of |
---|
408 | the velocity</td><td style="vertical-align: top;">m/s</td><td style="vertical-align: top;"></td></tr><tr><td style="width: 106px; vertical-align: top;"><span style="font-style: italic;">u*</span></td><td style="width: 196px; vertical-align: top;">(near surface) |
---|
409 | friction velocity</td><td style="vertical-align: top;">m/s</td><td style="vertical-align: top;">only horizontal cross section |
---|
410 | is allowed</td></tr><tr><td style="width: 106px; vertical-align: top;"><span style="font-style: italic;">v</span></td><td style="width: 196px; vertical-align: top;">v-component of |
---|
411 | the velocity</td><td style="vertical-align: top;">m/s</td><td style="vertical-align: top;"></td></tr><tr><td style="width: 106px; vertical-align: top;"><span style="font-style: italic;">vpt</span></td><td style="width: 196px; vertical-align: top;">virtual potential |
---|
412 | temperature</td><td style="vertical-align: top;">K</td><td style="vertical-align: top;">requires <a href="chapter_4.1.html#moisture">moisture</a> = <span style="font-style: italic;">.TRUE.</span></td></tr><tr><td style="width: 106px; vertical-align: top;"><span style="font-style: italic;">w</span></td><td style="width: 196px; vertical-align: top;">w-component of |
---|
413 | the velocity</td><td style="vertical-align: top;">m/s</td><td style="vertical-align: top;"></td></tr></tbody></table><br>Multiple |
---|
414 | quantities can be assigned, e.g. <span style="font-weight: bold;">data_output</span> |
---|
415 | = <span style="font-style: italic;">'e'</span>, <span style="font-style: italic;">'u'</span>, <span style="font-style: italic;">'w'</span>.<br><br>By |
---|
416 | assigning the pure strings from the above table, 3d volume data is |
---|
417 | output. Cross section data can be output by appending the string <span style="font-style: italic;">'_xy'</span>, <span style="font-style: italic;">'_xz'</span>, or <span style="font-style: italic;">'_yz'</span> to the |
---|
418 | respective quantities. Time averaged output is created by |
---|
419 | appending the string <span style="font-style: italic;">'_av' |
---|
420 | </span>(for |
---|
421 | cross section data, this string must be appended after the cross |
---|
422 | section string). Cross section data can also be (additionally) averaged |
---|
423 | along the direction normal to the respective section (see below). |
---|
424 | Assignments of quantities can be given in arbitrary |
---|
425 | order:<br><br>Example:<br><br><div style="margin-left: 40px;"><span style="font-weight: bold;">data_output</span> = <span style="font-style: italic;">'u'</span>, <span style="font-style: italic;">'pt_xz_av'</span>, <span style="font-style: italic;">'w_xy'</span>, <span style="font-style: italic;">'u_av'</span>.<br></div><br>This |
---|
426 | example will create the following output: instantaneous 3d volume data |
---|
427 | of u-velocity component (by default on file DATA_3D_NETCDF), temporally |
---|
428 | averaged 3d volume data of u-velocity component (by default on file |
---|
429 | DATA_3D_AV_NETCDF), instantaneous horizontal cross section data of |
---|
430 | w-velocity component (by default on file DATA_2D_XY_NETCDF), and |
---|
431 | temporally averaged vertical cross section data of potential |
---|
432 | temperature (by default on file DATA_2D_XZ_AV_NETCDF).<br><br>The |
---|
433 | user is allowed to extend the above list of quantities by defining his |
---|
434 | own output quantities (see the user-parameter <a href="chapter_4.3.html#data_output_user">data_output_user</a>).<br><br>The |
---|
435 | time interval of the output times is determined via <a href="#dt_data_output">dt_data_output</a>. |
---|
436 | This is valid for all types of output quantities by default. Individual |
---|
437 | time intervals for instantaneous (!) 3d and section data can |
---|
438 | be |
---|
439 | declared using <a href="#dt_do3d">dt_do3d</a>, <a href="#dt_do2d_xy">dt_do2d_xy</a>, <a href="#dt_do2d_xz">dt_do2d_xz</a>, and <a href="#dt_do2d_yz">dt_do2d_yz</a>.<br><br>Also, |
---|
440 | an individual time interval for output of temporally averaged data can |
---|
441 | be assigned using parameter <a href="#dt_data_output_av">dt_data_output_av</a>. |
---|
442 | This applies to both 3d volume and cross section data. The length of |
---|
443 | the averaging interval is controlled via parameter <a href="#averaging_interval">averaging_interval</a>.<br><br>The |
---|
444 | parameter <a href="#skip_time_data_output">skip_time_data_output</a> |
---|
445 | can be used to shift data output activities for a given time interval. |
---|
446 | Individual intervals can be set using <a href="#skip_time_do3d">skip_time_do3d</a>, |
---|
447 | <a href="#skip_time_do2d_xy">skip_time_do2d_xy</a>, <a href="#skip_time_do2d_xz">skip_time_do2d_xz</a>, <a href="#skip_time_do2d_yz">skip_time_do2d_yz</a>, and <a href="#skip_time_data_output_av">skip_time_data_output_av</a>.<br><p>With |
---|
448 | the parameter <a href="chapter_4.2.html#nz_do3d">nz_do3d</a> |
---|
449 | the output can be limited in the vertical direction up to a certain |
---|
450 | grid point.<br> </p> Cross sections extend through the |
---|
451 | total model |
---|
452 | domain. In the two horizontal directions all grid points with 0 |
---|
453 | <= i |
---|
454 | <= nx+1 and 0 <= j |
---|
455 | <= ny+1 are output so that in case of cyclic boundary conditions |
---|
456 | the |
---|
457 | complete total domain is represented. The location(s) of the cross |
---|
458 | sections can be defined with parameters <a href="#section_xy">section_xy</a>, |
---|
459 | <a href="#section_xz">section_xz</a>, and <a href="#section_yz">section_yz</a>. Assigning <span style="font-weight: bold;">section_..</span> = <span style="font-style: italic;">-1</span> |
---|
460 | causes the output data to be averaged along the direction |
---|
461 | normal to the respective section.<br><br><br><span style="font-weight: bold;">Output of user defined quantities:</span><br><br>Beside |
---|
462 | the standard quantities from the above list, the user can output any |
---|
463 | other quantities. These have to be defined and calculated within the |
---|
464 | user-defined code (see <a href="chapter_3.5.4.html">3.5.4</a>). |
---|
465 | They can be selected for output with the user-parameter <a href="chapter_4.3.html#data_output_user">data_output_user</a> |
---|
466 | for which the same rules apply as for <span style="font-weight: bold;">data_output</span>. |
---|
467 | Output of the user defined quantities (time interval, averaging, |
---|
468 | selection of cross sections, etc.) is controlled with the parameters |
---|
469 | listed above and data are written to the same file(s) as the standard |
---|
470 | quantities.<br><br><p style="font-weight: bold;">Output |
---|
471 | on parallel machines:</p><p> |
---|
472 | By default, with parallel runs, processors output only data |
---|
473 | of their respective subdomains into seperate local files (file names |
---|
474 | are |
---|
475 | constructed by appending the four digit processor ID, e.g. |
---|
476 | <filename>_0000, <filename>_0001, etc.). |
---|
477 | After PALM has |
---|
478 | finished, the contents of these individual |
---|
479 | files are sampled into one final file<span style="font-weight: bold;"></span> |
---|
480 | using the program <tt><font style="font-size: 11pt;" size="2">combine_plot_fields.x</font></tt> |
---|
481 | (to be started e.g. by a suitable OUTPUT command in the <span style="font-weight: bold;">mrun</span> |
---|
482 | configuration file).</p> <p>Alternatively, PALM is able to |
---|
483 | collect all grid points of a |
---|
484 | cross section on PE0 before output is done. In this case only |
---|
485 | one |
---|
486 | output file (DATA_2D_XY_NETCDF, etc.) is created and <tt><font style="font-size: 11pt;" size="2">combine_plot_fields.x</font></tt> |
---|
487 | does not have to be called. In case of very large numbers of horizontal |
---|
488 | gridpoints, sufficient |
---|
489 | memory is required on PE0. This method can be used by |
---|
490 | assigning <a href="chapter_4.2.html#data_output_2d_on_each_pe">data_output_2d_on_each_pe</a> |
---|
491 | = <i>.F.</i>.</p><p>3d volume data output is |
---|
492 | always handled seperately by each processor so that <span style="font-family: monospace;">combine_plot_fields.x</span> |
---|
493 | has to be called anyway after PALM has been finished.</p><p><br><span style="font-weight: bold;">Old formats:</span></p> |
---|
494 | <p>Beside |
---|
495 | the NetCDF format, 2d cross section data and 3d volume data |
---|
496 | can |
---|
497 | also be output, for historical reasons, in a different (binary) format |
---|
498 | using parameter <a href="#data_output_format">data_output_format</a>.</p><p>By |
---|
499 | assigning <span style="font-weight: bold;">data_output_format |
---|
500 | </span>= <span style="font-style: italic;">'avs'</span>, |
---|
501 | the 3d volume data is output to the local file <a href="chapter_3.4.html#PLOT3D_DATA">PLOT3D_DATA</a>. |
---|
502 | Output is in FORTRAN binary format readable by |
---|
503 | the plot software <span style="font-weight: bold;">AVS</span>. |
---|
504 | The order of data on the file follows the order used in the assignment |
---|
505 | for <b>data_output</b> (e.g. <b>data_output</b> |
---|
506 | = <span style="font-style: italic;">'p'</span>, <span style="font-style: italic;">'v'</span>,... |
---|
507 | means that the file starts with the pressure data, followed by the |
---|
508 | v-component of the velocity, etc.). Both instantaneous and time |
---|
509 | averaged data are written on this file! Additional to this file, PALM |
---|
510 | creates |
---|
511 | a second binary file (local name <a href="chapter_3.4.html#PLOT3D_COOR">PLOT3D_COOR</a>) |
---|
512 | with coordinate information needed by <span style="font-weight: bold;">AVS</span>. |
---|
513 | As third and |
---|
514 | fourth file two ASCII files are created (AVS-FLD-format, local name <a href="chapter_3.4.html#PLOT3D_FLD">PLOT3D_FLD</a> |
---|
515 | and <a href="chapter_3.4.html#PLOT3D_FLD_COOR">PLOT3D_FLD_COOR</a>), |
---|
516 | which describe the contents of the data file and/or coordinate file |
---|
517 | and are used by AVS. However, AVS expects the content description in |
---|
518 | one file. This needs the local file PLOT3D_FLD_COOR to be appended to |
---|
519 | the file |
---|
520 | PLOT3D_FLD (by suitable OUTPUT command in the <span style="font-weight: bold;">mrun</span> |
---|
521 | configuration file: “<span style="font-family: monospace;">cat |
---|
522 | PLOT3D_FLD_COOR >> PLOT3D_FLD</span>”) |
---|
523 | after PALM has |
---|
524 | finished. To reduce the amount of data, output to this file |
---|
525 | can be done |
---|
526 | in |
---|
527 | compressed form (see <a href="chapter_4.2.html#do3d_compress">do3d_compress</a>). |
---|
528 | Further details about plotting 3d volume data with <span style="font-weight: bold;">AVS </span>can be found in |
---|
529 | <a href="chapter_4.5.5.html">chapter |
---|
530 | 4.5.5</a>.</p>By assigning <span style="font-weight: bold;">data_output_format </span>= |
---|
531 | <span style="font-style: italic;">'iso2d'</span>, |
---|
532 | the cross section data is output to the local files <a href="chapter_3.4.html#PLOT2D_XY">PLOT2D_XY</a>, <a href="chapter_3.4.html#PLOT2D_XZ">PLOT2D_XZ</a>, and <a href="chapter_3.4.html#PLOT2D_YZ">PLOT2D_YZ</a>. |
---|
533 | Output is in FORTRAN binary format readable by |
---|
534 | the plot software <span style="font-weight: bold;">iso2d</span>. |
---|
535 | The order of data on the files follows the order used in the assignment |
---|
536 | for <b>data_output</b> (e.g. <b>data_output</b> |
---|
537 | = <span style="font-style: italic;">'p_xy'</span>, <span style="font-style: italic;">'v_xy_av'</span>,... |
---|
538 | means that the file containing the horizontal cross section data starts |
---|
539 | with the instantaneous pressure data, followed by the |
---|
540 | temporally averaged v-component of the velocity, etc.). Both |
---|
541 | instantaneous and time averaged data are written on this |
---|
542 | file!Additional to these binary files, PALM |
---|
543 | creates NAMELIST parameter files |
---|
544 | (local names <a href="chapter_3.4.html#PLOT2D_XY_GLOBAL">PLOT2D_XY_GLOBAL</a>, |
---|
545 | <a href="chapter_3.4.html#PLOT2D_XY_LOCAL">PLOT2D_XY_LOCAL</a>, |
---|
546 | <a href="chapter_3.4.html#PLOT2D_XZ_GLOBAL">PLOT2D_XZ_GLOBAL</a>, |
---|
547 | <a href="chapter_3.4.html#PLOT2D_XZ_LOCAL">PLOT2D_XZ_LOCAL</a>, |
---|
548 | <a href="chapter_3.4.html#PLOT2D_YZ_GLOBAL">PLOT2D_YZ_GLOBAL</a>, |
---|
549 | <a href="chapter_3.4.html#PLOT2D_YZ_LOCAL">PLOT2D_YZ_LOCAL</a>) |
---|
550 | which can be used as parameter input files for the plot software <a href="http://www.muk.uni-hannover.de/institut/software/iso2d_beschreibung.html">iso2d</a>. |
---|
551 | That needs local files with suffix _LOCAL to be appended to the |
---|
552 | respective files with suffix _GLOBAL (by |
---|
553 | suitable OUTPUT commands in the <span style="font-weight: bold;">mrun</span> |
---|
554 | configuration file, e.g.: “<span style="font-family: monospace;">cat |
---|
555 | PLOT2D_XY_LOCAL >> PLOT2D_XY_GLOBAL</span>”) |
---|
556 | after PALM has |
---|
557 | finished. Cross sections can be directly plotted with <span style="font-weight: bold;">iso2d</span> using the |
---|
558 | respective data and |
---|
559 | parameter file. The plot layout is steered via the parameter input |
---|
560 | file. |
---|
561 | The values of these <span style="font-weight: bold;">iso2d</span> |
---|
562 | parameters are determined by a set of mostly internal PALM parameters |
---|
563 | (exception: <a href="chapter_4.2.html#z_max_do2d">z_max_do2d</a>). |
---|
564 | All parameter values can be changed by editing the parameter input |
---|
565 | file. Further details about plotting 2d cross sections with <span style="font-weight: bold;">iso2d </span>can be found |
---|
566 | in <a href="chapter_4.5.4.html">chapter |
---|
567 | 4.5.4</a>.<br><br><span style="font-weight: bold;">Important:</span><br>There |
---|
568 | is no guarantee that iso2d- and avs-output will be available in future |
---|
569 | PALM versions (later than 3.0). </td> </tr> <tr> <td style="vertical-align: top;"><a name="data_output_format"></a><span style="font-weight: bold;">data_output_format</span><br> |
---|
570 | </td> <td style="vertical-align: top;">C * 10 (10) </td> |
---|
571 | <td style="vertical-align: top;"><span style="font-style: italic;">'netcdf'</span> </td> |
---|
572 | <td style="vertical-align: top;">Format of output data.<br><br>By |
---|
573 | default, all data (profiles, time |
---|
574 | series, spectra, particle data, cross sections, volume data) are output |
---|
575 | in NetCDF format (see chapter <a href="chapter_4.5.1.html">4.5.1</a>). |
---|
576 | Exception: restart data (local files <a href="chapter_3.4.html#BININ">BININ</a>, <a href="chapter_3.4.html#BINOUT">BINOUT</a>, <a href="chapter_3.4.html#PARTICLE_RESTART_DATA_IN">PARTICLE_RESTART_DATA_IN</a>, |
---|
577 | <a href="chapter_3.4.html#PARTICLE_RESTART_DATA_OUT">PARTICLE_RESTART_DATA_OUT</a>) |
---|
578 | are always output in FORTRAN binary format.<br><br>The |
---|
579 | numerical precision of the NetCDF output is determined with parameter <a href="#chapter_4.1.html#netcdf_precision">netcdf_precision</a>.<br><br>The |
---|
580 | maximum file size for NetCDF files is 2 GByte by default. Use the |
---|
581 | parameter <a href="#netcdf_64bit">netcdf_64bit</a> |
---|
582 | if larger files have to be created.<br><br>For historical |
---|
583 | reasons, other data formats are still available. Beside 'netcdf', <span style="font-weight: bold;">data_output_format</span> |
---|
584 | may be assigned the following values:<br><br><table style="text-align: left; width: 594px; height: 104px;" border="1" cellpadding="2" cellspacing="2"><tbody><tr><td style="vertical-align: top;"><span style="font-style: italic;">'profil'</span></td><td>output |
---|
585 | of profiles, time series and spectra in ASCII format to be |
---|
586 | read by the graphic software <span style="font-weight: bold;">profil |
---|
587 | </span>(see chapters <a href="chapter_4.5.2.html">4.5.2</a>, |
---|
588 | <a href="#chapter_4.5.3.html">4.5.3</a>)</td></tr><tr><td style="vertical-align: top;"><span style="font-style: italic;">'iso2d'</span></td><td>output |
---|
589 | of 2d cross-sections in FORTRAN binary format to be read by the graphic |
---|
590 | software <span style="font-weight: bold;">iso2d</span> |
---|
591 | (see chapter <a href="chapter_4.5.4.html">4.5.4</a>)</td></tr><tr><td style="vertical-align: top;"><span style="font-style: italic;">'avs'</span></td><td>output |
---|
592 | of 3d volume data in FORTRAN binary format to be read by the graphic |
---|
593 | software <span style="font-weight: bold;">AVS</span> |
---|
594 | (see chapter <a href="chapter_4.5.5.html">4.5.5</a>)</td></tr></tbody></table><br>Multiple |
---|
595 | values can be assigned to <span style="font-weight: bold;">data_output_format</span>, |
---|
596 | i.e. if the user wants to have both the "old" data format suitable for <span style="font-weight: bold;">iso2d</span> as well as |
---|
597 | cross section data in NetCDF format, then <span style="font-weight: bold;">data_output_format</span> = |
---|
598 | <span style="font-style: italic;">'iso2d'</span>, <span style="font-style: italic;">'netcdf'</span> has to be |
---|
599 | assigned.<br><br><span style="font-weight: bold;">Warning:</span> |
---|
600 | There is no guarantee that the "old" formats will be available in |
---|
601 | future PALM versions (beyond 3.0)!<br> </td> </tr> <tr> |
---|
602 | <td style="vertical-align: top;"> <p><a name="data_output_pr"></a><b>data_output_pr</b></p> |
---|
603 | </td> <td style="vertical-align: top;">C * |
---|
604 | 10 <br> |
---|
605 | (100)</td> <td style="vertical-align: top;"><i>100 |
---|
606 | * ' '</i></td> <td style="vertical-align: top;"> |
---|
607 | <p>Quantities for which vertical profiles (horizontally averaged) |
---|
608 | are to be output. </p> <p>By default vertical |
---|
609 | profile data is output to the local file <a href="chapter_3.4.html#DATA_1D_PR_NETCDF">DATA_1D_PR_NETCDF</a>. |
---|
610 | The file's format is NetCDF. Further details about processing |
---|
611 | NetCDF data are given in chapter <a href="chapter_4.5.1.html">4.5.1</a>.</p><p>For |
---|
612 | historical reasons, data can also be output in ASCII-format on local |
---|
613 | file <a href="chapter_3.4.html#PLOT1D_DATA">PLOT1D_DATA</a> |
---|
614 | which is readable by the graphic software <span style="font-weight: bold;">profil</span>. See |
---|
615 | parameter <a href="#data_output_format">data_output_format</a> |
---|
616 | for defining the format in which data shall be output.<br> </p> |
---|
617 | <p>For horizontally averaged vertical |
---|
618 | profiles always <span style="font-weight: bold;">all</span> |
---|
619 | vertical |
---|
620 | grid points (0 <= k <= nz+1) are output to file. Vertical |
---|
621 | profile data refers to the total domain but profiles for subdomains can |
---|
622 | also be output (see <a href="chapter_4.1.html#statistic_regions">statistic_regions</a>). |
---|
623 | </p> <p>The temporal interval of the output times of |
---|
624 | profiles is |
---|
625 | assigned via the parameter <a href="chapter_4.2.html#dt_dopr">dt_dopr</a>. |
---|
626 | Within the file <a href="chapter_3.4.html#PLOT1D_DATA">PLOT1D_DATA</a>, |
---|
627 | the profiles are ordered with respect to their |
---|
628 | output times.</p><p>Profiles can also be temporally |
---|
629 | averaged (see <a href="chapter_4.2.html#averaging_interval_pr">averaging_interval_pr</a>). <br> |
---|
630 | </p> <p>The following list shows the values which can be |
---|
631 | assigned to <span style="font-weight: bold;">data_output_pr</span>. |
---|
632 | The profile data is either defined on |
---|
633 | u-v-levels (variables marked in <font color="#ff6600">red</font>) |
---|
634 | or |
---|
635 | on w-levels (<font color="#33ff33">green</font>). |
---|
636 | According to this, |
---|
637 | the |
---|
638 | z-coordinates of the individual profiles vary. Beyond that, with a |
---|
639 | Prandtl layer switched on (<a href="chapter_4.1.html#prandtl_layer">prandtl_layer</a>) |
---|
640 | the lowest output |
---|
641 | level is z = zu(1) instead of z = zw(0) for profiles <i>w'' |
---|
642 | u'',w''v"</i>, <i>wu</i> and <i>wv</i> |
---|
643 | . <br> </p> <table style="text-align: left; width: 100%;" cellpadding="2" cellspacing="2"> <tbody> <tr> <td style="vertical-align: top;"><font color="#ff6600"><i>u</i></font></td> |
---|
644 | <td style="vertical-align: top;">u-component of the |
---|
645 | velocity (in m/s).</td> </tr> <tr> <td style="vertical-align: top;"><font color="#ff6600"><i>v</i></font></td> |
---|
646 | <td style="vertical-align: top;">v-component of the |
---|
647 | velocity (in m/s).</td> </tr> <tr> <td style="vertical-align: top;"><font color="#33ff33"><i>w</i></font></td> |
---|
648 | <td style="vertical-align: top;">w-component of the |
---|
649 | velocity (in m/s).</td> </tr> <tr> <td style="vertical-align: top;"><font color="#ff6600"><i>pt</i></font></td> |
---|
650 | <td style="vertical-align: top;">Potential temperature (in |
---|
651 | K).</td> </tr> <tr> <td style="vertical-align: top;"><font color="#ff6600"><i>vpt</i></font></td> |
---|
652 | <td style="vertical-align: top;">Virtual potential |
---|
653 | temperature (in K).</td> </tr> <tr> <td style="vertical-align: top;"><font color="#ff6600"><i>lpt</i></font></td> |
---|
654 | <td style="vertical-align: top;">Potential liquid water |
---|
655 | temperature (in K).</td> </tr> <tr> <td style="vertical-align: top;"><font color="#ff6600"><i>q</i></font></td> |
---|
656 | <td style="vertical-align: top;">Total water content |
---|
657 | (in kg/kg).</td> </tr> <tr> <td style="vertical-align: top;"><font color="#ff6600"><i>qv</i></font></td> |
---|
658 | <td style="vertical-align: top;">Specific humidity (in |
---|
659 | kg/kg).</td> </tr> <tr> <td style="vertical-align: top;"><font color="#ff6600"><i>ql</i></font></td> |
---|
660 | <td style="vertical-align: top;">Liquid water content |
---|
661 | (in kg/kg).</td> </tr> <tr> <td style="vertical-align: middle;"><font color="#ff6600">s</font></td> |
---|
662 | <td style="vertical-align: top;">Scalar concentration (in |
---|
663 | kg/m<sup>3</sup>).</td> </tr> <tr> <td style="vertical-align: middle;"><font color="#ff6600"><i>e</i></font></td> |
---|
664 | <td style="vertical-align: top;">Turbulent kinetic energy |
---|
665 | (TKE, subgrid-scale) (in m<sup>2</sup>/s<sup>2</sup>).</td> |
---|
666 | </tr> <tr> <td style="vertical-align: middle;"><font color="#ff6600"><i>e*</i></font></td> |
---|
667 | <td style="vertical-align: top;">Perturbation energy |
---|
668 | (resolved) (in m<sup>2</sup>/s<sup>2</sup>).</td> |
---|
669 | </tr> <tr> <td style="vertical-align: middle;"><font color="#ff6600"><i>km</i></font></td> |
---|
670 | <td style="vertical-align: top;">Eddy diffusivity for |
---|
671 | momentum (in m<sup>2</sup>/s).</td> </tr> <tr> |
---|
672 | <td style="vertical-align: middle;"><font color="#ff6600"><i>kh</i></font></td> |
---|
673 | <td style="vertical-align: top;">Eddy diffusivity for heat |
---|
674 | (in m<sup>2</sup>/s).</td> </tr> <tr> <td style="vertical-align: top;"><font color="#ff6600"><i>l</i></font></td> |
---|
675 | <td style="vertical-align: top;">Mixing length (in m).</td> |
---|
676 | </tr> <tr> <td style="vertical-align: middle;"><font color="#33ff33"><i>w"u"</i></font></td> |
---|
677 | <td style="vertical-align: top;">u-component of the |
---|
678 | subgrid-scale vertical momentum flux (in m<sup>2</sup>/s<sup>2</sup>).</td> |
---|
679 | </tr> <tr> <td style="vertical-align: middle;"><font color="#33ff33"><i>w*u*</i></font></td> |
---|
680 | <td style="vertical-align: top;">u-component of the |
---|
681 | resolved vertical momentum flux (in m<sup>2</sup>/s<sup>2</sup>).</td> |
---|
682 | </tr> <tr> <td style="vertical-align: middle;"><font color="#33ff33"><i>wu</i></font></td> |
---|
683 | <td style="vertical-align: top;">u-component of the total |
---|
684 | vertical momentum flux (<i>w"u"</i> + <i>w*u*</i>) |
---|
685 | (in m<sup>2</sup>/s<sup>2</sup>).</td> </tr> |
---|
686 | <tr> <td style="vertical-align: middle;"><font color="#33ff33"><i>w"v"</i></font></td> |
---|
687 | <td style="vertical-align: top;">v-component of the |
---|
688 | subgrid-scale vertical momentum flux (in m<sup>2</sup>/s<sup>2</sup>).</td> |
---|
689 | </tr> <tr> <td style="vertical-align: middle;"><font color="#33ff33"><i>w*v*</i></font></td> |
---|
690 | <td style="vertical-align: top;">v-component of the |
---|
691 | resolved vertical momentum flux (in m<sup>2</sup>/s<sup>2</sup>).</td> |
---|
692 | </tr> <tr> <td style="vertical-align: middle;"><font color="#33ff33"><i>wv</i></font></td> |
---|
693 | <td style="vertical-align: top;">v-component of the total |
---|
694 | vertical momentum flux (<i>w"v"</i> + <i>w*v*</i>) |
---|
695 | (in m<sup>2</sup>/s<sup>2</sup>).</td> </tr> |
---|
696 | <tr> <td style="vertical-align: top;"><font color="#33ff33"><i>w"pt"</i></font></td> |
---|
697 | <td style="vertical-align: top;">Subgrid-scale vertical |
---|
698 | sensible heat flux (in K m/s).</td> </tr> <tr> <td style="vertical-align: top;"><font color="#33ff33"><i>w*pt*</i></font></td> |
---|
699 | <td style="vertical-align: top;">Resolved vertical |
---|
700 | sensible |
---|
701 | heat flux (in K m/s).</td> </tr> <tr> <td style="vertical-align: top;"><font color="#33ff33"><i>wpt</i></font></td> |
---|
702 | <td style="vertical-align: top;">Total vertical sensible |
---|
703 | heat flux (<i>w"pt"</i> + <i>w*pt*</i>) |
---|
704 | (in K |
---|
705 | m/s).</td> </tr> <tr> <td style="vertical-align: top;"><font color="#33ff33"><i>w*pt*BC</i></font></td> |
---|
706 | <td style="vertical-align: top;">Subgrid-scale vertical |
---|
707 | sensible heat flux using the |
---|
708 | Bott-Chlond scheme (in K m/s).</td> </tr> <tr> <td style="vertical-align: top;"><font color="#33ff33"><i>wptBC</i></font></td> |
---|
709 | <td style="vertical-align: top;">Total vertical sensible |
---|
710 | heat flux using the Bott-Chlond scheme |
---|
711 | (<i>w"pt"</i> |
---|
712 | + <i>w*pt*BC</i>) (in K m/s).</td> </tr> <tr> |
---|
713 | <td style="vertical-align: top;"><font color="#33ff33"><i>w"vpt"</i></font></td> |
---|
714 | <td style="vertical-align: top;">Subgrid-scale vertical |
---|
715 | buoyancy flux (in K m/s).</td> </tr> <tr> <td style="vertical-align: top;"><font color="#33ff33"><i>w*pt*</i></font></td> |
---|
716 | <td style="vertical-align: top;">Resolved vertical |
---|
717 | buoyancy |
---|
718 | flux (in K m/s).</td> </tr> <tr> <td style="vertical-align: top;"><font color="#33ff33"><i>wvpt</i></font></td> |
---|
719 | <td style="vertical-align: top;">Total vertical buoyancy |
---|
720 | flux (w"vpt" + w*vpt*) (in K m/s).</td> </tr> <tr> <td style="vertical-align: top;"><font color="#33ff33"><i>w"q"</i></font></td> |
---|
721 | <td style="vertical-align: top;">Subgrid-scale vertical |
---|
722 | water flux (in kg/kg m/s).</td> </tr> <tr> <td style="vertical-align: top;"><font color="#33ff33"><i>w*q*</i></font></td> |
---|
723 | <td style="vertical-align: top;">Resolved vertical water |
---|
724 | flux (in kg/kg m/s).</td> </tr> <tr> <td style="vertical-align: top;"><font color="#33ff33"><i>wq</i></font></td> |
---|
725 | <td style="vertical-align: top;">Total vertical water flux |
---|
726 | (w"q" + w*q*) (in kg/kg m/s).</td> </tr> <tr> <td style="vertical-align: top;"><font color="#33ff33"><i>w"qv"</i></font></td> |
---|
727 | <td style="vertical-align: top;">Subgrid-scale vertical |
---|
728 | latent heat flux (in kg/kg m/s).</td> </tr> <tr> <td style="vertical-align: top;"><font color="#33ff33"><i>w*qv*</i></font></td> |
---|
729 | <td style="vertical-align: top;">Resolved vertical latent |
---|
730 | heat flux (in kg/kg m/s).</td> </tr> <tr> <td style="vertical-align: top;"><font color="#33ff33"><i>wqv</i></font></td> |
---|
731 | <td style="vertical-align: top;">Total vertical latent |
---|
732 | heat |
---|
733 | flux (w"qv" + w*qv*) (in kg/kg m/s).</td> </tr> <tr> |
---|
734 | <td style="vertical-align: middle;"><font color="#33ff33"><i>w"s"</i></font></td> |
---|
735 | <td style="vertical-align: top;">Subgrid-scale vertical |
---|
736 | scalar concentration flux (in kg/m<sup>3 </sup>m/s).</td> |
---|
737 | </tr> <tr> <td style="vertical-align: middle;"><font color="#33ff33"><i>w*s*</i></font></td> |
---|
738 | <td style="vertical-align: top;">Resolved vertical scalar |
---|
739 | concentration flux (in kg/m<sup>3</sup>)</td> </tr> |
---|
740 | <tr> <td style="vertical-align: middle;"><font color="#33ff33"><i>ws</i></font></td> |
---|
741 | <td style="vertical-align: top;">Total vertical scalar |
---|
742 | concentration flux (w"s" + w*s*) (in kg/m<sup>3 </sup>m/s).</td> |
---|
743 | </tr> <tr> <td style="vertical-align: top;"><font color="#33ff33"><i>w*e*</i></font></td> |
---|
744 | <td style="vertical-align: top;">Vertical flux of |
---|
745 | perturbation energy (resolved)</td> </tr> <tr> <td style="vertical-align: top;"><font color="#ff6600"><i>u*2</i></font></td> |
---|
746 | <td style="vertical-align: top;">Variance of the |
---|
747 | u-velocity |
---|
748 | component (resolved)</td> </tr> <tr> <td style="vertical-align: top;"><font color="#ff6600"><i>v*2</i></font></td> |
---|
749 | <td style="vertical-align: top;">Variance of the |
---|
750 | v-velocity |
---|
751 | component (resolved)</td> </tr> <tr> <td style="vertical-align: top;"><font color="#33ff33"><i>w*2</i></font></td> |
---|
752 | <td style="vertical-align: top;">Variance of the potential |
---|
753 | temperature (resolved)</td> </tr> <tr> <td style="vertical-align: top;"><font color="#ff6600"><i>pt*2</i></font></td> |
---|
754 | <td style="vertical-align: top;">Variance of the potential |
---|
755 | temperature (resolved)</td> </tr> <tr> <td style="vertical-align: top;"><font color="#33ff33"><i>w*3</i></font></td> |
---|
756 | <td style="vertical-align: top;">Third moment of the |
---|
757 | w-velocity component (resolved)</td> </tr> <tr> <td style="vertical-align: middle;"><font color="#33ff33"><i>Sw</i></font></td> |
---|
758 | <td style="vertical-align: top;">Skewness of the |
---|
759 | w-velocity |
---|
760 | component (resolved, S<sub>w</sub> |
---|
761 | = W<sup>3</sup>/(w<sup>2</sup>)<sup>1.5</sup>)</td> |
---|
762 | </tr> <tr> <td style="vertical-align: top;"><font color="#33ff33"><i>w*2pt*</i></font></td> |
---|
763 | <td style="vertical-align: top;">Third moment (resolved)</td> |
---|
764 | </tr> <tr> <td style="vertical-align: top;"><font color="#33ff33"><i>w*pt*2</i></font></td> |
---|
765 | <td style="vertical-align: top;">Third moment (resolved)</td> |
---|
766 | </tr> <tr> <td style="vertical-align: top;"><font color="#ff6666"><i>w*u*u*/dz</i></font></td> |
---|
767 | <td style="vertical-align: top;">Energy production by |
---|
768 | shear |
---|
769 | (resolved)</td> </tr> <tr> <td style="vertical-align: top;"><font color="#ff6666"><i>w*p*/dz</i></font></td> |
---|
770 | <td style="vertical-align: top;">Energy production by |
---|
771 | turbulent transport of pressure |
---|
772 | fluctuations (resolved)</td> </tr> <tr> <td style="vertical-align: top;"><font color="#ff6666"><i>w"e/dz</i></font></td> |
---|
773 | <td style="vertical-align: top;">Energy production by |
---|
774 | transport of resolved-scale TKE</td> </tr> </tbody> |
---|
775 | </table> <br>Beyond that, initial profiles (t=0) of some |
---|
776 | variables can be also be |
---|
777 | output (this output is only done once |
---|
778 | with the first plot output and not repeated with the profile output at |
---|
779 | later |
---|
780 | times). The names of these profiles result from the ones specified |
---|
781 | above leaded by a hash "#". Allowed values are:<br> <ul> |
---|
782 | <p><i>#u</i>, <i>#v</i>, <i>#pt</i>, |
---|
783 | <i>#km</i>, <i>#kh</i>, <i>#l</i></p> |
---|
784 | </ul> <p>These initial profiles have been either set by |
---|
785 | the user or |
---|
786 | have been calculated by a 1d-model prerun.<br> </p>In case |
---|
787 | of ASCII data output to local file PLOT1D_DATA, |
---|
788 | PALM additionally creates a NAMELIST parameter file (local name <a href="chapter_3.4.html#PLOT1D_PAR">PLOT1D_PAR</a>) |
---|
789 | which can be used as parameter input file for the plot software <a href="http://www.muk.uni-hannover.de/institut/software/profil_intro.html">profil</a>. |
---|
790 | Profiles can be directly plotted with <span style="font-weight: bold;">profil</span> |
---|
791 | using these two files. The |
---|
792 | plot layout is |
---|
793 | steered via the parameter input file. The values of these <span style="font-weight: bold;">profil</span>-parameters |
---|
794 | are determined by |
---|
795 | a set of PALM parameters (<a href="chapter_4.2.html#profile_columns">profile_columns</a>, |
---|
796 | <a href="chapter_4.2.html#profile_rows">profile_rows</a>, |
---|
797 | <a href="chapter_4.2.html#z_max_do1d">z_max_do1d</a>, |
---|
798 | <a href="chapter_4.2.html#cross_profiles">cross_profiles</a>, |
---|
799 | etc.) All parameter values can be changed by editing the parameter |
---|
800 | input |
---|
801 | file. <br><br>Further details about plotting vertical |
---|
802 | profiles with <span style="font-weight: bold;">profil </span>can |
---|
803 | be found in <a href="chapter_4.5.2.html">chapter |
---|
804 | 4.5.2</a></td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="data_output_2d_on_each_pe"></a><b>data_output_2d_on</b> |
---|
805 | <br> <b>_each_pe</b></p> </td> <td style="vertical-align: top;">L<br> </td> <td style="vertical-align: top;"><span style="font-style: italic;">.T.</span><br> </td> |
---|
806 | <td style="vertical-align: top;">Output 2d cross section |
---|
807 | data by one or |
---|
808 | all processors. <p>In runs with several processors, by |
---|
809 | default, each processor |
---|
810 | outputs cross section data of its subdomain into an individual |
---|
811 | file. After PALM |
---|
812 | has finished, the contents of these files have to be sampled into one |
---|
813 | file<span style="font-weight: bold;"></span> using |
---|
814 | the program <tt>combine_plot_fields.x</tt>. </p> |
---|
815 | <p>Alternatively, by assigning <b>data_output_2d_on_each_pe</b> |
---|
816 | = <i>.F.,</i> |
---|
817 | the respective data is gathered on PE0 and output is done directly |
---|
818 | into one file, so <tt>combine_plot_fields.x</tt> does not |
---|
819 | have to be |
---|
820 | called. However, in case of very large numbers of horizontal |
---|
821 | gridpoints, sufficient |
---|
822 | memory is required on PE0. </p> </td> </tr> |
---|
823 | <tr> <td style="vertical-align: top;"> <p><a name="disturbance_amplitude"></a><b>disturbance<br> |
---|
824 | _amplitude</b></p> </td> <td style="vertical-align: top;">R</td> <td style="vertical-align: top;"><i>0.25</i></td> |
---|
825 | <td style="vertical-align: top;"> <p>Maximum |
---|
826 | perturbation amplitude of the random perturbations |
---|
827 | imposed to the horizontal velocity field (in m/s). </p> |
---|
828 | <p>The parameter <a href="#create_disturbances">create_disturbances</a> |
---|
829 | describes how to impose random perturbations to the horizontal velocity |
---|
830 | field. Since the perturbation procedure includes two filter operations, |
---|
831 | the amplitude assigned by <b>disturbance_amplitude</b> is |
---|
832 | only an |
---|
833 | approximate value of the real magnitude of the perturbation.</p> </td> |
---|
834 | </tr> <tr> <td style="vertical-align: top;"><p><a name="disturbance_energy_limit"></a><b>disturbance_energy</b> |
---|
835 | <br> <b>_limit</b></p> </td> <td style="vertical-align: top;">R</td> <td style="vertical-align: top;"><i>0.01</i></td> |
---|
836 | <td style="vertical-align: top;"> <p lang="en-GB">Upper |
---|
837 | limit value of the perturbation energy of |
---|
838 | the velocity field used as a criterion for imposing random |
---|
839 | perturbations (in m<sup>2</sup>/s<sup>2</sup>). |
---|
840 | </p> <p><span lang="en-GB"><font face="Thorndale, serif">The parameter </font></span><a href="#create_disturbances"><span lang="en-GB"><font face="Thorndale, serif">create_disturbances</font></span></a><font face="Thorndale, serif"><span lang="en-GB"> |
---|
841 | describes how to impose |
---|
842 | random perturbations to the horizontal velocity field. The perturbation |
---|
843 | energy is defined as the volume average (over the total model domain) |
---|
844 | of the squares of the deviations of the velocity components from the |
---|
845 | mean flow (horizontal average). If the perturbation energy exceeds the |
---|
846 | assigned value, random perturbations to the fields of horizontal |
---|
847 | velocities are imposed no more. The value of this parameter usually |
---|
848 | must be determined by trial and error (it depends e.g. on the total |
---|
849 | number of grid points).</span> </font> </p> </td> |
---|
850 | </tr> <tr> <td style="vertical-align: top;"><p><a name="disturbance_level_b"></a><b>disturbance_level_b</b></p> |
---|
851 | </td> <td style="vertical-align: top;">R</td> |
---|
852 | <td style="vertical-align: top;"><i>zu(3)</i></td> |
---|
853 | <td style="vertical-align: top;"> <p lang="en-GB"><font face="Thorndale, serif"><font size="3">Lower |
---|
854 | limit of the vertical range for which random perturbations are to be |
---|
855 | imposed on the horizontal wind field (</font></font>in <font face="Thorndale, serif"><font size="3">m). |
---|
856 | </font></font> </p> <p><span lang="en-GB"><font face="Thorndale, serif">This |
---|
857 | parameter must hold the condition zu<i>(3)</i> <= <b>disturbance_level_b</b> |
---|
858 | <= <i>zu(</i></font></span><i><a href="chapter_4.1.html#nz"><span lang="en-GB"><font face="Thorndale, serif">nz-1</font></span></a><span lang="en-GB"><font face="Thorndale, serif">)</font></span></i><span lang="en-GB"><font face="Thorndale, serif">. |
---|
859 | Additionally, <b>disturbance_level_b</b> |
---|
860 | <= </font></span><a href="#disturbance_level_t"><span lang="en-GB"><font face="Thorndale, serif">disturbance_level_t</font></span></a> |
---|
861 | <span lang="en-GB"><font face="Thorndale, serif">must |
---|
862 | also hold. <br> </font></span></p> <p><span lang="en-GB"><font face="Thorndale, serif">The |
---|
863 | parameter </font></span><a href="#create_disturbances"><span lang="en-GB"><font face="Thorndale, serif">create_disturbances</font></span></a><font face="Thorndale, serif"><span lang="en-GB"> |
---|
864 | describes how to impose |
---|
865 | random perturbations to the horizontal velocity field</span></font><font face="Thorndale, serif"><span lang="en-GB">.</span> |
---|
866 | </font> </p> </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="disturbance_level_t"></a><b>disturbance_level_t</b></p> |
---|
867 | </td> <td style="vertical-align: top;">R</td> |
---|
868 | <td style="vertical-align: top;"><i>zu(nz/3)</i></td> |
---|
869 | <td style="vertical-align: top;"> <p lang="en-GB"><font face="Thorndale, serif"><font size="3">Upper |
---|
870 | limit of the vertical range for which random perturbations are to be |
---|
871 | imposed on the horizontal wind field (</font></font>in <font face="Thorndale, serif"><font size="3">m). |
---|
872 | </font></font> </p> <p><span lang="en-GB"><font face="Thorndale, serif">This |
---|
873 | parameter must hold the condition <b>disturbance_level_t</b> |
---|
874 | <= zu<i>(</i></font></span><i><a href="chapter_4.1.html#nz"><span lang="en-GB"><font face="Thorndale, serif">nz-1</font></span></a><span lang="en-GB"><font face="Thorndale, serif">)</font></span></i><span lang="en-GB"><font face="Thorndale, serif">. |
---|
875 | Additionally, </font></span><a href="#disturbance_level_b"><span lang="en-GB"><font face="Thorndale, serif">disturbance_level_b</font></span></a> |
---|
876 | <span lang="en-GB"><font face="Thorndale, serif"><= |
---|
877 | <b>disturbance_level_t</b> |
---|
878 | must also hold.<br> </font></span></p> <p><span lang="en-GB"><font face="Thorndale, serif">The |
---|
879 | parameter </font></span><a href="#create_disturbances"><span lang="en-GB"><font face="Thorndale, serif">create_disturbances</font></span></a><font face="Thorndale, serif"><span lang="en-GB"> |
---|
880 | describes how to impose |
---|
881 | random perturbations to the horizontal velocity field</span></font><font face="Thorndale, serif"><span lang="en-GB">.</span> |
---|
882 | </font> </p> </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="do2d_at_begin"></a><b>do2d_at_begin</b></p> |
---|
883 | </td> <td style="vertical-align: top;">L<br> </td> |
---|
884 | <td style="vertical-align: top;">.F.<br> </td> |
---|
885 | <td style="vertical-align: top;"> <p>Output of 2d |
---|
886 | cross section data at the beginning of a run. </p> <p>The |
---|
887 | temporal intervals of output times of 2d cross section data (see <a href="chapter_4.2.html#data_output">data_output</a>) |
---|
888 | are usually determined with parameters <a href="chapter_4.2.html#dt_do2d_xy">dt_do2d_xy</a>, <a href="chapter_4.2.html#dt_do2d_xz">dt_do2d_xz</a> |
---|
889 | and <a href="chapter_4.2.html#dt_do2d_yz">dt_do2d_yz</a>. |
---|
890 | By assigning <b>do2d_at_begin</b> = <i>.T.</i> |
---|
891 | an additional output |
---|
892 | will be made at the |
---|
893 | beginning of a run (thus at the time t = 0 or at the respective |
---|
894 | starting times of restart runs).</p> </td> </tr> <tr> |
---|
895 | <td style="vertical-align: top;"> <p><a name="do3d_at_begin"></a><b>do3d_at_begin</b></p> |
---|
896 | </td> <td style="vertical-align: top;">L<br> </td> |
---|
897 | <td style="vertical-align: top;">.F.<br> </td> |
---|
898 | <td style="vertical-align: top;">Output of 3d volume data |
---|
899 | at the beginning |
---|
900 | of a run.<br><br>The temporal intervals of output times of |
---|
901 | 3d volume data (see <a href="chapter_4.2.html#data_output">data_output</a>) |
---|
902 | is usually determined with parameter <a href="chapter_4.2.html#dt_do3d">dt_do3d</a>. |
---|
903 | By assigning <b>do3d_at_begin</b> = <i>.T.</i> |
---|
904 | an additional output |
---|
905 | will be made at the |
---|
906 | beginning of a run (thus at the time t = 0 or at the respective |
---|
907 | starting times of restart runs).</td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="do3d_compress"></a><b>do3d_compress</b></p> |
---|
908 | </td> <td style="vertical-align: top;">L<br> </td> |
---|
909 | <td style="vertical-align: top;">.F.<br> </td> |
---|
910 | <td style="vertical-align: top;"> <p>Output of data |
---|
911 | for 3d plots in compressed form. </p> <p>This |
---|
912 | parameter only applies for <a href="chapter_4.2.html#data_output_format">data_output_format</a> |
---|
913 | = <span style="font-style: italic;">'avs'</span>.</p><p>Output |
---|
914 | of 3d volume data may need huge amounts of disc storage |
---|
915 | (up to several Terabytes ore more). Data compression can serve to |
---|
916 | reduce this requirement. PALM is able to output 3d data in compressed |
---|
917 | form using 32-bit integers, if <span style="font-weight: bold;">do3d_compress</span> |
---|
918 | = <span style="font-style: italic;">.T.</span> is |
---|
919 | assigned. This |
---|
920 | yields a loss of accuracy, but the file size is clearly reduced. The |
---|
921 | parameter <a href="chapter_4.2.html#do3d_precision">do3d_precision</a> |
---|
922 | can be used to separately define the number of significant digits for |
---|
923 | each quantity.<br> </p> <p>So far compressed data |
---|
924 | output is only possible for Cray-T3E |
---|
925 | machines. Additional information for |
---|
926 | handling compressed data is given in <a href="chapter_4.5.6.html">chapter |
---|
927 | 4.5.6</a>.</p> </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="do3d_precision"></a><b>do3d_precision</b></p> |
---|
928 | </td> <td style="vertical-align: top;">C * |
---|
929 | 7 <br> |
---|
930 | (100)</td> <td style="vertical-align: top;">see |
---|
931 | right<br> </td> <td style="vertical-align: top;"> |
---|
932 | <p>Significant digits in case of compressed data |
---|
933 | output. </p> <p>This parameter only applies for |
---|
934 | <a href="chapter_4.2.html#data_output_format">data_output_format</a> |
---|
935 | = <span style="font-style: italic;">'avs'</span>.</p><p>In |
---|
936 | case that data compression is used for output of 3d data |
---|
937 | (see <a href="chapter_4.2.html#do3d_compress">do3d_compress</a>), |
---|
938 | this parameter determines the number of significant digits |
---|
939 | which are to be output.<br> </p> <p>Fewer digits |
---|
940 | clearly reduce the amount |
---|
941 | of data. Assignments have to be given separately for each individual |
---|
942 | quantity via a character string of the form <span style="font-style: italic;">'<quantity |
---|
943 | name><number of |
---|
944 | significant digits>'</span>, e.g. <span style="font-style: italic;">'pt2'</span>. |
---|
945 | Only those quantities listed in <a href="chapter_4.2.html#data_output">data_output</a> |
---|
946 | are admitted. Up to 9 significant digits are allowed (but large values |
---|
947 | are not very reasonable |
---|
948 | because they do not effect a significant compression).<br> </p> |
---|
949 | <p>The default assignment is <span style="font-weight: bold;">do3d_precision</span> |
---|
950 | = <span style="font-style: italic;">'u2'</span>, <span style="font-style: italic;">'v2'</span>, <span style="font-style: italic;">'w2'</span>, <span style="font-style: italic;">'p5'</span>, <span style="font-style: italic;">'pt2'</span>.</p> </td> |
---|
951 | </tr><tr> <td style="vertical-align: top;"> <p><a name="dt_laufparameter"></a><b>dt</b></p> |
---|
952 | </td> <td style="vertical-align: top;">R</td> |
---|
953 | <td style="vertical-align: top;"><i>variable</i></td> |
---|
954 | <td style="vertical-align: top;"> <p lang="en-GB"><font face="Thorndale, serif"><font size="3">Time |
---|
955 | step to be used by the 3d-model (</font></font>in <font face="Thorndale, serif"><font size="3">s). |
---|
956 | </font></font> </p> <p><span lang="en-GB"><font face="Thorndale, serif">This parameter</font></span> |
---|
957 | <font face="Thorndale, serif"><span lang="en-GB">is |
---|
958 | described in |
---|
959 | detail with the initialization parameters (see</span></font><span lang="en-GB"><font face="Thorndale, serif"> </font></span><a href="chapter_4.1.html#dt"><span lang="en-GB"><font face="Thorndale, serif">dt</font></span></a><font face="Thorndale, serif"><span lang="en-GB">). |
---|
960 | Additionally, it may be |
---|
961 | used as a run parameter and then applies to all restart runs (until it |
---|
962 | is changed again). A switch from a constant time step to a variable |
---|
963 | time step can be achieved with <b>dt</b> = <i>-1.0</i>.</span> |
---|
964 | </font> </p> </td> </tr> <tr> <td style="vertical-align: top;"><a name="dt_averaging_input"></a><span style="font-weight: bold;">dt_averaging_input</span><br> |
---|
965 | </td> <td style="vertical-align: top;">R<br> </td> |
---|
966 | <td style="vertical-align: top;"><span style="font-style: italic;">0.0</span><br> </td> |
---|
967 | <td style="vertical-align: top;">Temporal interval |
---|
968 | of data which are subject to temporal averaging (in s).<br><br>By |
---|
969 | default, data from each timestep within the interval defined by <a href="chapter_4.2.html#averaging_interval">averaging_interval</a> |
---|
970 | are used for calculating the temporal average. By choosing <span style="font-weight: bold;">dt_averaging_input</span> |
---|
971 | > <span lang="en-GB"><font face="Thorndale, serif"> </font></span><a href="chapter_4.1.html#dt"><span lang="en-GB"><font face="Thorndale, serif">dt</font></span></a><font face="Thorndale, serif"><span lang="en-GB"></span></font><span lang="en-GB"></span><span style="font-style: italic;"></span>, |
---|
972 | the number of time levels entering the average can be minimized. This |
---|
973 | reduces the CPU-time of a run but may worsen the quality of the |
---|
974 | average's statistics.<br><br><font face="Thorndale, serif"><span lang="en-GB">With |
---|
975 | variable time step (see <span style="font-weight: bold;">dt</span>), |
---|
976 | the number of time levels entering the average can vary from one |
---|
977 | averaging interval to the next (for a more detailed explanation see </span></font><font><a href="#averaging_interval"><span lang="en-GB"><font face="Thorndale, serif">averaging_interval</font></span></a>)</font><font face="Thorndale, serif"><span lang="en-GB">. It |
---|
978 | is approximately given by the quotient of <span style="font-weight: bold;">averaging_interval</span> / |
---|
979 | MAX(<span style="font-weight: bold;"> dt_averaging_input</span>, |
---|
980 | <span style="font-weight: bold;">dt</span>) (which |
---|
981 | gives a more or less exact value if a fixed timestep is used and if |
---|
982 | this is an integral divisor of <span style="font-weight: bold;">dt_averaging_input</span>).</span></font> |
---|
983 | <br><br><span style="font-weight: bold;">Example:</span><br>With |
---|
984 | an averaging interval of 100.0 s and <span style="font-weight: bold;">dt_averaging_input</span> = |
---|
985 | <span style="font-style: italic;">10.0</span>, |
---|
986 | the time levels entering the average have a (minimum) distance of 10.0 |
---|
987 | s (their distance may of course be larger if the current timestep is |
---|
988 | larger than 10.0 s), so the average is calculated from the data of |
---|
989 | (maximum) 10 time levels.<br><br><font face="Thorndale, serif"><span lang="en-GB">It |
---|
990 | is allowed |
---|
991 | to change <b>dt_averaging_input</b> during a job chain. If |
---|
992 | the last averaging |
---|
993 | interval of the run previous to the change could not be completed (i.e. |
---|
994 | has to be finished in the current run), the individual profiles and/or |
---|
995 | spectra entering the averaging are not uniformly distributed over the |
---|
996 | averaging interval.<br><br></span></font>Parameter <a href="#dt_averaging_input_pr">dt_averaging_input_pr</a> can |
---|
997 | be used to define a different temporal interval for |
---|
998 | vertical profile data and spectra.<br> </td> </tr> |
---|
999 | <tr> <td style="vertical-align: top;"> <p><a name="dt_averaging_input_pr"></a><b>dt_averaging_input_pr</b></p> |
---|
1000 | </td> <td style="vertical-align: top;">R</td> |
---|
1001 | <td style="vertical-align: top;"><span style="font-style: italic;">value of <a href="#dt_averaging_input">dt_<br>averaging_<br>input</a></span></td> |
---|
1002 | <td style="vertical-align: top;"> <p lang="en-GB">Temporal |
---|
1003 | interval of data which are subject to temporal averaging of <font face="Thorndale, serif"><font size="3">vertical |
---|
1004 | profiles and/or spectra (</font></font>in <font face="Thorndale, serif"><font size="3">s). |
---|
1005 | </font></font> </p> <p>By default, data from |
---|
1006 | each timestep within the interval defined by<font face="Thorndale, serif"><span lang="en-GB"> </span></font><a href="#averaging_interval_pr"><span lang="en-GB"><font face="Thorndale, serif">averaging_interval_pr</font></span></a><span lang="en-GB"><font face="Thorndale, serif">, </font></span><span lang="en-GB"><font face="Thorndale, serif">and </font></span><a href="#averaging_interval_sp"><span lang="en-GB"><font face="Thorndale, serif">averaging_interval_sp</font></span></a><span lang="en-GB"><font face="Thorndale, serif"> </font></span>are |
---|
1007 | used for calculating the temporal average. By choosing <span style="font-weight: bold;">dt_averaging_input_pr</span> |
---|
1008 | > <span lang="en-GB"><font face="Thorndale, serif"> </font></span><a href="chapter_4.1.html#dt"><span lang="en-GB"><font face="Thorndale, serif">dt</font></span></a><font face="Thorndale, serif"><span lang="en-GB"></span></font><span lang="en-GB"></span><span style="font-style: italic;"></span>, |
---|
1009 | the number of time levels entering the average can be minimized. This |
---|
1010 | reduces the CPU-time of a run but may worsen the quality of the |
---|
1011 | average's statistics. <span lang="en-GB"><font face="Thorndale, serif"><span style="font-weight: bold;"></span><span style="font-weight: bold;"></span></font></span><a href="chapter_4.1.html#dt"><span lang="en-GB"></span></a><font face="Thorndale, serif"><span lang="en-GB"></span></font><span lang="en-GB"></span><br> </p><p>For |
---|
1012 | more explanations see parameter <a href="#dt_averaging_input">dt_averaging_input</a>.<a href="chapter_4.1.html#dt"><span lang="en-GB"></span></a><font face="Thorndale, serif"><span lang="en-GB"></span></font></p></td> |
---|
1013 | </tr> <tr> <td style="vertical-align: top;"><a name="dt_data_output"></a><span style="font-weight: bold;">dt_data_output</span><br> |
---|
1014 | </td> <td style="vertical-align: top;">R<br> </td> |
---|
1015 | <td style="vertical-align: top;"><span style="font-style: italic;">9999999.9</span><br> |
---|
1016 | </td> <td style="vertical-align: top;"><p lang="en-GB"><font face="Thorndale"><font face="Thorndale, serif">Temporal interval</font> |
---|
1017 | at which data (3d volume data (instantaneous or time |
---|
1018 | averaged), |
---|
1019 | cross sections (instantaneous or time averaged), vertical profiles, |
---|
1020 | spectra) shall be output (</font>in <font face="Thorndale">s). </font></p> |
---|
1021 | <span lang="en-GB"><font face="Thorndale">If |
---|
1022 | data output is switched on (see </font></span><a href="chapter_4.2.html#data_output"><span lang="en-GB"><font face="Thorndale">data_output</font></span></a><span lang="en-GB"><font face="Thorndale">, <a href="#data_output_pr">data_output_pr</a>, <a href="#data_output_sp">data_output_sp</a>, and </font></span><a href="chapter_4.2.html#section_xy"><span lang="en-GB"><font face="Thorndale">section_xy</font></span></a><span lang="en-GB"><font face="Thorndale">), this |
---|
1023 | parameter can be used to |
---|
1024 | assign the temporal interval at which these data shall be |
---|
1025 | output. </font></span><span lang="en-GB"><font face="Thorndale">Output can be skipped at the beginning of a |
---|
1026 | simulation using parameter <a href="#skip_time_data_output">skip_time_data_output</a>, |
---|
1027 | which has zero value by default. </font></span><span lang="en-GB"><font face="Thorndale">Reference |
---|
1028 | time is the beginning of the simulation, i.e. output |
---|
1029 | takes place at times t = <b>skip_time_data_output + |
---|
1030 | dt_data_output</b>, <span style="font-weight: bold;">skip_time_data_output</span> |
---|
1031 | + 2*<b>dt_data_output</b>, <span style="font-weight: bold;">skip_time_data_output</span> |
---|
1032 | + 3*<b>dt_data_output</b>, |
---|
1033 | etc. Since output is only done at the discrete time levels given by |
---|
1034 | the timestep used, the actual output times can slightly |
---|
1035 | deviate |
---|
1036 | from these theoretical values</font></span><a href="chapter_4.2.html#dt_dopr_zeitpunkte"><span lang="en-GB"></span></a><span lang="en-GB"><font face="Thorndale">.<br><br>Individual temporal |
---|
1037 | intervals for the different output quantities can be assigned using |
---|
1038 | parameters <a href="#dt_do3d">dt_do3d</a>, <a href="#dt_do2d_xy">dt_do2d_xy</a>, <a href="dt_do2d_xz">dt_do2d_xz</a>, <a href="#dt_do2d_yz">dt_do2d_yz</a>, <a href="#dt_dopr">dt_dopr</a>, <a href="#dt_dosp">dt_dosp</a>, |
---|
1039 | and <a href="#dt_data_output_av">dt_data_output_av</a>.</font></span> |
---|
1040 | </td> </tr> <tr> <td style="vertical-align: top;"><a name="dt_data_output_av"></a><span style="font-weight: bold;">dt_data_output_av</span><br> |
---|
1041 | </td> <td style="vertical-align: top;">R<br> </td> |
---|
1042 | <td style="vertical-align: top;"><i>value of |
---|
1043 | <a href="chapter_4.2.html#dt_data_output">dt_data_<br>output</a></i> |
---|
1044 | </td> <td style="vertical-align: top;"><p lang="en-GB"><font face="Thorndale"><font face="Thorndale, serif">Temporal interval</font> |
---|
1045 | at which time averaged 3d volume data and/or 2d cross section data |
---|
1046 | shall be output (</font>in <font face="Thorndale">s). </font></p><span lang="en-GB"><font face="Thorndale">If data |
---|
1047 | output of time averaged 2d and 3d data is switched on (see </font></span><a href="chapter_4.2.html#data_output"><span lang="en-GB"><font face="Thorndale">data_output</font></span></a> <span lang="en-GB"><font face="Thorndale">and </font></span><a href="chapter_4.2.html#section_xy"><span lang="en-GB"><font face="Thorndale">section_xy</font></span></a><span lang="en-GB"><font face="Thorndale">), this |
---|
1048 | parameter can be used to |
---|
1049 | assign the temporal interval at which they shall be |
---|
1050 | output. </font></span><span lang="en-GB"><font face="Thorndale">Output can be skipped at the beginning of a |
---|
1051 | simulation using parameter <a href="#skip_time_data_output_av">skip_time_data_output_av</a>, |
---|
1052 | which has zero value by default. </font></span><span lang="en-GB"><font face="Thorndale">Reference |
---|
1053 | time is the beginning of the simulation, i.e. output |
---|
1054 | takes place at times t = <b>skip_time_data_output_av + |
---|
1055 | dt_data_output_av</b>, <span style="font-weight: bold;">skip_time_data_output_av</span> |
---|
1056 | + 2*<b>dt_data_output_av</b>, <span style="font-weight: bold;">skip_time_data_output_av</span> |
---|
1057 | + 3*<b>dt_data_output_av</b>, |
---|
1058 | etc. Since output is only done at the discrete time levels given by |
---|
1059 | the timestep used, the actual output times can slightly |
---|
1060 | deviate from |
---|
1061 | these theoretical values</font></span><a href="chapter_4.2.html#dt_dopr_zeitpunkte"><span lang="en-GB"></span></a><span lang="en-GB"><font face="Thorndale">.<br><br></font></span>The |
---|
1062 | length of the averaging interval is controlled via parameter <a href="chapter_4.2.html#averaging_interval">averaging_interval</a>.</td> |
---|
1063 | </tr><tr> <td style="vertical-align: top;"> <p><a name="dt_disturb"></a><b>dt_disturb</b></p> |
---|
1064 | </td> <td style="vertical-align: top;">R</td> |
---|
1065 | <td style="vertical-align: top;"><i>9999999.9</i></td> |
---|
1066 | <td style="vertical-align: top;"> <p lang="en-GB"><font face="Thorndale"><font face="Thorndale, serif">Temporal |
---|
1067 | interval</font> at which random |
---|
1068 | perturbations are to be imposed on the horizontal velocity field |
---|
1069 | (</font>in <font face="Thorndale">s). </font> |
---|
1070 | </p> <p><span lang="en-GB"><font face="Thorndale, serif">The parameter </font></span><a href="#create_disturbances"><span lang="en-GB"><font face="Thorndale, serif">create_disturbances</font></span></a><font face="Thorndale, serif"><span lang="en-GB"> |
---|
1071 | describes how to impose |
---|
1072 | random perturbations to the horizontal velocity field</span></font><font face="Thorndale, serif"><span lang="en-GB">.</span> |
---|
1073 | </font> </p> </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="dt_dopr"></a><b>dt_dopr</b></p> |
---|
1074 | </td> <td style="vertical-align: top;">R</td> |
---|
1075 | <td style="vertical-align: top;"><i>value of |
---|
1076 | <a href="#dt_data_output">dt_data_<br>output</a></i></td> |
---|
1077 | <td style="vertical-align: top;"> <p><span lang="en-GB"><font face="Thorndale">Temporal |
---|
1078 | interval at |
---|
1079 | which data of vertical profiles shall be output (to local |
---|
1080 | file <a href="chapter_3.4.html#DATA_1D_PR_NETCDF">DATA_1D_PR_NETCDF</a> |
---|
1081 | or/and </font></span><a href="chapter_3.4.html#PLOT1D_DATA"><span lang="en-GB"><font face="Thorndale">PLOT1D_DATA</font></span></a><span lang="en-GB"><font face="Thorndale">) (</font></span>in |
---|
1082 | <span lang="en-GB"><font face="Thorndale">s). |
---|
1083 | </font></span> </p> <p><span lang="en-GB"><font face="Thorndale">If output of |
---|
1084 | horizontally averaged vertical profiles is switched on (see </font></span><a href="chapter_4.2.html#data_output_pr"><span lang="en-GB"><font face="Thorndale">data_output_pr</font></span></a><span lang="en-GB"><font face="Thorndale">), </font></span><span lang="en-GB"><font face="Thorndale">this |
---|
1085 | parameter can be used to |
---|
1086 | assign the temporal interval at which profile data shall be output.</font></span><span lang="en-GB"><font face="Thorndale"> </font></span><span lang="en-GB"><font face="Thorndale">Output can |
---|
1087 | be skipped at the beginning of a simulation using parameter <a href="#skip_time_dopr">skip_time_dopr</a>, which has |
---|
1088 | zero value by default. </font></span><span lang="en-GB"></span><span lang="en-GB"><font face="Thorndale">Reference |
---|
1089 | time is the beginning |
---|
1090 | of the simulation, thus t = 0, </font></span><span lang="en-GB"><font face="Thorndale">i.e. output |
---|
1091 | takes place at times t = <b>skip_time_dopr + dt_dopr</b>, <span style="font-weight: bold;">skip_time_dopr</span> + 2*<b>dt_dopr</b>, |
---|
1092 | <span style="font-weight: bold;">skip_time_dopr</span> |
---|
1093 | + 3*<b>dt_dopr</b>, |
---|
1094 | etc.</font></span><span lang="en-GB"><font face="Thorndale"> Since |
---|
1095 | profiles can not be calculated for times lying within a time step |
---|
1096 | interval, the output times can deviate from these theoretical values. |
---|
1097 | If a time step ranges from t = 1799.8 to t = 1800.2, then in the |
---|
1098 | example above the output would take place at t = 1800.2. In general, |
---|
1099 | the output always lie between t = 1800.0 and t = 1800.0 + </font></span><a href="chapter_4.1.html#dt"><span lang="en-GB"><font face="Thorndale">dt</font></span></a><span lang="en-GB"><font face="Thorndale">. If the |
---|
1100 | model uses a variable time step, these |
---|
1101 | deviations from the theoretical output times will of course be |
---|
1102 | different for each output time.<br> </font></span></p> |
---|
1103 | <p><span lang="en-GB"><font face="Thorndale">In |
---|
1104 | order to |
---|
1105 | guarantee an output of profile data at the end of a simulation (see </font></span><font><a href="chapter_4.1.html#end_time"><span lang="en-GB"><font face="Thorndale">end_time</font></span></a></font><span lang="en-GB"><font face="Thorndale">) in any way</font></span><span lang="en-GB"><font face="Thorndale">, |
---|
1106 | <span style="font-weight: bold;">end_time</span> |
---|
1107 | should be equal or a little bit |
---|
1108 | larger than the respective theoretical output time. For example, if <b>dt_dopr</b> |
---|
1109 | = <i>900.0</i><span style="font-style: italic;"> |
---|
1110 | </span>and 3600.0 |
---|
1111 | seconds are to be simulated, then <b>end_time</b> |
---|
1112 | >= 3600.0 should be chosen.</font></span><a href="chapter_4.1.html#dt"><span lang="en-GB"></span></a><span lang="en-GB"><font face="Thorndale"><span style="font-weight: bold;"></span> </font></span> |
---|
1113 | </p> <p><span lang="en-GB"><font face="Thorndale">A selection of |
---|
1114 | profiles to be output can be done via parameter </font></span><a href="chapter_4.2.html#data_output_pr"><span lang="en-GB"><font face="Thorndale">data_output_pr</font></span></a><span lang="en-GB"><font face="Thorndale">. </font></span> |
---|
1115 | </p> </td> </tr> <tr> <td style="vertical-align: top;"><a name="dt_dopr_listing"></a><span style="font-weight: bold;">dt_dopr_listing</span><br> |
---|
1116 | </td> <td style="vertical-align: top;">R<br> </td> |
---|
1117 | <td style="vertical-align: top;"><i>9999999.9</i></td> |
---|
1118 | <td style="vertical-align: top;"> <p><span lang="en-GB"><font face="Thorndale, serif">Temporal |
---|
1119 | interval</font> at which data <font face="Thorndale">of |
---|
1120 | vertical |
---|
1121 | profiles shall be output (output for printouts, local file </font></span><a href="chapter_3.4.html#LIST_PROFIL"><span lang="en-GB"><font face="Thorndale">LIST_PROFIL</font></span></a><span lang="en-GB"><font face="Thorndale">) (</font></span>in |
---|
1122 | <span lang="en-GB"><font face="Thorndale">s). </font></span> |
---|
1123 | </p> <p>T<span lang="en-GB"></span><a href="chapter_4.2.html#pr1d"><span lang="en-GB"></span></a><span lang="en-GB"></span><span lang="en-GB"><font face="Thorndale">his |
---|
1124 | parameter can be used to |
---|
1125 | assign the temporal interval at which profile data shall be output.</font></span><span lang="en-GB"><font face="Thorndale"> Reference |
---|
1126 | time is the beginning |
---|
1127 | of the simulation, thus t = 0. For example if <b>dt_dopr_listing</b> |
---|
1128 | = 1800.0, |
---|
1129 | then output takes place at t = 1800.0, 3600.0, 5400.0, etc. Since |
---|
1130 | profiles can not be calculated for times lying within a time step |
---|
1131 | interval, the output times can deviate from these theoretical values. |
---|
1132 | If a time step ranges from t = 1799.8 to t = 1800.2, then in the |
---|
1133 | example above the output would take place at t = 1800.2. In general, |
---|
1134 | the output always lie between t = 1800.0 and t = 1800.0 + </font></span><a href="chapter_4.1.html#dt"><span lang="en-GB"><font face="Thorndale">dt</font></span></a> <span lang="en-GB"><font face="Thorndale">(numbers |
---|
1135 | are related to |
---|
1136 | the |
---|
1137 | example above). If the model uses a variable time step, these |
---|
1138 | deviations from the theoretical output times will of course be |
---|
1139 | different for each output time.<br> </font></span></p> |
---|
1140 | <p><span lang="en-GB"><font face="Thorndale">In |
---|
1141 | order to |
---|
1142 | guarantee an output of profile data at the end of a simulation (see </font></span><font><a href="chapter_4.1.html#end_time"><span lang="en-GB"><font face="Thorndale">end_time</font></span></a></font><span lang="en-GB"><font face="Thorndale">) in any way</font></span><span lang="en-GB"><font face="Thorndale">, |
---|
1143 | <span style="font-weight: bold;">end_time</span> |
---|
1144 | should be a little bit |
---|
1145 | larger than the respective theoretical output time. For example, if <b>dt_dopr_listing</b> |
---|
1146 | = <i>900.0</i><span style="font-style: italic;"> |
---|
1147 | </span>and 3600.0 |
---|
1148 | seconds are to be simulated, then it should be at least <b>end_time</b> |
---|
1149 | > 3600.0 + </font></span><a href="chapter_4.1.html#dt"><span lang="en-GB"><font face="Thorndale">dt</font></span></a><span lang="en-GB"><font face="Thorndale">. If |
---|
1150 | variable time steps are used |
---|
1151 | (which is the default), <span style="font-weight: bold;">dt</span> |
---|
1152 | should be properly estimated. </font></span> </p> |
---|
1153 | <p><span lang="en-GB"><font face="Thorndale">Data |
---|
1154 | and output |
---|
1155 | format of the file </font></span><a href="chapter_3.4.html#LIST_PROFIL"><span lang="en-GB"><font face="Thorndale">LIST_PROFIL</font></span></a> |
---|
1156 | <span lang="en-GB"><font face="Thorndale">is |
---|
1157 | internally fixed. In this file |
---|
1158 | the profiles of the most important model variables are arranged in |
---|
1159 | adjacent columns.</font></span> </p> </td> </tr> |
---|
1160 | <tr> <td style="vertical-align: top;"> <p><a name="dt_dots"></a><b>dt_dots</b></p> |
---|
1161 | </td> <td style="vertical-align: top;">R</td> |
---|
1162 | <td style="vertical-align: top;"><span style="font-style: italic;">see right</span></td> |
---|
1163 | <td style="vertical-align: top;"> <p lang="en-GB"><font face="Thorndale"><font face="Thorndale, serif">Temporal |
---|
1164 | interval</font> at which time series data shall be |
---|
1165 | output (</font>in <font face="Thorndale">s). </font> |
---|
1166 | </p> <p>The default interval for the output of timeseries |
---|
1167 | is calculated as shown below (this tries to minimize the number of |
---|
1168 | calls of <span style="font-family: Courier New,Courier,monospace;">flow_statistics</span>)</p><p style="font-family: Courier New,Courier,monospace;"> |
---|
1169 | IF ( <a href="#averaging_interval_pr">averaging_interval_pr</a> |
---|
1170 | == 0.0 ) THEN<br> |
---|
1171 | <span style="font-weight: bold;">dt_dots</span> = |
---|
1172 | MIN( <a href="#dt_run_control">dt_run_control</a>, <a href="#dt_dopr">dt_dopr</a> )<br> |
---|
1173 | ELSE<br> |
---|
1174 | <span style="font-weight: bold;">dt_dots</span> = |
---|
1175 | MIN( dt_run_control, <a href="#dt_averaging_input_pr">dt_averaging_input_pr</a> |
---|
1176 | )<br> |
---|
1177 | ENDIF</p><p>This parameter can be used to |
---|
1178 | assign the temporal interval at which data points shall be output. <span lang="en-GB"><font face="Thorndale">Reference |
---|
1179 | time is the beginning of |
---|
1180 | the simulation, i.e. output takes place at times t = <b>dt_dots</b>, |
---|
1181 | 2*<b>dt_dots</b>, 3*<b>dt_dots</b>, etc. The |
---|
1182 | actual output times can |
---|
1183 | deviate from these theoretical values (see </font></span><a href="#dt_dopr_zeitpunkte"><span lang="en-GB"><font face="Thorndale">dt_dopr</font></span></a><span lang="en-GB"><font face="Thorndale">). |
---|
1184 | Is <b>dt_dots</b> < </font></span><a href="chapter_4.1.html#dt"><span lang="en-GB"><font face="Thorndale">dt</font></span></a><span lang="en-GB"><font face="Thorndale">, then data |
---|
1185 | of the time series are |
---|
1186 | written after each time step (if this is requested it should be <b>dt_dots</b> |
---|
1187 | = <i>0</i>).</font></span></p><p><span lang="en-GB"><font face="Thorndale">The default |
---|
1188 | value of <span style="font-weight: bold;">dt_dots</span> |
---|
1189 | is calculated as follows:</font></span></p> |
---|
1190 | IF ( <a href="#averaging_interval_pr">averaging_interval_pr</a> |
---|
1191 | == 0.0 ) THEN<br> |
---|
1192 | <span style="font-weight: bold;">dt_dots</span> = |
---|
1193 | MIN( <a href="#dt_run_control">dt_run_control</a>, <a href="#dt_dopr">dt_dopr</a> )<br> |
---|
1194 | ELSE<br> |
---|
1195 | <span style="font-weight: bold;">dt_dots</span> = |
---|
1196 | MIN( <span style="font-weight: bold;">dt_run_control</span>, |
---|
1197 | <a href="#dt_averaging_input_pr">dt_averaging_input_pr</a> |
---|
1198 | )<br> |
---|
1199 | ENDIF<br><br>(which minimizes the number of calls of |
---|
1200 | routine flow_statistics).<br><p>By default time series data |
---|
1201 | is output to the local file <a href="chapter_3.4.html#DATA_1D_TS_NETCDF">DATA_1D_TS_NETCDF</a>. |
---|
1202 | Because of the default settings of <span style="font-weight: bold;">dt_dots</span>, |
---|
1203 | it will generally be created for each model run. The file's |
---|
1204 | format is NetCDF. Further details about processing NetCDF |
---|
1205 | data are given in chapter <a href="chapter_4.5.1.html">4.5.1</a>.</p>The |
---|
1206 | file contains the following timeseries quantities (the first column |
---|
1207 | gives the name of the quantities as used in the NetCDF file):<br><table style="text-align: left; width: 100%;" cellpadding="2" cellspacing="2"> <tbody> <tr> <td style="font-style: italic; vertical-align: middle;">E<br> |
---|
1208 | </td> <td style="vertical-align: top;">Total |
---|
1209 | kinetic energy of |
---|
1210 | the flow (in m<sup>2</sup>/s<sup>2</sup>) |
---|
1211 | (normalized with respect to the total number of grid points).</td> |
---|
1212 | </tr> <tr> <td style="font-style: italic; vertical-align: middle;">E*<br> |
---|
1213 | </td> <td style="vertical-align: top;">Perturbation |
---|
1214 | kinetic |
---|
1215 | energy of the flow (in m<sup>2</sup>/s<sup>2</sup>)<sup> |
---|
1216 | </sup>(normalized |
---|
1217 | with respect to the total number of grid |
---|
1218 | points)</td> </tr> <tr> <td style="vertical-align: top; font-style: italic;">dt<br> |
---|
1219 | </td> <td style="vertical-align: top;">Time step |
---|
1220 | size (in s).</td> </tr> <tr> <td style="vertical-align: top; font-style: italic;">u<sub>*</sub></td> |
---|
1221 | <td style="vertical-align: top;">Friction velocity (in |
---|
1222 | m/s) |
---|
1223 | (horizontal average).</td> </tr> <tr> <td style="vertical-align: top; font-style: italic;">w<sub>*</sub></td> |
---|
1224 | <td style="vertical-align: top;">Vertical velocity scale |
---|
1225 | of |
---|
1226 | the CBL (in m/s) (horizontal average)</td> </tr> <tr> |
---|
1227 | <td style="vertical-align: top; font-style: italic;">th<sub>*</sub></td> |
---|
1228 | <td style="vertical-align: top;">Temperature |
---|
1229 | scale (Prandtl layer), defined as <i>w"pt"0 |
---|
1230 | / </i><i>u<sub>*</sub></i> |
---|
1231 | (horizontal |
---|
1232 | average) (in K).</td> </tr> <tr> <td style="vertical-align: top; font-style: italic;">umax<br> |
---|
1233 | </td> <td style="vertical-align: top;">Maximum |
---|
1234 | u-component of the |
---|
1235 | velocity (in m/s).</td> </tr> <tr> <td style="vertical-align: top; font-style: italic;">vmax<br> |
---|
1236 | </td> <td style="vertical-align: top;">Maximum |
---|
1237 | v-component of the |
---|
1238 | velocity (in m/s).</td> </tr> <tr> <td style="vertical-align: top; font-style: italic;">wmax<br> |
---|
1239 | </td> <td style="vertical-align: top;">Maximum |
---|
1240 | w-component of the |
---|
1241 | velocity (in m/s).</td> </tr> <tr> <td style="vertical-align: top; font-style: italic;">div_old<br> |
---|
1242 | </td> <td style="vertical-align: top;">Divergence |
---|
1243 | of the velocity |
---|
1244 | field before the pressure |
---|
1245 | solver has been called (normalized with respect to the total number of |
---|
1246 | grid points) (in 1/s).</td> </tr> <tr> <td style="vertical-align: top; font-style: italic;">div_new</td> |
---|
1247 | <td style="vertical-align: top;">Divergence of the |
---|
1248 | velocity |
---|
1249 | field after the pressure |
---|
1250 | solver has been called (normalized with respect to the total number of |
---|
1251 | grid points) (in 1/s).</td> </tr> <tr> <td style="vertical-align: top; font-style: italic;">z_i_wpt</td> |
---|
1252 | <td style="vertical-align: top;">Height of the convective |
---|
1253 | boundary layer (horizontal average) |
---|
1254 | determined by the height of the minimum sensible heat flux (in m).</td> |
---|
1255 | </tr> <tr> <td style="vertical-align: top; font-style: italic;">z_i_pt</td> |
---|
1256 | <td style="vertical-align: top;">Height of the convective |
---|
1257 | boundary layer (horizontal average) |
---|
1258 | determined by the temperature profile (in m).</td> </tr> <tr> |
---|
1259 | <td style="vertical-align: top; font-style: italic;">w"pt"0</td> |
---|
1260 | <td style="vertical-align: top;">Subgrid-scale sensible |
---|
1261 | heat flux near the surface (horizontal |
---|
1262 | average) |
---|
1263 | between z = 0 and z = z<sub>p</sub> = zu(1) (there it |
---|
1264 | corresponds to |
---|
1265 | the total heat flux) (in K m/s).</td> </tr> <tr> <td style="vertical-align: top; font-style: italic;">w"pt"</td> |
---|
1266 | <td style="vertical-align: top;">Subgrid-scale heat flux |
---|
1267 | (horizontal average) for z = zw(1) (in K |
---|
1268 | m/s).</td> </tr> <tr> <td style="vertical-align: top; font-style: italic;">wpt</td> |
---|
1269 | <td style="vertical-align: top;">Total heat flux |
---|
1270 | (horizontal average) for z = zw(1) (in K m/s).</td> </tr> <tr> |
---|
1271 | <td style="vertical-align: top; font-style: italic;">pt(0)</td> |
---|
1272 | <td style="vertical-align: top;">Potential temperature at |
---|
1273 | the surface (horizontal average) (in K).</td> </tr> <tr> |
---|
1274 | <td style="vertical-align: top; font-style: italic;">pt(zp)</td> |
---|
1275 | <td style="vertical-align: top;">Potential temperature for |
---|
1276 | z = zu(1) (horizontal average) (in K).</td> </tr> <tr> |
---|
1277 | <td style="vertical-align: top; font-style: italic;">splptx</td> |
---|
1278 | <td style="vertical-align: top;">Percentage of grid points |
---|
1279 | using upstream scheme along x with |
---|
1280 | upstream-spline advection switched on.</td> </tr> <tr> |
---|
1281 | <td style="vertical-align: top; font-style: italic;">splpty</td> |
---|
1282 | <td style="vertical-align: top;">Percentage of grid points |
---|
1283 | using upstream scheme along y with |
---|
1284 | upstream-spline |
---|
1285 | advection switched on.</td> </tr> <tr> <td style="vertical-align: top; font-style: italic;">splptz</td> |
---|
1286 | <td style="vertical-align: top;">Percentage of grid points |
---|
1287 | using upstream scheme along z with |
---|
1288 | upstream-spline |
---|
1289 | advection switched on.<br> </td> </tr> <tr> <td style="vertical-align: top; font-style: italic;">L</td> |
---|
1290 | <td style="vertical-align: top;">Monin-Obukhov length.</td> |
---|
1291 | </tr> </tbody> </table><br>Additionally, the |
---|
1292 | user can add his own timeseries quantities to the file, by using the |
---|
1293 | user-interface subroutines<span style="font-family: Courier New,Courier,monospace;"> <a href="chapter_3.5.1.html#user_init">user_init</a> </span>and<span style="font-family: Courier New,Courier,monospace;"> <a href="chapter_3.5.1.html#user_statistics">user_statistics</a></span>. |
---|
1294 | These routines contain (as comment lines) a simple example how to do |
---|
1295 | this.<br><br>Time series data refers to the total |
---|
1296 | domain, but time series for subdomains can also be output (see <a href="chapter_4.1.html#statistic_regions">statistic_regions</a>). |
---|
1297 | However, the following time series always present the values of the |
---|
1298 | total model domain (even with output for subdomains): <i>umax</i>, |
---|
1299 | <i>vmax</i>, <i>wmax</i>, <i>div_old</i>, |
---|
1300 | <i>div_new</i>.</td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="dt_do2d_xy"></a><b>dt_do2d_xy</b></p> |
---|
1301 | </td> <td style="vertical-align: top;">R</td> |
---|
1302 | <td style="vertical-align: top;"><i>value of |
---|
1303 | <a href="chapter_4.2.html#dt_data_output">dt_data_<br>output</a></i></td> |
---|
1304 | <td style="vertical-align: top;"> <p lang="en-GB"><font face="Thorndale"><font face="Thorndale, serif">Temporal |
---|
1305 | interval</font> at which horizontal cross section data |
---|
1306 | shall be output (</font>in <font face="Thorndale">s). |
---|
1307 | </font> </p> <p><span lang="en-GB"><font face="Thorndale">If output of |
---|
1308 | horizontal cross sections is switched on (see </font></span><a href="#data_output"><span lang="en-GB"><font face="Thorndale">data_output</font></span></a> |
---|
1309 | <span lang="en-GB"><font face="Thorndale">and |
---|
1310 | </font></span><a href="#section_xy"><span lang="en-GB"><font face="Thorndale">section_xy</font></span></a><span lang="en-GB"><font face="Thorndale">), this |
---|
1311 | parameter can be used to |
---|
1312 | assign the temporal interval at which cross section data shall be |
---|
1313 | output. </font></span><span lang="en-GB"><font face="Thorndale">Output can be skipped at the beginning of a |
---|
1314 | simulation using parameter <a href="#skip_time_do2d_xy">skip_time_do2d_xy</a>, |
---|
1315 | which has zero value by default. </font></span><span lang="en-GB"><font face="Thorndale">Reference |
---|
1316 | time is the beginning of the simulation, i.e. output |
---|
1317 | takes place at times t = <b>skip_time_do2d_xy + dt_do2d_xy</b>, |
---|
1318 | <span style="font-weight: bold;">skip_time_do2d_xy</span> |
---|
1319 | + 2*<b>dt_do2d_xy</b>, <span style="font-weight: bold;">skip_time_do2d_xy</span> |
---|
1320 | + 3*<b>dt_do2d_xy</b>, |
---|
1321 | etc. The actual output times can deviate from these theoretical values |
---|
1322 | (see </font></span><a href="#dt_dopr_zeitpunkte"><span lang="en-GB"><font face="Thorndale">dt_dopr</font></span></a><span lang="en-GB"><font face="Thorndale">).<br> |
---|
1323 | </font></span></p> <p><span lang="en-GB"><font face="Thorndale">Parameter </font></span><a href="#do2d_at_begin"><span lang="en-GB"><font face="Thorndale">do2d_at_begin</font></span></a> |
---|
1324 | has to be used if an additional output is wanted at the start of a run <span lang="en-GB"><font face="Thorndale">(thus at |
---|
1325 | the time t = 0 or at the |
---|
1326 | respective starting times of restart runs).</font></span> </p> |
---|
1327 | </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="dt_do2d_xz"></a><b>dt_do2d_xz</b></p> |
---|
1328 | </td> <td style="vertical-align: top;">R</td> |
---|
1329 | <td style="vertical-align: top;"><i>value of |
---|
1330 | <a href="chapter_4.2.html#dt_data_output">dt_data_<br>output</a></i></td> |
---|
1331 | <td style="vertical-align: top;"> <p lang="en-GB"><font face="Thorndale"><font face="Thorndale, serif">Temporal |
---|
1332 | interval</font> at which vertical cross sections data |
---|
1333 | (xz) shall be output (</font>in <font face="Thorndale">s). |
---|
1334 | </font> </p> <p><span lang="en-GB"><font face="Thorndale">If output of |
---|
1335 | horizontal cross sections is switched on (see </font></span><a href="#data_output"><span lang="en-GB"><font face="Thorndale">data_output</font></span></a> |
---|
1336 | <span lang="en-GB"><font face="Thorndale">and |
---|
1337 | </font></span><a href="#section_xz"><span lang="en-GB"><font face="Thorndale">section_xz</font></span></a><span lang="en-GB"><font face="Thorndale">), |
---|
1338 | this parameter can be used to assign the temporal interval at which |
---|
1339 | cross section data shall be output. </font></span><span lang="en-GB"><font face="Thorndale">Output can |
---|
1340 | be skipped at the beginning of a simulation using parameter <a href="#skip_time_do2d_xz">skip_time_do2d_xz</a>, which |
---|
1341 | has zero value by default. </font></span><span lang="en-GB"></span><span lang="en-GB"><font face="Thorndale">Reference time is the beginning of |
---|
1342 | the simulation, i.e. output takes place at times t = <b>skip_time_do2d_xz |
---|
1343 | + dt_do2d_xz</b>, |
---|
1344 | <span style="font-weight: bold;">skip_time_do2d_xz</span> |
---|
1345 | + 2*<b>dt_do2d_xz</b>, <span style="font-weight: bold;">skip_time_do2d_xz</span> |
---|
1346 | + 3*<b>dt_do2d_xz</b>, etc. The actual output times |
---|
1347 | can deviate from these theoretical values (see </font></span><a href="#dt_dopr_zeitpunkte"><span lang="en-GB"><font face="Thorndale">dt_dopr</font></span></a><span lang="en-GB"><font face="Thorndale">).<br> |
---|
1348 | </font></span></p> <p><span lang="en-GB"><font face="Thorndale">Parameter </font></span><a href="#do2d_at_begin"><span lang="en-GB"><font face="Thorndale">do2d_at_begin</font></span></a> |
---|
1349 | has to be used if an additional output is wanted at the start of a run <span lang="en-GB"><font face="Thorndale">(thus at |
---|
1350 | the time t = 0 or at the |
---|
1351 | respective starting times of restart runs).</font></span> </p> |
---|
1352 | </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="dt_do2d_yz"></a><b>dt_do2d_yz</b></p> |
---|
1353 | </td> <td style="vertical-align: top;">R</td> |
---|
1354 | <td style="vertical-align: top;"><i>value of |
---|
1355 | <a href="chapter_4.2.html#dt_data_output">dt_data_<br>output</a></i></td> |
---|
1356 | <td style="vertical-align: top;"> <p lang="en-GB"><font face="Thorndale"><font face="Thorndale, serif">Temporal |
---|
1357 | interval</font> at which vertical cross section data |
---|
1358 | (yz) shall be output (</font>in s<font face="Thorndale">). |
---|
1359 | </font> </p> <p><span lang="en-GB"><font face="Thorndale">If output of |
---|
1360 | horizontal cross sections is switched on (see </font></span><a href="#data_output"><span lang="en-GB"><font face="Thorndale">data_output</font></span></a> |
---|
1361 | <span lang="en-GB"><font face="Thorndale">and |
---|
1362 | </font></span><a href="#section_yz"><span lang="en-GB"><font face="Thorndale">section_yz</font></span></a><span lang="en-GB"><font face="Thorndale">), |
---|
1363 | this parameter can be used to assign the temporal interval at which |
---|
1364 | cross section data shall be output. </font></span><span lang="en-GB"><font face="Thorndale">Output can |
---|
1365 | be skipped at the beginning of a simulation using parameter <a href="#skip_time_do2d_yz">skip_time_do2d_yz</a>, which |
---|
1366 | has zero value by default. </font></span><span lang="en-GB"></span><span lang="en-GB"></span><span lang="en-GB"><font face="Thorndale">Reference |
---|
1367 | time is the beginning of |
---|
1368 | the simulation, i.e. output takes place at times t = <b>skip_time_do2d_yz |
---|
1369 | + dt_do2d_yz</b>, |
---|
1370 | <span style="font-weight: bold;">skip_time_do2d_yz</span> |
---|
1371 | + 2*<b>dt_do2d_yz</b>, <span style="font-weight: bold;">skip_time_do2d_yz |
---|
1372 | </span>+ 3*<b>dt_do2d_yz</b>, etc. The actual output |
---|
1373 | times |
---|
1374 | can deviate from these theoretical values (see </font></span><a href="#dt_dopr_zeitpunkte"><span lang="en-GB"><font face="Thorndale">dt_dopr</font></span></a><span lang="en-GB"><font face="Thorndale">).<br> |
---|
1375 | </font></span></p> <p><span lang="en-GB"><font face="Thorndale">Parameter </font></span><a href="#do2d_at_begin"><span lang="en-GB"><font face="Thorndale">do2d_at_begin</font></span></a> |
---|
1376 | has to be used if an additional output is wanted at the start of a run <span lang="en-GB"><font face="Thorndale">(thus at |
---|
1377 | the time t = 0 or at the |
---|
1378 | respective starting times of restart runs).</font></span> </p> |
---|
1379 | </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="dt_do3d"></a><b>dt_do3d</b></p> |
---|
1380 | </td> <td style="vertical-align: top;">R</td> |
---|
1381 | <td style="vertical-align: top;"><i>value of |
---|
1382 | <a href="chapter_4.2.html#dt_data_output">dt_data_<br>output</a></i></td> |
---|
1383 | <td style="vertical-align: top;"> <p lang="en-GB"><font face="Thorndale"><font face="Thorndale, serif">Temporal |
---|
1384 | interval</font> at which 3d volume data shall be output (</font>in |
---|
1385 | <font face="Thorndale">s). </font> </p> |
---|
1386 | <p><span lang="en-GB"><font face="Thorndale">If |
---|
1387 | output of |
---|
1388 | 3d-volume data is switched on (see </font></span><font><a href="#data_output"><span lang="en-GB"><font face="Thorndale">data_output</font></span></a>)<span style="font-family: thorndale;">, this parameter can be used |
---|
1389 | to assign |
---|
1390 | th</span></font><span lang="en-GB"><font face="Thorndale">e temporal |
---|
1391 | interval at which 3d-data shall be output. </font></span><span lang="en-GB"><font face="Thorndale">Output can |
---|
1392 | be skipped at the beginning of a simulation using parameter <a href="#skip_time_do3d">skip_time_do3d</a>, which has |
---|
1393 | zero value by default. </font></span><span lang="en-GB"></span><span lang="en-GB"></span><span lang="en-GB"></span><span lang="en-GB"><font face="Thorndale">Reference |
---|
1394 | time is the |
---|
1395 | beginning of the simulation, i.e. output takes place at times t = <b>skip_time_do3d |
---|
1396 | + dt_do3d</b>, |
---|
1397 | <span style="font-weight: bold;">skip_time_do3d</span> |
---|
1398 | + 2*<b>dt_do3d</b>, <span style="font-weight: bold;">skip_time_do3d</span> |
---|
1399 | + 3*<b>dt_do3d</b>, etc. The actual output times can |
---|
1400 | deviate from these theoretical values (see </font></span><a href="#dt_dopr_zeitpunkte"><span lang="en-GB"><font face="Thorndale">dt_dopr</font></span></a><span lang="en-GB"><font face="Thorndale">). <br> |
---|
1401 | </font></span></p> <p><span lang="en-GB"><font face="Thorndale">Parameter </font></span><a href="#do3d_at_begin"><span lang="en-GB"><font face="Thorndale">do3d_at_begin</font></span></a> |
---|
1402 | has to be used if an additional output is wanted at the start of a run <span lang="en-GB"><font face="Thorndale">(thus at |
---|
1403 | the time t = 0 or at the |
---|
1404 | respective starting times of restart runs).</font></span> </p> |
---|
1405 | </td> </tr> <tr><td style="vertical-align: top;"><a name="dt_max"></a><span style="font-weight: bold;">dt_max</span></td><td style="vertical-align: top;">R</td><td style="vertical-align: top;"><span style="font-style: italic;">20.0</span></td><td>Maximum |
---|
1406 | allowed value of the timestep (in s).<br><br>By default, |
---|
1407 | the maximum timestep is restricted to be 20 s. This might be o.k. for |
---|
1408 | simulations of any kind of atmospheric turbulence but may have to be |
---|
1409 | changed for other situations.</td></tr><tr> <td style="vertical-align: top;"> <p><a name="dt_restart"></a><b>dt_restart</b></p> |
---|
1410 | </td> <td style="vertical-align: top;">R</td> |
---|
1411 | <td style="vertical-align: top;"><i>9999999.9</i></td> |
---|
1412 | <td style="vertical-align: top;"> <p lang="en-GB"><font face="Thorndale"><font face="Thorndale, serif">Temporal |
---|
1413 | interval</font> at which a new |
---|
1414 | restart run is to be carried out (</font>in <font face="Thorndale">s). </font> </p> <p><span lang="en-GB"><font face="Thorndale">For a |
---|
1415 | description |
---|
1416 | how to assign restart times manually see run time parameter </font></span><a href="#restart_time"><span lang="en-GB"><font face="Thorndale">restart_time</font></span></a><span lang="en-GB"><font face="Thorndale">. <span style="font-weight: bold;">dt_restart</span> |
---|
1417 | does not show any effect, if <span style="font-weight: bold;">restart_time</span> |
---|
1418 | has not been set.</font></span> </p> </td> </tr> |
---|
1419 | <tr> <td style="vertical-align: top;"> <p><a name="dt_run_control"></a><b>dt_run_control</b></p> |
---|
1420 | </td> <td style="vertical-align: top;">R</td> |
---|
1421 | <td style="vertical-align: top;"><i>60.0</i></td> |
---|
1422 | <td style="vertical-align: top;"> <p lang="en-GB"><font face="Thorndale"><font face="Thorndale, serif">Temporal |
---|
1423 | interval</font> at which run control |
---|
1424 | output is to be made (</font>in <font face="Thorndale">s). |
---|
1425 | </font> </p> <p><span lang="en-GB"><font face="Thorndale">Run control |
---|
1426 | information is output to the local ASCII-file </font></span><a href="chapter_3.4.html#RUN_CONTROL"><span lang="en-GB"><font face="Thorndale">RUN_CONTROL</font></span></a><span lang="en-GB"><font face="Thorndale">. At each |
---|
1427 | output time, one line |
---|
1428 | with information about the size of the time step, maximum speeds, total |
---|
1429 | kinetic energy etc. is written to this file. Reference time is the |
---|
1430 | beginning of the simulation, i.e. output takes place at times t = <b>dt_run_control</b>, |
---|
1431 | 2*<b>dt_run_control</b>, 3*<b>dt_run_control</b>, |
---|
1432 | etc., and always at |
---|
1433 | the beginning of a model run (thus at the time t = 0 or at the |
---|
1434 | respective starting times of restart runs). The actual output times can |
---|
1435 | deviate from these theoretical values (see </font></span><a href="#dt_dopr_zeitpunkte"><span lang="en-GB"><font face="Thorndale">dt_dopr</font></span></a><span lang="en-GB"><font face="Thorndale">).<br> |
---|
1436 | </font></span></p> <p><span lang="en-GB"><font face="Thorndale">Run control |
---|
1437 | information is output after each time step can be achieved via <b>dt_run_control</b> |
---|
1438 | = <i>0.0</i>.</font></span> </p> </td> |
---|
1439 | </tr> <tr> <td style="vertical-align: top;"><p><a name="end_time"></a><b>end_time</b></p> |
---|
1440 | </td> <td style="vertical-align: top;">R</td> |
---|
1441 | <td style="vertical-align: top;"><i>0.0</i></td> |
---|
1442 | <td style="vertical-align: top;"> <p lang="en-GB"><font face="Thorndale">Simulation time of the 3D |
---|
1443 | model (</font>in <font face="Thorndale">s). |
---|
1444 | </font> </p> <p><span lang="en-GB"><font face="Thorndale">The simulation time |
---|
1445 | is starting from the beginning of the initialization run (t = 0), not |
---|
1446 | starting from the beginning of the respective restart run.</font></span> |
---|
1447 | </p> </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="force_print_header"></a><b>force_print_header</b></p> |
---|
1448 | </td> <td style="vertical-align: top;">L</td> |
---|
1449 | <td style="vertical-align: top;"><i>.F.</i></td> |
---|
1450 | <td style="vertical-align: top;"> <p>Steering of |
---|
1451 | header output to the local file <a href="chapter_3.4.html#RUN_CONTROL">RUN_CONTROL</a>. |
---|
1452 | </p> <p>By default, informations about the model |
---|
1453 | parameters in use are |
---|
1454 | output to the beginning of file RUN_CONTROL for initial runs only |
---|
1455 | (these informations are identical to that which are output to the local |
---|
1456 | file <a href="chapter_3.4.html#HEADER">HEADER</a>). |
---|
1457 | With <b>force_print_header</b> = <i>.T.</i>, |
---|
1458 | these informations are |
---|
1459 | also output to <a href="chapter_3.4.html#RUN_CONTROL">RUN_CONTROL</a> |
---|
1460 | at restart runs.</p> </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="mg_cycles"></a><b>mg_cycles</b></p> |
---|
1461 | </td> <td style="vertical-align: top;">I</td> |
---|
1462 | <td style="vertical-align: top;"><i>-1</i></td> |
---|
1463 | <td style="vertical-align: top;"> <p>Number of |
---|
1464 | cycles to be used with the multi-grid scheme.<br> <br> |
---|
1465 | This parameter determines the number of cycles to be carried out in the |
---|
1466 | multi-grid method used for solving the Poisson equation for |
---|
1467 | perturbation pressure (see <a href="#psolver">psolver</a>). |
---|
1468 | The type of the cycles can be set with <a href="#cycle_mg">cycle_mg</a>.<br> |
---|
1469 | </p> <br>By default (<b>mg_cyles</b> = <i>- |
---|
1470 | 1</i>), the |
---|
1471 | number of cycles |
---|
1472 | depends on the requested accuracy of the scheme (see <a href="#residual_limit">residual_limit</a>) |
---|
1473 | and may vary from time step to time step. In this case, the CPU time |
---|
1474 | for a run will be difficult to estimate, since it heavily depends on |
---|
1475 | the total number of the cycles to be carried out.<br> <br> |
---|
1476 | By assigning <b>mg_cycles</b> a value (>=<span style="font-style: italic;">1</span>), the number of |
---|
1477 | cycles can be |
---|
1478 | fixed so that the CPU time can be clearly estimated. <br> <br> |
---|
1479 | <b>Note:</b> When using a fixed number of cycles, the user |
---|
1480 | must |
---|
1481 | examine the local file <a href="chapter_3.4.html#RUN_CONTROL">RUN_CONTROL</a> |
---|
1482 | regularly to check whether the divergence of the velocity field is |
---|
1483 | sufficiently reduced by the pressure solver. It should be reduced at |
---|
1484 | least by two orders of magnitude. For cyclic boundary conditions along |
---|
1485 | both horizontal directions (see <a href="chapter_4.1.html#bc_lr">bc_lr</a> |
---|
1486 | and <a href="chapter_4.1.html#bc_ns">bc_ns</a>) <span style="font-weight: bold;">mg_cycles</span> = <span style="font-style: italic;">2</span> is typically a |
---|
1487 | good choice, for |
---|
1488 | non-cyclic lateral boundary conditions <span style="font-weight: bold;">mg_cycles</span> |
---|
1489 | = <span style="font-style: italic;">4</span> may be |
---|
1490 | sufficient.</td> </tr> <tr> <td style="vertical-align: top;"><a name="mg_switch_to_pe0_level"></a><b>mg_switch_to_pe0_<br> |
---|
1491 | level</b></td> <td style="vertical-align: top;">I</td> |
---|
1492 | <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;">Grid |
---|
1493 | level at which data shall be gathered on PE0.<br> <br> |
---|
1494 | In case of a run using several PEs and the multigrid method for solving |
---|
1495 | the Poisson equation for perturbation pressure (see <a href="#psolver">psolver</a>), |
---|
1496 | the value of this parameter defines on which grid level the data are |
---|
1497 | gathered on PE0 in order to allow for a further coarsening of the grid. |
---|
1498 | The finest grid defines the largest grid level. By default, the |
---|
1499 | gathering level is determined automatically and displayed in file <a href="chapter_3.4.html#RUN_CONTROL">RUN_CONTROL</a>. |
---|
1500 | It is only possible to gather data from a level larger than the one |
---|
1501 | determined automatically. A test run may be neccessary to determine |
---|
1502 | this level.</td> </tr> <tr> <td style="vertical-align: top;"><a name="netcdf_64bit"></a><span style="font-weight: bold;">netcdf_64bit</span><br> |
---|
1503 | </td> <td style="vertical-align: top;">L<br> </td> |
---|
1504 | <td style="vertical-align: top;"><span style="font-style: italic;">.F.</span><br> </td> |
---|
1505 | <td style="vertical-align: top;">NetCDF files will have 64 |
---|
1506 | bit offset format.<br><br>By |
---|
1507 | default, the maximum file size of the NetCDF files opened by PALM is 2 |
---|
1508 | GByte. Using netcdf_64bit = .TRUE. allows file sizes larger than 2 |
---|
1509 | GByte.<br><br>The 64 bit offset format can be separately |
---|
1510 | switched off for those NetCDF files containing 3d volume date (<span style="font-family: Courier New,Courier,monospace;">DATA_3D_NETCDF</span>, |
---|
1511 | <span style="font-family: Courier New,Courier,monospace;">DATA_3D_AV_NETCDF</span>) |
---|
1512 | using <a href="#netcdf_64bit_3d">netcdf_64bit_3d</a>.<br><br><span style="font-weight: bold;">Warning:</span><br>Some |
---|
1513 | (PD or commercial) software may not support the 64 bit offset format.<br> |
---|
1514 | </td> </tr><tr><td style="vertical-align: top;"><a name="netcdf_64bit_3d"></a><span style="font-weight: bold;">netcdf_64bit_3d</span></td><td style="vertical-align: top;">L</td><td style="vertical-align: top;">.T.</td><td style="vertical-align: top;">NetCDF files containing 3d |
---|
1515 | volume data will have 64 bit offset format.<br><br>This |
---|
1516 | switch only comes into effect if <a href="#netcdf_64bit">netcdf_64bit</a> |
---|
1517 | = .TRUE.. It allows to switch off separately the 64 bit offset format |
---|
1518 | for those NetCDF files containing 3d volume data (<span style="font-family: Courier New,Courier,monospace;">DATA_3D_NETCDF</span>, |
---|
1519 | <span style="font-family: Courier New,Courier,monospace;">DATA_3D_AV_NETCDF</span>).</td></tr><tr> |
---|
1520 | <td style="vertical-align: top;"> <p><a name="ngsrb"></a><b>ngsrb</b></p> </td> |
---|
1521 | <td style="vertical-align: top;">I</td> <td style="vertical-align: top;"><i>2</i></td> |
---|
1522 | <td style="vertical-align: top;">Grid |
---|
1523 | level at which data shall be gathered on PE0.<br> <br> |
---|
1524 | In case of a run using several PEs and the multigrid method for solving |
---|
1525 | the Poisson equation for perturbation pressure (see <a href="#psolver">psolver</a>), |
---|
1526 | the value of this parameter defines on which grid level the data are |
---|
1527 | gathered on PE0 in order to allow for a further coarsening of the grid. |
---|
1528 | The finest grid defines the largest grid level. By default, the |
---|
1529 | gathering level is determined automatically and displayed in file <a href="chapter_3.4.html#RUN_CONTROL">RUN_CONTROL</a>. |
---|
1530 | It is only possible to gather data from a level larger than the one |
---|
1531 | determined automatically. A test run may be neccessary to determine |
---|
1532 | this level.</td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="normalizing_region"></a><b>normalizing_region</b></p> |
---|
1533 | </td> <td style="vertical-align: top;">I</td> |
---|
1534 | <td style="vertical-align: top;"><span style="font-style: italic;">0</span><br> </td> |
---|
1535 | <td style="vertical-align: top;"> <p>Determines the |
---|
1536 | subdomain from which the normalization |
---|
1537 | quantities are calculated. </p> <p>If output |
---|
1538 | data of the horizontally averaged vertical profiles |
---|
1539 | (see <a href="#data_output_pr">data_output_pr</a>) |
---|
1540 | is to be normalized (see <a href="#cross_normalized_x">cross_normalized_x</a>, |
---|
1541 | <a href="#cross_normalized_y">cross_normalized_y</a>), |
---|
1542 | the respective normalization quantities are by default calculated from |
---|
1543 | the averaged data of the total model domain (<b>normalizing_region</b> |
---|
1544 | = <i>0</i>) and are thus representative for the total |
---|
1545 | domain. Instead |
---|
1546 | of that, normalization quantities can also be calculated for a |
---|
1547 | subdomain. The wanted subdomain can be given with the parameter <span style="font-weight: bold;">normalizing_region</span>, |
---|
1548 | where <i>1</i> |
---|
1549 | <= <b>normalizing_region</b> <= <i>9 </i>must |
---|
1550 | hold. These |
---|
1551 | quantities are then used for normalizing of all profiles (even for that |
---|
1552 | of the total domain).</p> </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="nsor"></a><b>nsor</b></p> |
---|
1553 | </td> <td style="vertical-align: top;">I</td> |
---|
1554 | <td style="vertical-align: top;"><i>20</i></td> |
---|
1555 | <td style="vertical-align: top;"> <p>Number of |
---|
1556 | iterations to be used with the SOR-scheme. </p> <p>This |
---|
1557 | parameter is only effective if the SOR-scheme is selected |
---|
1558 | as pressure solver (<a href="#psolver">psolver</a> |
---|
1559 | = <span style="font-style: italic;">'sor'</span>). |
---|
1560 | The number of |
---|
1561 | iterations necessary for a sufficient convergence of the scheme depends |
---|
1562 | on the grid point numbers and is to be determined by appropriate test |
---|
1563 | runs (the default value will not at all be sufficient for larger grid |
---|
1564 | point numbers). The number of iterations used for the first call of the |
---|
1565 | SOR-scheme (t = 0) is determined via the parameter <a href="chapter_4.1.html#nsor_ini">nsor_ini</a>.</p> |
---|
1566 | </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="nz_do3d"></a><b>nz_do3d</b></p> |
---|
1567 | </td> <td style="vertical-align: top;">I</td> |
---|
1568 | <td style="vertical-align: top;"><i>nz+1</i></td> |
---|
1569 | <td style="vertical-align: top;"> Limits the output of 3d |
---|
1570 | volume data along the vertical direction (grid point index k).<br><br>By |
---|
1571 | default, data for all grid points along z are output. The parameter <span style="font-weight: bold;">nz_do3d</span> |
---|
1572 | can be used to limit the output up to a certain vertical grid point |
---|
1573 | (e.g. in order to reduce the amount of output data). It affects all |
---|
1574 | output of volume data ("normal" output to file, see <a href="#data_output">data_output</a>, as well as output |
---|
1575 | for <span style="font-weight: bold;">dvrp</span>-software, |
---|
1576 | see <a href="#mode_dvrp">mode_dvrp</a>).</td> |
---|
1577 | </tr> <tr> <td style="vertical-align: top;"><p><a name="omega_sor"></a><b>omega_sor</b></p> |
---|
1578 | </td> <td style="vertical-align: top;">R</td> |
---|
1579 | <td style="vertical-align: top;"><i>1.8</i></td> |
---|
1580 | <td style="vertical-align: top;"> <p>Convergence |
---|
1581 | factor to be used with the the SOR-scheme. </p> <p>If |
---|
1582 | the SOR-scheme is selected (<a href="#psolver">psolver</a> |
---|
1583 | = <span style="font-style: italic;">'sor'</span>), |
---|
1584 | this parameter |
---|
1585 | determines the value of the convergence factor, where <i>1.0</i> |
---|
1586 | <= <b>omega_sor</b> < <i>2.0 </i>. |
---|
1587 | The optimum value of <b>omega_sor</b> |
---|
1588 | depends on the number of grid points along the different directions in |
---|
1589 | space. For non-equidistant grids it can only be determined by |
---|
1590 | appropriate test runs.</p> </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="prandtl_number"></a><b>prandtl_number</b></p> |
---|
1591 | </td> <td style="vertical-align: top;">R</td> |
---|
1592 | <td style="vertical-align: top;"><i>1.0</i></td> |
---|
1593 | <td style="vertical-align: top;"> <p>Ratio of the |
---|
1594 | eddy diffusivities for momentum and heat (K<sub>m</sub>/K<sub>h</sub>). |
---|
1595 | </p> <p>For runs with constant eddy diffusivity (see <a href="chapter_4.1.html#km_constant">km_constant</a>), |
---|
1596 | this parameter can be used to assign the Prandtl number (ratio K<sub>m</sub> |
---|
1597 | / K<sub>h</sub>).</p> </td> </tr> <tr> |
---|
1598 | <td style="vertical-align: top;"> <p><a name="profile_columns"></a><b>profile_columns</b></p> |
---|
1599 | </td> <td style="vertical-align: top;">I</td> |
---|
1600 | <td style="vertical-align: top;"><i>3</i></td> |
---|
1601 | <td style="vertical-align: top;"> <p>Number of |
---|
1602 | coordinate systems to be plotted<span style="font-weight: bold;"></span> |
---|
1603 | in one row by <span style="font-weight: bold;">profil</span>. |
---|
1604 | </p> <p>This parameter only applies for <a href="chapter_4.2.html#data_output_format">data_output_format</a> |
---|
1605 | = <span style="font-style: italic;">'profil'</span>.</p><p>It |
---|
1606 | determines the layout of plots of |
---|
1607 | horizontally averaged profiles (<a href="#data_output_pr">data_output_pr</a>) |
---|
1608 | when plotted with the plot software <span style="font-weight: bold;">profil</span>. |
---|
1609 | Generally, the number and sequence of coordinate systems (panels) to be |
---|
1610 | plotted on one page are |
---|
1611 | determined by <a href="#cross_profiles">cross_profiles</a>. |
---|
1612 | <b>profile_columns</b> |
---|
1613 | determines how many panels are to be |
---|
1614 | arranged next to each other in one row (number of columns). The |
---|
1615 | respective number of rows on a page is assigned by <a href="#profile_rows">profile_rows</a>. |
---|
1616 | According to their order given by <a href="#data_output_pr">data_output_pr</a>, |
---|
1617 | the panels are arranged beginning in the top row from left to right and |
---|
1618 | then continued in the following row. If the number of panels cranz |
---|
1619 | > <b>profile_columns</b> * <b>profile_rows</b>, |
---|
1620 | the remaining |
---|
1621 | panels are drawn on an additional page. If cranz < <b>profile_columns</b>, |
---|
1622 | then <b>profile_columns</b> = cranz is automatically set. |
---|
1623 | If |
---|
1624 | row contains any panel, then the value of <b>profile_rows</b> |
---|
1625 | is reduced automatically.</p> </td> </tr> <tr> |
---|
1626 | <td style="vertical-align: top;"> <p><a name="profile_rows"></a><b>profile_rows</b></p> |
---|
1627 | </td> <td style="vertical-align: top;">I</td> |
---|
1628 | <td style="vertical-align: top;"><i>2</i></td> |
---|
1629 | <td style="vertical-align: top;"> <p>Number of rows |
---|
1630 | of coordinate systems to be plotted on one page |
---|
1631 | by <span style="font-weight: bold;">profil</span>. |
---|
1632 | </p> <p>This parameter only applies for <a href="chapter_4.2.html#data_output_format">data_output_format</a> |
---|
1633 | = <span style="font-style: italic;">'profil'</span>.</p><p>It |
---|
1634 | determines the layout of plots of horizontally averaged |
---|
1635 | profiles. See <a href="#profile_columns">profile_columns</a>.</p> |
---|
1636 | </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="psolver"></a><b>psolver</b></p> |
---|
1637 | </td> <td style="vertical-align: top;">C * 10</td> |
---|
1638 | <td style="vertical-align: top;"><i>'poisfft'</i></td> |
---|
1639 | <td style="vertical-align: top;"> <p>Scheme to be |
---|
1640 | used to solve the Poisson equation for the |
---|
1641 | perturbation pressure. </p> <br> |
---|
1642 | The user can choose between the following schemes:<br> <table style="text-align: left; width: 100%;" cellpadding="2" cellspacing="2"> <tbody> <tr> <td style="vertical-align: top;"><i>poisfft</i></td> |
---|
1643 | <td style="vertical-align: top;">Direct method using FFT |
---|
1644 | along x and y, solution of a |
---|
1645 | tridiagonal matrix along z, and backward |
---|
1646 | FFT (see Siano, institute reports, volume 54). The FFT routines to be |
---|
1647 | used can be determined via the initialization parameter <a href="chapter_4.1.html#fft_method">fft_method</a>.<br> |
---|
1648 | This solver is specially optimized for 1d domain decompositions. |
---|
1649 | Vectorization is optimized for domain decompositions along x only.</td> |
---|
1650 | </tr> <tr> <td style="vertical-align: top;"><p><i>poisfft_</i> |
---|
1651 | <br> <i>hybrid</i></p> |
---|
1652 | </td> <td style="vertical-align: top;">Direct |
---|
1653 | method using FFT |
---|
1654 | along x and y, solution of a |
---|
1655 | tridiagonal matrix along z, and backward |
---|
1656 | FFT (see Siano, institute reports, volume 54). The FFT routines to be |
---|
1657 | used can be determined via the initialization parameter <a href="chapter_4.1.html#fft_method">fft_method</a>.<br> |
---|
1658 | This solver is specially optimized for 1d domain decompositions. |
---|
1659 | Vectorization is optimized for domain decompositions along x only.</td> |
---|
1660 | </tr> <tr> <td style="vertical-align: top;"><i>multigrid</i></td> |
---|
1661 | <td style="vertical-align: top;"> <p>Multi-grid |
---|
1662 | scheme (see Uhlenbrock, diploma thesis). v- |
---|
1663 | and |
---|
1664 | w-cycles (see <a href="#cycle_mg">cycle_mg</a>) |
---|
1665 | are implemented. The convergence of the iterative scheme can be |
---|
1666 | steered by the number of v-/w-cycles to be carried out for each call of |
---|
1667 | the scheme (<a href="#mg_cycles">mg_cycles</a>) |
---|
1668 | and by the number of Gauss-Seidel iterations (see <a href="#ngsrb">ngsrb</a>) |
---|
1669 | to be carried out on each grid level. Instead the requested accuracy |
---|
1670 | can be given via <a href="#residual_limit">residual_limit</a>. |
---|
1671 | <span style="font-weight: bold;">This is the default!</span> |
---|
1672 | The |
---|
1673 | smaller this limit is, the more cycles have to be carried out in this |
---|
1674 | case and the number of cycles may vary from timestep to timestep.</p> |
---|
1675 | <br>If <a href="#mg_cycles">mg_cycles</a> |
---|
1676 | is set to its optimal value, the computing time of the |
---|
1677 | multi-grid scheme amounts approximately to that of the direct solver <span style="font-style: italic;">poisfft</span>, as long as |
---|
1678 | the number of |
---|
1679 | grid points in the three directions |
---|
1680 | of space corresponds to a power-of-two (2<sup>n</sup>) |
---|
1681 | where <i>n</i> >= 5 must hold. With large <i>n, |
---|
1682 | </i>the |
---|
1683 | multi-grid scheme can even be faster than the direct solver (although |
---|
1684 | its accuracy is several orders of magnitude worse, but this does not |
---|
1685 | affect the accuracy of the simulation). Nevertheless, the user should |
---|
1686 | always carry out some test runs in order to find out the optimum value |
---|
1687 | for <a href="#mg_cycles">mg_cycles</a>, |
---|
1688 | because the CPU time of a run very critically depends on this |
---|
1689 | parameter. <p>This scheme requires that the number of grid |
---|
1690 | points of |
---|
1691 | the |
---|
1692 | subdomains (or of the total domain, if only one PE is uesd) along each |
---|
1693 | of the directions can at least be devided once by 2 without rest.</p> |
---|
1694 | With parallel runs, starting from a certain grid level the |
---|
1695 | data of the subdomains are possibly gathered on PE0 in order to allow |
---|
1696 | for a further coarsening of the grid. The grid level for gathering can |
---|
1697 | be manually set by <a href="#mg_switch_to_pe0_level">mg_switch_to_pe0_level</a>.<br> |
---|
1698 | <p>Using this procedure requires the subdomains to be of |
---|
1699 | identical size (see <a href="chapter_4.1.html#grid_matching">grid_matching</a>).</p> |
---|
1700 | </td> </tr> <tr> <td style="vertical-align: top;"><i>sor</i></td> |
---|
1701 | <td style="vertical-align: top;">Successive over |
---|
1702 | relaxation |
---|
1703 | method (SOR). The convergence of |
---|
1704 | this |
---|
1705 | iterative scheme is steered with the parameters <a href="#omega_sor">omega_sor</a>, |
---|
1706 | <a href="chapter_4.1.html#nsor_ini">nsor_ini</a> |
---|
1707 | and <a href="chapter_4.1.html#nsor">nsor</a>. |
---|
1708 | <br>Compared to the direct method and the multi-grid method, this |
---|
1709 | scheme |
---|
1710 | needs substantially |
---|
1711 | more computing time. It should only be used for test runs, e.g. if |
---|
1712 | errors in the other pressure solver methods are assumed.</td> </tr> |
---|
1713 | </tbody> </table> <br>In order to speed-up |
---|
1714 | performance, the Poisson equation is by default |
---|
1715 | only solved at the last substep of a multistep Runge-Kutta scheme (see <a href="#call_psolver_at_all_substeps">call_psolver |
---|
1716 | at_all_substeps</a> and <a href="chapter_4.1.html#timestep_scheme">timestep_scheme</a>). |
---|
1717 | </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="rayleigh_damping_factor"></a><b>rayleigh_damping</b> |
---|
1718 | <br> <b>_factor</b></p> </td> <td style="vertical-align: top;">R</td> <td style="vertical-align: top;"><i>0.0 or</i><br> |
---|
1719 | <i>0.01</i></td> <td style="vertical-align: top;"> |
---|
1720 | <p>Factor for Rayleigh damping. </p> <p>A |
---|
1721 | so-called Rayleigh damping is applied to all prognostic |
---|
1722 | variables if a non-zero value is assigned to <b>rayleigh_damping_factor</b>. |
---|
1723 | If switched on, variables are forced towards the value of their |
---|
1724 | respective basic states (e.g. the geostrophic wind). The intensity of |
---|
1725 | damping is controlled by the value the <b>rayleigh_damping_factor</b> |
---|
1726 | is assigned to. |
---|
1727 | The damping starts weakly at a height defined by <a href="#rayleigh_damping_height">rayleigh_damping_height</a> |
---|
1728 | and rises according to a sin<sup>2</sup>-function to its |
---|
1729 | maximum value |
---|
1730 | at |
---|
1731 | the top boundary. </p> <p>This method |
---|
1732 | effectively damps gravity waves, caused by boundary layer convection, |
---|
1733 | which may spread out vertically in the inversion layer and which are |
---|
1734 | reflected at the top |
---|
1735 | boundary. This particularly happens with the upstream-spline scheme |
---|
1736 | switched on (see <a href="chapter_4.1.html#momentum_advec">momentum_advec</a> |
---|
1737 | or <a href="chapter_4.1.html#scalar_advec">scalar_advec</a>). |
---|
1738 | Therefore, with this scheme the Rayleigh damping is switched on (<b>rayleigh_damping_factor</b> |
---|
1739 | = <i>0.01</i>) by default. Otherwise it remains switched |
---|
1740 | off. </p> <p>The Rayleigh damping factor must |
---|
1741 | hold the condition <i>0.0</i> |
---|
1742 | <= <b>rayleigh_damping_factor</b> |
---|
1743 | <= <i>1.0</i>. Large values (close to <span style="font-style: italic;">1.0</span>) can cause |
---|
1744 | numerical instabilities.</p> </td> </tr> <tr> |
---|
1745 | <td style="vertical-align: top;"> <p><a name="rayleigh_damping_height"></a><b>rayleigh_damping</b> |
---|
1746 | <br> <b>_height</b></p> </td> <td style="vertical-align: top;">R</td> <td style="vertical-align: top;"> <p><i>2/3 *</i> |
---|
1747 | <br><span style="font-style: italic;">zu</span><i style="font-style: italic;">(nz)</i></p> |
---|
1748 | </td> <td style="vertical-align: top;"> <p>Height |
---|
1749 | where the Rayleigh damping starts (in m). </p> <p>With |
---|
1750 | Rayleigh damping switched on (see <a href="#rayleigh_damping_factor">rayleigh_damping_factor</a>), |
---|
1751 | this parameter determines the range where damping is applied. By |
---|
1752 | default, Rayleigh damping will be applied in the upper third of the |
---|
1753 | model |
---|
1754 | domain.</p> </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="residual_limit"></a><b>residual_limit</b></p> |
---|
1755 | </td> <td style="vertical-align: top;">R</td> |
---|
1756 | <td style="vertical-align: top;"><i>1.0E-6</i></td> |
---|
1757 | <td style="vertical-align: top;"> <p>Largest |
---|
1758 | residual permitted for the multi-grid scheme (in s<sup>-2</sup>m<sup>-3</sup>). |
---|
1759 | </p> <p>This is a parameter to steer the accuracy of the |
---|
1760 | multi-grid |
---|
1761 | scheme (see <a href="#psolver">psolver</a>). |
---|
1762 | The assigned cycle (v- or w-cycle, see <a href="#mg_cycles">mg_cycles</a>) |
---|
1763 | is passed through until the residual falls below the limit given by <span style="font-weight: bold;">residual_limit</span>. If |
---|
1764 | this |
---|
1765 | is not the case after 1000 cycles, the PALM aborts with a corresponding |
---|
1766 | error message.</p> <p>The reciprocal value of this |
---|
1767 | parameter can be interpreted as |
---|
1768 | a factor by the divergence of the provisional |
---|
1769 | velocity field is approximately reduced after the multi-grid scheme has |
---|
1770 | been applied (thus the default value causes a reduction of the |
---|
1771 | divergence by approx. 6 orders of magnitude). </p> </td> |
---|
1772 | </tr> <tr> <td style="vertical-align: top;"><p><a name="restart_time"></a><b>restart_time</b></p> |
---|
1773 | </td> <td style="vertical-align: top;">R</td> |
---|
1774 | <td style="vertical-align: top;"><i>9999999.9</i></td> |
---|
1775 | <td style="vertical-align: top;"> <p>Simulated time |
---|
1776 | after which a restart run is to be carried out |
---|
1777 | (in s). </p> <p>The simulated time refers to the |
---|
1778 | beginning of the |
---|
1779 | initial run (t = 0), not to the beginning of the respective |
---|
1780 | restart run. Restart runs can additionally be forced to be carried out |
---|
1781 | in regular intervals using the run time parameter <a href="#dt_restart">dt_restart</a>. </p> <p><span style="font-weight: bold;">Note:</span><br> |
---|
1782 | A successful operation of this parameter requires additional |
---|
1783 | modifications in the <span style="font-weight: bold;">mrun</span>-call |
---|
1784 | for the respective run (see <a href="chapter_3.3.html">chapter |
---|
1785 | 3.3</a>).<br> </p> <p>The choice of <b>restart_time</b> |
---|
1786 | or <b>dt_restart</b> does |
---|
1787 | not override the automatic start of restart runs in case that the job |
---|
1788 | runs out of CPU time. <br> </p> </td> </tr> |
---|
1789 | <tr> <td style="vertical-align: top;"> <p><a name="section_xy"></a><b>section_xy</b></p> |
---|
1790 | </td> <td style="vertical-align: top;">I(100)<br> |
---|
1791 | </td> <td style="vertical-align: top;"><span style="font-style: italic;">no section</span><br> |
---|
1792 | </td> <td style="vertical-align: top;"> <p lang="en-GB"><font face="Thorndale">Position |
---|
1793 | of cross section(s) for output of 2d horizontal cross |
---|
1794 | sections (grid point index k). </font> </p> <p><span lang="en-GB"><font face="Thorndale">If output |
---|
1795 | of |
---|
1796 | horizontal cross sections is selected (see </font></span><a href="chapter_4.2.html#data_output"><span lang="en-GB"><font face="Thorndale">data_output</font></span></a><span lang="en-GB"><font face="Thorndale">), this |
---|
1797 | parameter can be used to |
---|
1798 | define the position(s) of the cross section(s). Up to 100 positions of |
---|
1799 | cross sections can be selected by assigning <b>section_xy</b> |
---|
1800 | the |
---|
1801 | corresponding vertical grid point index/indices k of the requested |
---|
1802 | cross section(s). The exact location (height level) of the cross |
---|
1803 | section depends on the variable for which the output is made: zu(k) for |
---|
1804 | scalars and horizontal velocities, zw(k) for the vertical velocity. |
---|
1805 | Information about the exact location of the cross section is contained |
---|
1806 | in the NetCDF output file (if the default NetCDF output is switched on; |
---|
1807 | see <a href="#data_output_format">data_output_format</a>).</font></span></p><p><span lang="en-GB"><font face="Thorndale">Assigning <span style="font-weight: bold;">section_xy</span> = <span style="font-style: italic;">-1</span> |
---|
1808 | creates the output of horizontal cross sections averaged along z. In |
---|
1809 | the |
---|
1810 | NetCDF output file these (averaged) cross sections are given the |
---|
1811 | z-coordinate <span style="font-style: italic;">-1.0</span>.</font></span></p><p><span lang="en-GB"></span><span lang="en-GB"><font face="Thorndale">Assignments to <b>section_xy</b> |
---|
1812 | does not effect the output of horizontal cross sections of variable u<sub>*</sub> |
---|
1813 | and theta<sub>*</sub> and the liquid water path lwp*. For |
---|
1814 | these quantities always only one cross |
---|
1815 | section (for z=zu(1)) is output.</font></span></p><span lang="en-GB"><font face="Thorndale">In case of <span style="font-weight: bold;">data_output_format</span> = |
---|
1816 | <span style="font-style: italic;">'iso2d'</span> and |
---|
1817 | if several cross sections are selected (e.g. <b>section_xy</b> |
---|
1818 | = <i>1</i>, <i>10</i>, <i>15</i>), |
---|
1819 | then the respective data are |
---|
1820 | successively written to file. The output order follows the order given |
---|
1821 | by <b>section_xy</b>. </font></span></td> |
---|
1822 | </tr> <tr> <td style="vertical-align: top;"><p><a name="section_xz"></a><b>section_xz</b></p> |
---|
1823 | </td> <td style="vertical-align: top;">I(100)<br> |
---|
1824 | </td> <td style="vertical-align: top;"><span style="font-style: italic;">no section</span></td> |
---|
1825 | <td style="vertical-align: top;"> <p lang="en-GB"><font face="Thorndale">Position of cross section(s) |
---|
1826 | for output of 2d (xz) vertical cross sections (grid point |
---|
1827 | index j). </font> </p> <span lang="en-GB"><font face="Thorndale">If output of |
---|
1828 | vertical xz cross sections is selected (see </font></span><a href="chapter_4.2.html#data_output"><span lang="en-GB"><font face="Thorndale">data_output</font></span></a><span lang="en-GB"><font face="Thorndale">), this |
---|
1829 | parameter can be used to |
---|
1830 | define the position(s) of the cross section(s). Up to 100 positions of |
---|
1831 | cross sections can be selected by assigning <b>section_xz</b> |
---|
1832 | the |
---|
1833 | corresponding horizontal grid point index/indices j of the requested |
---|
1834 | cross section(s). The exact position (in y-direction) of the cross |
---|
1835 | section is given by j*</font></span><a href="chapter_4.1.html#dy"><span lang="en-GB"><font face="Thorndale">dy</font></span></a> <span lang="en-GB"><font face="Thorndale">or (j-0.5)*</font></span><a href="chapter_4.1.html#dy"><span lang="en-GB"><font face="Thorndale">dy</font></span></a><span lang="en-GB"><font face="Thorndale">, depending |
---|
1836 | on which grid the output quantity is defined. However, in |
---|
1837 | the NetCDF output file </font></span><span lang="en-GB"><font face="Thorndale">(if the |
---|
1838 | default NetCDF output is switched on; see <a href="chapter_4.2.html#data_output_format">data_output_format</a>) |
---|
1839 | no distinction is made between the quantities and j*<span style="font-weight: bold;">dy</span> is used for all |
---|
1840 | positions.<br><br>Assigning <span style="font-weight: bold;">section_xz</span> = <span style="font-style: italic;">-1</span> |
---|
1841 | creates the output of vertical cross sections averaged along y. In the |
---|
1842 | NetCDF output file these (averaged) cross sections are given the |
---|
1843 | y-coordinate <span style="font-style: italic;">-1.0</span>.<br></font></span><span lang="en-GB"><font face="Thorndale"><br></font></span><span lang="en-GB"><font face="Thorndale">In case of <span style="font-weight: bold;">data_output_format</span> = |
---|
1844 | <span style="font-style: italic;">'iso2d'</span> and |
---|
1845 | </font></span><span lang="en-GB"><font face="Thorndale">if several cross sections are |
---|
1846 | selected (e.g. <b>section_xz</b> = <i>0</i>, <i>12</i>, |
---|
1847 | <i>27</i>), |
---|
1848 | then the respective data are successively written to file. The output |
---|
1849 | order follows the order given by <b>section_xz</b>.</font></span></td> |
---|
1850 | </tr> <tr> <td style="vertical-align: top;"><p><a name="section_yz"></a><b>section_yz</b></p> |
---|
1851 | </td> <td style="vertical-align: top;">I(100)<br> |
---|
1852 | </td> <td style="vertical-align: top;"><span style="font-style: italic;">no section</span></td> |
---|
1853 | <td style="vertical-align: top;"> <p lang="en-GB"><font face="Thorndale">Position of cross section(s) |
---|
1854 | for output of 2d (yz) vertical cross sections (grid point |
---|
1855 | index i). </font> </p> <span lang="en-GB"><font face="Thorndale">If output of |
---|
1856 | vertical yz cross sections is selected (see </font></span><a href="chapter_4.2.html#data_output"><span lang="en-GB"><font face="Thorndale">data_output</font></span></a><span lang="en-GB"><font face="Thorndale">), |
---|
1857 | this parameter can be used to define the position(s) of the cross |
---|
1858 | section(s). Up to 100 positions of cross sections can be selected by |
---|
1859 | assigning <b>section_yz</b> the corresponding horizontal |
---|
1860 | grid point |
---|
1861 | index/indices i of the requested cross section(s). The exact position |
---|
1862 | (in x-direction) of the cross section is given by i*</font></span><a href="chapter_4.1.html#dx"><span lang="en-GB"><font face="Thorndale">dx</font></span></a> <span lang="en-GB"><font face="Thorndale">or (i-0.5)*</font></span><a href="chapter_4.1.html#dx"><span lang="en-GB"><font face="Thorndale">dx</font></span></a><span lang="en-GB"><font face="Thorndale">, depending |
---|
1863 | on which grid the output quantity is defined. </font></span><span lang="en-GB"><font face="Thorndale">However, in |
---|
1864 | the NetCDF output file </font></span><span lang="en-GB"><font face="Thorndale">(if the |
---|
1865 | default NetCDF output is switched on; see <a href="chapter_4.2.html#data_output_format">data_output_format</a>) |
---|
1866 | no distinction is made between the quantities and i*<span style="font-weight: bold;">dx</span> is used for all |
---|
1867 | positions.<br><br></font></span><span lang="en-GB"><font face="Thorndale">Assigning <span style="font-weight: bold;">section_yz</span> = <span style="font-style: italic;">-1</span> |
---|
1868 | creates the output of vertical cross sections averaged along x. In the |
---|
1869 | NetCDF output file these (averaged) cross sections are given the |
---|
1870 | x-coordinate <span style="font-style: italic;">-1.0</span>.</font></span><br><span lang="en-GB"></span><span lang="en-GB"><font face="Thorndale"> <br></font></span><span lang="en-GB"><font face="Thorndale">In case of <span style="font-weight: bold;">data_output_format</span> = |
---|
1871 | <span style="font-style: italic;">'iso2d'</span> and |
---|
1872 | </font></span><span lang="en-GB"><font face="Thorndale">if several cross sections are |
---|
1873 | selected (e.g. <b>section_yz</b> = <span style="font-style: italic;">3</span>, <span style="font-style: italic;">27</span>, 19), then the |
---|
1874 | respective data are successively written to file. The output order |
---|
1875 | follows the order given by <b>section_yz</b>.</font></span></td> |
---|
1876 | </tr> <tr> <td style="vertical-align: top;"><a name="skip_time_data_output"></a><span style="font-weight: bold;">skip_time_data_output</span><br> |
---|
1877 | </td> <td style="vertical-align: top;">R<br> </td> |
---|
1878 | <td style="vertical-align: top;"><span style="font-style: italic;">0.0</span><br> </td> |
---|
1879 | <td style="vertical-align: top;">No data output before |
---|
1880 | this interval has passed (in s).<br><br>This |
---|
1881 | parameter causes that data output activities are starting not before |
---|
1882 | this interval |
---|
1883 | (counting from the beginning of the simulation, t=0) has passed. By |
---|
1884 | default, this |
---|
1885 | applies for output of instantaneous 3d volume data, cross section data, |
---|
1886 | spectra and vertical profile data as well as for temporally averaged 2d |
---|
1887 | and 3d data. Individual intervals can be assigned using parameters <a href="#skip_time_do3d">skip_time_do3d</a>, <a href="#skip_time_do2d_xy">skip_time_do2d_xy</a>, <a href="#skip_time_do2d_xz">skip_time_do2d_xz</a>, <a href="#skip_time_do2d_yz">skip_time_do2d_yz</a>, <a href="#skip_time_dosp">skip_time_dosp</a>, <a href="#skip_time_dopr">skip_time_dopr</a>, and <a href="#skip_time_data_output_av">skip_time_data_output_av</a>.<br><br><span style="font-weight: bold;">Example:</span><br>If |
---|
1888 | the user has set <a href="#dt_data_output">dt_data_output</a> |
---|
1889 | = <span style="font-style: italic;">3600.0</span> |
---|
1890 | and <span style="font-weight: bold;">skip_time_data_output</span> |
---|
1891 | = <span style="font-style: italic;">1800.0</span>, |
---|
1892 | then the first output will be done at t = 5400 s.<br> </td> |
---|
1893 | </tr> <tr><td style="vertical-align: top;"><a name="skip_time_data_output_av"></a><span style="font-weight: bold;">skip_time_data_output_av</span></td><td style="vertical-align: top;">R</td><td style="vertical-align: top;"><span style="font-style: italic;">value of <a href="#skip_time_data_output">skip_time_<br>data_output</a></span></td><td style="vertical-align: top;">No output of temporally |
---|
1894 | averaged 2d/3d data before this interval has passed (in s).<br><br>This |
---|
1895 | parameter causes that data output activities are starting not before |
---|
1896 | this interval |
---|
1897 | (counting from the beginning of the simulation, t=0) has passed. <br><br><span style="font-weight: bold;">Example:</span><br>If |
---|
1898 | the user has set <a href="#dt_data_output_av">dt_data_output_av</a> |
---|
1899 | = <span style="font-style: italic;">3600.0</span> |
---|
1900 | and <span style="font-weight: bold;">skip_time_data_output_av</span> |
---|
1901 | = <span style="font-style: italic;">1800.0</span>, |
---|
1902 | then the first output will be done at t = 5400 s.</td></tr><tr> |
---|
1903 | <td style="vertical-align: top;"><a name="skip_time_dopr"></a><span style="font-weight: bold;">skip_time_dopr</span><br> |
---|
1904 | </td> <td style="vertical-align: top;">R<br> </td> |
---|
1905 | <td style="vertical-align: top;"><span style="font-style: italic;">value of <a href="chapter_4.2.html#skip_time_data_output">skip_time_<br>data_output</a></span> |
---|
1906 | </td> <td style="vertical-align: top;">No output of |
---|
1907 | vertical profile data before this interval has passed (in s).<br><br>This |
---|
1908 | parameter causes that data output activities are starting not before |
---|
1909 | this interval |
---|
1910 | (counting from the beginning of the simulation, t=0) has passed. <br><br><span style="font-weight: bold;">Example:</span><br>If |
---|
1911 | the user has set <a href="#dt_dopr">dt_dopr</a> = <span style="font-style: italic;">3600.0</span> and <span style="font-weight: bold;">skip_time_dopr</span> = <span style="font-style: italic;">1800.0</span>, then the |
---|
1912 | first output will be done at t = 5400 s. </td> </tr> <tr> |
---|
1913 | <td style="vertical-align: top;"><a name="skip_time_do2d_xy"></a><span style="font-weight: bold;">skip_time_do2d_xy</span><br> |
---|
1914 | </td> <td style="vertical-align: top;">R<br> </td> |
---|
1915 | <td style="vertical-align: top;"><span style="font-style: italic;">value of <a href="chapter_4.2.html#skip_time_data_output">skip_time_<br>data_output</a></span> |
---|
1916 | </td> <td style="vertical-align: top;">No output of |
---|
1917 | instantaneous horizontal cross section data before this interval has |
---|
1918 | passed (in s).<br><br>This |
---|
1919 | parameter causes that data output activities are starting not before |
---|
1920 | this interval |
---|
1921 | (counting from the beginning of the simulation, t=0) has passed. <br><br><span style="font-weight: bold;">Example:</span><br>If |
---|
1922 | the user has set <a href="#dt_do2d_xy">dt_do2d_xy</a> |
---|
1923 | = <span style="font-style: italic;">3600.0</span> |
---|
1924 | and <span style="font-weight: bold;">skip_time_do2d_xy</span> |
---|
1925 | = <span style="font-style: italic;">1800.0</span>, |
---|
1926 | then the first output will be done at t = 5400 s. </td> </tr> |
---|
1927 | <tr> <td style="vertical-align: top;"><a name="skip_time_do2d_xz"></a><span style="font-weight: bold;">skip_time_do2d_xz</span><br> |
---|
1928 | </td> <td style="vertical-align: top;">R<br> </td> |
---|
1929 | <td style="vertical-align: top;"><span style="font-style: italic;">value of <a href="chapter_4.2.html#skip_time_data_output">skip_time_<br>data_output</a></span> |
---|
1930 | </td> <td style="vertical-align: top;">No output of |
---|
1931 | instantaneous vertical (xz) cross section data before this interval has |
---|
1932 | passed (in s).<br><br>This |
---|
1933 | parameter causes that data output activities are starting not before |
---|
1934 | this interval |
---|
1935 | (counting from the beginning of the simulation, t=0) has passed. <br><br><span style="font-weight: bold;">Example:</span><br>If |
---|
1936 | the user has set <a href="#dt_do2d_xz">dt_do2d_xz</a> |
---|
1937 | = <span style="font-style: italic;">3600.0</span> |
---|
1938 | and <span style="font-weight: bold;">skip_time_do2d_xz</span> |
---|
1939 | = <span style="font-style: italic;">1800.0</span>, |
---|
1940 | then the first output will be done at t = 5400 s. </td> </tr> |
---|
1941 | <tr> <td style="vertical-align: top;"><a name="skip_time_do2d_yz"></a><span style="font-weight: bold;">skip_time_do2d_yz</span><br> |
---|
1942 | </td> <td style="vertical-align: top;">R<br> </td> |
---|
1943 | <td style="vertical-align: top;"><span style="font-style: italic;">value of <a href="chapter_4.2.html#skip_time_data_output">skip_time_<br>data_output</a></span> |
---|
1944 | </td> <td style="vertical-align: top;">No output of |
---|
1945 | instantaneous vertical (yz) cross section data before this interval has |
---|
1946 | passed (in s).<br><br>This |
---|
1947 | parameter causes that data output activities are starting not before |
---|
1948 | this interval |
---|
1949 | (counting from the beginning of the simulation, t=0) has passed. <br><br><span style="font-weight: bold;">Example:</span><br>If |
---|
1950 | the user has set <a href="#dt_do2d_yz">dt_do2d_yz</a> |
---|
1951 | = <span style="font-style: italic;">3600.0</span> |
---|
1952 | and <span style="font-weight: bold;">skip_time_do2d_yz</span> |
---|
1953 | = <span style="font-style: italic;">1800.0</span>, |
---|
1954 | then the first output will be done at t = 5400 s. </td> </tr> |
---|
1955 | <tr> <td style="vertical-align: top;"><a name="skip_time_do3d"></a><span style="font-weight: bold;">skip_time_do3d</span><br> |
---|
1956 | </td> <td style="vertical-align: top;">R<br> </td> |
---|
1957 | <td style="vertical-align: top;"><span style="font-style: italic;">value of <a href="chapter_4.2.html#skip_time_data_output">skip_time_<br>data_output</a></span> |
---|
1958 | </td> <td style="vertical-align: top;">No output of |
---|
1959 | instantaneous 3d volume data before this interval has passed (in s).<br><br>This |
---|
1960 | parameter causes that data output activities are starting not before |
---|
1961 | this interval |
---|
1962 | (counting from the beginning of the simulation, t=0) has passed. <br><br><span style="font-weight: bold;">Example:</span><br>If |
---|
1963 | the user has set <a href="#dt_do3d">dt_do3d</a> = <span style="font-style: italic;">3600.0</span> and <span style="font-weight: bold;">skip_time_do3d</span> = <span style="font-style: italic;">1800.0</span>, then the |
---|
1964 | first output will be done at t = 5400 s. </td> </tr> |
---|
1965 | <tr> <td style="vertical-align: top;"> <p><a name="termination_time_needed"></a><b>termination_time</b> |
---|
1966 | <br> <b>_needed</b></p> </td> <td style="vertical-align: top;">R<br> </td> <td style="vertical-align: top;"><span style="font-style: italic;">35.0</span><br> </td> |
---|
1967 | <td style="vertical-align: top;"> <p>CPU time |
---|
1968 | needed for terminal actions at the end of a run in |
---|
1969 | batch mode (in s).<br> </p> <p>If the environment |
---|
1970 | variable <b>write_binary </b>is |
---|
1971 | set <i>true</i> (see <a href="chapter_3.3.html">chapter |
---|
1972 | 3.3</a>), PALM checks the remaining CPU time of the job after |
---|
1973 | each |
---|
1974 | timestep. Time integration must not consume the CPU time completely, |
---|
1975 | since several actions still have to be carried out after time |
---|
1976 | integration has finished (e.g. writing of binary data for the restart |
---|
1977 | run, carrying out output commands, copying of local files to their |
---|
1978 | permanent destinations, etc.) which also takes some time. The maximum |
---|
1979 | possible time needed for these activities plus a reserve is to be given |
---|
1980 | with the parameter <b>termination_time_needed</b>. Among |
---|
1981 | other things, |
---|
1982 | it depends on |
---|
1983 | the number of grid points used. If its value is selected too small, |
---|
1984 | then the |
---|
1985 | respective job will be prematurely aborted by the queuing system, which |
---|
1986 | may result in a data loss and will possibly interrupt the job chain.<br> |
---|
1987 | </p> <p>An abort happens in any way, if the environment |
---|
1988 | variable <span style="font-weight: bold;">write_binary</span> |
---|
1989 | is not set to <span style="font-style: italic;">true</span> |
---|
1990 | and if moreover the job has |
---|
1991 | been assigned an insufficient CPU time by <b>mrun</b> |
---|
1992 | option <tt><tt>-t</tt></tt>. <i><br> |
---|
1993 | </i> </p> <p><span style="font-weight: bold;">Note:</span><br> |
---|
1994 | On the IBM computers of the HLRN the time used by the job <span style="font-weight: bold;">before</span> the start of |
---|
1995 | PALM |
---|
1996 | have also to be accounted for (e.g. for |
---|
1997 | compilation and copying of input files).</p> </td> </tr> |
---|
1998 | <tr> <td style="vertical-align: top;"> <p><a name="use_prior_plot1d_parameters"></a><b>use_prior_plot1d</b> |
---|
1999 | <br> <b>_parameters</b></p> </td> <td style="vertical-align: top;">L</td> <td style="vertical-align: top;"><i>.F.</i></td> |
---|
2000 | <td style="vertical-align: top;"> <p>Additional |
---|
2001 | plot of vertical profile data with <span style="font-weight: bold;">profil</span> |
---|
2002 | from preceding runs of the |
---|
2003 | job chain. </p> <p>This parameter only applies |
---|
2004 | for <a href="chapter_4.2.html#data_output_format">data_output_format</a> |
---|
2005 | = <span style="font-style: italic;">'profil'</span>.</p><p>By |
---|
2006 | default, plots of horizontally averaged vertical profiles |
---|
2007 | (see <a href="#data_output_pr">data_output_pr</a>) |
---|
2008 | only contain profiles of data produced by the model |
---|
2009 | run. If profiles of prior times (i.e. data of preceding jobs of a |
---|
2010 | job chain) shall be plotted additionally (e.g. for comparison |
---|
2011 | purposes), <b>use_prior_plot1d_parameters</b> = <i>.T</i>. |
---|
2012 | must be |
---|
2013 | set.<br> </p> <p>For further explanation see <a href="chapter_4.5.2.html">chapter |
---|
2014 | 4.5.2</a>.</p> </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="z_max_do1d"></a><b>z_max_do1d</b></p> |
---|
2015 | </td> <td style="vertical-align: top;">R</td> |
---|
2016 | <td style="vertical-align: top;"><i>zu(nzt+1) (model |
---|
2017 | top)</i></td> <td style="vertical-align: top;"> |
---|
2018 | <p>Height level up to which horizontally averaged profiles are to |
---|
2019 | be |
---|
2020 | plotted with <span style="font-weight: bold;">profil</span> |
---|
2021 | (in |
---|
2022 | m). </p> <p>This parameter only applies for |
---|
2023 | <a href="chapter_4.2.html#data_output_format">data_output_format</a> |
---|
2024 | = <span style="font-style: italic;">'profil'</span>.</p><p>It |
---|
2025 | affects plots of horizontally averaged profiles |
---|
2026 | (<a href="#data_output_pr">data_output_pr</a>) |
---|
2027 | when plotted with the plot software <span style="font-weight: bold;">profil</span>. |
---|
2028 | By default, profiles are plotted up to the top boundary. The height |
---|
2029 | level up to which profiles are plotted can be decreased by assigning <span style="font-weight: bold;">z_max_do1d</span> a smaller |
---|
2030 | value. |
---|
2031 | Nevertheless, <span style="font-weight: bold;">all</span> |
---|
2032 | vertical |
---|
2033 | grid points (0 <= k <= nz+1) are still output to file <a href="chapter_3.4.html#PLOT1D_DATA">PLOT1D_DATA</a>.</p> |
---|
2034 | <p>If a normalization for the vertical axis was selected (see <a href="#cross_normalized_y">cross_normalized_y</a>), <b>z_max_do1d</b> |
---|
2035 | has no effect. Instead, <a href="#z_max_do1d_normalized">z_max_do1d_normalized</a> |
---|
2036 | must be used.</p> </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="z_max_do1d_normalized"></a><b>z_max_do1d</b> |
---|
2037 | <br> <b>_normalized</b></p> </td> <td style="vertical-align: top;">R</td> <td style="vertical-align: top;"><i>determined by plot</i> |
---|
2038 | <br> <i>data</i> <br> </td> <td style="vertical-align: top;"> <p>Normalized height |
---|
2039 | level up to which horizontally averaged |
---|
2040 | profiles are to be plotted with <span style="font-weight: bold;">profil</span>. |
---|
2041 | </p> <p>This parameter only applies for <a href="chapter_4.2.html#data_output_format">data_output_format</a> |
---|
2042 | = <span style="font-style: italic;">'profil'</span>.</p><p>It |
---|
2043 | affects plots of horizontally averaged profiles |
---|
2044 | (<a href="#data_output_pr">data_output_pr</a>) |
---|
2045 | when plotted with the plot software <span style="font-weight: bold;">profil</span>, |
---|
2046 | if a normalization for the vertical axis is selected |
---|
2047 | (see <a href="#cross_normalized_y">cross_normalized_y</a>). |
---|
2048 | If e.g. the boundary layer height is used for normalization, then <b>z_max_do1d_normalized</b> |
---|
2049 | = <i>1.5</i> means that all profiles up to the height |
---|
2050 | level of z = |
---|
2051 | 1.5* z<sub>i </sub>are plotted.</p> </td> </tr> |
---|
2052 | <tr> <td style="vertical-align: top;"> <p><a name="z_max_do2d"></a><b>z_max_do2d</b></p> |
---|
2053 | </td> <td style="vertical-align: top;">R<br> </td> |
---|
2054 | <td style="vertical-align: top;"><span style="font-style: italic;">zu(nz)</span><br> </td> |
---|
2055 | <td style="vertical-align: top;"> <p>Height level |
---|
2056 | up to which 2d cross sections are to be plotted |
---|
2057 | with <span style="font-weight: bold;">iso2d</span> |
---|
2058 | (in m). </p> <p>This parameter only applies for |
---|
2059 | <a href="#data_output_format">data_output_format</a> |
---|
2060 | = <span style="font-style: italic;">'iso2d'</span>.</p><p>It |
---|
2061 | affects plots of 2d vertical cross |
---|
2062 | sections (<a href="#data_output">data_output</a>) |
---|
2063 | when plotted with <span style="font-weight: bold;">iso2d</span>. |
---|
2064 | By |
---|
2065 | default, vertical sections are plotted up to the top boundary. <span style="font-weight: bold;"></span>In contrast, with <b>z_max_do2d |
---|
2066 | </b>the |
---|
2067 | visualization within |
---|
2068 | the plot can be limited to a certain height level (0 <= z |
---|
2069 | <= <b>z_max_do2d</b>). |
---|
2070 | Nevertheless, <span style="font-weight: bold;">all</span> |
---|
2071 | grid points |
---|
2072 | of the complete cross section are still output to the local files <a href="chapter_3.4.html#PLOT2D_XZ">PLOT2D_XZ</a> |
---|
2073 | or <a href="chapter_3.4.html#PLOT2D_YZ">PLOT2D_YZ</a>. |
---|
2074 | The level up to which the section is visualized can later be changed by |
---|
2075 | manually editing the |
---|
2076 | file <a href="chapter_3.4.html#PLOT2D_XZ_GLOBAL">PLOT2D_XZ_GLOBAL</a> |
---|
2077 | or <a href="chapter_3.4.html#PLOT2D_YZ_GLOBAL">PLOT2D_YZ_GLOBAL</a> |
---|
2078 | (the respective <span style="font-weight: bold;">iso2d</span>-parameter |
---|
2079 | is <a href="http://www.muk.uni-hannover.de/institut/software/iso2d_beschreibung.html#YRIGHT">yright</a>).</p> |
---|
2080 | </td> </tr> </tbody></table><br> |
---|
2081 | <br><h3 style="line-height: 100%;"><a name="Paketparameter"></a>Package |
---|
2082 | parameters: </h3> |
---|
2083 | Package (<span style="font-weight: bold;">mrun</span> |
---|
2084 | option -p): <span style="font-weight: bold;"><a name="particles_package"></a>particles</span> |
---|
2085 | NAMELIST group name: <span style="font-weight: bold;">particles_par<br> |
---|
2086 | </span><table style="text-align: left; width: 100%;" border="1" cellpadding="2" cellspacing="2"> <tbody><tr> |
---|
2087 | <td style="vertical-align: top;"><font size="4"><b>Parameter |
---|
2088 | name</b></font></td> |
---|
2089 | <td style="vertical-align: top;"><font size="4"><b>Type</b></font></td> |
---|
2090 | <td style="vertical-align: top;"> <p><b><font size="4">Default</font></b> <br> <b><font size="4">value</font></b></p> </td> |
---|
2091 | <td style="vertical-align: top;"> <p><font size="4"><b>Explanation</b></font></p> |
---|
2092 | </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="dt_prel"></a><b>dt_prel</b></p> |
---|
2093 | </td> <td style="vertical-align: top;">R</td> |
---|
2094 | <td style="vertical-align: top;"><i>9999999.9</i></td> |
---|
2095 | <td style="vertical-align: top;"> <p><font face="Thorndale, serif"><span lang="en-GB">Temporal |
---|
2096 | interval at |
---|
2097 | which particles are to be released <span lang="en-GB">from |
---|
2098 | a particle |
---|
2099 | source </span>(</span></font>in <font face="Thorndale, serif"><span lang="en-GB">s).</span> |
---|
2100 | </font> </p> <p><span lang="en-GB"><font face="Thorndale, serif">By default |
---|
2101 | particles are released only at the beginning of a simulation |
---|
2102 | (t_init=0). The time of the first release (t_init) can be changed with |
---|
2103 | package parameter </font></span><span lang="en-GB"></span><font><a href="#particle_advection_start"><font face="Thorndale, serif">particle_advection_start</font></a>. |
---|
2104 | </font><span lang="en-GB"><font face="Thorndale, serif">The time of the last release can be |
---|
2105 | set with the package parameter <a href="#end_time_prel">end_time_prel</a>. |
---|
2106 | If <span style="font-weight: bold;">dt_prel</span> |
---|
2107 | has been set, additional |
---|
2108 | releases will be at t = t_init+<span style="font-weight: bold;">dt_prel</span>, |
---|
2109 | t_init+2*<span style="font-weight: bold;">dt_prel</span>, |
---|
2110 | t_init+3*<span style="font-weight: bold;">dt_prel</span>, |
---|
2111 | etc.. Actual release times |
---|
2112 | may slightly deviate from thesel values (</font></span><span lang="en-GB"><font face="Thorndale, serif">see |
---|
2113 | e.g. </font></span><a href="#dt_dopr"><span lang="en-GB"><font face="Thorndale, serif">dt_dopr</font></span></a><span lang="en-GB"><font face="Thorndale, serif">).</font></span></p> |
---|
2114 | <p><span lang="en-GB"><font face="Thorndale, serif"> The domain |
---|
2115 | of the particle <span lang="en-GB"><font face="Thorndale, serif">source </font></span>as |
---|
2116 | well as the distance of released particles |
---|
2117 | within this source </font></span><span lang="en-GB"><font face="Thorndale, serif">are determined via package |
---|
2118 | parameters </font></span><a href="#pst"><span lang="en-GB"><font face="Thorndale, serif">pst</font></span></a><span lang="en-GB"><font face="Thorndale, serif">, </font></span><a href="#psl"><span lang="en-GB"><font face="Thorndale, serif">psl</font></span></a><span lang="en-GB"><font face="Thorndale, serif">, </font></span><a href="#psr"><span lang="en-GB"><font face="Thorndale, serif">psr</font></span></a><span lang="en-GB"><font face="Thorndale, serif">, </font></span><a href="#pss"><span lang="en-GB"><font face="Thorndale, serif">pss</font></span></a><span lang="en-GB"><font face="Thorndale, serif">, </font></span><a href="#psn"><span lang="en-GB"><font face="Thorndale, serif">psn</font></span></a><span lang="en-GB"><font face="Thorndale, serif">, </font></span><a href="#psb"><span lang="en-GB"><font face="Thorndale, serif">psb</font></span></a><span lang="en-GB"><font face="Thorndale, serif">, </font></span><a href="#pdx"><span lang="en-GB"><font face="Thorndale, serif">pdx</font></span></a><span lang="en-GB"><font face="Thorndale, serif">, </font></span><a href="#pdy"><span lang="en-GB"><font face="Thorndale, serif">pdy</font></span></a> |
---|
2119 | <span lang="en-GB"><font face="Thorndale, serif">and |
---|
2120 | </font></span><a href="#pdz"><span lang="en-GB"><font face="Thorndale, serif">pdz</font></span></a><span lang="en-GB"><font face="Thorndale, serif">.</font></span><span lang="en-GB"><font face="Thorndale, serif"> By |
---|
2121 | default, one particle is released at all points defined by these |
---|
2122 | parameters. The package parameter <a href="#particles_per_point">particles_per_point</a> |
---|
2123 | can be used to start more than one particle per point.<br> |
---|
2124 | </font></span></p> <p><span lang="en-GB"><font face="Thorndale, serif">Up to 10 |
---|
2125 | different groups of particles can be released at the same time (see </font></span><a href="chapter_4.2.html#number_of_particle_groups"><span lang="en-GB"><font face="Thorndale, serif">number_of_particle_groups</font></span></a><span lang="en-GB"><font face="Thorndale, serif">) |
---|
2126 | where each group may have a different source. All particles belonging |
---|
2127 | to one group have the same density ratio and the same radius. All other |
---|
2128 | particle features (e.g. location of the source) are |
---|
2129 | identical for all groups of particles.</font></span></p>Subgrid |
---|
2130 | scale velocities can (optionally) be included for calculating the |
---|
2131 | particle advection, using the method of Weil et al. (2004, JAS, 61, |
---|
2132 | 2877-2887). This method is switched on by the package |
---|
2133 | parameter <a href="#use_sgs_for_particles">use_sgs_for_particles</a>. |
---|
2134 | This also forces the Euler/upstream method to be used for time |
---|
2135 | advancement of the TKE (see initialization parameter <a href="chapter_4.1.html#use_upstream_for_tke">use_upstream_for_tke</a>). |
---|
2136 | The minimum timestep during the sub-timesteps is controlled by package |
---|
2137 | parameter <a href="#dt_min_part">dt_min_part</a>. <p><span lang="en-GB"><font face="Thorndale, serif">By |
---|
2138 | default, particles are weightless and transported passively with the |
---|
2139 | resolved scale flow. Particles can be given a mass and thus an inertia |
---|
2140 | by assigning the |
---|
2141 | package parameter </font></span><a href="#density_ratio"><span lang="en-GB"><font face="Thorndale, serif">density_ratio</font></span></a><span lang="en-GB"><font face="Thorndale, serif"> a |
---|
2142 | non-zero value (it |
---|
2143 | defines the ratio of the density of the fluid and the density of the |
---|
2144 | particles). In these cases their </font></span><a href="#radius"><span lang="en-GB"><font face="Thorndale, serif">radius</font></span></a><span lang="en-GB"></span><span lang="en-GB"><font face="Thorndale, serif"> |
---|
2145 | must also be defined, which affects their flow resistance. </font></span><a href="#diameter"><span lang="en-GB"></span></a><span lang="en-GB"></span><span lang="en-GB"><font face="Thorndale, serif"> </font></span> </p> |
---|
2146 | <p><span lang="en-GB"><font face="Thorndale, serif">Boundary |
---|
2147 | conditions for the particle transport can be defined with package |
---|
2148 | parameters </font></span><a href="#bc_par_t"><span lang="en-GB"><font face="Thorndale, serif">bc_par_t</font></span></a><span lang="en-GB"><font face="Thorndale, serif">, </font></span><a href="#bc_par_lr"><span lang="en-GB"><font face="Thorndale, serif">bc_par_lr</font></span></a><span lang="en-GB"><font face="Thorndale, serif">, </font></span><a href="bc_par_ns"><span lang="en-GB"><font face="Thorndale, serif">bc_par_ns</font></span></a> |
---|
2149 | <span lang="en-GB"><font face="Thorndale, serif">and |
---|
2150 | </font></span><a href="#bc_par_b"><span lang="en-GB"><font face="Thorndale, serif">bc_par_b</font></span></a><span lang="en-GB"><font face="Thorndale, serif">.</font></span></p><span lang="en-GB"><font face="Thorndale, serif">Timeseries |
---|
2151 | of particle quantities in NetCDF format can be output to local file <a href="chapter_3.4.html#DATA_1D_PTS_NETCDF">DATA_1D_PTS_NETCDF</a> |
---|
2152 | by using package parameter <a href="#dt_dopts">dt_dopts</a>.<br></font></span><p>For |
---|
2153 | analysis, additional output of |
---|
2154 | particle |
---|
2155 | information in equidistant temporal intervals can be carried out using <a href="#dt_write_particle_data">dt_write_particle_data</a> |
---|
2156 | (file <a href="chapter_3.4.html#PARTICLE_DATA">PARTICLE_DATA</a>).<br> |
---|
2157 | </p> <p><span style="font-family: thorndale,serif;">Statistical |
---|
2158 | informations</span> (e.g. the total number of particles used, the |
---|
2159 | number of particles exchanged between the PEs, etc.) are output to the |
---|
2160 | local file <a href="chapter_3.4.html#PARTICLE_DATA">PARTICLE_INFOS</a>, |
---|
2161 | if switched on by the parameter <a href="#write_particle_statistics">write_particle_statistics</a>. |
---|
2162 | <br> </p> <p><span lang="en-GB"><font face="Thorndale, serif">If a job |
---|
2163 | chain is to be carried out, particle |
---|
2164 | informations </font></span><span lang="en-GB"><font face="Thorndale, serif">for the restart run (e.g. current |
---|
2165 | location of |
---|
2166 | all |
---|
2167 | particles at the end of the |
---|
2168 | run) is output to |
---|
2169 | the local file</font></span> <font><a href="chapter_4.2.html#dt_dvrp"><span lang="en-GB"></span></a></font><a href="chapter_3.4.html#PARTICLE_RESTART_DATA_OUT">PARTICLE_RESTART_DATA_OUT</a><font><a href="chapter_4.2.html#dt_dvrp"><span lang="en-GB"></span></a></font>, |
---|
2170 | <span lang="en-GB"><font face="Thorndale, serif">which |
---|
2171 | must be saved at the |
---|
2172 | end of the run <tt><span lang="en-GB"></span></tt>and |
---|
2173 | given as an |
---|
2174 | input file to the restart run |
---|
2175 | under local file name</font></span> <a href="chapter_3.4.html#PARTICLE_RESTART_DATA_IN">PARTICLE_RESTART_DATA_IN</a> |
---|
2176 | u<span lang="en-GB"><font face="Thorndale, serif">sing |
---|
2177 | respective file |
---|
2178 | connection statements in the <span style="font-weight: bold;">mrun</span> |
---|
2179 | configuration file. </font></span> <span lang="en-GB"></span></p><p><span lang="en-GB"></span><span lang="en-GB"><font face="Thorndale, serif">The output of |
---|
2180 | particles for visualization with the graphic software <span style="font-weight: bold;">dvrp</span> is steered by |
---|
2181 | the package |
---|
2182 | parameter </font></span><a href="chapter_4.2.html#dt_dvrp"><span lang="en-GB"><font face="Thorndale, serif">dt_dvrp</font></span></a><font face="Thorndale, serif"><span lang="en-GB">. |
---|
2183 | For visualization |
---|
2184 | purposes particles can be given a |
---|
2185 | diameter by the parameter <a href="chapter_4.2.html#dvrp_psize">dvrp_psize</a> |
---|
2186 | (this diameter only affects the visualization). All particles have the |
---|
2187 | same size. Alternatively, particles can be given an individual size and |
---|
2188 | a </span>color <span lang="en-GB">by modifying the |
---|
2189 | user-interface (subroutine</span></font> <span style="font-family: monospace;">user_init_particles</span>)<span lang="en-GB"><font face="Thorndale, serif">. |
---|
2190 | Particles can pull a |
---|
2191 | “tail” behind themselves to improve their |
---|
2192 | visualization. |
---|
2193 | This is steered via the parameter <a href="chapter_4.2.html#use_particle_tails">use_particle_tails</a>.</font></span><a href="chapter_4.2.html#maximum_number_of_tailpoints"><span lang="en-GB"></span></a></p> <span lang="en-GB"></span><p><b>So far, the |
---|
2194 | particle transport realized in PALM does only |
---|
2195 | work |
---|
2196 | duly in case of a constant vertical grid spacing!</b></p> </td> |
---|
2197 | </tr> <tr> <td style="vertical-align: top;"><p><a name="bc_par_b"></a><b>bc_par_b</b></p> |
---|
2198 | </td> <td style="vertical-align: top;">C*15</td> |
---|
2199 | <td style="vertical-align: top;"><i>´reflect´</i></td> |
---|
2200 | <td style="vertical-align: top;"> <p>Bottom |
---|
2201 | boundary condition for particle transport. </p> <p>By |
---|
2202 | default, particles are reflected at the bottom boundary. |
---|
2203 | Alternatively, a particle absorption can set by <b>bc_par_b</b> |
---|
2204 | = <i>´absorb´</i>.</p> </td> |
---|
2205 | </tr> <tr> <td style="vertical-align: top;"><p><a name="bc_par_lr"></a><b>bc_par_lr</b></p> |
---|
2206 | </td> <td style="vertical-align: top;">C*15</td> |
---|
2207 | <td style="vertical-align: top;"><i>´cyclic´</i></td> |
---|
2208 | <td style="vertical-align: top;"> <p>Lateral |
---|
2209 | boundary condition (x-direction) for particle |
---|
2210 | transport. </p> <p>By default, cyclic boundary conditions |
---|
2211 | are used along x. |
---|
2212 | Alternatively, reflection (<b>bc_par_lr</b> |
---|
2213 | = <i>´reflect´</i>) or absorption (<b>bc_par_lr</b> |
---|
2214 | = <i>´absorb´</i>) |
---|
2215 | can be set. <br> </p> <p>This lateral boundary |
---|
2216 | conditions should correspond to the |
---|
2217 | lateral boundary condition used for the flow (see <a href="chapter_4.1.html#bc_lr">bc_lr</a>).</p> </td> |
---|
2218 | </tr> <tr> <td style="vertical-align: top;"><p><a name="bc_par_ns"></a><b>bc_par_ns</b></p> |
---|
2219 | </td> <td style="vertical-align: top;">C*15</td> |
---|
2220 | <td style="vertical-align: top;"><i>´cyclic´</i></td> |
---|
2221 | <td style="vertical-align: top;"> <p>Lateral |
---|
2222 | boundary condition (y-direction) for particle |
---|
2223 | transport. </p> <p>By default, cyclic boundary conditions |
---|
2224 | are used along y. |
---|
2225 | Alternatively, reflection (<b>bc_par_ns</b> |
---|
2226 | = <i>´reflect´</i>) or absorption (<b>bc_par_ns</b> |
---|
2227 | = <i>´absorb´</i>) |
---|
2228 | can be set.<br> </p> |
---|
2229 | This lateral boundary conditions should correspond to the lateral |
---|
2230 | boundary condition used for the flow (see <a href="chapter_4.1.html#bc_ns">bc_ns</a>).</td> </tr> |
---|
2231 | <tr> <td style="vertical-align: top;"> <p><a name="bc_par_t"></a><b>bc_par_t</b></p> |
---|
2232 | </td> <td style="vertical-align: top;">C*15</td> |
---|
2233 | <td style="vertical-align: top;"><i>´absorb´</i></td> |
---|
2234 | <td style="vertical-align: top;"> <p>Top boundary |
---|
2235 | condition for particle transport. </p> <p>By default, |
---|
2236 | particles are absorbed at the top boundary. |
---|
2237 | Alternatively, a reflection condition can be set by <b>bc_par_t</b> |
---|
2238 | = <i>´reflect´</i>.</p> </td> |
---|
2239 | </tr> <tr> <td style="vertical-align: top;"><p><a name="density_ratio"></a><b>density_ratio</b></p> |
---|
2240 | </td> <td style="vertical-align: top;">R (10)</td> |
---|
2241 | <td style="vertical-align: top;"> <p><i>0.0, 9</i> |
---|
2242 | *<br> <i>9999999.9</i></p> </td> <td style="vertical-align: top;"> <p>Ratio of the density |
---|
2243 | of the fluid and the density of the |
---|
2244 | particles. </p> <p>With the default value<i> </i>the |
---|
2245 | particles are weightless and transported passively with the resolved |
---|
2246 | scale flow. |
---|
2247 | In case of <span style="font-weight: bold;">density_ratio</span> |
---|
2248 | /= |
---|
2249 | 0.0 particles have a mass and hence inertia so that their velocity |
---|
2250 | deviates more or less from the velocity of the surrounding flow. |
---|
2251 | Particle velocity is calculated analytically and depends on (besides |
---|
2252 | the density ratio and the current velocity difference between particles |
---|
2253 | and surrounding fluid) the |
---|
2254 | particle radius which is determined via <a href="#radius">radius</a> |
---|
2255 | as well as on the molecular viscosity (assumed as 1.461E-5 m<sup>2</sup>/s). |
---|
2256 | </p> <p>If <b>density_ratio</b> = <i>1.0</i>, |
---|
2257 | the particle density |
---|
2258 | corresponds to the density of the surrounding fluid and the particles |
---|
2259 | do not feel any buoyancy. Otherwise, particles will be accelerated |
---|
2260 | upwards (<b>density_ratio</b> > <i>1.0</i>) |
---|
2261 | or downwards (<b>density_ratio</b> < <i>1.0</i>).<br> |
---|
2262 | </p> <p>With several groups of particles (see <a href="chapter_4.2.html#number_of_particle_groups">number_of_particle_groups</a>), |
---|
2263 | each group can be assigned a different value. If the number of values |
---|
2264 | given for <span style="font-weight: bold;">density_ratio</span> |
---|
2265 | is less than the number of |
---|
2266 | groups defined by <span style="font-weight: bold;">number_of_particle_groups</span>, |
---|
2267 | then the last assigned value is used for all remaining groups. This |
---|
2268 | means that by default the particle density ratio for all groups will be |
---|
2269 | <span style="font-style: italic;">0.0</span>.</p> |
---|
2270 | </td> </tr> <tr><td align="left" valign="top"><a name="dt_dopts"></a><span style="font-weight: bold;">dt_dopts</span></td><td align="left" valign="top">R</td><td align="left" valign="top"><i>value of <a href="chapter_4.2.html#dt_data_output">dt_data_<br>output</a></i></td><td align="left" valign="top"><p lang="en-GB"><font face="Thorndale"><font face="Thorndale, serif">Temporal |
---|
2271 | interval</font> at which time series data of particle quantities |
---|
2272 | shall be output (</font>in <font face="Thorndale">s). </font></p> |
---|
2273 | <span lang="en-GB"><font face="Thorndale">If |
---|
2274 | particle advection is switched on (see</font></span><font><span style="font-family: thorndale;"> <a href="#dt_prel">dt_prel</a>) |
---|
2275 | this parameter can be used to assign |
---|
2276 | th</span></font><span lang="en-GB"><font face="Thorndale">e temporal |
---|
2277 | interval at which time series of particle quantities shall be output. |
---|
2278 | Output is written in NetCDF format on local file <a href="chapter_3.4.html#DATA_1D_PTS_NETCDF">DATA_1D_PTS_NETCDF</a>.<br><br>The |
---|
2279 | following list gives a short description of the quantities |
---|
2280 | available. Most quantities are averages over all particles. The |
---|
2281 | quantity name given in the first column is identical to the respective |
---|
2282 | name of the variable on the NetCDF file (see section <a href="chapter_4.5.1.html">4.5.1</a> for a general |
---|
2283 | description of the NetCDF files).<br><br>In case of using |
---|
2284 | more than one particle group (see <a href="#number_of_particle_groups">number_of_particle_groups</a>), |
---|
2285 | seperate time series are output for each of the groups. The long names |
---|
2286 | of the variables in the NetCDF file containing the respective |
---|
2287 | timeseries all end with the string</font><span style="font-style: italic; font-family: monospace;">' PG ##'</span><font face="Thorndale">, where ## is the number of the particle |
---|
2288 | group (<span style="font-style: italic;">01</span>, <span style="font-style: italic;">02</span>, etc.). <br> </font></span><table style="text-align: left; width: 631px; height: 652px;" border="1" cellpadding="2" cellspacing="2"><tbody><tr><td align="undefined" valign="undefined"><span style="color: rgb(255, 0, 0);">tnpt</span></td><td align="undefined" valign="undefined">total number of |
---|
2289 | particles</td></tr><tr><td align="undefined" valign="undefined"><span style="color: rgb(255, 0, 0);">x_</span></td><td align="undefined" valign="undefined">particle |
---|
2290 | x-coordinate with respect to the particle origin (in m)</td></tr><tr><td align="undefined" valign="undefined"><span style="color: rgb(255, 0, 0);">y_</span></td><td align="undefined" valign="undefined">particle |
---|
2291 | y-coordinate with respect to the particle origin (in m)</td></tr><tr><td align="undefined" valign="undefined"><span style="color: rgb(255, 0, 0);">z_</span></td><td align="undefined" valign="undefined">particle |
---|
2292 | z-coordinate with respect to the particle origin (in m)</td></tr><tr><td align="undefined" valign="undefined"><span style="color: rgb(255, 0, 0);">z_abs</span></td><td align="undefined" valign="undefined">absolute |
---|
2293 | particle z-coordinate (in m)</td></tr><tr><td align="undefined" valign="undefined"><span style="color: rgb(255, 0, 0);">u</span></td><td align="undefined" valign="undefined">u particle |
---|
2294 | velocity component (in m/s)</td></tr><tr><td align="undefined" valign="undefined"><span style="color: rgb(255, 0, 0);">v</span></td><td align="undefined" valign="undefined">v particle |
---|
2295 | velocity component (in m/s)</td></tr><tr><td align="undefined" valign="undefined"><span style="color: rgb(255, 0, 0);">w</span></td><td align="undefined" valign="undefined">w particle |
---|
2296 | velocity component (in m/s)</td></tr><tr><td align="undefined" valign="undefined"><span style="color: rgb(255, 0, 0);">u"</span></td><td align="undefined" valign="undefined">subgrid-scale u |
---|
2297 | particle velocity component (in m/s)</td></tr><tr><td align="undefined" valign="undefined"><span style="color: rgb(255, 0, 0);">v"</span></td><td align="undefined" valign="undefined">subgrid-scale v |
---|
2298 | particle velocity component (in m/s)</td></tr><tr><td align="undefined" valign="undefined"><span style="color: rgb(255, 0, 0);">w"</span></td><td align="undefined" valign="undefined">subgrid-scale w |
---|
2299 | particle velocity component (in m/s)</td></tr><tr><td align="undefined" valign="undefined"><span style="color: rgb(255, 0, 0);">npt_up</span></td><td align="undefined" valign="undefined">total number of |
---|
2300 | upward moving particles</td></tr><tr><td align="undefined" valign="undefined"><span style="color: rgb(255, 0, 0);">w_up</span></td><td align="undefined" valign="undefined">vertical |
---|
2301 | velocity of the upward moving particles (in m/s)</td></tr><tr><td align="undefined" valign="undefined"><span style="color: rgb(255, 0, 0);">w_down</span></td><td align="undefined" valign="undefined">vertical |
---|
2302 | velocity of the downward moving particles (in m/s)</td></tr><tr><td align="undefined" valign="undefined"><span style="color: rgb(255, 0, 0);">npt_max</span></td><td align="undefined" valign="undefined">maximum number |
---|
2303 | of particles in a subdomain (=tnpt for non-parallel runs)</td></tr><tr><td align="undefined" valign="undefined"><span style="color: rgb(255, 0, 0);">npt_min</span></td><td align="undefined" valign="undefined">minimum number |
---|
2304 | of particles in a subdomain (=tnpt for non-parallel runs)</td></tr><tr><td align="undefined" valign="undefined"><span style="color: rgb(255, 0, 0);">x*2</span></td><td align="undefined" valign="undefined">variance of the |
---|
2305 | particle x-coordinate with respect to <span style="color: rgb(255, 0, 0);">x_ </span>(in m<sup>2</sup>)</td></tr><tr><td align="undefined" valign="undefined"><span style="color: rgb(255, 0, 0);">y*2</span></td><td align="undefined" valign="undefined">variance of the |
---|
2306 | particle y-coordinate with respect to <span style="color: rgb(255, 0, 0);">y_</span> (in m<sup>2</sup>)</td></tr><tr><td align="undefined" valign="undefined"><span style="color: rgb(255, 0, 0);">z*2</span></td><td align="undefined" valign="undefined">variance of the |
---|
2307 | particle z-coordinate with respect to <span style="color: rgb(255, 0, 0);">z_</span> (in m<sup>2</sup>)</td></tr><tr><td align="undefined" valign="undefined"><span style="color: rgb(255, 0, 0);">u*2</span></td><td align="undefined" valign="undefined">variance of the |
---|
2308 | u particle velocity component with respect to <span style="color: rgb(255, 0, 0);">u </span>(in m<sup>2</sup>/s<sup>2</sup>)</td></tr><tr><td align="undefined" valign="undefined">v*2</td><td align="undefined" valign="undefined">variance of the |
---|
2309 | v particle velocity component with respect to <span style="color: rgb(255, 0, 0);">v </span>(in m<sup>2</sup>/s<sup>2</sup>)</td></tr><tr><td align="undefined" valign="undefined">w*2</td><td align="undefined" valign="undefined">variance of the |
---|
2310 | w particle velocity component with respect to <span style="color: rgb(255, 0, 0);">w </span>(in m<sup>2</sup>/s<sup>2</sup>)</td></tr><tr><td align="undefined" valign="undefined">u"2</td><td align="undefined" valign="undefined">variance of the |
---|
2311 | subgrid-scale u particle velocity component with respect to <span style="color: rgb(255, 0, 0);">u" </span>(in m<sup>2</sup>/s<sup>2</sup>)</td></tr><tr><td align="undefined" valign="undefined">v"2</td><td align="undefined" valign="undefined">variance of the |
---|
2312 | subgrid-scale v particle velocity component with respect to <span style="color: rgb(255, 0, 0);">v" </span>(in m<sup>2</sup>/s<sup>2</sup>)</td></tr><tr><td align="undefined" valign="undefined">w"2</td><td align="undefined" valign="undefined">variance of the |
---|
2313 | subgrid-scale w particle velocity component with respect to <span style="color: rgb(255, 0, 0);">w" </span>(in m<sup>2</sup>/s<sup>2</sup>)</td></tr><tr><td align="undefined" valign="undefined">npt*2</td><td align="undefined" valign="undefined">variance of the |
---|
2314 | number of particles with respect to the average number of particles per |
---|
2315 | subdomain</td></tr></tbody></table><span lang="en-GB"></span><span lang="en-GB"></span></td></tr><tr><td align="left" valign="top"><a name="dt_min_part"></a><span style="font-weight: bold;">dt_min_part</span></td><td align="left" valign="top">R</td><td align="left" valign="top"><span style="font-style: italic;">0.0002</span></td><td align="left" valign="top">Minimum value for the |
---|
2316 | particle timestep when SGS velocities are used (in s).<br><br>For |
---|
2317 | a further explanation see package parameter <a href="#use_sgs_for_particles">use_sgs_for_particles</a>.</td></tr><tr> |
---|
2318 | <td style="vertical-align: top;"> <p><a name="dt_write_particle_data"></a><b>dt_write_particle_</b> |
---|
2319 | <b>data</b></p> </td> <td style="vertical-align: top;">R<br> </td> <td style="vertical-align: top;"><i>9999999.9</i></td> |
---|
2320 | <td style="vertical-align: top;"> <p>Temporal |
---|
2321 | interval for output of particle data (in s). </p> <p>T<span lang="en-GB"></span><a href="#pr1d"><span lang="en-GB"></span></a><span lang="en-GB"></span><span lang="en-GB"><font face="Thorndale">his |
---|
2322 | parameter can be used to |
---|
2323 | assign the temporal interval at which particle data shall be output.</font></span> |
---|
2324 | Data are output to |
---|
2325 | the local file <a href="chapter_3.4.html#PARTICLE_DATA">PARTICLE_DATA</a>. |
---|
2326 | <span style="font-family: mon;">See the file description |
---|
2327 | for more |
---|
2328 | details about its format</span>. </p> <p>By |
---|
2329 | default, no particle data are output.</p> </td> </tr> |
---|
2330 | <tr> <td style="vertical-align: top;"> <p><a name="dvrp_psize"></a><b>dvrp_psize</b></p> |
---|
2331 | </td> <td style="vertical-align: top;">R<br> </td> |
---|
2332 | <td style="vertical-align: top;">0.2 * dx<br> </td> |
---|
2333 | <td style="vertical-align: top;"> <p>Diameter that |
---|
2334 | the particles is given in visualizations with |
---|
2335 | the <span style="font-weight: bold;">dvrp</span> |
---|
2336 | software (in |
---|
2337 | m). </p> <p>In case that particles are |
---|
2338 | visualized with the <span style="font-weight: bold;">dvrp</span> |
---|
2339 | software (see <a href="chapter_4.5.7.html">chapter |
---|
2340 | 4.5.7</a>), their size can be set by parameter <b>dvrp_psize</b>. |
---|
2341 | All |
---|
2342 | particles are displayed with this same size.<br> </p> <p>Alternatively, |
---|
2343 | the particle diameters can be set with the |
---|
2344 | user-interface in routine <span style="font-family: monospace;">user_init_particles</span> |
---|
2345 | (at the beginning of the simulation) and/or can be redefined after each |
---|
2346 | timestep in routine <tt>user<font style="font-size: 11pt;" size="2">_particle_attributes</font></tt> |
---|
2347 | (both routines can be found in file <tt><font style="font-size: 11pt;" size="2">user_interface.f90</font></tt><font style="font-size: 11pt;" size="2">)</font>. </p> |
---|
2348 | <p><b>Note:</b> This parameter determines exclusively |
---|
2349 | the size |
---|
2350 | under which particles appear in the <span style="font-weight: bold;">dvrp</span> |
---|
2351 | visualization. The flow relevant particle radius is determined via the |
---|
2352 | particle package parameter <a href="#radius">radius</a>!</p> |
---|
2353 | </td> </tr> <tr><td align="left" valign="top"><a name="end_time_prel"></a><span style="font-weight: bold;">end_time_prel</span></td><td align="left" valign="top">R</td><td align="left" valign="top"><span style="font-style: italic;">9999999.9</span></td><td align="left" valign="top">Time of the last release of |
---|
2354 | particles (in s).<br><br>See also <a href="#particle_advection_start">particle_advection_start</a>.</td></tr><tr> |
---|
2355 | <td style="vertical-align: top;"><span style="font-weight: bold;"><a name="initial_weighting_factor"></a>initial_weighting_factor</span></td> |
---|
2356 | <td style="vertical-align: top;">R<br> </td> <td style="vertical-align: top;"><span style="font-style: italic;">1.0</span><br> </td> |
---|
2357 | <td style="vertical-align: top;">Factor to define the real |
---|
2358 | number of initial droplets in a grid box.<br> <br> |
---|
2359 | In case of explicitly simulating cloud droplets (see <a href="chapter_4.1.html#cloud_droplets">cloud_droplets</a>), |
---|
2360 | the real number of initial droplets in a grid box is equal to the |
---|
2361 | initial number of droplets in this box (defined by the particle source |
---|
2362 | parameters <span lang="en-GB"><font face="Thorndale, serif"> </font></span><a href="chapter_4.2.html#pst"><span lang="en-GB"><font face="Thorndale, serif">pst</font></span></a><span lang="en-GB"><font face="Thorndale, serif">, </font></span><a href="chapter_4.2.html#psl"><span lang="en-GB"><font face="Thorndale, serif">psl</font></span></a><span lang="en-GB"><font face="Thorndale, serif">, </font></span><a href="chapter_4.2.html#psr"><span lang="en-GB"><font face="Thorndale, serif">psr</font></span></a><span lang="en-GB"><font face="Thorndale, serif">, </font></span><a href="chapter_4.2.html#pss"><span lang="en-GB"><font face="Thorndale, serif">pss</font></span></a><span lang="en-GB"><font face="Thorndale, serif">, </font></span><a href="chapter_4.2.html#psn"><span lang="en-GB"><font face="Thorndale, serif">psn</font></span></a><span lang="en-GB"><font face="Thorndale, serif">, </font></span><a href="chapter_4.2.html#psb"><span lang="en-GB"><font face="Thorndale, serif">psb</font></span></a><span lang="en-GB"><font face="Thorndale, serif">, </font></span><a href="chapter_4.2.html#pdx"><span lang="en-GB"><font face="Thorndale, serif">pdx</font></span></a><span lang="en-GB"><font face="Thorndale, serif">, </font></span><a href="chapter_4.2.html#pdy"><span lang="en-GB"><font face="Thorndale, serif">pdy</font></span></a> |
---|
2363 | <span lang="en-GB"><font face="Thorndale, serif">and |
---|
2364 | </font></span><a href="chapter_4.2.html#pdz"><span lang="en-GB"><font face="Thorndale, serif">pdz</font></span></a><span lang="en-GB"></span><span lang="en-GB"></span>) |
---|
2365 | times the <span style="font-weight: bold;">initial_weighting_factor</span>.</td> |
---|
2366 | </tr><tr> <td style="vertical-align: top;"> <p><a name="maximum_number_of_particles"></a><b>maximum_number_of_</b> |
---|
2367 | <br> <b>particles</b></p> </td> <td style="vertical-align: top;">I</td> <td style="vertical-align: top;"><i>1000</i></td> |
---|
2368 | <td style="vertical-align: top;"> <p>Maximum number |
---|
2369 | of particles (on a PE). </p> <p>This parameter |
---|
2370 | allows to set the number of particles for which |
---|
2371 | memory must be allocated at the beginning of the run. |
---|
2372 | If this memory becomes insufficient during the run, due to the |
---|
2373 | release of further particles (see <a href="#dt_prel">dt_prel</a>), |
---|
2374 | then more memory is automatically allocated.<br> </p> |
---|
2375 | For runs on several processors, <span style="font-weight: bold;">maximum_number_of_particles</span> |
---|
2376 | defines |
---|
2377 | the maximum number on each PE. This number must be larger than the |
---|
2378 | maximum number of particles initially released in a subdomain.</td> |
---|
2379 | </tr> <tr> <td style="vertical-align: top;"><p><a name="maximum_number_of_tailpoints"></a><b>maximum_number_of_</b> |
---|
2380 | <br> <b>tailpoints</b></p> </td> <td style="vertical-align: top;">I</td> <td style="vertical-align: top;"><i>100</i></td> |
---|
2381 | <td style="vertical-align: top;"> <p>Maximum number |
---|
2382 | of tailpoints that a particle tail can |
---|
2383 | have. </p> <p> <b>maximum_number_of_tailpoints</b> |
---|
2384 | sets the number of descrete points the tail consists of. A new point is |
---|
2385 | added to the particle tails after each time step. If the maximum number |
---|
2386 | of tail |
---|
2387 | points is reached after the corresponding number of timesteps, the |
---|
2388 | oldest respective tail points is deleted within the following |
---|
2389 | timestep. </p> <p>All particle tails have the |
---|
2390 | same number of points. The maximum |
---|
2391 | length of |
---|
2392 | these |
---|
2393 | tails is determined by the value of <b>maximum_number_of_tailpoints</b> |
---|
2394 | and by the minimum distance between each of the adjoining |
---|
2395 | tailpoints, which can be set by <a href="#minimum_tailpoint_distance">minimum_tailpoint_distance</a>. |
---|
2396 | Additionally, it can be determined that the temporal displacement |
---|
2397 | between the current position of the particle and the oldest point of |
---|
2398 | the tail may become not larger than a value to be assigned by <a href="#maximum_tailpoint_age">maximum_tailpoint_age</a>.</p> |
---|
2399 | </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="maximum_tailpoint_age"></a><b>maximum_tailpoint_</b> |
---|
2400 | <br> <b>age</b></p> </td> <td style="vertical-align: top;">R</td> <td style="vertical-align: top;">100000.0</td> <td style="vertical-align: top;"> <p>Maximum age that the |
---|
2401 | end point of a particle tail is allowed to have (in s). </p> |
---|
2402 | <p>If the temporal displacement between the oldest point of a |
---|
2403 | particle tail and the current position of the particle becomes larger |
---|
2404 | than the value given by <b>maximum_tailpoint_age</b>, this |
---|
2405 | oldest |
---|
2406 | point (which defines the end of the tail) is |
---|
2407 | removed. If this time is so small that the number of points defining |
---|
2408 | the particle tail do not exceed the value given by <b>maximum_number_of_tailpoints</b>, |
---|
2409 | then the length the particle tails is a measure for the distance the |
---|
2410 | particle travelled along during the time interval defined via <b>maximum_tailpoint_age</b>, |
---|
2411 | i.e. for the |
---|
2412 | particle velocity. Fast particles will have long tails, slow particles |
---|
2413 | shorter ones (note: this will not neccessarily hold if <a href="#minimum_tailpoint_distance">minimum_tailpoint_distance</a> |
---|
2414 | = <i>0.0</i>).</p> </td> </tr> <tr> |
---|
2415 | <td style="vertical-align: top;"> <p><a name="minimum_tailpoint_distance"></a><b>minimum_tailpoint_distance</b></p> |
---|
2416 | </td> <td style="vertical-align: top;">R</td> |
---|
2417 | <td style="vertical-align: top;"><i>0.0</i></td> |
---|
2418 | <td style="vertical-align: top;"> <p>Minimum |
---|
2419 | distance allowed between two adjacent points of a |
---|
2420 | particle tail (in m). </p> <p>In case of <b>minimum_tailpoint_distance</b> |
---|
2421 | > <i>0.0 </i>the |
---|
2422 | particle tail is extended by a new point only if the distance between |
---|
2423 | its current position and the most recent tail point exceed the |
---|
2424 | distance given via <b>minimum_tailpoint_distance</b>.<br> |
---|
2425 | </p> <p>If the length of the particle tails shall be |
---|
2426 | proportional to |
---|
2427 | the respective particle velocity, the parameter <a href="#maximum_tailpoint_age">maximum_tailpoint_age</a> |
---|
2428 | must also be set appropriately. </p> <b>Note:</b><br> |
---|
2429 | A suitable choice of <b>minimum_tailpoint_distance</b> |
---|
2430 | > <i>0.0</i> is recommended, because then the tail |
---|
2431 | coordinates of |
---|
2432 | slowly moving particles require less memory and can also be drawn |
---|
2433 | faster. The upper limit of <b>minimum_tailpoint_distance</b> |
---|
2434 | should be chosen in a way that the visualized particle |
---|
2435 | tails still appear as smooth lines. Example: with a model domain of |
---|
2436 | 1000 m and a monitor resolution of 1280 * 1024 pixels it |
---|
2437 | should be sufficient to set <b>minimum_tailpoint_distance</b> |
---|
2438 | = <i>5.0</i> |
---|
2439 | (m). </td> </tr> <tr> <td style="vertical-align: top;"><a name="number_of_particle_groups"></a><span style="font-weight: bold;">number_of_particle_groups</span><br> |
---|
2440 | </td> <td style="vertical-align: top;">I<br> </td> |
---|
2441 | <td style="vertical-align: top;">1<br> </td> <td style="vertical-align: top;">Number of particle groups to be |
---|
2442 | used.<br> <br> |
---|
2443 | Each particle group can be assigned its own source region (see <a href="#pdx">pdx</a>, <a href="#psl">psl</a>, |
---|
2444 | <a href="#psr">psr</a>, etc.), particle diameter (<a href="#radius">radius</a>) and particle density ratio (<a href="density_ratio">density_ratio</a>).<br><br>If |
---|
2445 | less values are given for <a href="#pdx">pdx</a>, <a href="#psl">psl</a>, |
---|
2446 | etc. than the number of particle groups, then the last value is used |
---|
2447 | for the remaining values (or the default value, if the user did not set |
---|
2448 | the parameter).<br> <br> |
---|
2449 | The maximum allowed number of particle groups is limited to <span style="font-style: italic;">10</span>.<br> </td> |
---|
2450 | </tr><tr><td align="left" valign="top"><a name="particles_per_point"></a><span style="font-weight: bold;">particles_per_point</span></td><td align="left" valign="top">I</td><td align="left" valign="top">1</td><td align="left" valign="top">Number of particles to be |
---|
2451 | started per point.<br><br>By default, one particle is |
---|
2452 | started at all points of the particle source, defined by <span style="font-family: Thorndale,serif;">the </span><span lang="en-GB"><font face="Thorndale, serif">package |
---|
2453 | parameters </font></span><a href="chapter_4.2.html#pst"><span lang="en-GB"><font face="Thorndale, serif">pst</font></span></a><span lang="en-GB"><font face="Thorndale, serif">, </font></span><a href="chapter_4.2.html#psl"><span lang="en-GB"><font face="Thorndale, serif">psl</font></span></a><span lang="en-GB"><font face="Thorndale, serif">, </font></span><a href="chapter_4.2.html#psr"><span lang="en-GB"><font face="Thorndale, serif">psr</font></span></a><span lang="en-GB"><font face="Thorndale, serif">, </font></span><a href="chapter_4.2.html#pss"><span lang="en-GB"><font face="Thorndale, serif">pss</font></span></a><span lang="en-GB"><font face="Thorndale, serif">, </font></span><a href="chapter_4.2.html#psn"><span lang="en-GB"><font face="Thorndale, serif">psn</font></span></a><span lang="en-GB"><font face="Thorndale, serif">, </font></span><a href="chapter_4.2.html#psb"><span lang="en-GB"><font face="Thorndale, serif">psb</font></span></a><span lang="en-GB"><font face="Thorndale, serif">, </font></span><a href="chapter_4.2.html#pdx"><span lang="en-GB"><font face="Thorndale, serif">pdx</font></span></a><span lang="en-GB"><font face="Thorndale, serif">, </font></span><a href="chapter_4.2.html#pdy"><span lang="en-GB"><font face="Thorndale, serif">pdy</font></span></a> |
---|
2454 | <span lang="en-GB"><font face="Thorndale, serif">and |
---|
2455 | </font></span><a href="chapter_4.2.html#pdz"><span lang="en-GB"><font face="Thorndale, serif">pdz</font></span></a><span lang="en-GB"><font face="Thorndale, serif">.</font></span><span lang="en-GB"></span></td></tr><tr> <td style="vertical-align: top;"> <p><a name="particle_advection_start"></a><b>particle_advection_</b> |
---|
2456 | <br> <b>start</b></p> </td> <td style="vertical-align: top;">R </td> <td style="vertical-align: top;">0.0 </td> <td style="vertical-align: top;"> <p>Time of the first |
---|
2457 | release of particles (in s). </p> <p>If particles are not |
---|
2458 | to be released at the beginning of the |
---|
2459 | run, the release time can be set via <b>particle_advection_start</b>.<br> |
---|
2460 | If particle transport is switched on in a restart run, then <a href="#read_particles_from_restartfile">read_particles_from_restartfile</a> |
---|
2461 | = <span style="font-style: italic;">.F.</span> is |
---|
2462 | also required.</p><p>See also <a href="#end_time_prel">end_time_prel</a>. |
---|
2463 | </p> </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="particle_maximum_age"></a><b>particle_maximum_age</b></p> |
---|
2464 | </td> <td style="vertical-align: top;">R </td> |
---|
2465 | <td style="vertical-align: top;"><i>9999999.9</i> |
---|
2466 | </td> <td style="vertical-align: top;"> <p>Maximum |
---|
2467 | allowed age of particles (in s). </p> <p>If the |
---|
2468 | age of a particle exceeds the time set by <b>particle_maximum_age</b>, |
---|
2469 | the particle as well as its tail is deleted.</p> </td> </tr> |
---|
2470 | <tr> <td style="vertical-align: top;"> <p><a name="pdx"></a><b>pdx</b></p> </td> |
---|
2471 | <td style="vertical-align: top;">R (10)<br> </td> |
---|
2472 | <td style="vertical-align: top;"><i>10 * dx</i> |
---|
2473 | </td> <td style="vertical-align: top;"> <p>Distance |
---|
2474 | along x between particles within a particle source |
---|
2475 | (in m). </p> <p>If the particle source shall be |
---|
2476 | confined to one grid point, |
---|
2477 | the distances given by <span style="font-weight: bold;">pdx</span>, |
---|
2478 | <a href="#pdy">pdy</a> |
---|
2479 | and <a href="#pdz">pdz</a> |
---|
2480 | must be set larger than the respective domain size or <a href="#psl">psl</a> |
---|
2481 | = <a href="#psr">psr</a> has to be set |
---|
2482 | alternatively.<br> |
---|
2483 | </p> <p><span style="font-weight: bold;">pdx</span> |
---|
2484 | can be assigned a different value for each particle group (see <a href="#number_of_particle_groups">number_of_particle_groups</a>).<br> |
---|
2485 | </p> </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="pdy"></a><b>pdy</b></p> |
---|
2486 | </td> <td style="vertical-align: top;">R (10)<br> |
---|
2487 | </td> <td style="vertical-align: top;"><i>10 |
---|
2488 | * dy</i> </td> <td style="vertical-align: top;">Distance |
---|
2489 | along y between |
---|
2490 | particles within a |
---|
2491 | particle source (in m). </td> </tr> <tr> |
---|
2492 | <td style="vertical-align: top;"> <p><a name="pdz"></a><b>pdz</b></p> </td> |
---|
2493 | <td style="vertical-align: top;">R (10)<br> |
---|
2494 | </td> <td style="vertical-align: top;"><i>10 |
---|
2495 | * ( zu(2) - zu(1) )</i> </td> <td style="vertical-align: top;">Distance along z between |
---|
2496 | particles within a particle source |
---|
2497 | (in m). </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="psb"></a><b>psb</b></p> |
---|
2498 | </td> <td style="vertical-align: top;">R (10)<br> |
---|
2499 | </td> <td style="vertical-align: top;"><i>10 |
---|
2500 | * zu(nz/2)</i> </td> <td style="vertical-align: top;">Bottom |
---|
2501 | edge of a particle |
---|
2502 | source (in m). </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="psl"></a><b>psl</b></p> |
---|
2503 | </td> <td style="vertical-align: top;">R (10)<br> |
---|
2504 | </td> <td style="vertical-align: top;"><i>10 |
---|
2505 | * 0.0</i> </td> <td style="vertical-align: top;">Left |
---|
2506 | edge of a particle source |
---|
2507 | (in m). </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="psn"></a><b>psn</b></p> |
---|
2508 | </td> <td style="vertical-align: top;">R (10)<br> |
---|
2509 | </td> <td style="vertical-align: top;"><i>10 |
---|
2510 | * (ny * dy)</i> </td> <td style="vertical-align: top;">Rear |
---|
2511 | (“north”) edge of a |
---|
2512 | particle source (in m). </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="psr"></a><b>psr</b></p> |
---|
2513 | </td> <td style="vertical-align: top;">R (10)<br> |
---|
2514 | </td> <td style="vertical-align: top;"><i>10 |
---|
2515 | * (nx * dx)</i> </td> <td style="vertical-align: top;">Right |
---|
2516 | edge of a particle |
---|
2517 | source (in m). </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="pss"></a><b>pss</b></p> |
---|
2518 | </td> <td style="vertical-align: top;">R (10)<br> |
---|
2519 | </td> <td style="vertical-align: top;"><i>10 |
---|
2520 | * 0.0</i> </td> <td style="vertical-align: top;">Front |
---|
2521 | (“south”) edge of a |
---|
2522 | particle source (in m). </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="pst"></a><b>pst</b></p> |
---|
2523 | </td> <td style="vertical-align: top;">R (10)<br> |
---|
2524 | </td> <td style="vertical-align: top;"><i>10 |
---|
2525 | * zu(nz/2)</i> </td> <td style="vertical-align: top;">Top |
---|
2526 | edge of a particle source |
---|
2527 | (in m). </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="radius"></a><b>radius</b></p> |
---|
2528 | </td> <td style="vertical-align: top;">R (10)</td> |
---|
2529 | <td style="vertical-align: top;"><i>0.0, 9</i>*<br> |
---|
2530 | <i>9999999.9</i></td> <td style="vertical-align: top;">Particle radius (in m).<br> |
---|
2531 | <br>The viscous friction (in case of a velocity difference |
---|
2532 | between |
---|
2533 | particles and surrounding fluid) depends on the particle radius which |
---|
2534 | must be assigned as soon as <a href="chapter_4.2.html#density_ratio">density_ratio</a> |
---|
2535 | /= <i>0.0</i>.<br> <br> |
---|
2536 | With several groups of particles (see <a href="#number_of_particle_groups">number_of_particle_groups</a>), |
---|
2537 | each group can be assigned a different value. If the number of values |
---|
2538 | given for <span style="font-weight: bold;">radius</span> |
---|
2539 | is less than the number of |
---|
2540 | groups defined by <span style="font-weight: bold;">number_of_particle_groups</span>, |
---|
2541 | then the last assigned value is used for all remaining groups. This |
---|
2542 | means that by default the particle radius for all groups will be <span style="font-style: italic;">0.0</span>.<br> </td> |
---|
2543 | </tr><tr> <td style="vertical-align: top;"> <p><a name="random_start_position"></a><b>random_start_position</b></p> |
---|
2544 | </td> <td style="vertical-align: top;">L<br> </td> |
---|
2545 | <td style="vertical-align: top;"><i>.F.</i> </td> |
---|
2546 | <td style="vertical-align: top;"> <p><span style="background: transparent none repeat scroll 0% 50%; -moz-background-clip: initial; -moz-background-origin: initial; -moz-background-inline-policy: initial;"><font color="#000000">Initial position of the</font></span> |
---|
2547 | particles is |
---|
2548 | varied randomly within certain limits. </p> <p>By |
---|
2549 | default, the initial positions of particles within the |
---|
2550 | source excatly correspond with the positions given by <a href="#psl">psl</a>, |
---|
2551 | <a href="#psr">psr</a>, <a href="#psn">psn</a>, |
---|
2552 | <a href="#pss">pss</a>, <a href="#psb">psb</a>, |
---|
2553 | <a href="#pst">pst</a>, <a href="#pdx">pdx</a>, |
---|
2554 | <a href="#pdy">pdy</a>, |
---|
2555 | and<a href="#pdz"> |
---|
2556 | pdz</a>. With <b>random_start_position</b> = <i>.T. |
---|
2557 | </i>the initial |
---|
2558 | positions of the particles are allowed to randomly vary from these |
---|
2559 | positions within certain limits. </p> <p><b>Very |
---|
2560 | important:<br> </b>In case of <b>random_start_position</b> |
---|
2561 | = <i>.T.</i>, the |
---|
2562 | random-number generators on the individual PEs no longer |
---|
2563 | run synchronously. If random disturbances are applied to the velocity |
---|
2564 | field |
---|
2565 | (see <a href="#create_disturbances">create_disturbances</a>), |
---|
2566 | <font color="#000000">then as consequence for parallel |
---|
2567 | runs the |
---|
2568 | realizations of the turbulent flow |
---|
2569 | fields will deviate between runs which used different numbers of PEs!</font></p> |
---|
2570 | </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="read_particles_from_restartfile"></a><b>read_particles_from_</b> |
---|
2571 | <br> <b>restartfile</b></p> </td> <td style="vertical-align: top;">L<br> </td> <td style="vertical-align: top;"><i>.T.</i> </td> |
---|
2572 | <td style="vertical-align: top;"> <p>Read particle |
---|
2573 | data from the previous run. </p> <p>By default, |
---|
2574 | with restart runs particle data is read |
---|
2575 | from file <a href="chapter_3.4.html#PARTICLE_RESTART_DATA_IN">PARTICLE_RESTART_DATA_IN</a>, |
---|
2576 | which is created by the preceding run. If this is not requested or if |
---|
2577 | in a restart run particle transport is switched on for the |
---|
2578 | first time (see <a href="#particle_advection_start">particle_advection_start</a>), |
---|
2579 | then <b>read_particles_from_restartfile</b> = <i>.F.</i> |
---|
2580 | is required.</p> </td> </tr> <tr> <td style="vertical-align: top;"><a name="skip_particles_for_tail"></a><span style="font-weight: bold;">skip_particles_for_tail</span><br> |
---|
2581 | </td> <td style="vertical-align: top;">I<br> </td> |
---|
2582 | <td style="vertical-align: top;"><span style="font-style: italic;">1</span><br> </td> |
---|
2583 | <td style="vertical-align: top;">Limit the number of |
---|
2584 | particle tails.<br> <br> |
---|
2585 | If particle tails are switched on (see <a href="#use_particle_tails">use_particle_tails</a>), |
---|
2586 | every particle is given a tail by default. <span style="font-weight: bold;">skip_particles_for_tail </span>can |
---|
2587 | be used to give only every n'th particle a tail.<br> <br> <span style="font-weight: bold;">Example:</span><br> <span style="font-weight: bold;">skip_particles_for_tail</span> |
---|
2588 | = <span style="font-style: italic;">10</span> means |
---|
2589 | that only every 10th particle will be given a tail.<br> </td> |
---|
2590 | </tr> <tr> <td style="vertical-align: top;"><a name="use_particle_tails"></a><span style="font-weight: bold;">use_particle_tails</span><br> |
---|
2591 | </td> <td style="vertical-align: top;">L<br> </td> |
---|
2592 | <td style="vertical-align: top;"><span style="font-style: italic;">.F.</span><br> </td> |
---|
2593 | <td style="vertical-align: top;">Give particles a tail.<br> |
---|
2594 | <br>A particle tail is defined by the path a particle has moved |
---|
2595 | along starting from some point of time in the past. It consists of a |
---|
2596 | set of descrete points in space which may e.g. be connected by a line |
---|
2597 | in order visualize how the particle has moved.<br> <br> |
---|
2598 | By default, particles have no tail. Parameter <a href="#skip_particles_for_tail">skip_particles_for_tail</a> |
---|
2599 | can be used to give only every n'th particle a tail.<br> <br> |
---|
2600 | The length of the tail is controlled by parameters <a href="#maximum_number_of_tailpoints">maximum_number_of_tailpoints</a>, <a href="#maximum_tailpoint_age">maximum_tailpoint_age</a>, |
---|
2601 | and <a href="#minimum_tailpoint_distance">minimum_tailpoint_distance</a>.<br> |
---|
2602 | </td> </tr><tr><td align="left" valign="top"><a name="use_sgs_for_particles"></a><span style="font-weight: bold;">use_sgs_for_particles</span></td><td align="left" valign="top">L</td><td align="left" valign="top"><span style="font-style: italic;">.F.</span></td><td align="left" valign="top">Use subgrid-scale |
---|
2603 | velocities for particle advection.<br><br>These |
---|
2604 | velocities are calculated from the resolved and subgrid-scale TKE using |
---|
2605 | the Monte-Carlo random-walk method described by Weil et al. (2004, JAS, |
---|
2606 | 61, |
---|
2607 | 2877-2887). When using this method, the timestep for the advancement of |
---|
2608 | the particles is limited by the so-called Lagrangian time scale. This |
---|
2609 | may be smaller than the current LES timestep so that several particle |
---|
2610 | (sub-) timesteps have to be carried out within one LES timestep. In |
---|
2611 | order to limit the number of sub-timesteps (and to limit the CPU-time), |
---|
2612 | the minimum value for the particle timestep is defined by the package |
---|
2613 | parameter <a href="#dt_min_part">dt_min_part</a>.<br><br>Setting |
---|
2614 | <span style="font-weight: bold;">use_sgs_for_particles</span> |
---|
2615 | = <span style="font-style: italic;">.TRUE.</span> |
---|
2616 | automatically forces <a href="chapter_4.1.html#use_upstream_for_tke">use_upstream_for_tke</a> |
---|
2617 | = <span style="font-style: italic;">.TRUE.</span>. |
---|
2618 | This inhibits the occurrence of large (artificial) spatial gradients of |
---|
2619 | the subgrid-scale TKE which otherwise would lead to wrong results for |
---|
2620 | the particle advection.</td></tr><tr> <td style="vertical-align: top;"> <p><a name="vertical_particle_advection"></a><b>vertical_particle_</b> |
---|
2621 | <br> <b>advection</b></p> </td> <td style="vertical-align: top;">L<br> </td> <td style="vertical-align: top;"><i>.T.</i> </td> |
---|
2622 | <td style="vertical-align: top;"> <p>Switch on/off |
---|
2623 | vertical particle transport. </p> <p>By default, |
---|
2624 | particles are transported along all three |
---|
2625 | directions in space. With <b>vertical_particle_advection</b> |
---|
2626 | = <i>.F., </i>the |
---|
2627 | particles will only be transported horizontally.</p> </td> |
---|
2628 | </tr> <tr> <td style="vertical-align: top;"><p><a name="write_particle_statistics"></a><b>write_particle_</b> |
---|
2629 | <br> <b>statistics</b></p> </td> <td style="vertical-align: top;">L<br> </td> <td style="vertical-align: top;"><i>.F.</i> </td> |
---|
2630 | <td style="vertical-align: top;"> <p>Switch on/off |
---|
2631 | output of particle informations.<br> </p> <p><br> |
---|
2632 | For <span style="font-weight: bold;">write_particle_statistics</span> |
---|
2633 | = <span style="font-style: italic;">.T.</span> s<span style="font-family: thorndale,serif;">tatistical |
---|
2634 | informations</span> (e.g. the total number of particles used, the |
---|
2635 | number of particles exchanged between the PEs, etc.) which may be used |
---|
2636 | for debugging are output to the |
---|
2637 | local file <a href="chapter_3.4.html#PARTICLE_DATA">PARTICLE_INFOS</a>. |
---|
2638 | </p> <p><b>Note:</b> For parallel runs files |
---|
2639 | may become very large |
---|
2640 | and performance of PALM may decrease.</p> </td> </tr> |
---|
2641 | </tbody></table><span style="font-weight: bold;"><br> |
---|
2642 | <br></span><span style="font-weight: bold;">Package |
---|
2643 | (<span style="font-weight: bold;">mrun</span> option |
---|
2644 | -p): <span style="font-weight: bold;"><a name="dvrp_graphics"></a>dvrp_graphics</span> |
---|
2645 | |
---|
2646 | NAMELIST group name: <span style="font-weight: bold;">dvrp_graphics_par<br> |
---|
2647 | <br></span></span><table style="text-align: left; width: 100%;" border="1" cellpadding="2" cellspacing="2"> <tbody> <tr> |
---|
2648 | <td style="vertical-align: top;"><font size="4"><b>Parameter |
---|
2649 | name</b></font></td> <td style="vertical-align: top;"><font size="4"><b>Type</b></font></td> |
---|
2650 | <td style="vertical-align: top;"> <p><b><font size="4">Default</font></b> <br> <b><font size="4">value</font></b></p> </td> |
---|
2651 | <td style="vertical-align: top;"><font size="4"><b>Explanation</b></font></td> |
---|
2652 | </tr> <tr> <td style="vertical-align: top;"><p><a name="dt_dvrp"></a><b>dt_dvrp</b></p> |
---|
2653 | </td> <td style="vertical-align: top;">R</td> |
---|
2654 | <td style="vertical-align: top;"><i>9999999.9</i></td> |
---|
2655 | <td style="vertical-align: top;"> <p>Temporal |
---|
2656 | interval of scenes to be displayed with the <span style="font-weight: bold;">dvrp</span> software (in |
---|
2657 | s). </p> <p>Isosurfaces, cross sections and |
---|
2658 | particles can be displayed |
---|
2659 | simultaneous. The display of particles requires that the particle |
---|
2660 | transport is switched on (see <a href="#dt_prel">dt_prel</a>). |
---|
2661 | Objects to be displayed have to be determined with <a href="#mode_dvrp">mode_dvrp</a>. </p> <p>If |
---|
2662 | output of scenes created by dvrp software is switched on |
---|
2663 | (see <a href="#mode_dvrp">mode_dvrp</a>), |
---|
2664 | this parameter can be used to assign the temporal interval at which |
---|
2665 | scenes are to be created (and the respective graphical data |
---|
2666 | is to |
---|
2667 | be output to the streaming server). <span lang="en-GB"><font face="Thorndale">Reference time is the beginning of |
---|
2668 | the simulation, i.e. output takes place at times t = <b>dt_dvrp</b>, |
---|
2669 | 2*<b>dt_dvrp</b>, 3*<b>dt_dvrp</b>, etc. The |
---|
2670 | actual output times can |
---|
2671 | deviate from these theoretical values (see </font></span><a href="#dt_dopr_zeitpunkte"><span lang="en-GB"><font face="Thorndale">dt_dopr</font></span></a><span lang="en-GB"><font face="Thorndale">). |
---|
2672 | Is <b>dt_dvrp</b> < </font></span><a href="chapter_4.1.html#dt"><span lang="en-GB"><font face="Thorndale">dt</font></span></a><span lang="en-GB"><font face="Thorndale">, then |
---|
2673 | scenes are created and |
---|
2674 | output after each time step (if this is requested it should be <b>dt_dvrp</b> |
---|
2675 | = <i>0</i>).</font></span> </p> </td> |
---|
2676 | </tr> <tr> <td style="vertical-align: top;"><p><a name="dvrp_directory"></a><b>dvrp_directory</b></p> |
---|
2677 | </td> <td style="vertical-align: top;">C*80</td> |
---|
2678 | <td style="vertical-align: top;"><i>'default'</i></td> |
---|
2679 | <td style="vertical-align: top;"> <p>Name of the |
---|
2680 | directory into which data created by the <span style="font-weight: bold;">dvrp</span> |
---|
2681 | software shall be saved. </p> <p>By default, |
---|
2682 | the directory name is generated from the user |
---|
2683 | name |
---|
2684 | (see package parameter <a href="#dvrp_username">dvrp_username</a>) |
---|
2685 | and the base file name (given as the argument of <span style="font-weight: bold;">mrun</span> option -d) as <span style="font-style: italic;">'<user |
---|
2686 | name>/<base file name>'</span>.</p> </td> |
---|
2687 | </tr> <tr> <td style="vertical-align: top;"><p><a name="dvrp_file"></a><b>dvrp_file</b></p> |
---|
2688 | </td> <td style="vertical-align: top;">C*80</td> |
---|
2689 | <td style="vertical-align: top;"><i>'default'</i></td> |
---|
2690 | <td style="vertical-align: top;"> <p>Name of the |
---|
2691 | file into which data created by the <span style="font-weight: bold;">dvrp</span> |
---|
2692 | software shall be output. </p> <p>This |
---|
2693 | parameter can be given a value only in case of <a href="#dvrp_output">dvrp_output</a> |
---|
2694 | = <span style="font-style: italic;">'local'</span><i> |
---|
2695 | </i>which |
---|
2696 | determines that the data created by <span style="font-weight: bold;">dvrp</span> |
---|
2697 | is output to a local file (on the machine where PALM is executed). |
---|
2698 | Apart from the default, it is only allowed to assign <span style="font-style: italic;">'/dev/null'</span> (which |
---|
2699 | means that no output is really stored). This can be used for special |
---|
2700 | runtime measurements of the <span style="font-weight: bold;">dvrp</span> |
---|
2701 | software.</p> </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="dvrp_host"></a><b>dvrp_host</b></p> |
---|
2702 | </td> <td style="vertical-align: top;">C*80</td> |
---|
2703 | <td style="vertical-align: top;"> <p><i>'origin.rvs.</i> |
---|
2704 | <br>u<i>ni- hanover.de'</i></p> </td> <td style="vertical-align: top;"> <p>Name of the computer |
---|
2705 | to which data created by the <span style="font-weight: bold;">dvrp</span> |
---|
2706 | software shall be |
---|
2707 | transferred. </p> <p>In case of <a href="#dvrp_output">dvrp_output</a> |
---|
2708 | = <span style="font-style: italic;">'rtsp'</span> |
---|
2709 | only the default |
---|
2710 | value is allowed (streaming server of |
---|
2711 | the RRZN). For <a href="#dvrp_output">dvrp_output</a> |
---|
2712 | = <span style="font-style: italic;">'local'</span><i> |
---|
2713 | </i>the |
---|
2714 | assigned value is ignored.</p> </td> </tr> <tr> |
---|
2715 | <td style="vertical-align: top;"> <p><a name="dvrp_output"></a><b>dvrp_output</b></p> |
---|
2716 | </td> <td style="vertical-align: top;">C*10</td> |
---|
2717 | <td style="vertical-align: top;"><i>'rtsp'</i></td> |
---|
2718 | <td style="vertical-align: top;"> <p>Output mode |
---|
2719 | for the <span style="font-weight: bold;">dvrp</span> |
---|
2720 | software. <br> <br> </p> |
---|
2721 | The following settings are allowed:<br> <br> <table style="text-align: left; width: 100%;" cellpadding="2" cellspacing="2"> <tbody> <tr> <td style="vertical-align: top;"><i>'rtsp'</i></td> |
---|
2722 | <td style="vertical-align: top;">Data created by the <span style="font-weight: bold;">dvrp</span> |
---|
2723 | software is transferred using |
---|
2724 | a special transmission protocol to a so-called streaming server, which |
---|
2725 | is able to continuously transfer visualization data with a |
---|
2726 | high transmission rate. <br> |
---|
2727 | Additionally, with this output mode a |
---|
2728 | set of files is generated automatically |
---|
2729 | within a directory on the streaming server (beside the visualization |
---|
2730 | data e.g. some html-files) which can be used to |
---|
2731 | visualize the data via an internet-browser plugin. Host |
---|
2732 | (streaming-server) and directory can be defined by the user with <a href="#dvrp_host">dvrp_host</a> |
---|
2733 | and <a href="#dvrp_directory">dvrp_directory</a>.</td> |
---|
2734 | </tr> <tr> <td style="vertical-align: top;"><i>'ftp'</i></td> |
---|
2735 | <td style="vertical-align: top;">Data created by the <span style="font-weight: bold;">dvrp</span> |
---|
2736 | software is transferred to the destination host (see <a href="#dvrp_host">dvrp_host</a> |
---|
2737 | and <a href="#dvrp_directory">dvrp_directory</a>) |
---|
2738 | using ftp.</td> </tr> <tr> <td style="vertical-align: top;"><i>'local'</i></td> |
---|
2739 | <td style="vertical-align: top;">Data created by the <span style="font-weight: bold;">dvrp</span> |
---|
2740 | software is output locally on a file defined by <a href="#dvrp_file">dvrp_file |
---|
2741 | </a>.</td> </tr> </tbody> </table> <br> |
---|
2742 | </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="dvrp_password"></a><b>dvrp_password</b></p> |
---|
2743 | </td> <td style="vertical-align: top;">C*80</td> |
---|
2744 | <td style="vertical-align: top;">'********'</td> <td style="vertical-align: top;"> <p>Password for the |
---|
2745 | computer to which data created by the <span style="font-weight: bold;">dvrp</span> software is to |
---|
2746 | be |
---|
2747 | transferred. </p> <p>Assigning a password is |
---|
2748 | only necessary in case of <a href="#dvrp_output">dvrp_output</a> |
---|
2749 | = <span style="font-style: italic;">'ftp'</span>. |
---|
2750 | For <a href="#dvrp_output">dvrp_output</a> |
---|
2751 | = <span style="font-style: italic;">'rtsp'</span><i> |
---|
2752 | </i>the default |
---|
2753 | value must not be changed!</p> </td> </tr> <tr> |
---|
2754 | <td style="vertical-align: top;"> <p><a name="dvrp_username"></a><b>dvrp_username</b></p> |
---|
2755 | </td> <td style="vertical-align: top;">C*80</td> |
---|
2756 | <td style="vertical-align: top;"><br> </td> <td style="vertical-align: top;"> <p>User name of a valid |
---|
2757 | account on the computer to which data |
---|
2758 | created by the <span style="font-weight: bold;">dvrp</span> |
---|
2759 | software |
---|
2760 | is to be |
---|
2761 | transferred. </p> <p>Assigning a value to this |
---|
2762 | parameter is required in case of <a href="#dvrp_output">dvrp_output</a> |
---|
2763 | = <span style="font-style: italic;">'rtsp'</span> |
---|
2764 | or <span style="font-style: italic;">'ftp'</span>.</p> |
---|
2765 | </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="mode_dvrp"></a><b>mode_dvrp</b></p> |
---|
2766 | </td> <td style="vertical-align: top;">C*20 |
---|
2767 | <br>(10)</td> <td style="vertical-align: top;"><i>10 |
---|
2768 | * ''</i></td> <td style="vertical-align: top;"> |
---|
2769 | <p>Graphical objects (isosurfaces, slicers, particles) which are |
---|
2770 | to be created by the <span style="font-weight: bold;">dvrp</span> |
---|
2771 | software. </p> <p>Several different objects can |
---|
2772 | be assigned simultaneously and |
---|
2773 | will be displayed in the same scene. Allowed values for <span style="font-weight: bold;">mode_dvrp</span> are <span style="font-style: italic;">'isosurface#'</span> |
---|
2774 | (isosurface), <span style="font-style: italic;">'slicer#'</span> |
---|
2775 | (cross sections), and <span style="font-style: italic;">'particles'</span>. |
---|
2776 | Within the strings the hash character ("#") has to be replaced by a |
---|
2777 | digit ≤9. Up to 10 objects |
---|
2778 | can be assigned at the same time, e.g. : </p> <blockquote><b>mode_dvrp</b> |
---|
2779 | = <span style="font-style: italic;">'isosurface2'</span><i>, |
---|
2780 | 'slicer1', |
---|
2781 | 'particles', 'slicer2'</i></blockquote> <p>In this |
---|
2782 | case one isosurface, two cross sections, and particles |
---|
2783 | will be created. The quantities for which an isosurface are to be |
---|
2784 | created have to be selected with |
---|
2785 | the parameter <a href="#data_output">data_output</a>, |
---|
2786 | those for cross sections with <a href="#data_output">data_output</a> |
---|
2787 | (<span style="font-weight: bold;">data_output</span> |
---|
2788 | also determines the |
---|
2789 | orientation of the cross section, thus xy, xz, or yz). Since for <span style="font-weight: bold;">data_output</span> and <span style="font-weight: bold;">data_output</span> lists of |
---|
2790 | variables may be |
---|
2791 | assigned, the digit at the end of the <span style="font-weight: bold;">mode_dvrp</span>-string |
---|
2792 | selects the quantity, which is given |
---|
2793 | at the respective position in the respective list (e.g. <span style="font-style: italic;">'isosurface2'</span> |
---|
2794 | selects the quantity |
---|
2795 | given in the second position of <span style="font-weight: bold;">data_output</span>). |
---|
2796 | If e.g. <span style="font-weight: bold;">data_output</span> |
---|
2797 | and <span style="font-weight: bold;">data_output</span> |
---|
2798 | are assigned as <b>data_output</b> = <span style="font-style: italic;">'u_xy'</span><i>, |
---|
2799 | 'w_xz', 'v_yz'</i> and <b>data_output</b> = <span style="font-style: italic;">'pt'</span><i>, |
---|
2800 | 'u', 'w' </i>are |
---|
2801 | indicated, then - assuming the above assignment of <span style="font-weight: bold;">mode_dvrp</span> - an |
---|
2802 | isosurface of u, a |
---|
2803 | horizontal cross section of u and |
---|
2804 | a vertical cross section (xz) of w is created. For locations of the |
---|
2805 | cross sections see <a href="#data_output">data_output</a>. |
---|
2806 | The theshold value for which the isosurface is |
---|
2807 | to be created can be defined with parameter <a href="#threshold">threshold</a>.<br> |
---|
2808 | </p> <p>The vertical extension of the displayed domain is |
---|
2809 | given by <a href="#nz_do3d">nz_do3d</a>.<br> </p> |
---|
2810 | <p>The vertical extension of the displayed domain is given by <a href="#nz_do3d">nz_do3d</a>. </p> <p><b>Assignments |
---|
2811 | of mode_dvrp must correspond to those of data_output |
---|
2812 | and |
---|
2813 | data_output! </b>If e.g. <b>data_output</b> = <span style="font-style: italic;">'pt_xy'</span> |
---|
2814 | and <b>data_output</b> |
---|
2815 | = 'w'<i> </i>was set, then only the digit "1" is allowed |
---|
2816 | for <b>mode_dvrp</b>, |
---|
2817 | thus <span style="font-style: italic;">'isosurface1'</span> |
---|
2818 | and/or <span style="font-style: italic;">'slicer1'</span><i>.</i> |
---|
2819 | </p> <p>Further details about using the <span style="font-weight: bold;">dvrp</span> software are |
---|
2820 | given in <a href="chapter_4.5.7.html">chapter |
---|
2821 | 4.5.7</a>.<br> </p> <b>Note:</b><br> |
---|
2822 | The declaration color charts to be |
---|
2823 | used still have to be given "manually" in subroutine <span style="font-family: monospace;">user_dvrp_coltab</span> |
---|
2824 | (file <tt><font style="font-size: 11pt;" size="2">user_interface.f90</font></tt>). |
---|
2825 | <br>A change of particle colors and/or particle diameters (e.g. |
---|
2826 | according |
---|
2827 | to the local characteristics of the flow field) to be used for the |
---|
2828 | visualization, must be carried out by adding respective code extensions |
---|
2829 | to <tt><font style="font-size: 11pt;" size="2">user_particle_attributes</font></tt> |
---|
2830 | (in file <tt><font style="font-size: 11pt;" size="2">user_interface.f90</font></tt>). </td> |
---|
2831 | </tr> <tr> <td style="vertical-align: top;"><a name="slicer_range_limits_dvrp"></a><span style="font-weight: bold;">slicer_range_limits_<br> |
---|
2832 | dvrp</span></td> <td style="vertical-align: top;">R(2,10)</td> |
---|
2833 | <td style="vertical-align: top;"><span style="font-style: italic;">10 |
---|
2834 | * (-1,1)</span></td> <td style="vertical-align: top;">Ranges |
---|
2835 | of values to which a color table has to be mapped (units of the |
---|
2836 | respective quantity).<br> <br> |
---|
2837 | In case that slicers have to be displayed (see <a href="#threshold">mode_dvrp</a>), |
---|
2838 | this parameter defines the ranges of values of the respective |
---|
2839 | quantities to which the colortable in use has to be mapped. If e.g. a |
---|
2840 | temperature slice shall be displayed, the colortable defines colors |
---|
2841 | from blue to red, and <span style="font-weight: bold;">slicer_range_limits_dvrp</span> |
---|
2842 | = 290.0, 305.0 then areas with temperature of 290 K are displayed in |
---|
2843 | dark blue and those with 305.0 are displayed in dark red. Temperatures |
---|
2844 | within these limits will be displayed by a continuous color gradient |
---|
2845 | from blue to red and Temperatures outside the limits will |
---|
2846 | be displayed either in dark blue or in dark red.<br> <br> |
---|
2847 | Up to ten different ranges can be assigned in case that more than one |
---|
2848 | slicer has to be displayed.<br> <br> |
---|
2849 | See <a href="#threshold">mode_dvrp</a> |
---|
2850 | for the declaration of color charts.</td> </tr> <tr> |
---|
2851 | <td style="vertical-align: top;"> <p><a name="superelevation"></a><b>superelevation</b></p> |
---|
2852 | </td> <td style="vertical-align: top;">R</td> |
---|
2853 | <td style="vertical-align: top;"><i>1.0</i></td> |
---|
2854 | <td style="vertical-align: top;"> <p>Superelevation |
---|
2855 | factor for the vertical coordinate. </p> <p>For |
---|
2856 | domains with unfavorable ratio between the vertical and |
---|
2857 | the horizontal size |
---|
2858 | (the vertical size is determined by <a href="#nz_do3d">nz_do3d</a>) |
---|
2859 | a <span style="font-weight: bold;">superelevation</span> |
---|
2860 | /= <span style="font-style: italic;">1.0</span> may |
---|
2861 | be used. If e.g. the |
---|
2862 | horizontal size is substantially larger |
---|
2863 | than the vertical size, a <span style="font-weight: bold;">superelevation</span> |
---|
2864 | much larger than <span style="font-style: italic;">1.0</span> |
---|
2865 | should |
---|
2866 | be used, since otherwise the domain appears as a |
---|
2867 | "flat disk" in the visualization and thus the vertical direction is |
---|
2868 | only very poorly resolved.</p> </td> </tr> <tr> |
---|
2869 | <td style="vertical-align: top;"> <p><a name="superelevation_x"></a><b>superelevation_x</b></p> |
---|
2870 | </td> <td style="vertical-align: top;">R<br> </td> |
---|
2871 | <td style="vertical-align: top; font-style: italic;">1.0<br> |
---|
2872 | </td> <td style="vertical-align: top;"> <p>Superelevation |
---|
2873 | factor for the horizontal (x) coordinate. </p> <p>This |
---|
2874 | parameter can be used to stretch the displayed domain |
---|
2875 | along the x-direction. See also <a href="#superelevation">superelevation</a>.</p> |
---|
2876 | </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="superelevation_y"></a><b>superelevation_y</b></p> |
---|
2877 | </td> <td style="vertical-align: top;">R<br> </td> |
---|
2878 | <td style="vertical-align: top; font-style: italic;">1.0<br> |
---|
2879 | </td> <td style="vertical-align: top;">Superelevation |
---|
2880 | factor for the |
---|
2881 | horizontal (y) coordinate. <p>This parameter can be |
---|
2882 | used to |
---|
2883 | stretch the displayed domain along the y-direction. See also <a href="#superelevation">superelevation</a>.</p> </td> |
---|
2884 | </tr> <tr> <td style="vertical-align: top;"><p><a name="threshold"></a><b>threshold</b></p> |
---|
2885 | </td> <td style="vertical-align: top;">R(10)<br> |
---|
2886 | </td> <td style="vertical-align: top; font-style: italic;">0.0<br> |
---|
2887 | </td> <td style="vertical-align: top;"> <p>Threshold |
---|
2888 | value for which an isosurface is to be created by |
---|
2889 | the <span style="font-weight: bold;">dvrp</span> |
---|
2890 | software. </p> <p>If the creation of |
---|
2891 | isosurfaces is switched on via |
---|
2892 | parameter <a href="#mode_dvrp">mode_dvrp</a>, |
---|
2893 | then the respective threshold value for which the isosurface is to be |
---|
2894 | created can be assigned via <b>threshold</b>. If several |
---|
2895 | isosurfaces |
---|
2896 | are given by <b>mode_dvrp</b>, then an individual |
---|
2897 | threshold value for |
---|
2898 | each isosurface can be assigned. The order of the threshold values |
---|
2899 | refers to the order of the isosurfaces given by <b>mode_dvrp</b>.</p> |
---|
2900 | </td> </tr> </tbody></table><span style="font-weight: bold;"><span style="font-weight: bold;"><br> |
---|
2901 | </span></span><span style="font-weight: bold;"><span style="font-weight: bold;">Package (<span style="font-weight: bold;">mrun</span> |
---|
2902 | option -p): <span style="font-weight: bold;"><a name="spectra"></a>spectra</span> |
---|
2903 | NAMELIST group name: <span style="font-weight: bold;">spectra_par<br> |
---|
2904 | <br></span></span></span><table style="text-align: left; width: 100%;" border="1" cellpadding="2" cellspacing="2"> <tbody> <tr> |
---|
2905 | <td style="vertical-align: top;"><font size="4"><b>Parameter |
---|
2906 | name</b></font></td> <td style="vertical-align: top;"><font size="4"><b>Type</b></font></td> |
---|
2907 | <td style="vertical-align: top;"> <p><b><font size="4">Default</font></b> <br> <b><font size="4">value</font></b></p> </td> |
---|
2908 | <td style="vertical-align: top;"><font size="4"><b>Explanation</b></font></td> |
---|
2909 | </tr> <tr> <td style="vertical-align: top;"><p><a name="averaging_interval_sp"></a><b>averaging_interval_sp</b></p> |
---|
2910 | </td> <td style="vertical-align: top;">R<br> </td> |
---|
2911 | <td style="vertical-align: top;"><span style="font-style: italic;">value of <a href="chapter_4.2.html#averaging_interval">averaging_<br> |
---|
2912 | interval</a></span> </td> <td style="vertical-align: top;"> <p>Averaging interval |
---|
2913 | for spectra output to local |
---|
2914 | file <font color="#000000"><font color="#000000"><a href="chapter_3.4.html#DATA_1D_SP_NETCDF">DATA_1D_SP_NETCDF</a> |
---|
2915 | </font></font>and/or <a href="chapter_3.4.html#PLOTSP_X_DATA">PLOTSP_X_DATA</a> |
---|
2916 | / <a href="chapter_3.4.html#PLOTSP_Y_DATA">PLOTSP_Y_DATA</a> |
---|
2917 | (in s). </p> <p>If |
---|
2918 | this parameter is given a non-zero value, temporally |
---|
2919 | averaged spectra data are output. By default, spectra data data are not |
---|
2920 | subject to temporal averaging. The interval length is limited by the |
---|
2921 | parameter <a href="#dt_dosp">dt_dosp</a>. In any |
---|
2922 | case <b>averaging_interval_sp</b> <= <b>dt_dosp |
---|
2923 | </b>must |
---|
2924 | hold.</p>If an interval is defined, then by default the average |
---|
2925 | is calculated |
---|
2926 | from the data values of all timesteps lying within this interval. The |
---|
2927 | number of time levels entering into the average can be reduced with the |
---|
2928 | parameter <a href="chapter_4.2.html#dt_averaging_input_pr">dt_averaging_input_pr</a>. |
---|
2929 | <p>If |
---|
2930 | an averaging interval can not be completed at the end of a run, it will |
---|
2931 | be finished at the beginning of the next restart run. Thus for restart |
---|
2932 | runs, averaging intervals do not |
---|
2933 | necessarily begin at the beginning of the run.</p></td> </tr> |
---|
2934 | <tr> <td style="vertical-align: top;"><b><a name="comp_spectra_level"></a>comp_spectra_level</b></td> |
---|
2935 | <td style="vertical-align: top;">I(10)</td> <td style="vertical-align: top;"><i>no level</i></td> |
---|
2936 | <td style="vertical-align: top;"> <p>Vertical level |
---|
2937 | for which horizontal spectra are to be |
---|
2938 | calculated and output (gridpoints).<br> </p> <br> |
---|
2939 | Spectra can be calculated for up to ten levels.</td> </tr> |
---|
2940 | <tr><td style="vertical-align: top;"><p><a name="data_output_sp"></a><b>data_output_sp</b></p></td><td style="vertical-align: top;">C*10 (10)</td><td style="vertical-align: top;"><i>10 * ' '</i></td><td style="vertical-align: top;"><p>Quantities for which |
---|
2941 | horizontal spectra are to be calculated |
---|
2942 | and output.</p> <p>Allowed values are: <b>data_output_sp</b> |
---|
2943 | = <span style="font-style: italic;">'u'</span>, <span style="font-style: italic;">'v'</span>, <span style="font-style: italic;">'w'</span>, <span style="font-style: italic;">'pt'</span>, <span style="font-style: italic;">'q'</span>.<br> </p> |
---|
2944 | <p>Spectra are calculated using the FFT-method defined by <a href="chapter_4.1.html#fft_method">fft_method</a>.</p> |
---|
2945 | <p>By default spectra data are output to the local file <a href="chapter_3.4.html#DATA_1D_SP_NETCDF">DATA_1D_SP_NETCDF</a>. |
---|
2946 | The file's format is NetCDF. Further details about processing |
---|
2947 | NetCDF data are given in chapter <a href="chapter_4.5.1.html">4.5.1</a>.</p><p>The |
---|
2948 | temporal interval of the output times of profiles is |
---|
2949 | assigned via the parameter <a href="chapter_4.2.html#dt_dosp">dt_dosp</a>. </p><p>The |
---|
2950 | vertical levels for which spectra are to be computed and output must be |
---|
2951 | given by parameter <font><a href="chapter_4.2.html#comp_spectra_level"><span lang="en-GB"><font face="Thorndale">comp_spectra_level</font></span></a></font>. |
---|
2952 | </p><span style="font-weight: bold;">Note:</span><br> |
---|
2953 | Beside <span style="font-weight: bold;">data_output_sp</span>, |
---|
2954 | values <span style="font-weight: bold;">must</span> |
---|
2955 | be given for each of the |
---|
2956 | parameters, <font><a href="chapter_4.2.html#comp_spectra_level"><span lang="en-GB"><font face="Thorndale">comp_spectra_level</font></span></a></font>, |
---|
2957 | and <font><a href="chapter_4.2.html#spectra_direction"><span lang="en-GB"><font face="Thorndale">spectra_direction</font></span></a></font>, |
---|
2958 | otherwise <span style="font-weight: bold;">no</span> |
---|
2959 | output will be |
---|
2960 | created!<br><br><br> |
---|
2961 | Calculation of spectra requires cyclic boundary conditions |
---|
2962 | along the respective directions (see <a href="chapter_4.1.html#bc_lr">bc_lr</a> |
---|
2963 | and <a href="chapter_4.1.html#bc_ns">bc_ns</a>).For |
---|
2964 | historical reasons, data can also be output in ASCII-format on local |
---|
2965 | files <a href="chapter_3.4.html#PLOTSP_X_DATA">PLOTSP_X_DATA</a> |
---|
2966 | and/or <a href="chapter_3.4.html#PLOTSP_Y_DATA">PLOTSP_Y_DATA</a> |
---|
2967 | (depending on the direction(s) along which spectra are to be |
---|
2968 | calculated; see <font><a href="chapter_4.2.html#spectra_direction"><span lang="en-GB"><font face="Thorndale">spectra_direction</font></span></a>),</font> |
---|
2969 | which are readable by the graphic software <span style="font-weight: bold;">profil</span>. See |
---|
2970 | parameter <a href="chapter_4.2.html#data_output_format">data_output_format</a> |
---|
2971 | for defining the format in which data shall be output. Within |
---|
2972 | these file, the spectra are ordered with respect to their |
---|
2973 | output times. Spectra can also be temporally averaged (see <a href="chapter_4.2.html#averaging_interval_sp">averaging_interval_sp</a> |
---|
2974 | ). <font><a href="chapter_4.2.html#spectra_direction"><span lang="en-GB"></span></a>Each data point of a |
---|
2975 | spectrum is output in a single line (1st column: |
---|
2976 | wavenumber, 2nd column: spectral coefficient). If spectra are to be |
---|
2977 | calculated and output for more than one height (see </font><font><a href="chapter_4.2.html#comp_spectra_level"><span lang="en-GB"><font face="Thorndale">comp_spectra_level</font></span></a></font><font>), |
---|
2978 | the spectral coefficients for the further heighs can be found in the |
---|
2979 | subsequent columns. </font>The order |
---|
2980 | of the data in the file follows the order used in the assignment for <b>data_output_sp</b> |
---|
2981 | (<b>data_output_sp</b> = <span style="font-style: italic;">'u'</span>, <span style="font-style: italic;">'v'</span>,… |
---|
2982 | means that the file starts with the spectra of the u-component, |
---|
2983 | followed by the v-component spectra, etc.). Additional to the files |
---|
2984 | PLOTSP_X_DATA and PLOTSP_Y_DATA which contain |
---|
2985 | the data, |
---|
2986 | PALM creates NAMELIST parameter files (local name <a href="chapter_3.4.html#PLOTSP_X_PAR">PLOTSP_X_PAR</a> |
---|
2987 | and <a href="chapter_3.4.html#PLOTSP_X_PAR">PLOTSP_Y_PAR</a>) |
---|
2988 | which can be used as parameter input file for the plot software <a href="http://www.muk.uni-hannover.de/institut/software/profil_intro.html">profil</a>. |
---|
2989 | Spectra can be directly plotted with <span style="font-weight: bold;">profil</span> |
---|
2990 | using the data and the corresponding parameter file. The |
---|
2991 | plot layout is |
---|
2992 | steered via the parameter input file. The vertical levels for which |
---|
2993 | spectra are to be plotted must be given by <font><a href="chapter_4.2.html#plot_spectra_level"><span lang="en-GB"><font face="Thorndale">plot_spectra_level</font></span></a></font><font><a href="chapter_4.2.html#comp_spectra_level"><span lang="en-GB"></span></a></font>. <span style="font-weight: bold;"></span>Otherwise, no |
---|
2994 | spectra |
---|
2995 | will appear on the plot, although data are available on file. All |
---|
2996 | parameter values can be changed by editing the parameter |
---|
2997 | input |
---|
2998 | file.<span style="font-weight: bold;"><br></span></td></tr><tr> |
---|
2999 | <td style="vertical-align: top;"> <p><a name="dt_dosp"></a><b>dt_dosp</b></p> |
---|
3000 | </td> <td style="vertical-align: top;">R</td> |
---|
3001 | <td style="vertical-align: top;"><i>value of |
---|
3002 | <a href="chapter_4.2.html#dt_data_output">dt_data_<br>output</a></i></td> |
---|
3003 | <td style="vertical-align: top;"> <p>Temporal |
---|
3004 | interval at which spectra data shall be output |
---|
3005 | (in s). </p> <p><span lang="en-GB"><font face="Thorndale">If output of |
---|
3006 | horizontal spectra is switched on (see </font></span><a href="#data_output_sp"><span lang="en-GB"><font face="Thorndale">data_output_sp</font></span></a><span lang="en-GB"><font face="Thorndale">), </font></span><span lang="en-GB"><font face="Thorndale">this |
---|
3007 | parameter can be used to |
---|
3008 | assign the temporal interval at which spectral data shall be |
---|
3009 | output. </font></span><span lang="en-GB"><font face="Thorndale">Output can be skipped at the beginning of a |
---|
3010 | simulation using parameter <a href="#skip_time_dosp">skip_time_dosp</a>, |
---|
3011 | which has zero value by default. </font></span><span lang="en-GB"></span><span lang="en-GB"><font face="Thorndale">Reference |
---|
3012 | time is the beginning of |
---|
3013 | the simulation, i.e. output takes place at times t = <span style="font-weight: bold;">skip_time_dosp</span> + <b>dt_dosp</b>, |
---|
3014 | <span style="font-weight: bold;">skip_time_dosp</span> |
---|
3015 | + 2*<b>dt_dosp</b>, skip_time_dosp + 3*<b>dt_dosp</b>, |
---|
3016 | etc. The actual output times can |
---|
3017 | deviate from these theoretical values (see </font></span><a href="#dt_dopr_zeitpunkte"><span lang="en-GB"><font face="Thorndale">dt_dopr</font></span></a><span lang="en-GB"><font face="Thorndale">). |
---|
3018 | If <b>dt_dosp</b> < </font></span><a href="chapter_4.1.html#dt"><span lang="en-GB"><font face="Thorndale">dt</font></span></a><span lang="en-GB"><font face="Thorndale">, then |
---|
3019 | spectral data are output |
---|
3020 | after each time step (if this is requested it should be <b>dt_dosp</b> |
---|
3021 | = <i>0</i>).</font></span> </p> </td> |
---|
3022 | </tr> <tr> <td style="vertical-align: top;"><p><a name="plot_spectra_level"></a><b>plot_spectra_level</b></p> |
---|
3023 | </td> <td style="vertical-align: top;">I(10)</td> |
---|
3024 | <td style="vertical-align: top;"><i>no level</i></td> |
---|
3025 | <td style="vertical-align: top;"> <p>Vertical |
---|
3026 | level(s) for which horizontal spectra are to be |
---|
3027 | plotted (in gridpoints). </p> <p>This parameter |
---|
3028 | only affects the display of spectra in plots |
---|
3029 | created with <span style="font-weight: bold;">profil</span>. |
---|
3030 | The |
---|
3031 | spectral data created and output to file are exclusively determined via |
---|
3032 | <font><a href="#comp_spectra_level"><span lang="en-GB"><font face="Thorndale">comp_spectra_level</font></span></a></font>.</p> |
---|
3033 | </td> </tr> <tr> <td style="vertical-align: top;"><a name="skip_time_dosp"></a><span style="font-weight: bold;">skip_time_dosp</span></td> |
---|
3034 | <td style="vertical-align: top;">R<br> </td> <td style="vertical-align: top;"><span style="font-style: italic;">value of <a href="chapter_4.2.html#skip_time_data_output">skip_time_<br>data_output</a></span> |
---|
3035 | </td> <td style="vertical-align: top;">No output of |
---|
3036 | spectra data before this interval has passed (in s).<br><br>This |
---|
3037 | parameter causes that data output activities are starting not before |
---|
3038 | this interval |
---|
3039 | (counting from the beginning of the simulation, t=0) has passed. <br><br><span style="font-weight: bold;">Example:</span><br>If |
---|
3040 | the user has set <a href="#dt_dosp">dt_dosp</a> = <span style="font-style: italic;">3600.0</span> and <span style="font-weight: bold;">skip_time_dosp</span> = <span style="font-style: italic;">1800.0</span>, then the |
---|
3041 | first output will be done at t = 5400 s. </td> </tr> |
---|
3042 | <tr> <td style="vertical-align: top;"> <p><a name="spectra_direction"></a><b>spectra_direction</b></p> |
---|
3043 | </td> <td style="vertical-align: top;">C*2 (10)</td> |
---|
3044 | <td style="vertical-align: top;"><i>10 * ' '</i></td> |
---|
3045 | <td style="vertical-align: top;"> <p>Direction(s) |
---|
3046 | along which spectra are to be calculated. </p> <p>Allowed |
---|
3047 | values are <span style="font-style: italic;">'x'</span>, |
---|
3048 | <span style="font-style: italic;">'y'</span> and <span style="font-style: italic;">'xy'</span>. For |
---|
3049 | every quantity given by <a href="#data_output_sp">data_output_sp</a> |
---|
3050 | a corresponding |
---|
3051 | direction<span style="font-weight: bold;"> </span>must |
---|
3052 | be assigned.<br> </p> <p>Calculation of spectra |
---|
3053 | requires cyclic boundary conditions |
---|
3054 | along the respective directions (see <a href="chapter_4.1.html#bc_lr">bc_lr</a> |
---|
3055 | and <a href="chapter_4.1.html#bc_ns">bc_ns</a>).</p> |
---|
3056 | </td> </tr> </tbody></table><span style="font-weight: bold;"><span style="font-weight: bold;"><span style="font-weight: bold;"><br> |
---|
3057 | </span></span></span><h3 style="line-height: 100%;"><br><a href="chapter_4.1.html"><img src="left.gif" name="Grafik1" align="bottom" border="2" height="32" width="32"></a><a href="index.html"><img src="up.gif" name="Grafik2" align="bottom" border="2" height="32" width="32"></a><a href="chapter_4.3.html"><img src="right.gif" name="Grafik3" align="bottom" border="2" height="32" width="32"></a></h3><span style="font-style: italic;">Last change:</span> |
---|
3058 | $Id: chapter_4.2.html 61 2007-03-12 05:42:06Z raasch $ <span style="font-weight: bold;"><span style="font-weight: bold;"><br> |
---|
3059 | </span></span><br></body></html> |
---|