1 | <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN"> |
---|
2 | <html><head> |
---|
3 | <meta content="text/html; charset=windows-1252" http-equiv="CONTENT-TYPE"><title>PALM chapter 4.4</title> |
---|
4 | <meta content="StarOffice 7 (Win32)" name="GENERATOR"> |
---|
5 | <meta content="Siegfried Raasch" name="AUTHOR"> <meta content="20041013;13430732" name="CREATED"> <meta content="20041117;11162734" name="CHANGED"> <meta content="parallel LES model" name="KEYWORDS"> <style> |
---|
6 | <!-- |
---|
7 | @page { size: 21cm 29.7cm } |
---|
8 | --> |
---|
9 | </style></head> |
---|
10 | <body style="direction: ltr;" lang="en-US"><h3 style="line-height: 100%;">4.4.2 A parameter set for ocean runs</h3> |
---|
11 | <p style="line-height: 100%;">... to be completed ...</p><p style="line-height: 100%;">In this chapter a brief, |
---|
12 | simple and |
---|
13 | complete parameter set is described, which can be used to carry out a |
---|
14 | model run. The presented example is available via <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/INSTALL/example_p3d">example |
---|
15 | file</a> and can be used (together with the <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/INSTALL/.mrun.config">configuration |
---|
16 | file</a> described in <a href="chapter_3.2.html">chapter |
---|
17 | 3.2)</a> for the execution of a simple model run. </p> |
---|
18 | <p style="line-height: 100%;">This run simulates a |
---|
19 | quasi-stationary, |
---|
20 | convective, atmospheric boundary layer with <font color="#000000">zero |
---|
21 | mean horizontal |
---|
22 | wind.</font> For evaluation purposes, cross sections and |
---|
23 | horizontally averaged vertical |
---|
24 | profiles of typical boundary layer variables |
---|
25 | are output at the end of the run. The run shall be carried out in |
---|
26 | batch mode on the IBM Regatta "hanni" of the HLRN.</p> |
---|
27 | <p style="line-height: 100%;">The parameter file necessary |
---|
28 | to carry |
---|
29 | out a run must be provided to the model as an input file under the |
---|
30 | local name <a href="chapter_3.4.html#PARIN">PARIN</a> |
---|
31 | and has the following contents:</p> |
---|
32 | <pre style="line-height: 100%;">&inipar <a href="chapter_4.1.html#nx">nx</a> = <span style="font-style: italic;">39</span>, <a href="chapter_4.1.html#ny">ny</a> = <span style="font-style: italic;">39</span>, <a href="chapter_4.1.html#nz">nz</a> = <span style="font-style: italic;">40</span>,<br> <a href="chapter_4.1.html#dx">dx</a> = <span style="font-style: italic;">50.0</span>, <a href="chapter_4.1.html#dy">dy</a> = <span style="font-style: italic;">50.0</span>, <a href="chapter_4.1.html#dz">dz</a> = <span style="font-style: italic;">50.0</span>,<br> <a href="chapter_4.1.html#dz_stretch_level">dz_stretch_level</a> = <span style="font-style: italic;">1200.0</span>,<br> <a href="chapter_4.1.html#fft_method">fft_method</a> = <span style="font-style: italic;">'temperton-algorithm'</span>,<br> <a href="chapter_4.1.html#initializing_actions">initializing_actions</a> = <span style="font-style: italic;">'set_constant_profiles'</span>,<br> <a href="chapter_4.1.html#ug_surface">ug_surface</a> = <span style="font-style: italic;">0.0</span>, <a href="chapter_4.1.html#vg_surface">vg_surface</a> = <span style="font-style: italic;">0.0</span>,<br> <a href="chapter_4.1.html#pt_vertical_gradient">pt_vertical_gradient</a> = <span style="font-style: italic;">0.0</span>, <span style="font-style: italic;">1.0</span>,<br> <a href="chapter_4.1.html#pt_vertical_gradient_level">pt_vertical_gradient_level</a> = <span style="font-style: italic;">0.0</span>, <span style="font-style: italic;">800.0</span>,<br> <a href="chapter_4.1.html#surface_heatflux">surface_heatflux</a> = <span style="font-style: italic;">0.1</span>, <a href="chapter_4.1.html#bc_pt_b">bc_pt_b</a> = <span style="font-style: italic;">'neumann'</span>,/<br><br>&d3par <a href="chapter_4.2.html#end_time">end_time</a> = <span style="font-style: italic;">3600.0</span>,<br> <a href="chapter_4.2.html#create_disturbances">create_disturbances</a> = <span style="font-style: italic;">.T.</span>,<br> <a href="chapter_4.2.html#dt_disturb">dt_disturb</a> = <span style="font-style: italic;">150.0</span>, <a href="chapter_4.2.html#disturbance_energy_limit">disturbance_energy_limit</a> = <span style="font-style: italic;">0.01</span>,<br> <a href="chapter_4.2.html#dt_run_control">dt_run_control</a> = <span style="font-style: italic;">0.0</span>,<br> <a href="chapter_4.2.html#data_output">data_output</a> = <span style="font-style: italic;">'w_xy'</span>, <span style="font-style: italic;">'w_xz'</span>, <span style="font-style: italic;">'w_xz_av'</span>, <span style="font-style: italic;">'pt_xy'</span>, <span style="font-style: italic;">'pt_xz'</span>,<br> <a href="chapter_4.2.html#dt_data_output">dt_data_output</a> = <span style="font-style: italic;">900.0</span>,<br> <a href="chapter_4.2.html#dt_data_output_av">dt_data_output_av</a> = <span style="font-style: italic;">1800.0</span>,<br> <a href="chapter_4.2.html#averaging_interval">averaging_interval</a> = <span style="font-style: italic;">900.0</span>,<br> <a href="chapter_4.2.html#dt_averaging_input">dt_averaging_input</a> = <span style="font-style: italic;">10.0</span>,<br> <a href="chapter_4.2.html#section_xy">section_xy</a> = <span style="font-style: italic;">2</span>, <span style="font-style: italic;">10</span>, <a href="chapter_4.2.html#section_xz">section_xz</a> = <span style="font-style: italic;">20</span>,<br> <a href="chapter_4.2.html#data_output_2d_on_each_pe">data_output_2d_on_each_pe</a> = <span style="font-style: italic;">.F.</span>,<br> <a href="chapter_4.2.html#dt_dopr">dt_dopr</a> = <span style="font-style: italic;">900.0</span>, <a href="chapter_4.2.html#averaging_interval_pr">averaging_interval_pr</a> = <span style="font-style: italic;">600.0</span>,<br> <a href="chapter_4.2.html#dt_averaging_input_pr">dt_averaging_input_pr</a> = <span style="font-style: italic;">10.0</span>,<br> <a href="chapter_4.2.html#data_output_pr">data_output_pr</a> = <span style="font-style: italic;">'#pt'</span>, <span style="font-style: italic;">'w”pt”'</span>, <span style="font-style: italic;">'w*pt*'</span>, <span style="font-style: italic;">'wpt'</span>, <span style="font-style: italic;">'w*2'</span>, <span style="font-style: italic;">'pt*2'</span>,<br> <a href="chapter_4.2.html#cross_profiles">cross_profiles</a> = <span style="font-style: italic;">' pt '</span>, <span style="font-style: italic;">' w"pt" w*pt* wpt '</span>, <span style="font-style: italic;">' w*2 '</span>, <span style="font-style: italic;">' pt*2 '</span>,<br> <a href="chapter_4.2.html#cross_xtext">cross_xtext</a> = <span style="font-style: italic;">'pot. temperature in K'</span>,<br> <span style="font-style: italic;">'heat flux in K ms>->1'</span>,<br> <span style="font-style: italic;">'velocity variance in m>2s>->2'</span>,<br> <span style="font-style: italic;">'temperature variance in K>2'</span>,<br> <a href="chapter_4.2.html#z_max_do1d">z_max_do1d</a> = <span style="font-style: italic;">1500.0</span>, /</pre><p style="line-height: 100%;"><br><br></p> |
---|
33 | <p style="line-height: 100%;">The initialization |
---|
34 | parameters (<tt><font style="font-size: 10pt;" size="2">&inipar</font></tt>) |
---|
35 | are located at the beginning of the file. For analysis of a |
---|
36 | convective boundary layer of approx. 1000 m thickness the horizontal |
---|
37 | size of the model domain should amount to at least 2 km x 2 km. In |
---|
38 | order to resolve the convective structures a grid spacing of <b>dx</b> |
---|
39 | = |
---|
40 | <b>dy</b> = <b>dz</b> = <i>50 m</i> |
---|
41 | is enough, since the typical |
---|
42 | diameter of convective plumes is more than 100 m. Thereby the |
---|
43 | upper array index in the two horizontal directions needs to be <b>nx</b> |
---|
44 | = <b>ny</b> = <i>39</i>. <font color="#000000">Since in |
---|
45 | each case the lower array index has the value 0, 40 grid points are |
---|
46 | used along both horizontal directions.</font> In the vertical |
---|
47 | direction |
---|
48 | the domain must be high enough to include the entrainment processes at |
---|
49 | the top of the boundary layer as well as the propagation of gravity |
---|
50 | waves, which were stimulated by |
---|
51 | the convection. However, in the stably stratified region the grid |
---|
52 | resolution has not necessarily to be as high as within the boundary |
---|
53 | layer. This can be obtained by a vertical stretching of the grid |
---|
54 | starting |
---|
55 | from 1200 m via <b>dz_stretch_level</b> = <i>1200.0 |
---|
56 | m.</i> This saves |
---|
57 | grid points and computing time. <font color="#800000">T</font><font color="#000000">he |
---|
58 | upper boundary of the model is located at (see </font><a href="chapter_4.1.html#dz_stretch_factor"><font color="#000000">dz_stretch_factor</font></a><font color="#000000">) |
---|
59 | … m (computed by the model)</font>.</p><p style="line-height: 100%;">Fast Fourier transformations are |
---|
60 | calculated using the Temperton-algorithm, which -on the IBM Regatta- is |
---|
61 | faster than the default system-specific algorithm (from IBM essl |
---|
62 | library).</p><p style="line-height: 100%;">The |
---|
63 | initial profiles for |
---|
64 | wind and temperature can be assigned via <b>initializing_actions</b> |
---|
65 | = <span style="font-style: italic;">'set_constant_profiles'</span>. |
---|
66 | The wind speed, constant with |
---|
67 | height, amounts to <b>ug_surface</b> = <b>vg_surface</b> |
---|
68 | = <i>0.0 m/s</i>. In order |
---|
69 | to allow for a fast onset of convection, a neutral stratified layer up |
---|
70 | to z |
---|
71 | = 800 m capped by an inversion with dtheta/dz = 1K/100 m is given: |
---|
72 | <b>pt_vertical_gradient</b> = <i>0.0, 1.0</i>, |
---|
73 | <b>pt_vertical_gradient_level</b> = <i>0.0, 800.0.</i> |
---|
74 | The surface |
---|
75 | temperature, which by default amounts to 300 K, provides the fixed |
---|
76 | point for the temperature profile (see <a href="chapter_4.1.html#pt_surface">pt_surface</a>). |
---|
77 | Convection is driven by a given, near-surface sensible heat flux via <b>surface_heatflux</b> |
---|
78 | = <i>0.1 K m/s.</i> A given surface sensible heta flux |
---|
79 | requires the |
---|
80 | bottom boundary condition for potential temperature to be <b>bc_pt_b</b> |
---|
81 | = |
---|
82 | <span style="font-style: italic;">'neumann'</span> . |
---|
83 | Thus |
---|
84 | all initialization parameters are determined. These can not be |
---|
85 | changed during the run (also not for restart runs). </p> |
---|
86 | <p style="line-height: 100%;">Now the run parameters (<tt><font style="font-size: 10pt;" size="2">&d3par</font></tt>) |
---|
87 | must be specified. To produce a quasi stationary boundary layer the |
---|
88 | simulated time should be at least one hour, i.e. <b>end_time</b> |
---|
89 | = <i>3600 |
---|
90 | s.</i> To stimulate convection, the initially homogeneous (zero) |
---|
91 | wind |
---|
92 | field must be disturbed (<b>create_disturbances</b> = <i>.T.</i>). |
---|
93 | These perturbations should be repeated in a temporal interval of |
---|
94 | <b>dt_disturb</b> = <i>150.0 s</i> until the |
---|
95 | energy of the |
---|
96 | perturbations exceeds the value <b>disturbance_energy_limit</b> |
---|
97 | = 0.<i>01 |
---|
98 | m<sup>2</sup>/s<sup>2</sup></i>. After |
---|
99 | each time step run time |
---|
100 | informations (e.g. size of the timestep, maximum velocities, etc.) are |
---|
101 | to be written to the local file <a href="chapter_3.4.html#RUN_CONTROL">RUN_CONTROL</a> |
---|
102 | (<b>dt_run_control</b> = <i>0.0 s</i>).</p><p style="line-height: 100%;">Instantaneous cross section data |
---|
103 | of vertical velocity (<span style="font-style: italic;">w</span>) |
---|
104 | and potential temperature (<span style="font-style: italic;">pt</span>) |
---|
105 | are to be output for horizontal (<span style="font-style: italic;">xy</span>) |
---|
106 | and vertical (<span style="font-style: italic;">xz</span>) |
---|
107 | cross sections, and additionally, time averaged (<span style="font-style: italic;">av</span>) vertical cross |
---|
108 | section data are to be output for the vertical velocity: <span style="font-weight: bold;">data_output</span> = <span style="font-style: italic;">'w_xy'</span>, <span style="font-style: italic;">'w_xz'</span>, <span style="font-style: italic;">'w_xz_av'</span>, <span style="font-style: italic;">'pt_xy'</span>, <span style="font-style: italic;">'pt_xz'</span>. Output of |
---|
109 | instantaneous (time averaged) data is done after each 900 (1800)s: <span style="font-weight: bold;">dt_data_output</span> = <span style="font-style: italic;">900.0</span>, <span style="font-weight: bold;">dt_data_output_av</span> = <span style="font-style: italic;">1800.0</span>. The |
---|
110 | averaged data are time averaged over the last 900.0 s, where the |
---|
111 | temporal interval of data entering the average is 10 s: <span style="font-weight: bold;">averaging_interval</span> = |
---|
112 | <span style="font-style: italic;">900.0</span>, <span style="font-weight: bold;">dt_averaging_input</span> = |
---|
113 | <span style="font-style: italic;">10.0</span>. |
---|
114 | Horizontal cross sections are output for vertical levels with grid |
---|
115 | index k=2 and k=10, vertical cross sections are output for index j=20: <span style="font-weight: bold;">section_xy</span> = <span style="font-style: italic;">2</span>, <span style="font-style: italic;">10</span>, <span style="font-weight: bold;">section_xz</span> = <span style="font-style: italic;">20</span>. For runs on |
---|
116 | more than one processor, cross section data are collected and output on |
---|
117 | PE0: <span style="font-weight: bold;">data_output_2d_on_each_pe</span> |
---|
118 | = <span style="font-style: italic;">.F.</span>.</p><p style="line-height: 100%;">Output |
---|
119 | of vertical profiles is to be done after each 900 s. The profiles shall |
---|
120 | be temporally averaged<font color="#000000"> over the last |
---|
121 | <font color="#000000">600 </font>seconds, </font>whereby |
---|
122 | the temporal interval of the profiles entering the average has to be |
---|
123 | 10 s: <b>dt_dopr</b> = <i>900.0 s</i>, <b>averaging_interval_pr</b> |
---|
124 | = |
---|
125 | <i>600.0 s</i>, <b>dt_averaging_input_pr</b> = |
---|
126 | <i>10.0 s.</i> The temperature |
---|
127 | profile including the initial temperature profile (therefore <span style="font-style: italic;">'#pt'</span>), |
---|
128 | the subgrid scale, resolved and total vertical sensible heat flux as |
---|
129 | well as the variances of the vertical velocity and the potential |
---|
130 | temperature are to be output: <b>data_output_pr</b> |
---|
131 | = <span style="font-style: italic;">'#pt'</span><i>, |
---|
132 | 'w"pt”', |
---|
133 | 'w*pt*', 'wpt', 'w*2', 'pt*2'</i>.</p><p style="line-height: 100%;">If the data output format for |
---|
134 | graphic software <span style="font-weight: bold;">profil</span> |
---|
135 | is selected (see <a href="chapter_4.2.html#data_output_format">data_output_format</a>), |
---|
136 | the temperature |
---|
137 | profile and the individual variances are to be drawn into independent |
---|
138 | coordinate systems, and in contrast to this all heat flux profiles are |
---|
139 | to |
---|
140 | be |
---|
141 | drawn into the same system: <b>cross_profiles</b> = <span style="font-style: italic;">'pt'</span><i>, |
---|
142 | 'w"pt"w*pt*wpt', 'w*2', 'pt*2'</i>. The legend of the x |
---|
143 | axes of these systems is set to <b>cross_xtext</b>= <i>'pot. |
---|
144 | temperature in K', 'heat flux in K ms>->1', 'velocity |
---|
145 | variance |
---|
146 | in m>2s>->2', 'temperature variance in K>2'</i>. |
---|
147 | The profiles are to be drawn up to a height level of <b>z_max_do1d</b> |
---|
148 | = |
---|
149 | <i>1500.0 m</i>. </p> |
---|
150 | <p style="line-height: 100%;">Before starting the mo<font color="#000000">del |
---|
151 | on the parallel computer, the number of processing elements must be |
---|
152 | specified.</font> Since relatively few grid points are used for |
---|
153 | this run, choosing of e.g. 8 PEs is sufficient. By default, a 1d domain |
---|
154 | decomposition along x is used on the IBM-Regatta, which means that a |
---|
155 | virtual processor topology (grid) of 8*1 (x*y) is used. (<span style="font-weight: bold;">Note:</span> the user may |
---|
156 | adjust this |
---|
157 | default domain decomposition with the help of the parameters <a href="chapter_4.1.html#npex">npex</a> |
---|
158 | and <a href="chapter_4.1.html#npey">npey</a>). |
---|
159 | </p><p style="line-height: 100%;">Provided that the |
---|
160 | parameters |
---|
161 | file described above are set within the file </p> |
---|
162 | <ul> <pre style="margin-bottom: 0.5cm; line-height: 100%;"><font style="font-size: 10pt;" size="2">~/palm/current_version/JOBS/example/INPUT/example_p3d</font></pre></ul><p style="line-height: 100%;">and that the conditions |
---|
163 | mentioned in the |
---|
164 | first sections of <a href="chapter_3.2.html">chapter |
---|
165 | 3.2</a> are met, the model run can be started with the command </p> |
---|
166 | <p style="line-height: 100%;"><font face="Cumberland, monospace"><font style="font-size: 10pt;" size="2">mrun |
---|
167 | -d example -h ibmh -K parallel -X 8 -T 8 -t 1800 -q cdev -r |
---|
168 | “d3# xy# xz# pr#”</font></font></p> |
---|
169 | <p style="line-height: 100%;">The output files will appear |
---|
170 | in the |
---|
171 | directories </p> |
---|
172 | <blockquote style="line-height: 100%;"><tt><font style="font-size: 10pt;" size="2">~/palm/current_version/JOBS/example/MONITORING</font></tt><font style="font-size: 10pt;" size="2"><br> </font><tt><font style="font-size: 10pt;" size="2">~/palm/current_version/JOBS/example/OUTPUT |
---|
173 | ,</font></tt></blockquote> |
---|
174 | <p style="line-height: 100%;">while the job protocol will |
---|
175 | appear in |
---|
176 | directory <font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace">~/</font></font><tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace">job_queue</font></font></tt>. |
---|
177 | <br> |
---|
178 | </p> |
---|
179 | <hr><p style="line-height: 100%;"><br><font color="#000080"><font color="#000080"><a href="chapter_4.4.1.html"><font color="#000080"><img style="border: 2px solid ; width: 32px; height: 32px;" alt="" name="Grafik1" src="left.gif"></font></a><a href="index.html"><font color="#000080"><img name="Grafik2" src="up.gif" align="bottom" border="2" height="32" width="32"></font></a><a href="chapter_4.5.html"><font color="#000080"><img style="border: 2px solid ; width: 32px; height: 32px;" alt="" name="Grafik3" src="right.gif"></font></a></font></font></p><p style="line-height: 100%;"><i>Last change: |
---|
180 | </i>$Id: chapter_4.4.2.html 97 2007-06-21 08:23:15Z maronga $ |
---|
181 | <br> <br> |
---|
182 | </p> |
---|
183 | </body></html> |
---|