Home

Context Navigation

source: palm/trunk/DOC/app/chapter_3.5.4.html @ 402

Last change on this file since 402 was 343, checked in by maronga, 15 years ago
adjustments for lcxt4 and ibmy, allow user 2d xy cross section output at z=nzb+1
Property svn:keywords set to `Id`
File size: 46.4 KB

Line
1	<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
2	<html><head>
3	<meta http-equiv="CONTENT-TYPE" content="text/html; charset=windows-1252"><title>PALM chapter 3.5.4</title> <meta name="GENERATOR" content="StarOffice 7 (Win32)"> <meta name="AUTHOR" content="Siegfried Raasch"> <meta name="CREATED" content="20040802;14050943"> <meta name="CHANGED" content="20041117;12180008"> <meta name="KEYWORDS" content="parallel LES model"> <style>
4	<!--
5	@page { size: 21cm 29.7cm }
6	-->
7	</style></head>
8	<body style="direction: ltr;" lang="en-US"><h3 style="line-height: 100%;"><font size="4">3.5.4
9	User-defined output quantities<br>
10	</font></h3>A very typical request of users is the
11	calculation and
12	output of
13	quantities which are not part of PALM's standard output. The basic user
14	interface includes a number of subroutines which allow the calculation
15	of user-defined quantities and output of these quantities as 1.
16	(horizontally averaged) <a href="#vertical_profiles">vertical
17	profiles</a>, 2. <a href="#time_series">time
18	series</a>, 3. <a href="#2D_and_3D">2d cross
19	section</a> or <a href="#2D_and_3D">3d volume data</a>,
20	4. <a href="#dvrp">dvrp</a> objects, and 5. <a href="#spectra">spectra</a>. The respective
21	subroutines
22	contain sample code lines (written as comment lines) for defining,
23	calculating and
24	output of such quantities.<br><br>Output times, averaging
25	intervals, etc. are steered by the same variables as used for the
26	standard PALM output quantities, e.g. <a href="chapter_4.2.html#dt_data_output">dt_data_output</a>.<br><br>The
27	following five parts of this chapter explains step-by-step how
28	to modify/extend the
29	respective default user interface subroutines in order to generate the respective
30	output:<br><ol><li>(horizontally averaged) <a href="chapter_3.5.4.html#vertical_profiles">vertical profiles</a>, </li><li><a href="chapter_3.5.4.html#time_series">time series</a>, </li><li><a href="chapter_3.5.4.html#2D_and_3D">2d cross
31	section</a> or <a href="chapter_3.5.4.html#2D_and_3D">3d
32	volume data</a>, </li><li><a href="chapter_3.5.4.html#dvrp">dvrp</a> objects,
33	and </li><li><a href="chapter_3.5.4.html#spectra">spectra</a>.</li></ol><br><h4><a name="vertical_profiles"></a>1. Output of user-defined
34	vertical profiles</h4>This example shows the output of the
35	quantity "turbulent resolved-scale horizontal momentum flux" (uv). If
36	more than one user-defined
37	quantity shall be output, the following steps have to be carried out in
38	the
39	same way for each of the quantities.<br><br><ol><li>The
40	quantity has to be given a unique string identifier, e.g. <span style="font-style: italic;">'uv'</span>.
41	This identifier must be different from the identifiers used for the
42	PALM standard output (see list in description of parameter <a href="chapter_4.2.html#data_output_pr">data_output_pr</a>).
43	To switch on output of this quantity, the user has to assign the string
44	identifier to the parameter <a href="chapter_4.3.html#data_output_pr_user">data_output_pr_user</a>,
45	eg.:<br><br><span style="font-family: monospace;">
46	data_output_pr_user</span> = <span style="font-style: italic;">'uv'</span>,<span style="font-style: italic;"></span><span style="font-style: italic;"></span><span style="font-style: italic;"></span><span style="font-style: italic;"></span><span style="font-style: italic;"></span><span style="font-style: italic;"></span><br><br></li><li>For
47	the
48	quantity, an identification number, a physical unit, and the vertical
49	grid on which it is defined (u- or w-grid), has to be assigned
50	(subroutine <a href="chapter_3.5.1.html#user_check_data_output_pr"><span style="font-family: monospace;">user_check_data_output_pr</span></a>):<span style="font-family: monospace;"><br><br>
51	CASE (
52	'uv' )<br></span><span style="font-family: monospace;">
53	index = pr_palm + <span style="color: rgb(255, 0, 0);">1<span style="color: rgb(0, 0, 0);">
54
55
56
57
58	! identification number</span></span><br>
59	dopr_index(var_count) = index<br>
60	dopr_unit(var_count)  = '<span style="color: rgb(255, 0, 0);">m2/s2</span>'
61
62
63	!
64	physical unit<br>
65	hom(:,2,index,:)      = SPREAD( <span style="color: rgb(255, 0, 0);">zu</span>, 2,
66	statistic_regions+1 )    ! vertical grid<br></span><br>Here
67	only the those parts in <span style="color: rgb(255, 0, 0);">red</span>
68	color have to be given by the user appropriately.<br>The
69	identification number (<span style="font-family: Courier New,Courier,monospace;">index</span>)
70	must be within the range [<span style="font-family: Courier New,Courier,monospace;">
71	pr_palm+1 , pr_palm+max_pr_user</span> ], where <span style="font-family: Courier New,Courier,monospace;">max_pr_user</span>
72	is the number of user-defined profiles as given by parameter <a href="chapter_4.3.html#data_output_pr_user">data_output_pr_user</a>
73	in the respective PALM run. The physical unit has to be given
74	with respect to the NetCDF conventions. If no unit is given,
75	PALM will abort. The vertical grid has to be either <span style="font-family: Courier New,Courier,monospace;">zu</span>
76	(u-grid) or <span style="font-family: Courier New,Courier,monospace;">zw</span>
77	(w-grid).<br><br></li><li>The quantity has to
78	be calculated for all gridpoints (subroutine <a href="chapter_3.5.1.html#user_statistics">user_statistics</a>):<br><br><span style="font-family: monospace;">    !$OMP
79	DO</span><br style="font-family: monospace;"><span style="font-family: monospace;">
80	DO  i = nxl, nxr</span><br style="font-family: monospace;"><span style="font-family: monospace;">
81	DO  j = nys, nyn</span><br style="font-family: monospace;"><span style="font-family: monospace;">
82	DO  k = nzb_s_inner(j,i)+1, nzt</span><br style="font-family: monospace;"><span style="font-family: monospace;"></span><span style="font-family: monospace;"></span><span style="font-family: monospace;"></span><span style="font-family: monospace;">
83	sums_l(k,pr_palm+<span style="color: rgb(255, 0, 0);">1</span>,tn)
84	= sums_l(k,pr_palm+<span style="color: rgb(255, 0, 0);">1</span>,tn)
85	+
86	&</span><br style="font-family: monospace;"><span style="font-family: monospace;">
87	<span style="color: rgb(255, 0, 0);">( 0.5 * ( u(k,j,i) +
88	u(k,j,i+1) ) - hom(k,1,1,sr) ) * &</span></span><br style="font-family: monospace; color: rgb(255, 0, 0);"><span style="font-family: monospace; color: rgb(255, 0, 0);">
89	( 0.5 * ( v(k,j,i) + v(k,j+1,i) ) - hom(k,1,2,sr) )   &</span><br style="font-family: monospace;"><span style="font-family: monospace;">
90	* rmask(j,i,sr)<br>
91	ENDDO<br>
92	ENDDO<br>    ENDDO<br></span><br style="font-family: monospace;">Once again, only those parts
93	in <span style="color: rgb(255, 0, 0);">red</span>
94	have to be adjusted by the user.<br>The
95	turbulent resolved-scale momentum flux uv is defined as the product
96	of the deviations of the horizontal velocities from their respective
97	horizontally averaged mean values. These mean values are stored in
98	array <span style="font-family: Courier New,Courier,monospace;">hom(..,1,1,sr)</span>
99	and <span style="font-family: Courier New,Courier,monospace;">hom(..,1,2,sr)</span>
100	for the u- and v-component, respectively. Since due to the staggered
101	grid, <span style="font-family: Courier New,Courier,monospace;">u</span>
102	and <span style="font-family: Courier New,Courier,monospace;">v</span>
103	are not defined at the same gridpoints, they have to be interpolated
104	appropriately (here to the center of the gridbox). The result of the
105	calculation is stored in array <span style="font-family: Courier New,Courier,monospace;">sums_l</span>.
106	The second index of this array is the identification number of the
107	profile which must match the one given in the previous step 2.</li></ol><br><h4><a name="time_series"></a>2. Output of user-defined
108	timeseries</h4>This example shows the output of two time series
109	for the absolut extremal values of the horizontal velocities <span style="font-style: italic;">u</span> and <span style="font-style: italic;">v</span>. If more than one
110	user-defined
111	quantity shall be output, the following steps have to be carried out in
112	the
113	same way for each of the quantities.<br><br><ol><li>For
114	each time series quantity you have to give a label and a unit
115	(subroutine <a href="chapter_3.5.1.html#user_init">user_init</a>),
116	which will be used for the NetCDF file. They must not contain more than
117	seven characters. The value of <code>dots_num</code> has
118	to be increased by the number of new time series quantities. Its old
119	value has to be stored in <code>dots_num_palm</code>.<br><br><code>
120	dots_label(dots_num<span style="color: rgb(255, 0, 0);">+1</span>)
121	= '<span style="color: rgb(255, 0, 0);">abs_umx</span>'<br>
122	dots_unit(dots_num<span style="color: rgb(255, 0, 0);">+1</span>)
123	= '<span style="color: rgb(255, 0, 0);">m/s</span>'<br>
124	dots_label(dots_num<span style="color: rgb(255, 0, 0);">+2</span>)
125	= '<span style="color: rgb(255, 0, 0);">abs_vmx</span>'<br>
126	dots_unit(dots_num<span style="color: rgb(255, 0, 0);">+2</span>)
127	= '<span style="color: rgb(255, 0, 0);">m/s</span>'<br><br>
128	dots_num_palm = dots_num<br>
129	dots_num = dots_num <span style="color: rgb(255, 0, 0);">+
130	2</span><br><br></code>Only those parts in <span style="color: rgb(255, 0, 0);">red</span> have to be
131	adjusted by the user.<br><br></li><li>These
132	quantities are calculated and output in subroutine <a href="chapter_3.5.1.html#user_statistics">user_statistics</a>
133	for every statistic region <code>sr</code> defined by the
134	user, but at least for the region <span style="font-style: italic;">"total domain"</span>:<br> <br><code>
135	ts_value(dots_num_palm+1,sr) = ABS( u_max )<br>
136	ts_value(dots_num_palm+2,sr) = ABS( v_max )<br></code></li></ol> <br><h4><span style="font-weight: bold;"><a name="2D_and_3D"></a>3.
137	Output of user-defined 2d cross sections or 3d volume data</span></h4><br>This
138	example shows the output of the
139	quantity "square of the u-component" (<span style="font-weight: bold;">Note:</span> this quantity
140	could of course easily be calculated from the u-component by
141	postprocessing the PALM output so that calculation within PALM is not
142	necessarily required). If more than one user-defined
143	quantity shall be output, the following steps have to be carried out in
144	the
145	same way for each of the quantities.<br><br><ol><li>The
146	quantity has to be given a unique string identifier, e.g. <span style="font-style: italic;">'u2'</span>.
147	This identifier must be different from the identifiers used for the
148	PALM standard output (see list in description of parameter <a href="chapter_4.2.html#data_output">data_output</a>).
149	To switch on output of this quantity, the user has to assign the string
150	identifier to the parameter <a href="chapter_4.3.html#data_output_user">data_output_user</a>,
151	eg.:<br><br><span style="font-family: monospace;">
152	data_output_user</span> = <span style="font-style: italic;">'u2'</span>, <span style="font-style: italic;"> 'u2_xy_av'</span><br><br>The
153	pure string <span style="font-style: italic;">'u2'</span>
154	switches on the output of instantaneous 3d volume data. Output of cross
155	section data and time averaged data is switched on by additionally
156	appending the strings <span style="font-style: italic;">'_xy'</span>,
157	<span style="font-style: italic;">'_xz'</span>, <span style="font-style: italic;">'_yz'</span>, and/or <span style="font-style: italic;">'_av'</span> (for a
158	detailed explanation see parameter <a href="chapter_4.2.html#data_output">data_output</a>).<br><br></li><li>In
159	order to store the quantities' grid point data within PALM, a 3d data
160	array has to be declared in module <a href="chapter_3.5.1.html#user"><span style="font-family: monospace;">user</span></a>:<br><br><span style="font-family: monospace;">    REAL,
161	DIMENSION(:,:,:), ALLOCATABLE ::  u2, u2_av</span><br><br>The
162	second array <span style="font-family: monospace;">u2_av</span>
163	is needed in case that output of time averaged data is requested. It is
164	used to store the sum of the data of the respective time levels over
165	which the average has to be carried out.<br><br><br></li><li>The
166	data array has to be allocated in subroutine <a href="chapter_3.5.1.html#user_init"><span style="font-family: monospace;">user_init</span></a>:<br><br><span style="font-family: monospace;">
167	ALLOCATE( u2(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )</span><br><br>In
168	case that output of time averaged data is requested, the array
169	containing the sum has possibly to be read from the restart file (local
170	filename <a href="chapter_3.4.html#BININ">BININ</a>)
171	by executing the following code in <span style="font-family: monospace;">user_init</span>:<br><br><span style="font-family: monospace;">
172	IF ( initializing_actions == 'read_restart_data' )
173	THEN</span><br style="font-family: monospace;"><span style="font-family: monospace;">
174	READ ( 13 )  field_chr</span><br style="font-family: monospace;"><span style="font-family: monospace;">
175	DO  WHILE ( TRIM( field_chr ) /= '*** end
176	user ***' )</span><br style="font-family: monospace;"><br style="font-family: monospace;"><span style="font-family: monospace;">
177	SELECT CASE ( TRIM( field_chr ) )</span><br style="font-family: monospace;"><br style="font-family: monospace;"><span style="font-family: monospace;">
178	CASE ( 'u2_av' )</span><br style="font-family: monospace;"><span style="font-family: monospace;">
179	ALLOCATE(
180	u2_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )</span><br style="font-family: monospace;"><span style="font-family: monospace;">
181	READ ( 13 )  u2_av</span><br style="font-family: monospace;"><br style="font-family: monospace;"><span style="font-family: monospace;">
182	CASE DEFAULT</span><br style="font-family: monospace;"><span style="font-family: monospace;">
183	PRINT*,
184	'+++ user_init: unknown
185	variable named "', &</span><br style="font-family: monospace;"><span style="font-family: monospace;">
186
187
188	TRIM(
189	field_chr ), '" found in'</span><br style="font-family: monospace;"><span style="font-family: monospace;">
190	PRINT*,
191	'
192	data from prior run on PE ', myid</span><br style="font-family: monospace;"><span style="font-family: monospace;">
193
194	CALL local_stop</span><br style="font-family: monospace;"><br style="font-family: monospace;"><span style="font-family: monospace;">
195	END SELECT<br>
196	ENDDO</span><br style="font-family: monospace;"><span style="font-family: monospace;">
197	ENDIF</span><br style="font-family: monospace;"><br><br></li><li>The
198	quantity has to be given a unit (subroutine <a href="chapter_3.5.1.html#user_check_data_output"><span style="font-family: monospace;">user_check_data_output</span></a>):<br><br><span style="font-family: monospace;">    CASE (
199	'u2' )</span><br style="font-family: monospace;"><span style="font-family: monospace;">
200	unit = 'm2/s2'</span><br> <br>Otherwise,
201	PALM will abort.<br><br><br></li><li>The
202	vertical grid on which the quantity is defined (given by the levels
203	'zu' or 'zw', on which the u- or w-component of the velocity are
204	defined) has to be specified for the NetCDF output files in subroutine <a href="chapter_3.5.1.html#user_define_netcdf_grid"><span style="font-family: monospace;">user_define_netcdf_grid</span></a>:<br><br><span style="font-family: monospace;">
205	CASE ( 'u2', 'u2_xy', 'u2_xz', 'u2_yz' )</span><br style="font-family: monospace;"><span style="font-family: monospace;">
206	grid = 'zu'</span><br> <br>As
207	the example shows, this grid has to be defined for the 3d volume data
208	as well as for all of the three cross sections.<br><span style="font-family: monospace;"><br><br></span></li><li>After
209	each timestep, the quantity has to be calculated at all gridpoints and
210	to be stored. This has to be done in subroutine <a href="chapter_3.5.1#user_actions"><span style="font-family: monospace;">user_actions</span></a>
211	at location 'after_integration':<br><br><span style="font-family: monospace;">    CASE
212	( 'after_integration' )</span><br style="font-family: monospace;"><span style="font-family: monospace;">!</span><br style="font-family: monospace;"><span style="font-family: monospace;">!--
213	Enter actions to be done after every time integration (before</span><br style="font-family: monospace;"><span style="font-family: monospace;">!--
214	data output)</span><br style="font-family: monospace;"><span style="font-family: monospace;">!--
215	Sample for user-defined output:</span><br style="font-family: monospace;"><span style="font-family: monospace;">
216	DO  i = nxl-1, nxr+1</span><br style="font-family: monospace;"><span style="font-family: monospace;">
217	DO  j = nys-1, nyn+1</span><br style="font-family: monospace;"><span style="font-family: monospace;">
218	DO  k = nzb, nzt+1</span><br style="font-family: monospace;"><span style="font-family: monospace;">
219	u2(k,j,i) = u(k,j,i)**2</span><br style="font-family: monospace;"><span style="font-family: monospace;">
220	ENDDO</span><br style="font-family: monospace;"><span style="font-family: monospace;">
221	ENDDO</span><br style="font-family: monospace;"><span style="font-family: monospace;">
222	ENDDO</span><br><br> <br></li><li>In
223	case that output of time-averaged data is requested, the sum- and
224	average-operations as well as the allocation of the sum-array have to
225	be carried out in subroutine <a href="chapter_3.5.1.html#user_3d_data_averaging"><span style="font-family: monospace;">user_3d_data_averaging</span></a>:<br><br><span style="font-family: monospace;">
226	IF ( mode == 'allocate' )  THEN<br>
227	...</span><br style="font-family: monospace;"><span style="font-family: monospace;">
228	CASE ( 'u2' )</span><br style="font-family: monospace;"><span style="font-family: monospace;">
229	IF ( .NOT. ALLOCATED( u2_av ) )  THEN</span><br style="font-family: monospace;"><span style="font-family: monospace;">
230	ALLOCATE( u2_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )</span><br style="font-family: monospace;"><span style="font-family: monospace;">
231	ENDIF</span><br style="font-family: monospace;"><span style="font-family: monospace;">
232	u2_av = 0.0<br>    ...</span><br style="font-family: monospace;"><span style="font-family: monospace;">
233	ELSEIF ( mode == 'sum' )  THEN<br>
234	...</span><br style="font-family: monospace;"><span style="font-family: monospace;">       CASE
235	( 'u2' )</span><br style="font-family: monospace;"><span style="font-family: monospace;">
236	DO  i = nxl-1, nxr+1</span><br style="font-family: monospace;"><span style="font-family: monospace;">
237	DO  j = nys-1, nyn+1</span><br style="font-family: monospace;"><span style="font-family: monospace;">
238	DO  k = nzb, nzt+1</span><br style="font-family: monospace;"><span style="font-family: monospace;">
239	u2_av(k,j,i) = u2_av(k,j,i) + u2(k,j,i)</span><br style="font-family: monospace;"><span style="font-family: monospace;">
240	ENDDO</span><br style="font-family: monospace;"><span style="font-family: monospace;">
241	ENDDO</span><br style="font-family: monospace;"><span style="font-family: monospace;">
242	ENDDO<br>    ...</span><br style="font-family: monospace;"><span style="font-family: monospace;">
243	ELSEIF ( mode == 'average' )  THEN<br>
244	...</span><br style="font-family: monospace;"><span style="font-family: monospace;">       CASE
245	( 'u2' )</span><br style="font-family: monospace;"><span style="font-family: monospace;">
246	DO  i = nxl-1, nxr+1</span><br style="font-family: monospace;"><span style="font-family: monospace;">
247	DO  j = nys-1, nyn+1</span><br style="font-family: monospace;"><span style="font-family: monospace;">
248	DO  k = nzb, nzt+1</span><br style="font-family: monospace;"><span style="font-family: monospace;">
249	u2_av(k,j,i) = u2_av(k,j,i) / REAL( average_count_3d )</span><br style="font-family: monospace;"><span style="font-family: monospace;">
250	ENDDO</span><br style="font-family: monospace;"><span style="font-family: monospace;">
251	ENDDO</span><br style="font-family: monospace;"><span style="font-family: monospace;">
252	ENDDO</span><br><br> </li><li>For
253	output of 2d cross sections, the gridpoint data of the quantity has to
254	be resorted to array <span style="font-family: monospace;">local_pf</span>
255	in subroutine <a href="chapter_3.5.1.html#user_data_output_2d"><span style="font-family: monospace;">user_data_output_2d</span></a>.
256	Also the vertical grid, on which the quantity is defined, has to be set
257	again:<br><br><span style="font-family: monospace;">    CASE
258	( 'u2_xy', 'u2_xz', 'u2_yz' )</span><br style="font-family: monospace;"><span style="font-family: monospace;">
259	IF ( av == 0 )  THEN</span><br style="font-family: monospace;"><span style="font-family: monospace;">
260	DO  i = nxl-1, nxr+1</span><br style="font-family: monospace;"><span style="font-family: monospace;">
261	DO  j = nys-1, nyn+1</span><br style="font-family: monospace;"><span style="font-family: monospace;">
262	DO  k = nzb, nzt+1</span><br style="font-family: monospace;"><span style="font-family: monospace;">
263	local_pf(i,j,k) = u2(k,j,i)</span><br style="font-family: monospace;"><span style="font-family: monospace;">
264	ENDDO</span><br style="font-family: monospace;"><span style="font-family: monospace;">
265	ENDDO</span><br style="font-family: monospace;"><span style="font-family: monospace;">
266	ENDDO</span><br style="font-family: monospace;"><span style="font-family: monospace;">
267	ELSE</span><br style="font-family: monospace;"><span style="font-family: monospace;">
268	DO  i = nxl-1, nxr+1</span><br style="font-family: monospace;"><span style="font-family: monospace;">
269	DO  j = nys-1, nyn+1</span><br style="font-family: monospace;"><span style="font-family: monospace;">
270	DO  k = nzb, nzt+1</span><br style="font-family: monospace;"><span style="font-family: monospace;">
271	local_pf(i,j,k) = u2_av(k,j,i)</span><br style="font-family: monospace;"><span style="font-family: monospace;">
272	ENDDO</span><br style="font-family: monospace;"><span style="font-family: monospace;">
273	ENDDO</span><br style="font-family: monospace;"><span style="font-family: monospace;">
274	ENDDO</span><br style="font-family: monospace;"><span style="font-family: monospace;">
275	ENDIF</span><br style="font-family: monospace;"><br style="font-family: monospace;"><span style="font-family: monospace;">
276	grid = 'zu'</span><br> <br>The <span style="font-family: monospace;">ELSE</span> case is
277	only needed in case that output of time-averaged data is requested.<br><br>
278	As a special case, xy cross section output can also be defined at one single level at height <span style="font-family: monospace;">k=nzb+1</span> on the u-grid. This features is useful for output of surface data (e.g. heat fluxes). In this case, the corresponding 2d data has to be resorted to the array <span style="font-family: monospace;">local_pf(i,j,nzb+1).</span> In addition to this, the grid in <a href="chapter_3.5.1.html#user_define_netcdf_grid"><span style="font-family: monospace;">user_define_netcdf_grid</span></a> as well as in <a href="chapter_3.5.1.html#user_data_output_2d"><span style="font-family: monospace;">user_data_output_2d</span></a> must be set to <span style="font-family: monospace;">grid = 'zu1'</span>.<br><br>
279	<span style="font-family: monospace;">    CASE
280	( 'u2_xy' )</span><br style="font-family: monospace;"><span style="font-family: monospace;">
281	IF ( av == 0 )  THEN</span><br style="font-family: monospace;"><span style="font-family: monospace;">
282	DO  i = nxl-1, nxr+1</span><br style="font-family: monospace;"><span style="font-family: monospace;">
283	DO  j = nys-1, nyn+1</span><br style="font-family: monospace;"><span style="font-family: monospace;">
284	local_pf(i,j,nzb+1) = u2(j,i)</span><br style="font-family: monospace;"><span style="font-family: monospace;">
285	ENDDO</span><br style="font-family: monospace;"><span style="font-family: monospace;">
286	ENDDO</span><br style="font-family: monospace;"><span style="font-family: monospace;">
287	ELSE</span><br style="font-family: monospace;"><span style="font-family: monospace;">
288	DO  i = nxl-1, nxr+1</span><br style="font-family: monospace;"><span style="font-family: monospace;">
289	DO  j = nys-1, nyn+1</span><br style="font-family: monospace;"><span style="font-family: monospace;">
290	local_pf(i,j,nzb+1) = u2_av(j,i)</span><br style="font-family: monospace;"><span style="font-family: monospace;">
291	ENDDO</span><br style="font-family: monospace;"><span style="font-family: monospace;">
292	ENDDO</span><br style="font-family: monospace;"><span style="font-family: monospace;">
293	ENDIF</span><br style="font-family: monospace;"><br style="font-family: monospace;"><span style="font-family: monospace;">
294	grid = 'zu1'</span><br><span style="font-family: monospace;">
295	two_d = .TRUE.</span><br><br>
296	Note that <span style="font-family: monospace;">two_d = .TRUE.</span> is necessary for output of a 2d data slice.<br><br>
297	<br></li><li>For
298	output of 3d volume data, the gridpoint data of the quantity has to be
299	resorted to array <span style="font-family: monospace;">local_pf</span>
300	in subroutine <a href="chapter_3.5.1.html#user_data_output_3d"><span style="font-family: monospace;">user_data_output_3d</span></a>.:<br><br><span style="font-family: monospace;">    CASE
301	( 'u2' )</span><br style="font-family: monospace;"><span style="font-family: monospace;">
302	IF ( av == 0 )  THEN</span><br style="font-family: monospace;"><span style="font-family: monospace;">
303	DO  i = nxl-1, nxr+1</span><br style="font-family: monospace;"><span style="font-family: monospace;">
304	DO  j = nys-1, nyn+1</span><br style="font-family: monospace;"><span style="font-family: monospace;">
305	DO  k = nzb, nz_do</span><br style="font-family: monospace;"><span style="font-family: monospace;">
306	local_pf(i,j,k) = u2(k,j,i)</span><br style="font-family: monospace;"><span style="font-family: monospace;">
307	ENDDO</span><br style="font-family: monospace;"><span style="font-family: monospace;">
308	ENDDO</span><br style="font-family: monospace;"><span style="font-family: monospace;">
309	ENDDO</span><br style="font-family: monospace;"><span style="font-family: monospace;">
310	ELSE</span><br style="font-family: monospace;"><span style="font-family: monospace;">
311	DO  i = nxl-1, nxr+1</span><br style="font-family: monospace;"><span style="font-family: monospace;">
312	DO  j = nys-1, nyn+1</span><br style="font-family: monospace;"><span style="font-family: monospace;">
313	DO  k = nzb, nz_do</span><br style="font-family: monospace;"><span style="font-family: monospace;">
314	local_pf(i,j,k) = u2_av(k,j,i)</span><br style="font-family: monospace;"><span style="font-family: monospace;">
315	ENDDO</span><br style="font-family: monospace;"><span style="font-family: monospace;">
316	ENDDO</span><br style="font-family: monospace;"><span style="font-family: monospace;">
317	ENDDO</span><br style="font-family: monospace;"><span style="font-family: monospace;">
318	ENDIF</span><br><br>The <span style="font-family: monospace;">ELSE</span> case is
319	only needed in case that output of time-averaged data is requested.<br><br><br></li><li>In
320	case of job chains, the sum array has to be written to the (binary)
321	restart file (local filename <a href="chapter_3.4.html#BINOUT">BINOUT</a>)
322	in subroutine <a href="chapter_3.5.1.html#user_last_actions"><span style="font-family: monospace;">user_last_actions</span></a>:<br><br><span style="font-family: monospace;">
323	IF ( ALLOCATED( u2_av ) )  THEN</span><br style="font-family: monospace;"><span style="font-family: monospace;">
324	WRITE ( 14 )
325	'u2_av
326	';  WRITE ( 14 )  u2_av</span><br style="font-family: monospace;"><span style="font-family: monospace;">
327	ENDIF</span><br><br>Otherwise,
328	the calculated
329	time-average may be wrong. In the restart run, this quantity has to be
330	read from the restart file by including the following code in
331	subroutine <a href="chapter_3.5.1.html#user_read_restart_data"><span style="font-family: monospace;">user_read_restart_data</span></a>:<br><br><span style="font-family: monospace;">    IF ( initializing_actions == 'read_restart_data' )  THEN<br>       READ ( 13 )  field_char<br>       DO  WHILE ( TRIM( field_char ) /= '* end user *' )<br><br>          SELECT CASE ( TRIM( field_char ) )<br><br>             CASE ( 'u2_av' )<br>                IF ( .NOT. ALLOCATED( u2_av ) ) THEN<br>
332	ALLOCATE( u2_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )<br>                ENDIF<br>                READ ( 13 )  tmp_3d<br>
333	u2_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &<br>
334	tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)<br><br>             CASE DEFAULT<br>
335	PRINT*, '+++ user_init: unknown variable named "', &<br>
336	TRIM( field_char ), '" found in'<br>
337	PRINT*,
338	'
339	data from prior run on PE ', myid<br>                CALL local_stop<br><br>          END SELECT<br><br>          READ ( 13 )  field_char<br><br>       ENDDO<br>    ENDIF<br></span></li></ol><h4><a name="dvrp"></a>4. Output of user-defined DVRP objects</h4>This
340	example shows the output of the
341	quantity "square of the u-component", <span style="font-style: italic;">u<sup>2</sup></span>.
342	If more than one user-defined
343	quantity shall be output, the following steps have to be carried out in
344	the
345	same way for each of the quantities. First, steps 1 - 6
346	of <a href="chapter_3.5.4.html#2D_and_3D"><span style="font-weight: bold;">part 3.</span> </a> <a href="chapter_3.5.4.html#2D_and_3D">(2d cross
347	section or </a><a href="chapter_3.5.4.html#2D_and_3D">3d
348	volume data)</a> are required. Second, the gridpoint
349	data
350	of the quantity has to
351	be resorted to array <span style="font-family: monospace;">local_pf</span>
352	in subroutine <a href="chapter_3.5.1.html#user_data_output_dvrp"><span style="font-family: monospace;">user_data_output_dvrp</span></a>
353	as follows:<br><br><code>
354	CASE ( '<span style="color: rgb(255, 0, 0);">u2</span>',
355	'<span style="color: rgb(255, 0, 0);">u2</span>_xy',
356	'<span style="color: rgb(255, 0, 0);">u2</span>_xz',
357	'<span style="color: rgb(255, 0, 0);">u2</span>_yz'
358	)<br>
359	DO  i = nxl, nxr+1<br>
360	DO  j = nys, nyn+1<br>
361	DO  k = nzb, nz_do3d<br>
362	local_pf(i,j,k) = <span style="color: rgb(255, 0, 0);">u2</span>(k,j,i)<br>
363	ENDDO<br>
364	ENDDO<br>
365	ENDDO</code><br><br>Only those parts in <span style="color: rgb(255, 0, 0);">red</span> have to be
366	adjusted by the user.<br><br>After performing these steps,
367	the user-defined quantity <code>'u2'</code> can be
368	selected like standard model quantities by the <a href="chapter_4.2.html#dvrp_graphics">dvrp_graphics</a>
369	package parameter <a href="chapter_4.2.html#mode_dvrp">mode_dvrp</a>.<br><br><h4><a name="spectra"></a>5. Output of user-defined spectra</h4>This
370	example shows the output of the
371	quantity "turbulent resolved-scale horizontal momentum flux" (uv). If
372	more than one user-defined
373	quantity shall be output, the following steps have to be carried out in
374	the
375	same way for each of the quantities.<br><br><ol><li>The
376	calculation of user-defined spectra is closely linked with the
377	calculation of <a href="chapter_3.5.4.html#vertical_profiles"><span style="font-weight: bold;">part 1.</span>
378	(user-defined
379	</a><a href="chapter_3.5.4.html#vertical_profiles">vertical
380	profiles</a><a href="chapter_3.5.4.html#vertical_profiles">)</a><span style="font-weight: bold;"></span> and <a href="chapter_3.5.4.html#2D_and_3D"><span style="font-weight: bold;">part 3.</span> (user-defined</a><a href="chapter_3.5.4.html#2D_and_3D"> </a><a href="chapter_3.5.4.html#2D_and_3D">3d volume data</a><a href="chapter_3.5.4.html#2D_and_3D">)</a><span style="font-weight: bold;"></span>. Therefore, the
381	following prerequisites apply
382	for each user-defined
383	spectra quantity:</li></ol><ol style="list-style-type: upper-alpha;"><ol style="list-style-type: upper-alpha;"><li>From <a href="chapter_3.5.4.html#vertical_profiles"><span style="font-weight: bold;">part 1.</span>
384	(user-defined
385	</a><a href="chapter_3.5.4.html#vertical_profiles">vertical
386	profiles</a><a href="chapter_3.5.4.html#vertical_profiles">)</a>
387	steps 2 and 3. <br>See the sample code (as comment lines) for <span style="font-style: italic;">'uv'</span> and <code>ustvst</code>,
388	respectively. (Actual output of vertical profiles - step 1 - is not
389	required.)</li><li>From <a href="chapter_3.5.4.html#2D_and_3D"><span style="font-weight: bold;">part 3.</span> (user-defined</a><a href="chapter_3.5.4.html#2D_and_3D"> </a><a href="chapter_3.5.4.html#2D_and_3D">3d volume data</a><a href="chapter_3.5.4.html#2D_and_3D">)</a> steps 2, 3,
390	4, 5, and 6. <br>See the sample code (as comment lines) for <span style="font-style: italic;">'uv'</span> and <code>ustvst</code>,
391	respectively. (Actual output of 3d volume data - step 1 - is not
392	required.)</li><li>The quantity has to be given a unique
393	string identifier, e.g. <span style="font-style: italic;">'uv'</span>.
394	This identifier must be different from the identifiers used for the
395	PALM standard output (see list in description of package parameter <a href="chapter_4.2.html#data_output_sp">data_output_sp</a>).
396	To switch on output of this quantity, the user has to assign the string
397	identifier to the package parameter <a href="chapter_4.2.html#data_output_sp">data_output_sp</a>,
398	eg.:<br><br><span style="font-family: monospace;">
399	data_output_sp</span> =  <span style="font-style: italic;">'uv'</span><br><br></li></ol>A.
400	and B. as prerequisites for C. require a naming convention of identical
401	identifiers, e.g. data_output_pr_user = 'uv', data_output_user =
402	'uv' and data_output_sp = 'uv'. This naming convention applies only
403	in case of user-defined spectra.</ol><ol start="2"><li>Edit
404	the subroutine <a href="chapter_3.5.1.html#user_spectra">user_spectra</a>
405	(only those parts in <span style="color: rgb(255, 0, 0);">red</span>)
406	as follows:<br><br><code>
407	IF ( mode == 'preprocess' )  THEN<br><br>
408	SELECT CASE ( TRIM( data_output_sp(m) ) )<br>
409	<br>
410	CASE ( 'u', 'v', 'w', 'pt', 'q' )<br>!--
411	Not allowed here since these are the standard quantities used in <br>!--
412	preprocess_spectra.<br>
413	<br><span style="color: rgb(255, 0, 0); background-color: rgb(255, 255, 255);">
414	CASE ( 'uv' )</span><br style="color: rgb(255, 0, 0); background-color: rgb(255, 255, 255);"><span style="color: rgb(255, 0, 0); background-color: rgb(255, 255, 255);">
415	pr = pr_palm+1</span><br style="color: rgb(255, 0, 0); background-color: rgb(255, 255, 255);"><span style="color: rgb(255, 0, 0); background-color: rgb(255, 255, 255);">
416	d(nzb+1:nzt,nys:nyn,nxl:nxr) = ustvst(nzb+1:nzt,nys:nyn,nxl:nxr)</span><br style="color: rgb(255, 0, 0); background-color: rgb(255, 255, 255);">
417	<br>
418	CASE DEFAULT<br>
419	PRINT*, '+++ user_spectra/preprocess: Spectra of ', &<br>
420	TRIM( data_output_sp(m) ), ' can not be calculated'<br>
421	<br>
422	END SELECT<br><br>    ELSEIF (
423	mode == 'data_output' )  THEN<br><br>
424	SELECT CASE ( TRIM( data_output_sp(m) ) )<br><br>
425	CASE ( 'u', 'v', 'w', 'pt', 'q' )<br>!--
426	Not allowed here since these are the standard quantities used in <br>!--
427	data_output_spectra.<br><br><span style="color: rgb(255, 0, 0);">
428	CASE ( 'uv' )</span><br style="color: rgb(255, 0, 0);"><span style="color: rgb(255, 0, 0);">
429	pr = 6</span><br style="color: rgb(255, 0, 0);"><br>
430	CASE DEFAULT<br>
431	PRINT*, '+++ user_spectra/data_output: Spectra of ', &<br>
432	TRIM( data_output_sp(m) ), ' are not defined'<br><br>
433	END SELECT<br><br>    ENDIF</code><br></li></ol>Note
434	that spectra output is an <a href="chapter_3.7.html">optional
435	software package</a> (see <a href="chapter_3.7.html">chapter
436	3.7</a>). Therefore user-defined spectra also require the package
437	activation via <code>mrun -p spectra</code>.<hr>
438	<p style="line-height: 100%;"><br><font color="#000080"><font color="#000080"><a href="chapter_3.5.3.html"><font color="#000080"><img src="left.gif" name="Grafik1" align="bottom" border="2" height="32" width="32"></font></a><a href="index.html"><font color="#000080"><img src="up.gif" name="Grafik2" align="bottom" border="2" height="32" width="32"></font></a><a href="chapter_3.5.5.html"><font color="#000080"><img style="border: 2px solid ; width: 32px; height: 32px;" alt="" src="right.gif" name="Grafik3"></font></a></font></font></p><p style="line-height: 100%;"><i>Last
439	change: </i> $Id: chapter_3.5.4.html 136 2007-11-26
440	02:47:32Z letzel $</p>
441	</body></html>

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |