Subversion Repositories lagranto.wrf

Rev

Go to most recent revision | Details | Last modification | View Log | RSS feed

Rev Author Line No. Line
2 michaesp 1
%!PS-Adobe-3.0
2
%%Creator: groff version 1.21
3
%%CreationDate: Tue Apr  2 16:15:02 2013
4
%%DocumentNeededResources: font Times-Roman
5
%%+ font Times-Bold
6
%%+ font Times-Italic
7
%%DocumentSuppliedResources: procset grops 1.21 0
8
%%Pages: 26
9
%%PageOrder: Ascend
10
%%DocumentMedia: Default 612 792 0 () ()
11
%%Orientation: Portrait
12
%%EndComments
13
%%BeginDefaults
14
%%PageMedia: Default
15
%%EndDefaults
16
%%BeginProlog
17
%%BeginResource: procset grops 1.21 0
18
%!PS-Adobe-3.0 Resource-ProcSet
19
/setpacking where{
20
pop
21
currentpacking
22
true setpacking
23
}if
24
/grops 120 dict dup begin
25
/SC 32 def
26
/A/show load def
27
/B{0 SC 3 -1 roll widthshow}bind def
28
/C{0 exch ashow}bind def
29
/D{0 exch 0 SC 5 2 roll awidthshow}bind def
30
/E{0 rmoveto show}bind def
31
/F{0 rmoveto 0 SC 3 -1 roll widthshow}bind def
32
/G{0 rmoveto 0 exch ashow}bind def
33
/H{0 rmoveto 0 exch 0 SC 5 2 roll awidthshow}bind def
34
/I{0 exch rmoveto show}bind def
35
/J{0 exch rmoveto 0 SC 3 -1 roll widthshow}bind def
36
/K{0 exch rmoveto 0 exch ashow}bind def
37
/L{0 exch rmoveto 0 exch 0 SC 5 2 roll awidthshow}bind def
38
/M{rmoveto show}bind def
39
/N{rmoveto 0 SC 3 -1 roll widthshow}bind def
40
/O{rmoveto 0 exch ashow}bind def
41
/P{rmoveto 0 exch 0 SC 5 2 roll awidthshow}bind def
42
/Q{moveto show}bind def
43
/R{moveto 0 SC 3 -1 roll widthshow}bind def
44
/S{moveto 0 exch ashow}bind def
45
/T{moveto 0 exch 0 SC 5 2 roll awidthshow}bind def
46
/SF{
47
findfont exch
48
[exch dup 0 exch 0 exch neg 0 0]makefont
49
dup setfont
50
[exch/setfont cvx]cvx bind def
51
}bind def
52
/MF{
53
findfont
54
[5 2 roll
55
 
56
neg 0 0]makefont
57
dup setfont
58
[exch/setfont cvx]cvx bind def
59
}bind def
60
/level0 0 def
61
/RES 0 def
62
/PL 0 def
63
/LS 0 def
64
/MANUAL{
65
statusdict begin/manualfeed true store end
66
}bind def
67
/PLG{
68
gsave newpath clippath pathbbox grestore
69
exch pop add exch pop
70
}bind def
71
/BP{
72
/level0 save def
73
1 setlinecap
74
1 setlinejoin
75
DEFS/BPhook known{DEFS begin BPhook end}if
76
72 RES div dup scale
77
LS{
78
90 rotate
79
}{
80
 
81
}ifelse
82
1 -1 scale
83
}bind def
84
/EP{
85
level0 restore
86
showpage
87
}def
88
/DA{
89
newpath arcn stroke
90
}bind def
91
/SN{
92
transform
93
.25 sub exch .25 sub exch
94
round .25 add exch round .25 add exch
95
itransform
96
}bind def
97
/DL{
98
SN
99
moveto
100
SN
101
lineto stroke
102
}bind def
103
/DC{
104
newpath 0 360 arc closepath
105
}bind def
106
/TM matrix def
107
/DE{
108
TM currentmatrix pop
109
translate scale newpath 0 0 .5 0 360 arc closepath
110
TM setmatrix
111
}bind def
112
/RC/rcurveto load def
113
/RL/rlineto load def
114
/ST/stroke load def
115
/MT/moveto load def
116
/CL/closepath load def
117
/Fr{
118
setrgbcolor fill
119
}bind def
120
/setcmykcolor where{
121
pop
122
/Fk{
123
setcmykcolor fill
124
}bind def
125
}if
126
/Fg{
127
setgray fill
128
}bind def
129
/FL/fill load def
130
/LW/setlinewidth load def
131
/Cr/setrgbcolor load def
132
/setcmykcolor where{
133
pop
134
/Ck/setcmykcolor load def
135
}if
136
/Cg/setgray load def
137
/RE{
138
findfont
139
dup maxlength 1 index/FontName known not{1 add}if dict begin
140
{
141
1 index/FID ne
142
2 index/UniqueID ne
143
and
144
{def}{pop pop}ifelse
145
}forall
146
/Encoding exch def
147
dup/FontName exch def
148
currentdict end definefont pop
149
}bind def
150
/DEFS 0 def
151
/EBEGIN{
152
moveto
153
DEFS begin
154
}bind def
155
/EEND/end load def
156
/CNT 0 def
157
/level1 0 def
158
/PBEGIN{
159
/level1 save def
160
translate
161
div 3 1 roll div exch scale
162
neg exch neg exch translate
163
 
164
 
165
1 setlinewidth
166
 
167
10 setmiterlimit
168
[]0 setdash
169
/setstrokeadjust where{
170
pop
171
false setstrokeadjust
172
}if
173
/setoverprint where{
174
pop
175
false setoverprint
176
}if
177
newpath
178
/CNT countdictstack def
179
userdict begin
180
/showpage{}def
181
/setpagedevice{}def
182
mark
183
}bind def
184
/PEND{
185
cleartomark
186
countdictstack CNT sub{end}repeat
187
level1 restore
188
}bind def
189
end def
190
/setpacking where{
191
pop
192
setpacking
193
}if
194
%%EndResource
195
%%EndProlog
196
%%BeginSetup
197
%%BeginFeature: *PageSize Default
198
<< /PageSize [ 612 792 ] /ImagingBBox null >> setpagedevice
199
%%EndFeature
200
%%IncludeResource: font Times-Roman
201
%%IncludeResource: font Times-Bold
202
%%IncludeResource: font Times-Italic
203
grops begin/DEFS 1 dict def DEFS begin/u{.001 mul}bind def end/RES 72
204
def/PL 792 def/LS false def/ENC0[/asciicircum/asciitilde/Scaron/Zcaron
205
/scaron/zcaron/Ydieresis/trademark/quotesingle/Euro/.notdef/.notdef
206
/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef
207
/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef/.notdef
208
/.notdef/.notdef/space/exclam/quotedbl/numbersign/dollar/percent
209
/ampersand/quoteright/parenleft/parenright/asterisk/plus/comma/hyphen
210
/period/slash/zero/one/two/three/four/five/six/seven/eight/nine/colon
211
/semicolon/less/equal/greater/question/at/A/B/C/D/E/F/G/H/I/J/K/L/M/N/O
212
/P/Q/R/S/T/U/V/W/X/Y/Z/bracketleft/backslash/bracketright/circumflex
213
/underscore/quoteleft/a/b/c/d/e/f/g/h/i/j/k/l/m/n/o/p/q/r/s/t/u/v/w/x/y
214
/z/braceleft/bar/braceright/tilde/.notdef/quotesinglbase/guillemotleft
215
/guillemotright/bullet/florin/fraction/perthousand/dagger/daggerdbl
216
/endash/emdash/ff/fi/fl/ffi/ffl/dotlessi/dotlessj/grave/hungarumlaut
217
/dotaccent/breve/caron/ring/ogonek/quotedblleft/quotedblright/oe/lslash
218
/quotedblbase/OE/Lslash/.notdef/exclamdown/cent/sterling/currency/yen
219
/brokenbar/section/dieresis/copyright/ordfeminine/guilsinglleft
220
/logicalnot/minus/registered/macron/degree/plusminus/twosuperior
221
/threesuperior/acute/mu/paragraph/periodcentered/cedilla/onesuperior
222
/ordmasculine/guilsinglright/onequarter/onehalf/threequarters
223
/questiondown/Agrave/Aacute/Acircumflex/Atilde/Adieresis/Aring/AE
224
/Ccedilla/Egrave/Eacute/Ecircumflex/Edieresis/Igrave/Iacute/Icircumflex
225
/Idieresis/Eth/Ntilde/Ograve/Oacute/Ocircumflex/Otilde/Odieresis
226
/multiply/Oslash/Ugrave/Uacute/Ucircumflex/Udieresis/Yacute/Thorn
227
/germandbls/agrave/aacute/acircumflex/atilde/adieresis/aring/ae/ccedilla
228
/egrave/eacute/ecircumflex/edieresis/igrave/iacute/icircumflex/idieresis
229
/eth/ntilde/ograve/oacute/ocircumflex/otilde/odieresis/divide/oslash
230
/ugrave/uacute/ucircumflex/udieresis/yacute/thorn/ydieresis]def
231
/Times-Italic@0 ENC0/Times-Italic RE/Times-Bold@0 ENC0/Times-Bold RE
232
/Times-Roman@0 ENC0/Times-Roman RE
233
%%EndSetup
234
%%Page: 1 1
235
%%BeginPageSetup
236
BP
237
%%EndPageSetup
238
/F0 10/Times-Roman@0 SF 407.76(caltra\(\) caltra\(\))72 48 R/F1 10.95
239
/Times-Bold@0 SF -.219(NA)72 84 S(ME).219 E/F2 10/Times-Bold@0 SF
240
(caltra -)108 96 Q F0(calculate air parcel trajectories)2.5 E F1
241
(SYNOPSIS)72 112.8 Q F2(caltra)108 124.8 Q/F3 10/Times-Italic@0 SF .37
242
(startdate enddate)2.84 F 2.09(start\214le \214lename)3.02 F F0([)2.68 E
243
F3(optional ar)2.73 E(guments)-.37 E F0(])2.77 E F1(DESCRIPTION)72 141.6
244
Q F0 .201(Calculate trajectories for the air parcels starting at the po\
245
sitions speci\214ed in)108 153.6 R F3(start\214le)3.041 E(.)-.15 E F0
246
.201(The trajectories co)5.201 F -.15(ve)-.15 G(r).15 E .017
247
(the time period from)108 165.6 R F3(startdate)2.857 E F0(to)2.697 E F3
248
(enddate)2.707 E F0 .017(and the trajectories are sa)2.697 F -.15(ve)-.2
249
G 2.517(di).15 G 2.517(nt)-2.517 G .017(he output \214le)-2.517 F F3
250
(\214lename)4.427 E(.)-.15 E F0 -.15(Fo)5.018 G(rw).15 E(ard)-.1 E
251
(and backw)108 177.6 Q(ard trajectories can be calculated according to \
252
the order of the start and end date.)-.1 E F1 -.81(PA)72 194.4 S
253
(RAMETERS).81 E F3(startdate)108.34 206.4 Q F0 .97(start time of the ai\
254
r parcels in the format YYYYMMDD_HH\(MM\) \(e.g. 20100101_00)183 206.4 R
255
.502(or 20100101_0030 for 1 January 2010, 00 UTC and 00:30 UTC\). Note \
256
that the minutes)183 218.4 R(\(MM\) are optional.)183 230.4 Q F3
257
(enddate)108.19 247.2 Q F0 1.059
258
(end time of the air parcels \(same format as the)183 247.2 R F3
259
(startdate)3.898 E F0 1.058(\). If the end time is after the)3.738 F
260
.319(start time, forw)183 259.2 R .319(ard trajectories are calculated;\
261
 otherwise, i.e. for end date before the start)-.1 F(date, backw)183
262
271.2 Q(ard trajectories result.)-.1 E F3(start\214le)108.34 288 Q F0
263
.442(\214le with the starting positions of the trajectories \(possibly \
264
created with)183 288 R F2(cr)2.941 E(eate_strartf)-.18 E F0(\).)2.941 E
265
(Dif)183 300 Q .227
266
(ferent formats for the "start\214le" are supported \(see)-.25 F F2 -.18
267
(re)2.728 G -.25(fo).18 G(rmat).25 E F0 .228
268
(for details\). If no format)2.728 F 3.82
269
(speci\214er \(appendix .[1234]\) is gi)183 312 R -.15(ve)-.25 G 3.82
270
(n, a simple \(longitude,latitude,pressure\) list is).15 F -.15(ex)183
271
324 S(pected.).15 E F3(\214lename)109.91 340.8 Q F0 .353
272
(output trajectory \214le with trajectories. Dif)183 340.8 R .354
273
(ferent formats are supported \(see)-.25 F F2 -.18(re)2.854 G -.25(fo)
274
.18 G(rmat).25 E F0(for)2.854 E(details\).)183 352.8 Q F1(OPTION)72
275
369.6 Q(AL)-.219 E F3(-j)108.01 381.6 Q F0 .584
276
(Jumping \215ag: if a trajectory crosses the lo)183 381.6 R .584
277
(wer boundary)-.25 F 3.084(,i)-.65 G 3.084(ti)-3.084 G 3.084(sr)-3.084 G
278
.583(aised a little and hence is)-3.084 F(allo)183 393.6 Q .023
279
(wed to mo)-.25 F .323 -.15(ve o)-.15 H .023
280
(n. Otherwise, i.e. no "-j" \215ag set, the trajectory w).15 F .023
281
(ould stick at the same)-.1 F(position. The def)183 405.6 Q
282
(ault is that "-j" is)-.1 E F2(not)2.5 E F0(set.)2.5 E F3(-i hour)108.01
283
422.4 Q(s)-.1 E F0 1.404(time increments \(in hours\) for input P and)
284
183 422.4 R 3.903<538c>6.403 G 1.403(les. If not e)-3.903 F 1.403
285
(xplicitely speci\214ed, this is)-.15 F
286
(determined from the P and S \214les in the current directory)183 434.4
287
Q(.)-.65 E F3(-t min)108.01 451.2 Q F0 1.383
288
(time step \(in minutes\) for trajcetory calculation. Per def)183 451.2
289
R 1.383(ault, the time step is 1/12 the)-.1 F .028
290
(time increment of the input \214les. F)183 463.2 R .027
291
(or instance, 6-h input P and S \214les result in a time step)-.15 F
292
1.383(of 5 min. The time step must be consisten with the output interv)
293
183 475.2 R 1.384(al \(see ne)-.25 F 1.384(xt optional)-.15 F
294
(parameter "-o"\).)183 487.2 Q F3(-o hour)108.01 504 Q(s)-.1 E F0 .757
295
(Output interv)183 504 R .757
296
(al \(in minutes\) of the air parcel positions. Per def)-.25 F .757
297
(ault it is the same as the)-.1 F 1.295(time increment between the inpu\
298
t P and S \214les \(see option "-i"\). Note that the output)183 516 R
299
(interv)183 528 Q .495(al must be a multiple of the time step for traje\
300
ctory calculation \(see optional ar)-.25 F(gu-)-.18 E(ment "-t"\).)183
301
540 Q F3(-p)108.01 556.8 Q F0 .629(Periodicity \215ag. If set, a period\
302
ic domain is assumed in zonal direction. Per def)183 556.8 R .63
303
(ault, the)-.1 F(\215ag is)183 568.8 Q F2(not)2.5 E F0(set.)2.5 E F3(-c)
304
108.01 585.6 Q(hang)-.15 E(et)-.1 E F0 .202(\215ag whether the times of\
305
 the P and S \214les should be changed or not before a calculation;)183
306
585.6 R(the def)183 597.6 Q(ault is that the times are)-.1 E F2(not)2.5
307
E F0(changed.)2.5 E F3(-noclean)108.01 614.4 Q F0 .7
308
(\215ag whether parameter and criterion \214les should be k)183 614.4 R
309
.7(ept; this is particularly helpful for)-.1 F(deb)183 626.4 Q(ugging.)
310
-.2 E F3(-timec)108.01 643.2 Q(hec)-.15 E(k)-.2 E F0
311
(enforce a time check on the data \214le)183 643.2 Q F1(EXAMPLES)72 660
312
Q F2([1] caltra 19891020_00 19891020_18 startf OUT)108 672 Q(.1)-.9 E F0
313
.336(Calculate forw)133 684 R .335(ard trajectories from 20/10/1989 00 \
314
UTC to 20/10/1989 18 UTC. The starting starting)-.1 F .45
315
(positions are gi)133 696 R -.15(ve)-.25 G 2.95(no).15 G 2.95(nt)-2.95 G
316
.451(he \214le "startf" as a list of \(longitude,latitude,pressure\) v)
317
-2.95 F .451(alues. The output tra-)-.25 F
318
(jectories are written to the \214le "OUT)133 708 Q
319
(.1", where the appendix 1 denotes ASCII format..B)-.74 E(1)535 768 Q 0
320
Cg EP
321
%%Page: 2 2
322
%%BeginPageSetup
323
BP
324
%%EndPageSetup
325
/F0 10/Times-Roman@0 SF 407.76(caltra\(\) caltra\(\))72 48 R/F1 10
326
/Times-Bold@0 SF([2] caltra 19891020_18 19891020_00 startf OUT)108 84 Q
327
(.1)-.9 E F0(As in e)133 96 Q 2.5(xample 1,)-.15 F -.2(bu)2.5 G 2.5(tb)
328
.2 G(ackw)-2.5 E
329
(ard trajectories from 20/10/1989 18 UTC to 20/10/1989 00 UTC.)-.1 E F1
330
([3] caltra 19891020_00 19891020_18 startf OUT -j)108 112.8 Q F0 .17
331
(As in e)133 124.8 R .17(xample [1], b)-.15 F .169
332
(ut with jumping \215ag set: if a trajectory crosses the lo)-.2 F .169
333
(wer boundary \(topograph)-.25 F(y\),)-.05 E
334
(it is raised a little and then is allo)133 136.8 Q(wed to mo)-.25 E .3
335
-.15(ve o)-.15 H(n.).15 E F1
336
([4] caltra 19891020_00 19891020_18 startf OUT -j -o 15 -t 15)108 153.6
337
Q F0 .375(As in e)133 165.6 R .375(xample [3], b)-.15 F .376
338
(ut the output interv)-.2 F .376
339
(al is set to 15 min with the optional ar)-.25 F .376
340
(gument "-o". Note that)-.18 F .859(the output interv)133 177.6 R .858
341
(al \(15 min\) must be a multiple of the time step, which is here set e)
342
-.25 F .858(xplicitely to 15)-.15 F(min with "-t".)133 189.6 Q F1([5] c\
343
altra 19891020_0130 19891020_1730 startf1 OUT -j -o 15 -t 15 -changet)
344
108 206.4 Q F0 .342(Start from non-analysis time 01:30 UTC to non-analy\
345
sis time 17:30 UTC. Furthermore, the times on)133 218.4 R
346
(the prinmary netCDF \214les are changed accordingly)133 230.4 Q(.)-.65
347
E/F2 10.95/Times-Bold@0 SF -.548(AU)72 247.2 S(THOR).548 E F0
348
(Written by Michael Sprenger and Heini W)108 259.2 Q
349
(ernli \(January 2011\))-.8 E(2)535 768 Q 0 Cg EP
350
%%Page: 1 3
351
%%BeginPageSetup
352
BP
353
%%EndPageSetup
354
/F0 10/Times-Roman@0 SF 394.4(datelist\(\) datelist\(\))72 48 R/F1 10.95
355
/Times-Bold@0 SF -.219(NA)72 84 S(ME).219 E/F2 10/Times-Bold@0 SF
356
(datelist - handling of datelists)108 96 Q F1(SYNOPSIS)72 112.8 Q F2
357
(datelist)108 124.8 Q/F3 10/Times-Italic@0 SF .56(\214lename mode)4.41 F
358
2.5([p)2.97 G(ar)-2.5 E(ameter)-.15 E 2.5(s])-.1 G F1(DESCRIPTION)72
359
141.6 Q F0 3.007(Ad)108 153.6 S .507(ate list is a \214le of dates in t\
360
he format {YYYYMMDD_HH}, e.g. 19900101_00 for 00 UTC, 1 January)-3.007 F
361
(1990. This command of)108 165.6 Q(fers se)-.25 E -.15(ve)-.25 G(ral w)
362
.15 E(ays ho)-.1 E 2.5(wt)-.25 G 2.5(oc)-2.5 G
363
(reate date lists and to w)-2.5 E(ork with them.)-.1 E F1 -.81(PA)72
364
182.4 S(RAMETERS).81 E F3(\214lename)109.91 194.4 Q F0 .527
365
(name of the date list \214le. If the name)158 194.4 R F2 -.37('s)3.027
366
G(tdout').37 E F0(or)3.027 E F2 -.37('s)3.027 G(cr).37 E(een')-.18 E F0
367
.527(is gi)3.027 F -.15(ve)-.25 G .528(n, the output will be directed)
368
.15 F(to standardoutput: no \214le will be created.)158 206.4 Q F3(mode)
369
108.38 223.2 Q F0(one of se)158 223.2 Q -.15(ve)-.25 G
370
(ral modes \(see belo).15 E(w\).)-.25 E F1(CREA)72 240 Q(TING D)-1.04 E
371
-1.04(AT)-.383 G 2.738(EL)1.04 G(ISTS)-2.738 E F2(-cr)108 252 Q
372
(eate stardate enddate)-.18 E F0 .487(create a datelist from startdate \
373
\(in format {YYYYMMDD_HHMM}\) to enddate; the time interv)133 264 R .486
374
(al is)-.25 F .667(per def)133 276 R .668(ault 6 h \(see option -interv)
375
-.1 F .668
376
(al\). If the start and end date do not match with the analysis times,)
377
-.25 F 2.949(the date list will contain the enclosing analysis times: f\
378
or instance, for)133 288 R F2(-cr)5.449 E 2.949(eate 20100201_04)-.18 F
379
(20100201_19)133 300 Q F0 5.17(the date list will contain the follo)7.67
380
F 5.17(wing dates: 20100201_00, 20100201_06,)-.25 F
381
(20100201_12, 20100201_18, 20100202_00.)133 312 Q F2(-indir dir)108
382
328.8 Q(name)-.15 E F0 .83
383
(search for dates \(in format {YYYYMMDD_HH}\) in the direcory gi)133
384
340.8 R -.15(ve)-.25 G 3.33(nw).15 G .83(ith {dirname} - the dates)-3.33
385
F(are written in ascending order to the datelist \214le and repeating d\
386
ates are remo)133 352.8 Q -.15(ve)-.15 G(d.).15 E F2(-inter)108 369.6 Q
387
-.1(va)-.1 G 2.5(lv).1 G(alue)-2.6 E F0(change the interv)133 381.6 Q
388
(al to {v)-.25 E(alue} hours, instead of the def)-.25 E(ault 6 hours.)
389
-.1 E F1(INFO ABOUT D)72 398.4 Q -1.04(AT)-.383 G 2.738(EL)1.04 G(ISTS)
390
-2.738 E F2(-ndates)108 410.4 Q F0
391
(write the number of dates in the list)133 422.4 Q F2(-timerange)108
392
439.2 Q F0
393
(write the time range {last date} - {\214rst date} [in hours].)133 451.2
394
Q F2(-isin date)108 468 Q F0
395
(check whether the date is in the list \(1\) or not \(0\).)133 480 Q F1
396
(STEPPING THR)72 496.8 Q(OUGH D)-.329 E -1.04(AT)-.383 G 2.738(EL)1.04 G
397
(ISTS)-2.738 E F2(-\214rst)108 508.8 Q F0
398
(write the \214rst date of the date list)133 520.8 Q F2(-last)108 537.6
399
Q F0(write the last date of the date list)133 537.6 Q F2(-next date)108
400
554.4 Q F0 .547(\214nd the date {date} in the list and write the)133
401
566.4 R F2(next)3.047 E F0 .547(date to screen; if no ne)3.047 F .548
402
(xt date is in the list, i.e. the)-.15 F
403
(end of the list is reached, 'nil' will be returned.)133 578.4 Q F2(-pr)
404
108 595.2 Q .3 -.15(ev d)-.18 H(ate).15 E F0 .329
405
(\214nd the date {date} in the list and write the)133 607.2 R F2(pr)
406
2.829 E -.15(ev)-.18 G(ious).15 E F0 .328(date to screen; if no pre)
407
2.829 F .328(vious date is in the list,)-.25 F(i.e. the be)133 619.2 Q
408
(ginning of the list is reached, 'nil' will be returned.)-.15 E F1(COMP)
409
72 636 Q(ARING D)-.81 E -1.04(AT)-.383 G 2.738(EL)1.04 G(ISTS)-2.738 E
410
F2(-o)108 648 Q -.1(ve)-.1 G(rlap \214le1 \214le1).1 E F0
411
(determine the o)133 660 Q -.15(ve)-.15 G(rlap of tw).15 E 2.5(od)-.1 G
412
(ate lists.)-2.5 E F2(-onlyin1 \214le1 \214le1)108 676.8 Q F0
413
(determine the dates which only occur in date list 1, b)133 688.8 Q
414
(ut not in datelist 2)-.2 E F2(-onlyin2 \214le1 \214le1)108 705.6 Q F0
415
(determine the dates which only occur in date list 2, b)133 717.6 Q
416
(ut not in datelist 1)-.2 E(1)535 768 Q 0 Cg EP
417
%%Page: 2 4
418
%%BeginPageSetup
419
BP
420
%%EndPageSetup
421
/F0 10/Times-Roman@0 SF 394.4(datelist\(\) datelist\(\))72 48 R/F1 10.95
422
/Times-Bold@0 SF(Examples)72 84 Q/F2 10/Times-Bold@0 SF
423
([1] datelist scr)108 96 Q(een -indir ./)-.18 E F0
424
(look for dates in the current directory and write them to screen)133
425
108 Q F2([2] datelist dates -cr)108 124.8 Q
426
(eate 19890101_00 20110101_00 -inter)-.18 E -.1(va)-.1 G 2.5(l2).1 G F0
427
.589(creates dates from 00 UTC, 1 January 1989 to 00 UTC, 1 January 201\
428
1 with an interv)133 136.8 R .59(al of 2 hours.)-.25 F
429
(The output will be written to the \214le ')133 148.8 Q(dates'.)-.5 E F1
430
-.548(AU)72 165.6 S(THOR).548 E F0
431
(Written by Michael Sprenger and Heini W)108 177.6 Q
432
(ernli \(January 2011\))-.8 E(2)535 768 Q 0 Cg EP
433
%%Page: 1 5
434
%%BeginPageSetup
435
BP
436
%%EndPageSetup
437
/F0 10/Times-Roman@0 SF(dif)72 48 Q 371.62(ference\(\) dif)-.25 F
438
(ference\(\))-.25 E/F1 10.95/Times-Bold@0 SF -.219(NA)72 84 S(ME).219 E
439
/F2 10/Times-Bold@0 SF(differ)108 96 Q(ence - calculate the differ)-.18
440
E(ence of tw)-.18 E 2.5(ot)-.1 G(rajecrtory \214les)-2.5 E F1(SYNOPSIS)
441
72 112.8 Q F2(differ)108 124.8 Q(ence)-.18 E/F3 10/Times-Italic@0 SF .01
442
(inp\214l1 inp\214le2)2.51 F 2.09(out\214le \214eld)2.75 F 2.5([-)3.56 G
443
(single|-max ])-2.5 E F1(DESCRIPTION)72 141.6 Q F0 .521
444
(Calculate the dif)108 153.6 R .521(ference of tw)-.25 F 3.021(ot)-.1 G
445
.521
446
(rajectory \214les "inp\214le1" and "inp\214le2" and write it to a ne)
447
-3.021 F 3.02(wt)-.25 G .52(rajectory \214le)-3.02 F .311
448
("out\214le". The dif)108 165.6 R .311
449
(ference is calculated of the \214eld "\214eld", which must be a)-.25 F
450
-.25(va)-.2 G .311(ilable on both input \214les. If the).25 F .036
451
(\214eld "LA)108 177.6 R .036
452
(TLON" is speci\214ed, the spherical distance between the tw)-1.11 F
453
2.536(ot)-.1 G .035(rajectories is tak)-2.536 F .035
454
(en. Furthermore, with)-.1 F("-single" the dif)108 189.6 Q(ference is w\
455
ritten at all times, with "-max" only the time of maximum dif)-.25 E
456
(ference is written.)-.25 E F1 -.81(PA)72 206.4 S(RAMETERS).81 E F3
457
(inp\214le1)108.01 218.4 Q F0(\214rst trajectory \214le)183 218.4 Q F3
458
(inp\214le2)108.01 235.2 Q F0(second trajectory \214le)183 235.2 Q F3
459
(out\214le)108.23 252 Q F0 .406(output trajectory \214le - note that th\
460
is is not a standard trajectory \214le which can be further)183 252 R
461
.266(processed! The position of both trajectories are listed and also t\
462
he \214eld for both trajecto-)183 264 R(ries and their dif)183 276 Q
463
(ference.)-.25 E F3(\214eld)109.91 292.8 Q F0 .37
464
(Name of \214eld for which the dif)183 292.8 R .37
465
(ference should be calculated. It must be a)-.25 F -.25(va)-.2 G .37
466
(ilable in both).25 F .583(trajectory \214les. IF "LA)183 304.8 R .583
467
(TLON" is speci\214ed, the spherical distance between the trajectory)
468
-1.11 F(positions is calculated.)183 316.8 Q F1(OPTION)72 333.6 Q(AL P)
469
-.219 E(ARAMETERS)-.81 E F3(-single)108.01 345.6 Q F0(the dif)183 345.6
470
Q(ference is written to "out\214le" for all trajectory times)-.25 E F3
471
(-max)108.01 362.4 Q F0 1.1(the dif)183 362.4 R 1.101(ference is writte\
472
n to "out\214le" only for the trajectory time with maximum dif)-.25 F
473
(fer)-.25 E(-)-.2 E(ence.)183 374.4 Q F1(EXAMPLES)72 391.2 Q F2
474
([1] differ)108 403.2 Q(ence tra1 tra2 out LA)-.18 E(TLON -single)-.95 E
475
F0(gi)133 415.2 Q -.15(ve)-.25 G 3.295(st).15 G .795
476
(he speherical distance \(LA)-3.295 F .795
477
(TLON\) between trajectories in "tra1" and "tra2". The dif)-1.11 F .794
478
(ference is)-.25 F(written for all trajectory times.)133 427.2 Q F2
479
([2] differ)108 444 Q(ence tra1 tra2 out TH -max)-.18 E F0(gi)133 456 Q
480
-.15(ve)-.25 G 2.88(st).15 G .38(he dif)-2.88 F .381(ference of potenti\
481
al temperature \(TH\) between trajectory \214le "ra1" and "tra2" and wr\
482
ites)-.25 F(the maximum dif)133 468 Q
483
(ference \("-max"\) to the output \214le "out".)-.25 E F1 -.548(AU)72
484
484.8 S(THOR).548 E F0(Written by Michael Sprenger and Heini W)108 496.8
485
Q(ernli \(January 2011\))-.8 E(1)535 768 Q 0 Cg EP
486
%%Page: 1 6
487
%%BeginPageSetup
488
BP
489
%%EndPageSetup
490
/F0 10/Times-Roman@0 SF -.15(ex)72 48 S 398.06(tract\(\) e).15 F
491
(xtract\(\))-.15 E/F1 10.95/Times-Bold@0 SF(extract)72 84 Q/F2 10
492
/Times-Bold@0 SF(extract - extract columns, times, single trajectories \
493
and starting positions)108 96 Q F1(SYNOPSIS)72 112.8 Q F2(extract)108
494
124.8 Q/F3 10/Times-Italic@0 SF(inptr)2.51 E 2.99(ao)-.15 G(uttr)-2.99 E
495
3.14(am)-.15 G(ode)-3.14 E F1(DESCRIPTION)72 141.6 Q F0 .146(Extract co\
496
lumns, times, single trajectories or starting positions from an input t\
497
rajectory \214le)108 153.6 R F3(inptr)2.656 E(a)-.15 E F0 .146
498
(and write)2.906 F .914(output to a ne)108 165.6 R 3.414(wt)-.25 G .914
499
(rajectory \214le)-3.414 F F3(outtr)3.644 E(a.)-.15 E F0 .914(The dif)
500
5.914 F .914(ferent e)-.25 F .915(xtraction modes are speci\214ed with)
501
-.15 F F3(mode)3.795 E(.)-.15 E F0 .915(Note: the)5.915 F .093
502
(time, longitude, latitude and pressure need not be e)108 177.6 R .093
503
(xtracted beacuse the)-.15 F 2.593(ya)-.15 G .093(re an inte)-2.593 F
504
.093(gral part of e)-.15 F -.15(ve)-.25 G .092(ry trajec-).15 F
505
(tory \214le - the)108 189.6 Q 2.5(ya)-.15 G(re e)-2.5 E
506
(xtracted by def)-.15 E(ault.)-.1 E F1(EXTRA)72 206.4 Q(CTION MODE)-.602
507
E F2(-v)108 218.4 Q(ar)-.1 E F0 -.15(ex)133 218.4 S 1.022(tract columns\
508
 of a trajectory \214le; the columns can be listed by name \(e.g. -v).15
509
F 1.022(ar TH PV RH\) or a)-.25 F .254
510
(range of columns can be speci\214ed by the "to" operator \(e.g. -v)133
511
230.4 R .254(ar TH to PV\). The tw)-.25 F 2.754(om)-.1 G .254
512
(odes can also)-2.754 F(be combined: "-v)133 242.4 Q(ar TH to PV RH" e)
513
-.25 E(xtracts all columns between TH and PV)-.15 E 2.5(,a)-1.29 G
514
(nd additionally RH.)-2.5 E F2(-time)108 259.2 Q F0 -.15(ex)133 271.2 S
515
.587(tract trajectory times; the times can be gi).15 F -.15(ve)-.25 G
516
3.087(na).15 G 3.088(sal)-3.087 G .588
517
(ist of times \(e.g. -time 6 12\) or as a time range)-3.088 F
518
(\(e.g. -time 6 to 18\).)133 283.2 Q F2(-tra)108 300 Q F0 -.15(ex)133
519
300 S .258(tract single trajectories; the inde).15 F 2.758(xo)-.15 G
520
2.758(ft)-2.758 G .257
521
(he trajectories can be speci\214ed as a list \(e.g. -tra 10 12 14\) or)
522
-2.758 F(as a range of trajectories \(e.g. -tra 10 to 20\).)133 312 Q F2
523
(-startf)108 328.8 Q F0 -.15(ex)133 340.8 S .461(tract list \(longitude\
524
, latitude, pressure\) of starting positions ofvthe trajectory \214le \
525
\(corresponding to).15 F(time 0\).)133 352.8 Q F2(-index)108 369.6 Q F0
526
-.15(ex)133 381.6 S 1.381(tract single trajectories - the trajectory in\
527
dices \(from 1 to #trajectories\) are gi).15 F -.15(ve)-.25 G 3.88(no)
528
.15 G 3.88(na\214)-3.88 G 1.38(le \(e.g.)-3.88 F(-inde)133 393.6 Q 2.5
529
<788c>-.15 G(lename\).)-2.5 E F2(-boolean)108 410.4 Q F0 -.15(ex)133
530
422.4 S(tract single trajectories - the trajectory are speci\214ed on a\
531
 boolean \(0/1\) \214le \(e.g. -boolean \214lename\).).15 E F2(-patter)
532
108 439.2 Q(n)-.15 E F0 -.15(ex)133 451.2 S 1.82
533
(tract all trajectories which match the pattern gi).15 F -.15(ve)-.25 G
534
1.82(n; the pattern is a list of numbers. It is then).15 F(check)133
535
463.2 Q
536
(ed whether these numbers occur in a trajectory \(all at one time\).)-.1
537
E F1(EXAMPLES)72 480 Q F2([1] extract inptra outtra -time 6 to 36 72)108
538
492 Q F0 .547(read input trajectory \214le "inptra", e)133 504 R .547
539
(xtract times 6 to 36 and additionally time 72, and write output to)-.15
540
F(trajectory \214le "outtra".)133 516 Q F2
541
([2] extract inptra outtra -index ind\214le)108 532.8 Q F0 .637(reads i\
542
nput trajectories from "inptra" and write all trajectories to "outtra" \
543
which are listed in the \214le)133 544.8 R("ind\214le". In "in\214le" t\
544
he indices of selected trajectories are listed line-by-line.)133 556.8 Q
545
F2([3] extract inptra outtra -patter)108 573.6 Q 2.5(n-)-.15 G(999.99)
546
-2.5 E F0 -.15(ex)133 585.6 S(tracts all trajectories which ha).15 E .3
547
-.15(ve a m)-.2 H(issing data \(-999.99\) entry).15 E(.)-.65 E F2
548
([4] extract inptra outtra -patter)108 602.4 Q 2.5(n0)-.15 G
549
(.00 -44.25 -28.47 140)-2.5 E F0 -.15(ex)133 614.4 S .207
550
(tracts all trajectories which ha).15 F .507 -.15(ve t)-.2 H .207
551
(he numbers '0.00).15 F 2.707(-44.25 -28.47)7.707 F .206
552
(140' in their list. This options)7.707 F(is con)133 626.4 Q -.15(ve)-.4
553
G(nient to search for speci\214c times and positions, and then to see t\
554
he complete trajectory).15 E(.)-.65 E F1 -.548(AU)72 643.2 S(THOR).548 E
555
F0(Written by Michael Sprenger and Heini W)108 655.2 Q
556
(ernli \(January 2011\))-.8 E(1)535 768 Q 0 Cg EP
557
%%Page: 1 7
558
%%BeginPageSetup
559
BP
560
%%EndPageSetup
561
/F0 10/Times-Roman@0 SF(gettidif)72 48 Q 388.24(f\(\) gettidif)-.25 F
562
(f\(\))-.25 E/F1 10.95/Times-Bold@0 SF -.219(NA)72 84 S(ME).219 E/F2 10
563
/Times-Bold@0 SF(gettidiff - get the differ)108 96 Q
564
(ence \(in hours\) between tw)-.18 E 2.5(od)-.1 G(ates)-2.5 E F1
565
(SYNOPSIS)72 112.8 Q F2(gettidiff)108 124.8 Q/F3 10/Times-Italic@0 SF
566
.35(date1 date2)2.85 F F1(DESCRIPTION)72 141.6 Q F0 .68
567
(Get the time dif)108 153.6 R .68(ference \(in hours\) between tw)-.25 F
568
3.181(od)-.1 G .681(ates \( in form YYYYMMDD_HH\(MM\), i.e. the minutes)
569
-3.181 F
570
(are optional\): date1 - date2. The output is in format HH or HH.MM)108
571
165.6 Q F1 -.81(PA)72 182.4 S(RAMETERS).81 E F2(date1)108 194.4 Q F0
572
(\214rst date \( in form YYYYMMDD_HH\(MM\) \))183 194.4 Q F2(date2)108
573
211.2 Q F0(second date \( in form YYYYMMDD_HH\(MM\) \))183 211.2 Q F1
574
(EXAMPLES)72 228 Q F2([1] gettidiff 20110102_18 20110205_00)108 240 Q F0
575
(gi)133 252 Q -.15(ve)-.25 G 2.5(st).15 G(he time dif)-2.5 E
576
(ference between 20110102_18 and 20110205_00; the result is -798 h.)-.25
577
E F1 -.548(AU)72 268.8 S(THOR).548 E F0
578
(Written by Michael Sprenger and Heini W)108 280.8 Q
579
(ernli \(January 2011\))-.8 E(1)535 768 Q 0 Cg EP
580
%%Page: 1 8
581
%%BeginPageSetup
582
BP
583
%%EndPageSetup
584
/F0 10/Times-Roman@0 SF(getv)72 48 Q 394.92(ars\(\) getv)-.25 F(ars\(\))
585
-.25 E/F1 10.95/Times-Bold@0 SF -.219(NA)72 84 S(ME).219 E/F2 10
586
/Times-Bold@0 SF(getv)108 96 Q
587
(ars - get a list of a \214elds on a netCDF \214le)-.1 E F1(SYNOPSIS)72
588
112.8 Q F2(getv)108 124.8 Q(ars)-.1 E/F3 10/Times-Italic@0 SF
589
(\214lename)4.41 E F1(DESCRIPTION)72 141.6 Q F0(Get a list of all \214e\
590
lds on a netCDF \214le - the names of the \214elds are listed line-by-l\
591
ine.)108 153.6 Q F1 -.81(PA)72 170.4 S(RAMETERS).81 E F2(\214lename)108
592
182.4 Q F0(name of the netCDF \214le \(e.g. P20110102_00\).)183 182.4 Q
593
F1(EXAMPLES)72 199.2 Q F2([1] getv)108 211.2 Q(ars P20110102_00)-.1 E F0
594
(gi)133 223.2 Q -.15(ve)-.25 G 2.938(sal).15 G .438
595
(ist of all \214elds on the netCDF \214le P20110102_00; the results is:)
596
-2.938 F .438(time, Q, L)5.438 F .438(WC, IWC, T)-.74 F 2.937(,U)-.74 G
597
(,)-2.937 E 2.58 -1.29(V, O)133 235.2 T(MEGA, PS, SLP)1.29 E F1 -.548
598
(AU)72 252 S(THOR).548 E F0(Written by Michael Sprenger and Heini W)108
599
264 Q(ernli \(January 2011\))-.8 E(1)535 768 Q 0 Cg EP
600
%%Page: 1 9
601
%%BeginPageSetup
602
BP
603
%%EndPageSetup
604
/F0 10/Times-Roman@0 SF 345.54(Lagrantohelp\(\) Lagrantohelp\(\))72 48 R
605
/F1 10.95/Times-Bold@0 SF -.219(NA)72 84 S(ME).219 E/F2 10/Times-Bold@0
606
SF(Lagranto -)108 96 Q F0
607
(calculate trajectories for ECMWF analyses and forecasts)2.5 E F1
608
(SYNOPSIS)72 112.8 Q F2 .204(lagrantohelp [topic|tutorial|r)108 124.8 R
609
(efer)-.18 E(ence|futur)-.18 E .204(e] -)-.18 F F0(sho)2.704 E 2.704(wt)
610
-.25 G .205(his man page or the one corresponding to a speci\214c)-2.704
611
F .853(topic. All)108 136.8 R/F3 10/Times-Italic@0 SF(underlined)3.353 E
612
F0 .853(names in this document ha)3.353 F 1.152 -.15(ve t)-.2 H .852
613
(heir o).15 F .852(wn man page. If "tutorial" is chosen, a detailed)-.25
614
F 1.248(tutorial for Lagranto is opened. F)108 148.8 R 1.249
615
(or "reference" a PDF document with all man pages is opened. Finally)
616
-.15 F(,)-.65 E
617
(with option "future" plans for additional features will be listed.)108
618
160.8 Q F1(DESCRIPTION)72 177.6 Q F0 1.24(Lagranto is a softw)108 189.6
619
R 1.239(are tool \(UNIX shell-scripts and F)-.1 F 1.239
620
(ortran programs\) to analyze Lagrangian aspects of)-.15 F 1.4(atmosphe\
621
ric phenomena. It requires a time-series of 3-dimensional wind \214elds\
622
 \(and if necessary further)108 201.6 R -.25(va)108 213.6 S
623
(riables\) on netCDF \214les.).25 E F1(CORE PR)72 230.4 Q(OGRAMS)-.329 E
624
F3(cr)108.2 242.4 Q(eate_startf)-.37 E F0(create starting \214les)183
625
242.4 Q F3(caltr)108.2 259.2 Q(a)-.15 E F0
626
(calculate trajectory positions)183 259.2 Q F3(tr)108.13 276 Q(ace)-.15
627
E F0(trace meteorological \214elds along trajectories)183 276 Q F3
628
(select)108.34 292.8 Q F0(select trajectories based on se)183 292.8 Q
629
-.15(ve)-.25 G(ral criteria).15 E F3(density)108.35 309.6 Q F0
630
(create netCDF \214les with trajectory densities)183 309.6 Q F1(TRAJECT)
631
72 326.4 Q(OR)-.197 E 2.738(YT)-.383 G(OOLS)-2.935 E F3 -.2(ex)108 338.4
632
S(tr).2 E(act)-.15 E F0 -.15(ex)183 338.4 S
633
(tract single trajectories, times or columns from a trajectory \214le)
634
.15 E F3(list2lsl)108.09 355.2 Q F0(transform a \(longitude,latitude,pr\
635
essure\) list into a trajectory \214le)183 355.2 Q F3(lsl2list)108.09
636
372 Q F0(transform a trajectory \214le into a \(longitude,latitude,pres\
637
sure\) list)183 372 Q F3(mer)108.38 388.8 Q -.1(ge)-.37 G(tr).1 E(a)-.15
638
E F0(mer)183 388.8 Q(ge tw)-.18 E 2.5(ot)-.1 G(rajectory \214les)-2.5 E
639
F3 -.37(re)108 405.6 S(format).37 E F0
640
(change the format of a trajectory \214le)183 405.6 Q F3(timer)108.13
641
422.4 Q(es)-.37 E F0
642
(change the time resolution of a trajectory \214le with interpolation)
643
183 422.4 Q F3(tr)108.13 439.2 Q(ainfo)-.15 E F0
644
(write some information about a trajectory \214le)183 439.2 Q F3(dif)
645
108.35 456 Q(fer)-.18 E(ence)-.37 E F0(calculate the dif)183 456 Q
646
(ference between tw)-.25 E 2.5(ot)-.1 G(rajectory \214les)-2.5 E F1
647
(NETCDF T)72 472.8 Q(OOLS)-.197 E F3 -.15(ch)108 484.8 S(ang).15 E(ecst)
648
-.1 E F0(change the constants \214lename on a netCDF \214le)183 484.8 Q
649
F3 -.15(ch)108 501.6 S(ang).15 E(et)-.1 E F0
650
(change the time on a netCDF \214le)183 501.6 Q F3 -.1(ge)108 518.4 S
651
(tmima).1 E F0(get the minimum and maximum v)183 518.4 Q
652
(alue of a \214eld on a netCDF \214le)-.25 E F3 -.1(ge)108 535.2 S(tvar)
653
.1 E(s)-.1 E F0(get a list of \214elds on a netCDF \214le)183 535.2 Q F1
654
(GENERAL T)72 552 Q(OOLS)-.197 E F3(ne)108.36 564 Q(wtime)-.15 E F0
655
(get a ne)183 564 Q 2.5(wd)-.25 G(ate from an initial date and an of)
656
-2.5 E(fset in hours)-.25 E F3 -.1(ge)108 580.8 S(ttidif).1 E(f)-.18 E
657
F0(get the time dif)183 580.8 Q(ference between tw)-.25 E 2.5(od)-.1 G
658
(ates)-2.5 E F1(VISU)72 597.6 Q(ALISA)-.657 E(TION)-1.04 E F0 1.927
659
(Lagranto comes with tw)108 609.6 R 4.427(op)-.1 G 1.927
660
(ackages for visualisation of trajectories: one based on Matlab \(www)
661
-4.427 F(.math-)-.65 E -.1(wo)108 621.6 S 1.158
662
(rks.com\) and one based on NCL \(www).1 F(.ncl.ucar)-.65 E 1.159
663
(.edu\). Examples are pro)-.55 F 1.159
664
(vided in folder "\214gure" of the)-.15 F(Lagranto installation.)108
665
633.6 Q F1 -.548(AU)72 650.4 S(THOR).548 E F0
666
(Written by Michael Sprenger and Heini W)108 662.4 Q
667
(ernli \(January 2011\))-.8 E(1)535 768 Q 0 Cg EP
668
%%Page: 1 10
669
%%BeginPageSetup
670
BP
671
%%EndPageSetup
672
/F0 10/Times-Roman@0 SF 398.82(list2lsl\(\) list2lsl\(\))72 48 R/F1
673
10.95/Times-Bold@0 SF -.219(NA)72 84 S(ME).219 E/F2 10/Times-Bold@0 SF
674
(list2lsl - con)108 96 Q -.1(ve)-.4 G
675
(rt a list of longitude, latitude and pr).1 E(essur)-.18 E 2.5(ei)-.18 G
676
(nto a trajectory \214le.)-2.5 E F1(SYNOPSIS)72 112.8 Q F2(list2lsl)108
677
124.8 Q/F3 10/Times-Italic@0 SF .41(inp\214le out\214le)2.51 F 2.5([o)
678
2.97 G(ptional ar)-2.5 E(guments ])-.37 E F1(DESCRIPTION)72 141.6 Q F0
679
(Con)108 153.6 Q -.15(ve)-.4 G(rt an input \214le "inp\214le" wqith).15
680
E 2.5(al)5 G
681
(ist of longitude, latitude and pressure into a trajectory \214le.)-2.5
682
E F1 -.81(PA)72 170.4 S(RAMETERS).81 E F3(inp\214le)108.01 182.4 Q F0
683
(input \214le \(a list of longitude, latitude and pressure\))183 182.4 Q
684
F3(out\214le)108.23 199.2 Q F0(output trajectory \214le)183 199.2 Q F1
685
(OPTION)72 216 Q(AL P)-.219 E(ARAMETERS)-.81 E F2(-r)108 228 Q(ef r)-.18
686
E(efdate)-.18 E F0(set the reference date of the trajectory \214le to "\
687
refdate", e.g. "-ref 20110102_00".)133 240 Q F2(-time v)108 256.8 Q
688
(alue)-.1 E F0(set the time of the trajectory \214le to "v)133 268.8 Q
689
(alue" \(in format HH.MM\): e.g. "-time 6" sets the time to 6 h.)-.25 E
690
F1(EXAMPLES)72 285.6 Q F2([1] list2lsl list\214le tra\214le -r)108 297.6
691
Q(efdate 20110101_00 -time 6)-.18 E F0(con)133 309.6 Q -.15(ve)-.4 G
692
.457(rts the list "list\214le" into a trajectory \214le "tra\214le", wh\
693
ere the reference date is 20110101_00 and).15 F(the time is 6 h.)133
694
321.6 Q F1 -.548(AU)72 338.4 S(THOR).548 E F0
695
(Written by Michael Sprenger and Heini W)108 350.4 Q
696
(ernli \(January 2011\))-.8 E(1)535 768 Q 0 Cg EP
697
%%Page: 1 11
698
%%BeginPageSetup
699
BP
700
%%EndPageSetup
701
/F0 10/Times-Roman@0 SF 398.82(lsl2list\(\) lsl2list\(\))72 48 R/F1
702
10.95/Times-Bold@0 SF -.219(NA)72 84 S(ME).219 E/F2 10/Times-Bold@0 SF
703
(lsl2list - con)108 96 Q -.1(ve)-.4 G
704
(rt a trajectory \214le into a list of longitude, latituden and pr).1 E
705
(essur)-.18 E(e)-.18 E F1(SYNOPSIS)72 112.8 Q F2(lsl2list)108 124.8 Q/F3
706
10/Times-Italic@0 SF .41(inp\214le out\214le)2.51 F F1(DESCRIPTION)72
707
141.6 Q F0(Con)108 153.6 Q -.15(ve)-.4 G .195(rt an input trajectory \
708
\214le "inp\214le" into a list of longitude, latitude and pressure; all\
709
 times are included).15 F(in this list.)108 165.6 Q F1 -.81(PA)72 182.4
710
S(RAMETERS).81 E F3(inp\214le)108.01 194.4 Q F0(input trajectory \214le)
711
183 194.4 Q F3(out\214le)108.23 211.2 Q F0
712
(output \214le \(a list of longitude, latitude and pressure\))183 211.2
713
Q F1(EXAMPLES)72 228 Q F2(lsl2list tra\214le list\214le)108 240 Q F0
714
(con)133 252 Q -.15(ve)-.4 G(rts all entries of the trajectory \214le "\
715
tra\214le" into a list\214le with longitude, latitude and pressure.).15
716
E F1 -.548(AU)72 268.8 S(THOR).548 E F0
717
(Written by Michael Sprenger and Heini W)108 280.8 Q
718
(ernli \(January 2011\))-.8 E(1)535 768 Q 0 Cg EP
719
%%Page: 1 12
720
%%BeginPageSetup
721
BP
722
%%EndPageSetup
723
/F0 10/Times-Roman@0 SF(mer)72 48 Q 381.46(getra\(\) mer)-.18 F
724
(getra\(\))-.18 E/F1 10.95/Times-Bold@0 SF -.219(NA)72 84 S(ME).219 E/F2
725
10/Times-Bold@0 SF(mer)108 96 Q(getra -)-.1 E F0(combine tw)2.5 E 2.5
726
(ot)-.1 G(rajectory \214les)-2.5 E F1(SYNOPSIS)72 112.8 Q F2(mer)108
727
124.8 Q(getra)-.1 E/F3 10/Times-Italic@0 SF .01(inp\214le1 inp\214le2)
728
2.51 F(out\214le)2.75 E F1(DESCRIPTION)72 141.6 Q F0 .077(Combines tw)
729
108 153.6 R 2.577(oi)-.1 G .077
730
(nput trajectory \214les "inp\214le1" and "inp\214le2" and writes a ne)
731
-2.577 F 2.576(wc)-.25 G .076(ombined trajectory \214le "out-)-2.576 F
732
(\214le". Three dif)108 165.6 Q
733
(ferent modes for combination are supported \(see belo)-.25 E(w\).)-.25
734
E F1 -.81(PA)72 182.4 S(RAMETERS).81 E F3(inp\214le1)108.01 194.4 Q F0
735
(\214rst input trajectory \214le)183 194.4 Q F3(inp\214le2)108.01 211.2
736
Q F0(second input trajectory \214le)183 211.2 Q F3(out\214le)108.23 228
737
Q F0(output trajectory \214le)183 228 Q F1(MODES FOR COMBIN)72 244.8 Q
738
-1.04(AT)-.219 G(ION)1.04 E F0 .128(Three dif)108 256.8 R .128
739
(ferent combination modes are supported. The)-.25 F 2.629(ya)-.15 G .129
740
(re tested for sequentially and the \214rst mode which)-2.629 F
741
(is accepted, based on the criteria listed belo)108 268.8 Q 1.3 -.65
742
(w, i)-.25 H 2.5(sp).65 G(erformed. The three modes are:)-2.5 E F2 2.5
743
(-c)108 285.6 S(olumn)-2.5 E F0 .271(The input \214les ha)133 297.6 R
744
.571 -.15(ve t)-.2 H .271(he same trajectory times, longitudes, latitud\
745
es and pressures, and both \214les con-).15 F 1.19
746
(tain the same number of trajectories. The ne)133 309.6 R 3.691(wt)-.25
747
G 1.191(rajectory \214le contains all columns from both input)-3.691 F
748
.035(\214les, duplicate columns only included once, e.g. if "inp\214le1\
749
" has columns "time,lon,lat,p,TH,PV" and)133 321.6 R
750
("inp\214le2" columns "time,lon,lat,p,PV)133 333.6 Q
751
(,Q", the output \214le will ha)-1.29 E .3 -.15(ve ")-.2 H
752
(time,lon,lat,p,TH,PV).15 E(,Q".)-1.29 E F2 2.5(-a)108 350.4 S(ppend)
753
-2.5 E F0 .39(The input \214les ha)133 362.4 R .69 -.15(ve t)-.2 H .39
754
(he same trajectory times and the same columns, b).15 F .391(ut dif)-.2
755
F .391(fer in the number of tra-)-.25 F(jectories. Then the second \214\
756
le is appended to the \214rst one and written as a ne)133 374.4 Q 2.5
757
(wt)-.25 G(rajectory \214le.)-2.5 E F2 2.5(-t)108 391.2 S(imes)-2.5 E F0
758
.077(The input \214les dif)133 403.2 R .077(fer only in the time v)-.25
759
F .076(alues, e.g. the \214rst \214le might contain times -96 h to 0 h \
760
and the)-.25 F .399(second the times from 0 h to 96 hours. Then the ne)
761
133 415.2 R 2.9(wc)-.25 G .4
762
(ombined \214le will contain the time -96 h to 96)-2.9 F 1.238(h. Note \
763
that duplicate times are eliminated. Furthermore, the output times are \
764
sorted in increasing)133 427.2 R(order)133 439.2 Q(.)-.55 E F1(EXAMPLES)
765
72 456 Q F2([1] mer)108 468 Q(getra \214le1 \214le2 new\214le)-.1 E F0
766
1.283(combines the tw)133 480 R 3.783(ot)-.1 G 1.283
767
(rajectory \214les "\214le1" and "\214le2" and writes a ne)-3.783 F
768
3.783(wt)-.25 G 1.283(rajectory \214le; the mode for)-3.783 F
769
(combination is automatically detected.)133 492 Q F1 -.548(AU)72 508.8 S
770
(THOR).548 E F0(Written by Michael Sprenger and Heini W)108 520.8 Q
771
(ernli \(January 2011\))-.8 E(1)535 768 Q 0 Cg EP
772
%%Page: 1 13
773
%%BeginPageSetup
774
BP
775
%%EndPageSetup
776
/F0 10/Times-Roman@0 SF(ne)72 48 Q 383.8(wtime\(\) ne)-.25 F(wtime\(\))
777
-.25 E/F1 10.95/Times-Bold@0 SF -.219(NA)72 84 S(ME).219 E/F2 10
778
/Times-Bold@0 SF(newtime -)108 96 Q F0(get a ne)2.5 E 2.5(wd)-.25 G
779
(ate string from an initial date and an of)-2.5 E(fset in hours)-.25 E
780
F1(SYNOPSIS)72 112.8 Q F2(newtime)108 124.8 Q/F3 10/Times-Italic@0 SF
781
.41(date of)2.85 F(fset)-.18 E F1(DESCRIPTION)72 141.6 Q F0(Gi)108 153.6
782
Q -.15(ve)-.25 G 2.568(na).15 G 2.568(ni)-2.568 G .067
783
(nitial date string in format "YYYYMMDD_HH" and an of)-2.568 F .067
784
(fset in hours, create a ne)-.25 F 2.567(wd)-.25 G .067(ate string -)
785
-2.567 F(agin in the form "YYYYMMDD_HH".)108 165.6 Q F1 -.81(PA)72 182.4
786
S(RAMETERS).81 E F3(date)108.35 194.4 Q F0 1.035(initial date in form "\
787
YYYYMMDD_HH", e.g. 20110102_18 for 18 UTC, 2nd January)183 194.4 R
788
(2011.)183 206.4 Q F3(of)108.23 223.2 Q(fset)-.18 E F0(Of)183 223.2 Q
789
(fset in hours; the of)-.25 E(fset can be positi)-.25 E .3 -.15(ve \()
790
-.25 H(into the future\) or ne).15 E -.05(ga)-.15 G(ti).05 E .3 -.15
791
(ve \()-.25 H(into the past\).).15 E F1(EXAMPLES)72 240 Q F2
792
([1] newtime 20110101_00 45)108 252 Q F0(gi)133 264 Q -.15(ve)-.25 G 2.5
793
(s2).15 G(0110102_21, the datestring 45 h after 20110101_00.)-2.5 E F1
794
-.548(AU)72 280.8 S(THOR).548 E F0
795
(Written by Michael Sprenger and Heini W)108 292.8 Q
796
(ernli \(January 2011\))-.8 E(1)535 768 Q 0 Cg EP
797
%%Page: 1 14
798
%%BeginPageSetup
799
BP
800
%%EndPageSetup
801
/F0 10/Times-Roman@0 SF 383.32(reformat\(\) reformat\(\))72 48 R/F1
802
10.95/Times-Bold@0 SF -.219(NA)72 84 S(ME).219 E/F2 10/Times-Bold@0 SF
803
-.18(re)108 96 S -.25(fo).18 G(rmat - con).25 E -.1(ve)-.4 G
804
(rt trajectory \214les between differ).1 E(ent f)-.18 E(ormats)-.25 E F1
805
(SYNOPSIS)72 112.8 Q F2 -.18(re)108 124.8 S -.25(fo).18 G(rmat).25 E/F3
806
10/Times-Italic@0 SF .41(inp\214le out\214le)2.51 F F1(DESCRIPTION)72
807
141.6 Q F0(Con)108 153.6 Q -.15(ve)-.4 G .568
808
(rt a Lagranto trajectory \214le).15 F F3(inp\214le)3.078 E F0 .568
809
(from one format to another format; the ne)3.248 F 3.068<778c>-.25 G
810
.567(le is written to)-3.068 F F3(out-)3.297 E(\214le)108 165.6 Q(.)-.15
811
E F0 .278(The formats are speci\214ed with an appendix to the \214lenam\
812
e, e.g. "trajectory)5.277 F .278(.1" speci\214es format 1. If no)-.65 F
813
(appendix is gi)108 177.6 Q -.15(ve)-.25 G(n, format 1 is chosen.).15 E
814
F1 -.81(PA)72 194.4 S(RAMETERS).81 E F3(inp\214le)108.01 206.4 Q F0
815
(input trajectory \214le)183 206.4 Q F3(out\214le)108.23 223.2 Q F0
816
(output trajectory \214le \(can be the same as inp\214le\).)183 223.2 Q
817
F1(FORMA)72 240 Q(TS)-1.04 E F0 -.15(Fo)108 252 S(rmats must be speci\
818
\214ed with an appendix to the \214lename. If no appendix is gi).15 E
819
-.15(ve)-.25 G(n, format 1 is chosen.).15 E F2(.1)108 268.8 Q F0 .949(A\
820
SCII \214le; the trajectories are sorted according to their starting po\
821
sitions. Dif)133 268.8 R .948(ferent trajectories are)-.25 F .12(separa\
822
ted by a blank line. This format is the only one supported by the Matla\
823
b and NCL visualisation)133 280.8 R(scripts \(see)133 292.8 Q F2
824
(lagrantohelp/visualisation\).)2.5 E(.2)108 309.6 Q F0 .971(ASCII \214l\
825
e; the trajectories are sorted according to their times. Dif)133 309.6 R
826
.97(ferent times are separated with a)-.25 F(blank line.)133 321.6 Q F2
827
(.3)108 338.4 Q F0(unformatted F)133 338.4 Q
828
(ortran; this is the most ef)-.15 E(\214cient format, b)-.25 E
829
(ut least portable one.)-.2 E F2(.4)108 355.2 Q F0
830
(netCDF; portable and compact data format.)133 355.2 Q F1(EXAMPLES)72
831
372 Q F2([1] r)108 384 Q(ef)-.18 E(ormat trainp.1 trainp.4)-.25 E F0
832
(Con)133 396 Q -.15(ve)-.4 G(rt the input trajectory \214le in format 1\
833
 \(ASCII\) to format 4 \(netCDF\).).15 E F1 -.548(AU)72 412.8 S(THOR)
834
.548 E F0(Written by Michael Sprenger and Heini W)108 424.8 Q
835
(ernli \(January 2011\))-.8 E(1)535 768 Q 0 Cg EP
836
%%Page: 1 15
837
%%BeginPageSetup
838
BP
839
%%EndPageSetup
840
/F0 10/Times-Roman@0 SF 406.64(select\(\) select\(\))72 48 R/F1 10.95
841
/Times-Bold@0 SF -.219(NA)72 84 S(ME).219 E/F2 10/Times-Bold@0 SF
842
(select -)108 96 Q F0(select trajectories)2.5 E F1(SYNOPSIS)72 112.8 Q
843
F2(select)108 124.8 Q/F3 10/Times-Italic@0 SF(inptr)2.51 E 2.99(ao)-.15
844
G(uttr)-2.99 E 2.96(ac)-.15 G(riterion)-2.96 E F1(DESCRIPTION)72 141.6 Q
845
F0 .535(Select trajectories from the input trajectories in "inptra" bas\
846
ed upon meteorological conditions speci\214ed in)108 153.6 R .7
847
("criterion" \(to be described belo)108 165.6 R .7
848
(w\). The selected trajectories are then written to a ne)-.25 F 3.2(wt)
849
-.25 G .7(rajectory \214le "out-)-3.2 F(tra".)108 177.6 Q F1 -.81(PA)72
850
194.4 S(RAMETERS).81 E F3(inptr)108.01 206.4 Q(a)-.15 E F0
851
(input trajectory \214le)183 206.4 Q F3(outtr)108.23 223.2 Q(a)-.15 E F0
852
(output trajectory \214le)183 223.2 Q F3(criterion)108.2 240 Q F0 .234(\
853
speci\214cation of the selection criterion; the speci\214cation is eith\
854
er an e)183 240 R .234(xplicit criterion or a)-.15 F(\214le containing \
855
the speci\214cation. Each selection criterion has the follo)183 252 Q
856
(wing form:)-.25 E F2(COMMAND)183 264 Q F0(:)2.5 E F2(FIELD)2.5 E F0(:)
857
2.5 E F2(ARGUMENTS)2.5 E F0(:)2.5 E F2 2.5([T)2.5 G(IME ].)-2.5 E F0(Se)
858
183 276 Q -.15(ve)-.25 G 1.105(ral selection criteria can be combined w\
859
ith logical operator & \(AND\) and | \(OR\),).15 F(the AND ha)183 288 Q
860
(ving higher priority than the OR.)-.2 E F1(OPTION)72 304.8 Q
861
(AL ARGUMENTS)-.219 E F3(-noclean)108.01 316.8 Q F0 -.1(ke)133 328.8 S
862
(ep temporary \214les for deb).1 E(ugging.)-.2 E F3(-boolean)108.01
863
345.6 Q F0 .286(Write a boolean list \(0/1\) instead of a trajectory \
864
\214le; the \214le has #trajectories entries, each line corre-)133 357.6
865
R(sponding to an input trajectory /1=trajectory selected, 0=not selecte\
866
d\).)133 369.6 Q F3(-inde)108.01 386.4 Q(x)-.2 E F0 .903(Write an inde)
867
133 398.4 R 3.403(xl)-.15 G .903
868
(ist instead of a trajectory \214le: the inde)-3.403 F 3.403(xo)-.15 G
869
3.403(fa)-3.403 G .902(ll selected trajectories is written to the)-3.403
870
F(output \214le - the inde)133 410.4 Q 2.5(xr)-.15 G(anges from 1 to)
871
-2.5 E(#trajectories.)5 E F3(-count)108.01 427.2 Q F0
872
(Write only the number of selected trajectories to the output \214le.)
873
133 439.2 Q F3(-startf)108.01 456 Q F0(Write only the starting position\
874
s of selected trajectories to the output \214le.)133 468 Q F3(-r)108.01
875
484.8 Q -.4(eg)-.37 G(ionf \214lename).4 E F0 .295(change the re)133
876
496.8 R .295(gion \214le from its def)-.15 F .295(ault v)-.1 F .295
877
(alue "re)-.25 F .296(gionf" to a ne)-.15 F 2.796<778c>-.25 G .296
878
(le name: the syntax is "-re)-2.796 F .296(gionf \214le-)-.15 F(name".)
879
133 508.8 Q F1(COMMAND)72 525.6 Q F2 2.5(-G)108 537.6 S(T)-2.5 E F0
880
(greater than: e.g.)133 549.6 Q F2(GT)5 E(:PV)-.74 E(:2)-.92 E F0
881
(selects trajectories with \214rst potential v)2.5 E
882
(orticity \(PV\) lar)-.2 E(ger than 2 PVU.)-.18 E F2 2.5(-L)108 566.4 S
883
(T)-3.42 E F0(less than: e.g.)133 566.4 Q F2 -1.4 -.92(LT :)5 H(RH:70)
884
.92 E F0(selects trajectories with \214rst relati)2.5 E .3 -.15(ve h)
885
-.25 H(umidity \(RH\) belo).15 E 2.5(w7)-.25 G 2.5(0%)-2.5 G(.)-2.5 E F2
886
2.5(-I)108 583.2 S(N)-2.5 E F0(within: e.g.)133 583.2 Q F2(IN:lon:30,40)
887
5 E F0(selects trajectories with \214rst longitude between 30 and 40 de)
888
2.5 E(g.)-.15 E F2 2.5(-O)108 600 S(UT)-2.5 E F0 .578(outside: e.g.)133
889
612 R F2(OUT)5.578 E(:lat:-30,30)-.74 E F0 .578
890
(selects trajectories with \214rst latitude outside -30 and 30 de)3.078
891
F 3.077(g-n)-.15 G -.15(eg)-3.077 G(lect-).15 E
892
(ing an equatoriual/subtropical band.)133 624 Q F2 2.5(-E)108 640.8 S(Q)
893
-2.5 E F0(equal: e.g.)133 652.8 Q F2(EQ:p:460)5 E F0
894
(selects trajectories with \214rst pressure equal to 460 hP)2.5 E(a.)
895
-.15 E F2 2.5(-T)108 669.6 S -.3(RU)-2.5 G(E).3 E F0 .35
896
(check whether v)133 681.6 R .351(alue is dif)-.25 F .351
897
(ferent from zero \(logical TR)-.25 F .351(UE\): e.g.)-.4 F F2(TR)5.351
898
E(UE:CYCL::ALL\(ANY\))-.3 E F0(checks)2.851 E .244
899
(whether the trajectory passes through a c)133 693.6 R .244
900
(ycclone, which is mark)-.15 F .243(ed as 0/1 \214eld, at an)-.1 F 2.743
901
(yt)-.15 G .243(ime. Note that)-2.743 F(the command 'TR)133 705.6 Q
902
(UE' has no ar)-.4 E(guments of its o)-.18 E(wn!)-.25 E(1)535 768 Q 0 Cg
903
EP
904
%%Page: 2 16
905
%%BeginPageSetup
906
BP
907
%%EndPageSetup
908
/F0 10/Times-Roman@0 SF 406.64(select\(\) select\(\))72 48 R/F1 10
909
/Times-Bold@0 SF 2.5(-F)108 84 S(ALSE)-3.4 E F0 1.925(check whether v)
910
133 96 R 1.925(alue is equal to zero \(logical F)-.25 F 1.926
911
(ALSE\): e.g.)-.74 F F1 -.9(FA)6.926 G(LSE:CYCL::ALL\(ALL\)).9 E F0
912
(checks)4.426 E .9(whether the trajectory ne)133 108 R -.15(ve)-.25 G
913
3.4(rp).15 G .9(asses through a c)-3.4 F .9(ycclone, which is mark)-.15
914
F .9(ed as 0/1 \214eld. Note that the)-.1 F(command 'TR)133 120 Q
915
(UE' has no ar)-.4 E(guments of its o)-.18 E(wn!)-.25 E F1 2.5(-A)108
916
136.8 S(LL, ANY)-2.5 E 2.5(,N)-.92 G(ONE)-2.5 E F0 .167(these are speci\
917
al commands which only apply for the TRIGGER \214eld. Further e)133
918
148.8 R .168(xplanations are gi)-.15 F -.15(ve)-.25 G(n).15 E(belo)133
919
160.8 Q 2.5(wi)-.25 G 2.5(ns)-2.5 G(ection TRIGGER.)-2.5 E/F2 10.95
920
/Times-Bold@0 SF(FIELD)72 177.6 Q F1 2.5(-V)108 189.6 S(ALUE)-3.85 E F0
921
(tak)133 201.6 Q 3.587(ev)-.1 G 1.087(alue of the \214eld: e.g.)-3.837 F
922
F1(GT)6.087 E(:PV\(V)-.74 E(ALUE\):2)-1.35 E F0 1.087
923
(selects the trajectories with \214rst potential v)3.587 F(orticity)-.2
924
E 1.073(\(PV\) v)133 213.6 R 1.073
925
(alue greater than 2 PVU. This selection criterion is equi)-.25 F -.25
926
(va)-.25 G 1.074(lent to).25 F F1(GT)3.574 E(:PV)-.74 E(:2,)-.92 E F0
927
1.074(i.e. the V)3.574 F(ALUE)-1.35 E(ar)133 225.6 Q(gument is tak)-.18
928
E(en as def)-.1 E(ault.)-.1 E F1 2.5(-M)108 242.4 S(EAN)-2.5 E F0(tak)
929
133 254.4 Q 4.505(et)-.1 G 2.005(he mean o)-4.505 F -.15(ve)-.15 G 4.504
930
(rt).15 G 2.004(he selected times: e.g.)-4.504 F F1(GT)7.004 E
931
(:RH\(MEAN\):70:ALL)-.74 E F0 2.004(selects all trajectories for)4.504 F
932
(which the mean relati)133 266.4 Q .3 -.15(ve h)-.25 H(umidity \(RH\) o)
933
.15 E -.15(ve)-.15 G 2.5(ra).15 G
934
(ll times \(ALL\) is greater than 70 %.)-2.5 E F1 2.5(-V)108 283.2 S(AR)
935
-3.85 E F0(tak)133 295.2 Q 4.685(et)-.1 G 2.185(he v)-4.685 F 2.185
936
(ariance o)-.25 F -.15(ve)-.15 G 4.685(rt).15 G 2.185
937
(he selected times: e.g.)-4.685 F F1(GT)7.185 E(:lat\(V)-.74 E
938
(AR\):10:ALL)-1.35 E F0 2.185(selects all trajectories for)4.685 F
939
(which the v)133 307.2 Q(ariance of latitude \(lat\) o)-.25 E -.15(ve)
940
-.15 G 2.5(ra).15 G(ll times \(ALL\) is greater than 10.)-2.5 E F1 2.5
941
(-M)108 324 S(IN)-2.5 E F0(tak)133 336 Q 2.693(et)-.1 G .193
942
(he minimum of the selected times: e.g.)-2.693 F F1 -1.4 -.92(LT :)5.192
943
H(p\(MIN\):300:ALL).92 E F0 .192(select all trajectories which ha)2.692
944
F -.15(ve)-.2 G 2.5(am)133 348 S(inium pressure \(p\) less than 300 hP)
945
-2.5 E 2.5(ao)-.15 G -.15(ve)-2.65 G 2.5(ra).15 G(ll times \(ALL\).)-2.5
946
E F1 2.5(-M)108 364.8 S(AX)-2.5 E F0(tak)133 376.8 Q 4.045(et)-.1 G
947
1.545(he maximum of the selected times: e.g.)-4.045 F F1 -1.4 -.92(LT :)
948
6.545 H(p\(MAX\):300:ALL).92 E F0 1.545(select all trajectories which)
949
4.045 F(ha)133 388.8 Q .3 -.15(ve a m)-.2 H
950
(aximum pressure \(p\) less than 300 hP).15 E 2.5(ao)-.15 G -.15(ve)
951
-2.65 G 2.5(ra).15 G(ll times \(ALL\).)-2.5 E F1 2.5(-S)108 405.6 S(UM)
952
-2.5 E F0(tak)133 417.6 Q 3.174(et)-.1 G .674(he sum o)-3.174 F -.15(ve)
953
-.15 G 3.174(rt).15 G .674(he selected times: e.g.)-3.174 F F1(GT)5.674
954
E(:LHR\(SUM\):2:ALL)-.74 E F0 .674(selects all trajectories for which)
955
3.174 F(the sum o)133 429.6 Q -.15(ve)-.15 G 2.5(ra).15 G
956
(ll latent heating rates \(LHR\) o)-2.5 E -.15(ve)-.15 G 2.5(ra).15 G
957
(ll times \(ALL\) is greater than 2 K.)-2.5 E F1 2.5(-C)108 446.4 S
958
(HANGE)-2.5 E F0(tak)133 458.4 Q 4.364(ec)-.1 G 1.864(hange between tw)
959
-4.364 F 4.364(ot)-.1 G 1.864(imes: e.g.)-4.364 F F1(GT)6.865 E
960
(:p\(CHANGE\):600:FIRST)-.74 E(,LAST)-.74 E F0 1.865
961
(selects all trajectories)4.365 F .633(wich ha)133 470.4 R .933 -.15
962
(ve a p)-.2 H .632(ressure dif).15 F .632
963
(ference |p\(FIRST\)-p\(LAST\)| greater than 600 hP)-.25 F 3.132(ab)-.15
964
G .632(etween the \214rst and last)-3.132 F .8(time or vice v)133 482.4
965
R .8(ersa. Note that the change can be positi)-.15 F 1.1 -.15(ve o)-.25
966
H 3.3(rn).15 G -2.25 -.15(eg a)-3.3 H(ti).15 E -.15(ve)-.25 G 3.3(,i).15
967
G .8(.e. it is not clear whether it is)-3.3 F(ascent or descent.)133
968
494.4 Q F1 2.5(-D)108 511.2 S(IFF)-2.5 E F0(tak)133 523.2 Q 2.844(ed)-.1
969
G(if)-2.844 E .344(ference between tw)-.25 F 2.844(ot)-.1 G .344
970
(imes: e.g.)-2.844 F F1(GT)5.344 E(:p\(DIFF\):600:FIRST)-.74 E(,LAST)
971
-.74 E F0 .344(selects all trajectories wich)2.844 F(ha)133 535.2 Q .878
972
-.15(ve a p)-.2 H .578(ressure dif).15 F .579
973
(ference p\(FIRST\)-p\(LAST\) greater than 600 hP)-.25 F 3.079(ab)-.15 G
974
.579(etween the \214rst and last time -)-3.079 F .854
975
(corresponding to an ascending air stream. Correspondingly)133 547.2 R
976
F1(GT)3.353 E(:p\(DIFF\):600:LAST)-.74 E(,FIRST)-.74 E F0 .853
977
(\214nds a)3.353 F(descending air stream.)133 559.2 Q F2(TIME MODE)72
978
576 Q F0 1.394(The command are applied to a set of trajectory times; if\
979
 no time is speci\214ed, the the command is only)108 588 R .509
980
(applied to the \214rst time. Most generally)108 600 R 3.009(,t)-.65 G
981
.509(he time mode consists of tw)-3.009 F 3.009(op)-.1 G .509
982
(arts: time list\( time mode\), where)-3.009 F
983
(the \214rst speci\214es a list of times and the second ho)108 612 Q 2.5
984
(wt)-.25 G 2.5(oa)-2.5 G(pply the criterion to the selected times.)-2.5
985
E F1 2.5(-F)108 628.8 S(IRST)-2.5 E F0 .542(\214rst time: e.g.)133 640.8
986
R F1(IN:lat:-20,20:FIRST)5.542 E F0 .543
987
(selects all trajectories with \214rst latitude between 20 S to 20 N,)
988
3.042 F(i.e. wich start in an equatorial band.)133 652.8 Q F1 2.5(-L)108
989
669.6 S(AST)-2.5 E F0 .006(last time: e.g.)133 681.6 R F1
990
(IN:lat:-20,20:LAST)5.006 E F0 .006
991
(selects all trajectories with last latitude between 20 S to 20 N, i.e.)
992
2.506 F(which end in an equatorial band.)133 693.6 Q(2)535 768 Q 0 Cg EP
993
%%Page: 3 17
994
%%BeginPageSetup
995
BP
996
%%EndPageSetup
997
/F0 10/Times-Roman@0 SF 406.64(select\(\) select\(\))72 48 R/F1 10
998
/Times-Bold@0 SF 2.5(-T)108 84 S(1,T2,T3)-2.5 E F0 1.447(an e)133 96 R
999
1.448(xplict list of times: e.g.)-.15 F F1(IN:lat:-20,20:6,12)6.448 E F0
1000
1.448(selects all trajectories which are in the equatorial)3.948 F(band\
1001
 at times 6 h and 12 h. The criterion must apply at both times \(see be\
1002
lo)133 108 Q 2.5(wA)-.25 G(LL, ANY)-2.5 E 2.5(,N)-1.29 G(ONE\).)-2.5 E
1003
F1 2.5(-T)108 124.8 S 2.5(1t)-2.5 G 2.5(oT)-2.5 G(2)-2.5 E F0 2.805(at)
1004
133 136.8 S .305(ime range: e.g.)-2.805 F F1 .305(IN:lat:-20,20:6 to 18)
1005
5.305 F F0 .304
1006
(selects all trajectories which are in the equatorial band from)2.804 F
1007
3.874(6ht)133 148.8 S 3.874(o1)-3.874 G 3.874(8h)-3.874 G 6.374(.T)
1008
-3.874 G 1.375
1009
(he criterion must apply at all times between 6 h and 18 h \(see belo)
1010
-6.374 F 3.875(wA)-.25 G 1.375(LL, ANY)-3.875 F(,)-1.29 E(NONE\).)133
1011
160.8 Q F1 2.5(-A)108 177.6 S(LL)-2.5 E F0 1.535(all times: e.g.)133
1012
189.6 R F1(IN:lat:-20,20:ALL)6.535 E F0 1.535
1013
(selects all trajectories which stay at all times in the equatorial)
1014
4.035 F 1.084(band. This time mode is the same as)133 201.6 R F1
1015
(ALL\(ALL\),)3.584 E F0 1.084
1016
(i.e. all times are selected and the criterion must)3.584 F .373
1017
(apply to all times. W)133 213.6 R(ith)-.4 E F1
1018
(IN:lat:-20,20:12-24\(ALL\))2.873 E F0 .373
1019
(the criterion must apply for all times between 12)2.873 F 2.5(ha)133
1020
225.6 S(nd 24 h.)-2.5 E F1 2.5(-A)108 242.4 S(NY)-2.5 E F0(an)133 254.4
1021
Q 2.669(yt)-.15 G .169(imes: e.g.)-2.669 F F1(IN:lat:-20,20:ALL\(ANY\))
1022
5.169 E F0 .169(selects all trajectories which stay at an)2.669 F 2.67
1023
(yt)-.15 G .17(imes in the equa-)-2.67 F .526(torial band. Note that wi\
1024
th the \214rst "ALL" the times are selected, i.e. all times in this cas\
1025
e, and with)133 266.4 R 1.296(the second "ANY" it is speci\214ed that t\
1026
he criterion must only apply to at least one of the selected)133 278.4 R
1027
(times.)133 290.4 Q F1 2.5(-N)108 307.2 S(ONE)-2.5 E F0 .169
1028
(at no time: e.g.)133 319.2 R F1(IN:lat:-20,20:ALL\(NONE\))5.169 E F0
1029
.168(selects all trajectories which ne)2.669 F -.15(ve)-.25 G 2.668(rs)
1030
.15 G .168(tay in the equatorial)-2.668 F(band.)133 331.2 Q F1(OUT)6.086
1031
E(:lat:-20,20:FIRST\(NONE\))-.74 E F0 1.086
1032
(selects the trajectories which are not outside the equatorial)3.586 F
1033
(time at the \214rst time: the)133 343.2 Q 2.5(ym)-.15 G(ust be inside.)
1034
-2.5 E F1 2.5(-T)108 360 S(RIGGER)-2.5 E F0 .559
1035
(the trajectory is automatically selected, b)133 372 R .558
1036
(ut the trigger column is updated. A selection of trajectories)-.2 F
1037
.348(might then e accomplished based on this trigger column: e.g)133 384
1038
R F1(GT)2.849 E(:p:700:1\(TRIGGER\))-.74 E F0 -.1(wo)2.849 G .349
1039
(uld set the).1 F 3.432(trigger 1 for all trajectory times where the pr\
1040
essure \(p\) is greater than 700 hP)133 396 R 3.432(a. Similarly)-.15 F
1041
(,)-.65 E F1(GT)133 408 Q .162(:p:800:1\(TRIGGER\) & GT)-.74 F
1042
(:lat:50:2\(TRIGGER\))-.74 E F0 -.1(wo)2.662 G .163
1043
(uld set the trigger 1 for all times where the).1 F .463
1044
(pressure \(p\) is lar)133 420 R .463(ger than 800 hP)-.18 F 2.963(aa)
1045
-.15 G .463
1046
(nd set the trigger 2 for all times where the latitude \(lat\) is lar)
1047
-2.963 F(ger)-.18 E .749(than 50 de)133 432 R .749(grees north. Note, t\
1048
hat both eents might apply: then the rsulting trigger is 3. The trigger\
1049
s)-.15 F .128(are internally sa)133 444 R -.15(ve)-.2 G 2.628(da).15 G
1050
2.628(st)-2.628 G .128(he bits of an inte)-2.628 F .128(ger v)-.15 F
1051
.127(ariable, i.e. trigger 1 corresponds to v)-.25 F .127
1052
(alue 1=2^0, trigger)-.25 F 2.5(2t)133 456 S 2.5(o2)-2.5 G
1053
(=2^1, trigger 3 to 4=2^2...)-2.5 E/F2 10.95/Times-Bold@0 SF
1054
(LOGICAL OPERA)72 472.8 Q -.197(TO)-1.04 G(RS).197 E F1 2.5(-&)108 484.8
1055
S F0 .112(logical and: e.g.)133 484.8 R F1(GT)5.112 E .112
1056
(:lat:34:FIRST & GT)-.74 F(:lon:50:FIRST)-.74 E F0 .112
1057
(selects the trajectories to the north of 34 N)2.612 F
1058
(and to the east of 50 E at \214rst time. Se)133 496.8 Q -.15(ve)-.25 G
1059
(ral selection criteria can be combined with '&'.).15 E F1 2.5(-|)108
1060
513.6 S F0 .199(logical or: e.g.)133 513.6 R F1(GT)5.199 E .199
1061
(:lat:34:FIRST | GT)-.74 F(:lon:50:FIRST)-.74 E F0 .198
1062
(selects the trajectories to the north of 34 N or)2.699 F .536
1063
(to the east of 50 E at \214rst time. Se)133 525.6 R -.15(ve)-.25 G .537
1064
(ral selection criteria can be combined with '|'. Note that locical).15
1065
F .188(OR has a lo)133 537.6 R .188
1066
(wer priority than logical AND, i.e. in an e)-.25 F .187(xpression lik)
1067
-.15 F 2.687(eT)-.1 G 2.687(1|T)-2.687 G 2.687(2&T)-2.687 G 2.687<338c>
1068
-2.687 G .187(rst the e)-2.687 F(xpression)-.15 E(T2 & T3 is e)133 549.6
1069
Q -.25(va)-.25 G(luated and only then logically OR-combined with T1.).25
1070
E F2(IMPLICIT FIELDS)72 566.4 Q F0 .448(Implicit v)108 578.4 R .448
1071
(ariables can be used in the selection criteria, although the)-.25 F
1072
2.948(yd)-.15 G 2.948(on)-2.948 G .448(ot e)-2.948 F .449
1073
(xplicitely appear as a column)-.15 F(in the trajectory \214le. The)108
1074
590.4 Q 2.5(ya)-.15 G(re calculated on-the-\215y during the selection.)
1075
-2.5 E F1 2.5(-D)108 607.2 S(IST)-2.5 E F0 .31
1076
(length of the trajectory \(in km\), inte)133 619.2 R .31
1077
(grated along great circle sections between the trajectory v)-.15 F
1078
(ertices:)-.15 E(e.g.)133 631.2 Q F1(GT)5.283 E(:DIST)-.74 E(:1000:LAST)
1079
-.74 E F0 .283
1080
(tests whether the total path length of the trajectory \(DIST\) e)2.783
1081
F .284(xceeds 1000)-.15 F(km.)133 643.2 Q F1 2.5(-D)108 660 S(IST0)-2.5
1082
E F0 .555
1083
(distance of the trajectory \(in km\) from its starting position: e.g.)
1084
133 672 R F1(GT)5.555 E(:DIST0:1000:18,24\(ANY\))-.74 E F0(tests)3.055 E
1085
.379(whether the air parcel is more than 1000 km a)133 684 R -.1(wa)-.15
1086
G 2.879(yf).1 G .379
1087
(rom its starting position, either at time 18 h or at)-2.879 F
1088
(time 24 h.)133 696 Q(3)535 768 Q 0 Cg EP
1089
%%Page: 4 18
1090
%%BeginPageSetup
1091
BP
1092
%%EndPageSetup
1093
/F0 10/Times-Roman@0 SF 406.64(select\(\) select\(\))72 48 R/F1 10
1094
/Times-Bold@0 SF 2.5(-I)108 84 S(NPOL)-2.5 E(YGON)-.92 E F0 .86
1095
(speci\214cation of a polygon: e.g.)133 96 R F1(TR)5.86 E(UE:INPOL)-.3 E
1096
(YGON:\214lename:ALL\(ANY\))-.92 E F0 .86(selects all trajectories)3.36
1097
F .123(which pass through the polygon speci\214ed in "\214lename". The \
1098
polygon is speci\214ed as described in)133 108 R F1(cr)2.624 E(e-)-.18 E
1099
(ate_startf)133 120 Q F0 .77(\(see comments there\). W)3.27 F(ith)-.4 E
1100
F1 -.9(FA)3.27 G(LSE:INPOL).9 E(YGON:\214lename:ALL\(ALL\))-.92 E F0 .77
1101
(all trajectories)3.27 F .627(are selected which ne)133 132 R -.15(ve)
1102
-.25 G 3.127(rp).15 G .628(ass through the polygon. Note that for e)
1103
-3.127 F -.15(ve)-.25 G .628(ry call of).15 F F1(select)3.128 E F0 .628
1104
(only one poly-)3.128 F(gon can be used in the criteria!)133 144 Q F1
1105
2.5(-I)108 160.8 S(NBO)-2.5 E(X)-.4 E F0 20.002
1106
(speci\214cation of a longitude/latitude rectangle <lonw)133 172.8 R
1107
20.001(,lone,lats,latn>: e.g.)-.65 F F1(TR)133 184.8 Q(UE:INBO)-.3 E
1108
(X:20,40,30,60:ALL\(ANY\))-.4 E F0 .1
1109
(selects all trajectories which pass through the longitude/lat-)2.6 F
1110
.158(itude rectangle with lo)133 196.8 R(wer)-.25 E .157
1111
(-left corner at 20 E / 30 N and the upper)-.2 F .157
1112
(-right corner at 40 E / 60 N. Corre-)-.2 F(spondingly)133 208.8 Q 2.767
1113
(,w)-.65 G(ith)-2.767 E F1(TR)2.767 E(UE:INBO)-.3 E
1114
(X:20,40,30,60:ALL\(NONE\))-.4 E F0 .268
1115
(the trajectories are selected which ne)2.767 F -.15(ve)-.25 G(r).15 E
1116
(pass through the rectangle.)133 220.8 Q F1 2.5(-I)108 237.6 S(NCIRCLE)
1117
-2.5 E F0 .651
1118
(speci\214cation of circle around a speci\214ed point: e.g.)133 249.6 R
1119
F1(TR)5.651 E(UE:INCIRCLE:40,50,500:LAST)-.3 E F0 .651(select all)3.151
1120
F 1.077(trajectories which ha)133 261.6 R 1.377 -.15(ve t)-.2 H 1.077(h\
1121
eir \214nal position \(LAST\) in the circle centred at 40 E / 50 N and \
1122
with a).15 F(radius of 500 km.)133 273.6 Q F1 2.5(-I)108 290.4 S
1123
(NREGION)-2.5 E F0 .412(speci\214cation of tar)133 302.4 R .412(get re)
1124
-.18 F .412(gions in a re)-.15 F .411(gion \214le \(def)-.15 F .411
1125
(ault "re)-.1 F .411(gionf"\) - please consider the documentation)-.15 F
1126
.521(of "create_startf" for details concerning the format of the re)133
1127
314.4 R .521(gion \214le. As an e)-.15 F .521(xample, if a re)-.15 F
1128
.522(gion 1 is)-.15 F 1.39(de\214ned on the re)133 326.4 R 1.39
1129
(gion \214le, the criterion)-.15 F F1(TR)3.89 E(UE:INREGION:1:18)-.3 E
1130
F0 1.39(selects all trejactories which are)3.89 F(within re)133 338.4 Q
1131
(gion 1 at time 18 h.)-.15 E/F2 10.95/Times-Bold@0 SF(SPECIAL CRITERIA)
1132
72 355.2 Q F0 .624(Special criteria are and can be implemented into)108
1133
367.2 R F1(select.)3.124 E F0 .624
1134
(The call to the special criteria is of the follo)5.624 F(wing)-.25 E
1135
(form:)108 379.2 Q F1(SPECIAL:command:parameters,)2.738 E F0 .237
1136
(where "command" is a command string \(e.g. WCB\) and "param-)2.738 F
1137
(eters" is a list of parameter v)108 391.2 Q(alues.)-.25 E F1 2.5(-W)108
1138
408 S(CB)-2.5 E F0 1.195(identi\214cation of W)133 420 R 1.195(arm Con)
1139
-.8 F -.15(vey)-.4 G 1.196(or Belts \(WCB\): e.g.).15 F F1
1140
(SPECIAL:WCB:300,0,24)6.196 E F0 1.196(selects trajectories)3.696 F
1141
2.239(which ascend more than 300 hP)133 432 R 4.738(ab)-.15 G 2.238
1142
(etween time 0 and 24 h. Note, the ascent is determined as)-4.738 F .628
1143
(min{p\(0...24\)}-p\(0\), i.e. the \214rst pressure p\(0\) is \214x)133
1144
444 R .629(ed whereas the lo)-.15 F .629(west pressure can occur at an)
1145
-.25 F(y)-.15 E(time between 0...24 h.)133 456 Q F2(TRIGGER FIELD)72
1146
472.8 Q F0 3.324(At)108 484.8 S .824(rigger \(or \215ag\) \214eld can b\
1147
e de\214ned in select. This special column of the trajectory \214le all\
1148
o)-3.324 F .823(ws to mark)-.25 F .934(speci\214ed e)108 496.8 R -.15
1149
(ve)-.25 G .934(nts for each trajectory and its times. As an e).15 F
1150
.934(xample, you w)-.15 F .934(ould lik)-.1 F 3.434(et)-.1 G 3.434(os)
1151
-3.434 G .935(elect all trajectories)-3.434 F .205(which are belo)108
1152
508.8 R 2.705(w7)-.25 G .205(00 hP)-2.705 F 2.705(aa)-.15 G 2.705(tac)
1153
-2.705 G .205(ertain time and ha)-2.705 F .505 -.15(ve r)-.2 H(elati).15
1154
E .505 -.15(ve h)-.25 H .205(umidity abo).15 F .505 -.15(ve 8)-.15 H
1155
2.705(0%).15 G F1 .205(at the same times:)B F0(Then)2.705 E .749(you co\
1156
uld set a \214rst trigger for the pressure criterion and a second one f\
1157
or the relati)108 520.8 R 1.05 -.15(ve h)-.25 H(umidity).15 E 3.25(,a)
1158
-.65 G .75(nd then)-3.25 F(you w)108 532.8 Q
1159
(ould check for the simultaneous occurence of the tw)-.1 E 2.5(ot)-.1 G
1160
(riggers. More speci\214cally)-2.5 E(,)-.65 E F1(GT)108 549.6 Q
1161
(:p:700:1\(TRIGGER\) & GT)-.74 E(:RH:80:2\(TRIGGER\))-.74 E F0 .256
1162
(this will de\214ne the triggers for the tw)133 561.6 R 2.756(oe)-.1 G
1163
-.15(ve)-3.006 G .256(nts. Note that both e).15 F -.15(ve)-.25 G .255
1164
(nts might be ful\214lled, in which case).15 F .082
1165
(both triggers are set. The v)133 573.6 R .083(alue of the trigger is: \
1166
1 if only the \214rst criterion is ful\214lled \(binary 01\); 2 if)-.25
1167
F(only the second is ful\214lled \(binary 10\) and 3 if both are ful\
1168
\214lled \(binary 11\).)133 585.6 Q F1(ALL:TRIGGER:1,2:ALL\(ANY\))108
1169
602.4 Q F0 .462(this will select all trajectories for which both trigge\
1170
rs 1 and 2 are set - the)133 614.4 R 2.961(yc)-.15 G .461
1171
(an be set at an)-2.961 F 2.961(yt)-.15 G .461(ime of)-2.961 F .519
1172
(all the trajectory times; if, for instance, you w)133 626.4 R .519
1173
(ould lik)-.1 F 3.019(et)-.1 G 3.019(ot)-3.019 G .519
1174
(est whether the tw)-3.019 F 3.019(ot)-.1 G .519(riggers are set at the)
1175
-3.019 F(tw)133 638.4 Q 2.5(ot)-.1 G(imes 6 and 12, the command w)-2.5 E
1176
(ould be)-.1 E F1(ALL:TRIGGER:1,2:6,12\(ALL\).)2.5 E
1177
(ALL:TRIGGER:1,2:ALL\(ANY\))108 655.2 Q F0
1178
(will only check whether one of the tw)133 667.2 Q 2.5(ot)-.1 G
1179
(riggers 1 and 2 is set)-2.5 E F1(NONE:TRIGGER:1,2:ALL\(ANY\))108 684 Q
1180
F0(will check whether none of the tw)133 696 Q 2.5(ot)-.1 G
1181
(riggers 1 and 2 is set)-2.5 E(4)535 768 Q 0 Cg EP
1182
%%Page: 5 19
1183
%%BeginPageSetup
1184
BP
1185
%%EndPageSetup
1186
/F0 10/Times-Roman@0 SF 406.64(select\(\) select\(\))72 48 R/F1 10.95
1187
/Times-Bold@0 SF(IMPLEMENTING COMPLEX CRITERIA)72 84 Q F0(Ne)108 96 Q
1188
3.969(ws)-.25 G 1.469(pecial criteria can easily be implemented into th\
1189
e code - to this aim the follo)-3.969 F 1.468(wing steps must be)-.25 F
1190
(tak)108 108 Q(en:)-.1 E/F2 10/Times-Bold@0 SF([1] <special.f>)108 124.8
1191
Q F0 .957(in directory "${LA)133 136.8 R(GRANT)-.4 E .957
1192
(O}/select/" must be modi\214ed. T)-.18 F(ak)-.8 E 3.458(et)-.1 G .958
1193
(he e)-3.458 F .958(xample "WCB" to see ho)-.15 F 3.458(wt)-.25 G(he)
1194
-3.458 E -.15(Fo)133 148.8 S(rtran code must be changed.).15 E F2
1195
([2] <install.sh select>)108 165.6 Q F0 2.458(must be in)133 177.6 R -.2
1196
(vo)-.4 G -.1(ke).2 G 4.958(dt).1 G 4.958(or)-4.958 G 2.458
1197
(ecompile the program "select". F)-4.958 F 2.457
1198
(or successful compilation, the e)-.15 F -.15(xe)-.15 G(cutable).15 E
1199
("select" will be listed.)133 189.6 Q -.8(Ty)108 206.4 S(pe <select -sp\
1200
ecial> to get a list of all special commands and a detailed description\
1201
.).8 E F1(EXAMPLES)72 223.2 Q F2([1] select trainp traout 'GT)108 235.2
1202
Q(:PV)-.74 E(:2:LAST')-.92 E F0 1.762(selects all trajectories with PV>\
1203
2 PVU for the last time step. The input trajectories are gi)133 247.2 R
1204
-.15(ve)-.25 G 4.263(ni).15 G(n)-4.263 E .229
1205
(trainp, the selected ones are written to traout. If the tw)133 259.2 R
1206
2.728<6f8c>-.1 G .228(lenames are the same, the input \214le is o)-2.728
1207
F -.15(ve)-.15 G -.2(r-).15 G(written.)133 271.2 Q F2
1208
([2] select trainp traout 'IN:lat:-20,20:6,12')108 288 Q F0 .701(select\
1209
s all trajectories which are in the equatorial band at times 6 h and)
1210
135.5 300 R .701(12 h.)5.701 F 3.202(The criterion must)5.701 F
1211
(apply at both times.)133 312 Q F2([3] select trainp traout 'GT)108
1212
328.8 Q(:lat:34:FIRST & GT)-.74 E(:lon:50:FIRST')-.74 E F0 2.763
1213
(selects the)133 340.8 R .263(trajectories to the north of 34 N and to \
1214
the east of 50 E at \214rst time. Se)2.763 F -.15(ve)-.25 G .262
1215
(ral selection cri-).15 F(teria can be combined with '&'.)133 352.8 Q F2
1216
([4] select trainp traout 'L)108 369.6 Q -.74(T:)-.92 G
1217
(p\(MIN\):300:ALL).74 E(')-1.1 E F0(select all trajectories which ha)133
1218
381.6 Q .3 -.15(ve a m)-.2 H(inium pressure \(p\) less than 300 hP).15 E
1219
2.5(ao)-.15 G -.15(ve)-2.65 G 2.5(ra).15 G(ll times \(ALL\).)-2.5 E F1
1220
-.548(AU)72 398.4 S(THOR).548 E F0
1221
(Written by Michael Sprenger and Heini W)108 410.4 Q
1222
(ernli \(January 2011\))-.8 E(5)535 768 Q 0 Cg EP
1223
%%Page: 1 20
1224
%%BeginPageSetup
1225
BP
1226
%%EndPageSetup
1227
/F0 10/Times-Roman@0 SF 393.3(timeres\(\) timeres\(\))72 48 R/F1 10.95
1228
/Times-Bold@0 SF -.219(NA)72 84 S(ME).219 E/F2 10/Times-Bold@0 SF(timer)
1229
108 96 Q(es - change the time r)-.18 E(esolution of a trajectory)-.18 E
1230
F1(SYNOPSIS)72 112.8 Q F2(timer)108 124.8 Q(es)-.18 E/F3 10
1231
/Times-Italic@0 SF .41(inp\214le out\214le)2.51 F(-[h|min] value)2.69 E
1232
([-cubic|-linear])2.97 E F1(DESCRIPTION)72 141.6 Q F0 .294(Change the t\
1233
ime resolution of an input trajectory \214le "inp\214le" through interp\
1234
olation and write a ne)108 153.6 R 2.794(wt)-.25 G(rajec-)-2.794 E .664
1235
(tory \214le "out\214le". The ne)108 165.6 R 3.164(wt)-.25 G .664
1236
(ime resolution "v)-3.164 F .664(alue" is gi)-.25 F -.15(ve)-.25 G 3.164
1237
(ne).15 G .664(ither in hours "-h" or in minutes "-min". The)-3.164 F
1238
1.182(interpolation is performed on the trajectory \214le either in lin\
1239
ear mode \("-linear"\) or in cubic spline mode)108 177.6 R
1240
(\("-cubic"\). The def)108 189.6 Q(ault is "-cubic".)-.1 E F1 -.81(PA)72
1241
206.4 S(RAMETERS).81 E F3(inp\214le)108.01 218.4 Q F0
1242
(input trajectory \214le)183 218.4 Q F3(out\214le)108.23 235.2 Q F0
1243
(output trajectory \214le \(can be the same as inp\214le\).)183 235.2 Q
1244
F3(-h value)108.01 252 Q F0(ne)183 252 Q 2.5(wt)-.25 G
1245
(ime resolution in hours \(e.g. "-h 1"\).)-2.5 E F3(-min value)108.01
1246
268.8 Q F0(ne)183 268.8 Q 2.5(wt)-.25 G
1247
(ime resolution in minutes \(e.g. "-min 15"\).)-2.5 E F3(-linear)108.01
1248
285.6 Q F0 .423(linear interpolation between tw)183 285.6 R 2.923(ot)-.1
1249
G .422(rajectory times; note that this mode conserv)-2.923 F .422
1250
(es the sign)-.15 F 2.511(between tw)183 297.6 R 5.011(ot)-.1 G 2.512
1251
(rajectory times - possibly of important for speci\214c humidity)-5.011
1252
F 5.012(,r)-.65 G(elati)-5.012 E -.15(ve)-.25 G(humidity)183 309.6 Q
1253
(,...)-.65 E F3(-linear)108.01 326.4 Q F0 1.294
1254
(subic spline interpolation between tw)183 326.4 R 3.794(ot)-.1 G 1.293
1255
(rajectory times; note that this mode conserv)-3.794 F(es)-.15 E
1256
(can change the sign between tw)183 338.4 Q 2.5(ot)-.1 G
1257
(rajectory times!)-2.5 E F1(EXAMPLES)72 355.2 Q F2([1] timer)108 367.2 Q
1258
(es tra\214le tra\214le -min 15 -linear)-.18 E F0 .008(changes the time\
1259
 resolution to 15 minutes, using linear interpolation, and o)133 379.2 R
1260
-.15(ve)-.15 G .009(rwrites the old trajectory).15 F(\214le.)133 391.2 Q
1261
F1 -.548(AU)72 408 S(THOR).548 E F0
1262
(Written by Michael Sprenger and Heini W)108 420 Q(ernli \(January 2011)
1263
-.8 E(1)535 768 Q 0 Cg EP
1264
%%Page: 1 21
1265
%%BeginPageSetup
1266
BP
1267
%%EndPageSetup
1268
/F0 10/Times-Roman@0 SF 407.76(tracal\(\) tracal\(\))72 48 R/F1 10.95
1269
/Times-Bold@0 SF -.219(NA)72 84 S(ME).219 E/F2 10/Times-Bold@0 SF
1270
(tracal - simple calculations with trajectory \214les)108 96 Q F1
1271
(SYNOPSIS)72 112.8 Q F2(tracal)108 124.8 Q/F3 10/Times-Italic@0 SF .41
1272
(inp\214le out\214le)2.51 F -.2(ex)2.68 G(pr).2 E(ession)-.37 E F1
1273
(DESCRIPTION)72 141.6 Q F0 .398
1274
(Simple calculation on the input trajectory \214le)108 153.6 R F3
1275
(inp\214le)2.907 E F0 2.897(;t)3.077 G .397
1276
(he output is written to the output trajectory \214le)-2.897 F F3(out-)
1277
3.127 E(\214le)108 165.6 Q F0(and the e)2.68 E
1278
(xpression for the calculation is gi)-.15 E -.15(ve)-.25 G 2.5(ni).15 G
1279
(n)-2.5 E F3 -.2(ex)2.5 G(pr).2 E(ession.)-.37 E F1 -.81(PA)72 182.4 S
1280
(RAMETERS).81 E F3(inp\214le)108.01 194.4 Q F0(input trajectory \214le)
1281
183 194.4 Q F3(out\214le)108.23 211.2 Q F0
1282
(output trajectory \214le \(can be the same as inp\214le\).)183 211.2 Q
1283
F3 -.2(ex)108 228 S(pr).2 E(ession)-.37 E F0 .5(arithmetic e)183 228 R
1284
.5(xpression, e.g. 'DIFF=PS-p' to get dif)-.15 F .5
1285
(ference between surf)-.25 F .5(ace pressure \(PS\))-.1 F .832
1286
(and the pressure height of the trajectory \(p\). Note that the v)183
1287
240 R .831(ariable names must e)-.25 F(xactly)-.15 E .362
1288
(match to the column names in the trajectory \214le.)183 252 R .363
1289
(The mathematical e)5.362 F .363(xpression can con-)-.15 F .045
1290
(tain pre)183 264 R .045(viously de\214ned v)-.25 F .045(ariables, numb\
1291
ers \(e.g., 1.5 or 1.5e-6\), arithmetic operators \(+, -,)-.25 F .178(*\
1292
, /, ^\), and a number of mathematical functions \(sin, cos, tan, sqrt,\
1293
 e)183 276 R .179(xp, log, ln, abs, ang,)-.15 F .236
1294
(real, imag, conjg, comple)183 288 R .236(x\). The e)-.15 F .235
1295
(xpression can also use nested le)-.15 F -.15(ve)-.25 G .235
1296
(ls of parentheses for).15 F .385(grouping. There are tw)183 300 R 2.885
1297
(op)-.1 G .385(rede\214ned v)-2.885 F .385(ariables that are a)-.25 F
1298
-.25(va)-.2 G .385(ilable to the user \342 the constant).25 F
1299
(pi=3.14159265358979 and the imaginary unit i.)183 312 Q F1(EXAMPLES)72
1300
328.8 Q F2([1] tracal inp out 'agl=z-zb')108 340.8 Q F0 .004
1301
(caculates the dif)133 352.8 R .004
1302
(ference between air parcel height \(z\) and surf)-.25 F .003
1303
(ace height \(zb\). The result is sa)-.1 F -.15(ve)-.2 G 2.503(di).15 G
1304
2.503(na)-2.503 G(n)-2.503 E(ne)133 364.8 Q 2.5(wc)-.25 G
1305
(olumn \(agl\).)-2.5 E F1 -.548(AU)72 381.6 S(THOR).548 E F0 .608
1306
(Written by Michael Sprenger and Heini W)108 393.6 R .609
1307
(ernli \(August 2012\); the e)-.8 F -.25(va)-.25 G .609
1308
(luation of the arithmetic e).25 F(xpression)-.15 E
1309
(is based on the string manipulation routines by 'Dr)108 405.6 Q 2.5(.G)
1310
-.55 G(eor)-2.5 E(ge Benthien' \(http://gbenthien.net/inde)-.18 E
1311
(x.html\).)-.15 E(1)535 768 Q 0 Cg EP
1312
%%Page: 1 22
1313
%%BeginPageSetup
1314
BP
1315
%%EndPageSetup
1316
/F0 10/Times-Roman@0 SF 413.32(trace\(\) trace\(\))72 48 R/F1 10.95
1317
/Times-Bold@0 SF -.219(NA)72 84 S(ME).219 E/F2 10/Times-Bold@0 SF
1318
(trace - trace meteor)108 96 Q(ological \214elds along trajectories)-.18
1319
E F1(SYNOPSIS)72 112.8 Q F2(trace)108 124.8 Q/F3 10/Times-Italic@0 SF
1320
.41(inp\214le out\214le)2.51 F F0([)2.68 E F3(optional ar)2.73 E
1321
(guments)-.37 E F0(])2.77 E F1(DESCRIPTION)72 141.6 Q F0 -.35(Tr)108
1322
153.6 S .91(ace meteorological \214elds along the trajectories gi).35 F
1323
-.15(ve)-.25 G 3.409(ni).15 G 3.409(nt)-3.409 G .909(he input \214le)
1324
-3.409 F F3(inp\214le)3.419 E F0 .909(and write a ne)3.589 F 3.409(wt)
1325
-.25 G(rajectory)-3.409 E(\214le)108 165.6 Q F3 .327(out\214le tr)2.877
1326
F .147(acing \214le)-.15 F F0(\(def)2.647 E .147(ault: trace)-.1 F -.25
1327
(va)-.25 G .147(rs\). P).25 F .147(artly the)-.15 F 2.647(yc)-.15 G .147
1328
(an be computed "online" \(see ONLINE CALCULA-)-2.647 F(TIONS belo)108
1329
177.6 Q(w\), normally the)-.25 E 2.5(ya)-.15 G(re a)-2.5 E -.25(va)-.2 G
1330
(ilbale on the primary and secondary P and S \214les.).25 E F1 -.81(PA)
1331
72 194.4 S(RAMETERS).81 E F3(inp\214le)108.01 206.4 Q F0
1332
(input trajectory \214le; the appendix determines the format \(see)183
1333
206.4 Q F2 -.18(re)2.5 G -.25(fo).18 G(rmat).25 E F0(for details\).)2.5
1334
E F3(out\214le)108.23 223.2 Q F0
1335
(output trajectory \214le; the appendix determines the format \(see)183
1336
223.2 Q F2 -.18(re)2.5 G -.25(fo).18 G(rmat).25 E F0(for details\).)2.5
1337
E F1(TRA)72 240 Q(CING FILE)-.602 E F0 .809(Normally the meteorological\
1338
 \214elds for tracing are listed in a \214le with name)108 252 R F2
1339
(trace)3.308 E -.1(va)-.15 G .808(rs .).1 F F0 .808(The name of the)
1340
5.808 F .317(tracing \214le can be changed with the optional ar)108 264
1341
R .317(gument "-v" \(see belo)-.18 F .317
1342
(w\). The format of the tracing \214le is as)-.25 F(follo)108 276 Q(ws:)
1343
-.25 E -.15(Fo)108 292.8 S(rmat).15 E F3 .77(\214eld[:shift] scale)
1344
134.91 304.8 R 1.49(computation pr)2.88 F(e\214x)-.37 E F0
1345
(Shifts \(optional\))108 321.6 Q F2 3.633<2d8c>133 333.6 S
1346
(eld:+100km[lat])-3.633 E F0 3.633(-g)3.633 G 1.133(et \214eld at traje\
1347
ctory position + 100 km shifted to north. A shift to south is)-3.633 F
1348
(obtained with \214eld:-100km[lat].)133 345.6 Q F2 4.003<2d8c>133 357.6
1349
S(eld:+100km[lon])-4.003 E F0 4.003(-g)4.003 G 1.503(et \214eld at traj\
1350
ectory position + 100 km shifted to east. A shift to west is)-4.003 F
1351
(obtained with \214eld:-100km[lon].)133 369.6 Q F2 2.507<2d8c>133 381.6
1352
S(eld:+2[dlat])-2.507 E F0 2.507(-g)2.507 G .006(et \214eld at trajecto\
1353
ry position + 2 grid spacings dlat shifted to north. A shift to south)
1354
-2.507 F(is obtained with \214eld:-2[dlat].)133 393.6 Q F2 2.711<2d8c>
1355
133 405.6 S(eld:+2[dlon])-2.711 E F0 2.711(-g)2.711 G .211(et \214eld a\
1356
t trajectory position + 2 grid spacings dlon shifted to east. A shift t\
1357
o west)-2.711 F(is obtained with \214eld:-2[dlon].)133 417.6 Q F2 2.683
1358
<2d8c>133 429.6 S(eld:+50hP)-2.683 E(a)-.1 E F0 2.683(-g)2.683 G .183
1359
(et \214eld at trajectory position + 50 hP)-2.683 F 2.683(as)-.15 G .183
1360
(hifted in v)-2.683 F .182(ertical. A shift to lo)-.15 F .182
1361
(wer pressures)-.25 F(is obtained with \214eld:-50hP)133 441.6 Q(a.)-.15
1362
E F2 2.966<2d8c>133 453.6 S(eld:+1dp)-2.966 E F0 2.966(-g)2.966 G .466
1363
(et \214eld at trajectory position + 1 grid spacing DP shifted in v)
1364
-2.966 F .467(ertical. A shift to lo)-.15 F(wer)-.25 E
1365
(pressures is obtained with \214eld:-1dp. Note that DP is not \214x)133
1366
465.6 Q(ed b)-.15 E(ut v)-.2 E(aries with height.)-.25 E F2 2.993<2d8c>
1367
133 477.6 S(eld:+6h)-2.993 E F0 2.993(-g)2.993 G .493
1368
(et \214eld at trajectory position, b)-2.993 F .493
1369
(ut 6 h in the future. Shifts to the past are poeeible with)-.2 F(\214e\
1370
ld:-6h. In addition to hours \(h\), the time shift can be speci\214ed i\
1371
n minutes \(min\).)133 489.6 Q(Examples)108 506.4 Q F2 2.664(-T)133
1372
518.4 S 2.664(H1)-2.664 G 2.664(.0S:)-2.664 G F0 .164(trace potential t\
1373
emperature \(TH\), scale it with 1 \(no scaling\); it is a)B -.25(va)-.2
1374
G .164(ilable on the S \214le).25 F(\(no computation is needed: 0\).)133
1375
530.4 Q F2 2.515(-Q1)133 542.4 S .015(000. 0 P :)-2.515 F F0 .014
1376
(trace speci\214c humidity \(Q\), scale it with 1000 to ha)2.515 F .314
1377
-.15(ve g)-.2 H .014(/kg; it is a).15 F -.25(va)-.2 G .014
1378
(ilable on the P \214le).25 F(\(no computation is needed: 0\).)133 554.4
1379
Q F2 2.556(-R)133 566.4 S 2.556(H1)-2.556 G -5.055 2.556(.1 * :)-2.556 H
1380
F0 .057(trace relati).001 F .357 -.15(ve h)-.25 H .057
1381
(umidity \(RH\), no scaling is needed \(1.\); relati).15 F .357 -.15
1382
(ve h)-.25 H .057(umidity is not a).15 F -.25(va)-.2 G(ilable).25 E
1383
(on either P or S \214le and must be computed \(1\).)133 578.4 Q F2
1384
3.549(-T)133 590.4 S(H:100hP)-3.549 E 3.549(a1)-.1 G 3.549(.0S:)-3.549 G
1385
F0 1.049(As in the \214rst e)B 1.048(xample, b)-.15 F 1.048(ut no)-.2 F
1386
3.548(wt)-.25 G 1.048(he potential temperature is tak)-3.548 F 1.048
1387
(en 100 hP)-.1 F(a)-.15 E(belo)133 602.4 Q 2.5(wt)-.25 G
1388
(he air parcel position.)-2.5 E F1(OPTION)72 619.2 Q(AL ARGUMENTS)-.219
1389
E F3(-i hour)108.01 631.2 Q(s)-.1 E F0 .859
1390
(time increments \(in hours\) for input P and)183 631.2 R 3.359<538c>
1391
5.859 G .86(les. If not e)-3.359 F .86(xplicitely speci\214ed,)-.15 F
1392
(this is determined from the P and S \214les i n the current directory)
1393
183 643.2 Q(.)-.65 E F3(-v var\214le)108.01 660 Q F0 2.998
1394
(Change the name of the tracing \214le from its def)183 660 R 2.997
1395
(ault v)-.1 F 2.997(alue "trace)-.25 F -.25(va)-.25 G 2.997(rs" to).25 F
1396
("v)183 672 Q(ar\214le".)-.25 E F3(-f \214eld scale)108.01 688.8 Q F0
1397
-.35(Tr)183 688.8 S .658(ace \214eld \(with scaling scale\) along the t\
1398
rajectories; the computation \215ag and).35 F 1.337(the pre\214x for th\
1399
e data \214le is automatically set. This options allo)183 700.8 R 1.337
1400
(ws the quick)-.25 F
1401
(tracing of a \214eld, without specifying a tracing \214le.)183 712.8 Q
1402
(1)535 768 Q 0 Cg EP
1403
%%Page: 2 23
1404
%%BeginPageSetup
1405
BP
1406
%%EndPageSetup
1407
/F0 10/Times-Roman@0 SF 413.32(trace\(\) trace\(\))72 48 R/F1 10
1408
/Times-Italic@0 SF(-c)108.01 84 Q(hang)-.15 E(et)-.1 E F0 1.037(\215ag \
1409
whether the times of the P and S \214les should be changed or not befor\
1410
e a)183 84 R(calculation; the def)183 96 Q(ault is that the times are)
1411
-.1 E/F2 10/Times-Bold@0 SF(not)2.5 E F0(changed.)2.5 E F1(-noclean)
1412
108.01 112.8 Q F0 1.645
1413
(\215ag whether parameter and criterion \214les should be k)183 112.8 R
1414
1.644(ept; this is particularly)-.1 F(helpful for deb)183 124.8 Q
1415
(ugging.)-.2 E F1(-timec)108.01 141.6 Q(hec)-.15 E(k)-.2 E F0
1416
(enforce a time check for the data \214le)183 141.6 Q/F3 10.95
1417
/Times-Bold@0 SF(SPECIAL INTERPOLA)72 158.4 Q(TION MODES)-1.04 E F1
1418
(-near)108.01 170.4 Q(est)-.37 E F0 1.391
1419
(Do no interpolation between grid points; just tak)183 170.4 R 3.891(et)
1420
-.1 G 1.391(he nearest neighbor! This)-3.891 F .132
1421
(option is useful, if a discrete input \214eld is gi)183 182.4 R -.15
1422
(ve)-.25 G 2.632(n\().15 G .132(e.g. labels\), where interpolated)-2.632
1423
F -.25(va)183 194.4 S(lues are meaningless.).25 E F1(-cir)108.01 211.2 Q
1424
(cle_avg r)-.37 E(adius)-.15 E F0 .284(calculate area-weighted a)183
1425
223.2 R -.15(ve)-.2 G .284(rage o).15 F -.15(ve)-.15 G 2.784(ra).15 G
1426
.284(ll grid points within a circle of the speci-)-2.784 F .124
1427
(\214ed radius [km]. e.g. -circle_a)183 235.2 R .123
1428
(vg 200; note that the tracing of \214elds within a cir)-.2 F(-)-.2 E
1429
(cle is quite slo)183 247.2 Q -.65(w.)-.25 G F1(-cir)108.01 264 Q
1430
(cle_max r)-.37 E(adius)-.15 E F0(calculate maximum within a circle of \
1431
a radius [km]. e.g. -circle_max 200)183 276 Q F1(-cir)108.01 292.8 Q
1432
(cle_min r)-.37 E(adius)-.15 E F0(calculate minimum within a circle of \
1433
a radius [km]. e.g. -circle_min 200)183 304.8 Q F1(-clustering)108.01
1434
321.6 Q F0 .846
1435
(special mode to trace strat/trop label; the labels are attrib)183 321.6
1436
R .846(uted according to the)-.2 F 1.182
1437
(program \(tropopause\), which clusters the atmosphere into \214v)183
1438
333.6 R 3.682(ed)-.15 G 1.182(istinct classes)-3.682 F .103(according t\
1439
o the de\214nition of the tropopause: PV \(2 PVU\) and potential temper)
1440
183 345.6 R(-)-.2 E 1.397
1441
(ature \(380 K\). The clustering mode is a re\214ned v)183 357.6 R 1.397
1442
(ersion of the nearest mode,)-.15 F
1443
(where all surrounding eight grid points v)183 369.6 Q
1444
(ote for the \214nal v)-.2 E(alue.)-.25 E F3(ONLINE CALCULA)72 386.4 Q
1445
(TIONS)-1.04 E F0 .073(If the computation \215ag in the tracing \214le \
1446
is set to 1, a meteorological \214eld is calculated based upon)108 398.4
1447
R .106(the already traced \214elds and/or based on the \214elds on the \
1448
primary and secondary P and S \214les. The)108 410.4 R(follo)108 422.4 Q
1449
(wing \214elds are implemented for online calculations:)-.25 E F2 2.5
1450
(-T)108 439.2 S(H)-2.5 E F0(potential temperature \(in K\).)133 451.2 Q
1451
F2 2.5(-R)108 468 S(HO)-2.5 E F0(density \(in kg/m^-3\).)133 480 Q F2
1452
2.5(-R)108 496.8 S(H)-2.5 E F0(relati)133 508.8 Q .3 -.15(ve h)-.25 H
1453
(umidity \(in %\).).15 E F2 2.5(-T)108 525.6 S(HE)-2.5 E F0(equi)133
1454
537.6 Q -.25(va)-.25 G(lent-potential temperature \(in K\).).25 E F2 2.5
1455
(-L)108 554.4 S(HR)-2.5 E F0
1456
(latent heating rate \(K per input time step, typically K/6h\).)133
1457
566.4 Q F2 2.5(-D)108 583.2 S([U)-2.5 E(,V)-.5 E(,T)-1.29 E(,TH]DX)-.74
1458
E F0 .884(horizontal deri)133 595.2 R -.25(va)-.25 G(ti).25 E 1.184 -.15
1459
(ve d)-.25 H([U,V).15 E(,T)-1.29 E .884
1460
(,TH]/dx in west-east direction along pressure surf)-.74 F .884
1461
(aces - zonal)-.1 F 2.836(distance in m. U=zonal wind component \(m/s\)\
1462
, V=meridional wind component \(m/s\),)133 607.2 R(T=temperature \(de)
1463
133 619.2 Q 2.5(gCo)-.15 G 2.5(rK)-2.5 G
1464
(\), TH=potential temperature \(K\).)-2.5 E F2 2.5(-D)108 636 S([U)-2.5
1465
E(,V)-.5 E(,T)-1.29 E(,TH]D)-.74 E(Y)-.4 E F0 3.741(horizontal deri)133
1466
648 R -.25(va)-.25 G(ti).25 E 4.041 -.15(ve d)-.25 H([U,V).15 E(,T)-1.29
1467
E 3.741(,TH]/dy in south-north direction along pressure surf)-.74 F
1468
(aces)-.1 E(-meridional distance in m.)133 660 Q F2 2.5(-D)108 676.8 S
1469
([U)-2.5 E(,V)-.5 E(,T)-1.29 E(,TH]DP)-.74 E F0 -.15(ve)133 688.8 S
1470
(rtical deri).15 E -.25(va)-.25 G(ti).25 E .3 -.15(ve d)-.25 H([U,V).15
1471
E(,T)-1.29 E(,TH]/dp - pressure p in P)-.74 E(a.)-.15 E F2 2.5(-N)108
1472
705.6 S(SQ)-2.5 E F0(squared Brunt-V)133 717.6 Q
1473
(aisala frequence \(in m^-2\).)-1.11 E(2)535 768 Q 0 Cg EP
1474
%%Page: 3 24
1475
%%BeginPageSetup
1476
BP
1477
%%EndPageSetup
1478
/F0 10/Times-Roman@0 SF 413.32(trace\(\) trace\(\))72 48 R/F1 10
1479
/Times-Bold@0 SF 2.5(-R)108 84 S(EL)-2.5 E -.45(VO)-.92 G -.4(RT).45 G
1480
F0(relati)133 96 Q .3 -.15(ve v)-.25 H(orticity \(in s^-1\) - REL)-.05 E
1481
-.4(VO)-1 G 1.2 -.6(RT = D).4 H(VDX - DUD).2 E -1.29(Y.)-.55 G F1 2.5
1482
(-A)108 112.8 S(BSV)-2.5 E(OR)-.45 E(T)-.4 E F0 .083(absolute v)133
1483
124.8 R .083(orticity \(in s^-1\) - ABSV)-.2 F(OR)-.4 E 2.583(T=D)-.6 G
1484
.083(VDX - DUD)-2.983 F 2.582(Y+F)-.55 G 2.582(,Fb)-3.382 G .082
1485
(eing the Coriolis parame-)-2.582 F(ter)133 136.8 Q(.)-.55 E F1 2.5(-D)
1486
108 153.6 S(IV)-2.5 E F0(horizontal di)133 165.6 Q -.15(ve)-.25 G -.18
1487
(rg).15 G(ence of the v).18 E
1488
(elocity \214eld \(in s^-1\) - DIV = DUDX + D)-.15 E(VD)-.4 E -1.29(Y.)
1489
-.55 G F1 2.5(-D)108 182.4 S(EF)-2.5 E F0 .474
1490
(horizontal deformation of the v)133 194.4 R .474
1491
(elocity \214eld \(in s^-1\) - DEF = SQR)-.15 F .475(T\( \( D)-.6 F .475
1492
(VDX + DUD)-.4 F 2.975(Y\))-.55 G(^2)-2.975 E 2.5(+\()133 206.4 S
1493
(DUDX-D)-2.5 E(VD)-.4 E(Y\)^2 \).)-.55 E F1 2.5(-P)108 223.2 S(V)-2.5 E
1494
F0 .485(Ertel potential v)133 223.2 R .485
1495
(orticity \(in PVU\) - PV = g * \( ABSV)-.2 F(OR)-.4 E 2.984(T*D)-.6 G
1496
.484(THDP + DUDP * DTHD)-2.984 F 2.984(Y-)-.55 G -.4(DV)133 235.2 S
1497
(DP * DTHDX \).).4 E F1 2.5(-R)108 252 S(I)-2.5 E F0
1498
(Richardson number - RI = NSQ / \(DUDP^2 + D)133 252 Q(VDP^2 \).)-.4 E
1499
F1 2.5(-T)108 268.8 S(I)-2.5 E F0(tub)133 268.8 Q .379
1500
(ulence indicator according to Ellrod & Knapp - TI = DEF * SQR)-.2 F .38
1501
(T\( DUDP^2 + D)-.6 F(VDP^2)-.4 E 2.5(\)*\(R)133 280.8 S(HO * G\).)-2.5
1502
E F1 2.5(-D)108 297.6 S(IR)-2.5 E F0 2.16(wind direction relati)133
1503
309.6 R 2.46 -.15(ve t)-.25 H 4.659(oz).15 G 2.159(onal \215o)-4.659 F
1504
2.159(w: \(U,V\)=\(1,1\) -> 45 de)-.25 F 2.159
1505
(g; \(U,V\)=\(1,-1\) -> -45 de)-.15 F(g;)-.15 E .411
1506
(\(U,V\)=\(-1,-1\) -> -135 de)133 321.6 R .411
1507
(g; \(U,V\)=\(-1,1\) -> 135 de)-.15 F .411(g. A westerly \215o)-.15 F
1508
2.912(wh)-.25 G .412(as 0 de)-2.912 F .412(g, a southerly)-.15 F<8d6f>
1509
133 333.6 Q 2.5(w9)-.25 G 2.5(0d)-2.5 G -.15(eg)-2.5 G 2.5(,a).15 G
1510
(nd a northerly one -90 de)-2.5 E(g.)-.15 E F1 2.5(-D)108 350.4 S(IST0)
1511
-2.5 E F0(spherical distance \(in km\) from starting position.)133 362.4
1512
Q F1 2.5(-D)108 379.2 S(IST)-2.5 E F0 .33
1513
(length of the trajectory \(in km\): inte)133 391.2 R .33
1514
(grated along great circle sections between the trajectory)-.15 F -.15
1515
(ve)133 403.2 S(rtices.).15 E F1 2.5(-H)108 420 S(EAD)-2.5 E F0 4.261
1516
(heading of the trajectory: \(DX,D)133 432 R 4.261(Y\)=\(1,1\) -> 45 de)
1517
-.55 F 4.261(g; \(DX,D)-.15 F 4.261(Y\)=\(1,-1\) -> -45 de)-.55 F(g;)
1518
-.15 E(\(DX,D)133 444 Q 1.363(Y\)=\(-1,-1\) -> -135 de)-.55 F 1.363
1519
(g; \(DX,D)-.15 F 1.362(Y\)=\(-1,1\) -> 135 de)-.55 F 1.362
1520
(g. A path increment to east has)-.15 F(heading of 0 de)133 456 Q
1521
(g; to the north 90 de)-.15 E(g; to the south -90 de)-.15 E
1522
(g; and to the west -180 de)-.15 E(g.)-.15 E/F2 10.95/Times-Bold@0 SF
1523
(EXAMPLES)72 472.8 Q F1([1] trace TRAJECT)108 484.8 Q(OR)-.18 E -.92(Y.)
1524
-.35 G 2.5(1T).92 G(RAJECT)-2.5 E(OR)-.18 E -.92(Y.)-.35 G 2.5(1-).92 G
1525
(changet)-2.5 E F0 .315(Read the trajectory \214le TRAJECT)133 496.8 R
1526
(OR)-.18 E -1.29(Y.)-.65 G .315
1527
(1, trace all \214elds in the \214le "trace)1.29 F -.25(va)-.25 G .315
1528
(rs" along the tra-).25 F .934(jectories and o)133 508.8 R -.15(ve)-.15
1529
G .934(rwrite the e).15 F .933
1530
(xisting trajectory \214le. In preparation, all times on the P and S)
1531
-.15 F(\214les are changed prior to the tracing.)133 520.8 Q F1
1532
([2] trace INPTRA.1 OUTTRA.1 -f PV 1.)108 537.6 Q F0 -.35(Tr)133 549.6 S
1533
1.22(ace PV \(with scaling f).35 F 1.22
1534
(actor 1.\) along the trajectories in trajectory \214le "INPTRA.1" and)
1535
-.1 F(write a ne)133 561.6 Q 2.5(wt)-.25 G(rajectory \214le "OUTTRA.1".)
1536
-2.5 E F1([3] trace INPTRA.1 OUTTRA.1 -f PV)108 578.4 Q(:-100HP)-.92 E
1537
2.5(A1)-.74 G(.)-2.5 E F0 .122(As in e)133 590.4 R .121(xample [2], b)
1538
-.15 F .121(ut the PV is tak)-.2 F .121(en at a position 100 hP)-.1 F
1539
2.621(ah)-.15 G .121(igher \(lo)-2.621 F .121(wer pressure\) than the)
1540
-.25 F(air parcel')133 602.4 Q 2.5(sp)-.55 G(osition.)-2.5 E F1
1541
([4] trace INPTRA.1 OUTTRA.1 -f DIST0 1.)108 619.2 Q F0(Get the spheric\
1542
al distance \(in km\) of the air parcel from its starting position.)133
1543
631.2 Q F2 -.548(AU)72 648 S(THOR).548 E F0
1544
(Written by Michael Sprenger and Heini W)108 660 Q
1545
(ernli \(January 2011\).)-.8 E(3)535 768 Q 0 Cg EP
1546
%%Page: 1 25
1547
%%BeginPageSetup
1548
BP
1549
%%EndPageSetup
1550
/F0 10/Times-Roman@0 SF 398.86(trainfo\(\) trainfo\(\))72 48 R/F1 10.95
1551
/Times-Bold@0 SF -.219(NA)72 84 S(ME).219 E/F2 10/Times-Bold@0 SF
1552
(trainf)108 96 Q 2.5(o-w)-.25 G(rite meta-inf)-2.5 E(ormation f)-.25 E
1553
(or a trajectory \214le)-.25 E F1(SYNOPSIS)72 112.8 Q F2(trainf)108
1554
124.8 Q(o)-.25 E/F3 10/Times-Italic@0 SF(tr)2.63 E(a\214le)-.15 E F0([)
1555
2.68 E F3(option)2.73 E F0(].)2.74 E F1(DESCRIPTION)72 141.6 Q F0 .768
1556
(Write meta-information for a trajectory \214le)108 153.6 R F3(tr)3.399
1557
E(a\214le)-.15 E F0 .769(to screen. If no option is gi)3.449 F -.15(ve)
1558
-.25 G .769(n, all meta-information is).15 F
1559
(written, otherwise the speci\214c piece of information is passed with)
1560
108 165.6 Q F3(option.)2.73 E F1 -.81(PA)72 182.4 S(RAMETERS).81 E F3
1561
(tr)108.13 194.4 Q(a\214le)-.15 E F0
1562
(name of the input trajectory \214le)133 206.4 Q F1(OPTION)72 223.2 Q
1563
(AL ARGUMENTS)-.219 E F2 2.5(-d)108 235.2 S(im)-2.5 E F0
1564
(dimensions of the trajectory \214le: #tra, #ntimes, #ncolumns.)133
1565
247.2 Q F2 2.5(-n)108 264 S(tra)-2.5 E F0(number of trajectories.)133
1566
276 Q F2 2.5(-n)108 292.8 S(tim)-2.5 E F0(number of times.)133 304.8 Q
1567
F2 2.5(-n)108 321.6 S(col)-2.5 E F0
1568
(number of columns \(including time, longitude, latitude, pressure\).)
1569
133 333.6 Q F2 2.5(-v)108 350.4 S(ars)-2.6 E F0
1570
(list of \214eld names \(columns\) on the trajectory \214le.)133 362.4 Q
1571
F2 2.5(-r)108 379.2 S(efdate)-2.68 E F0
1572
(reference date in format \(YYYYMMDD_HHMM\).)133 391.2 Q F2 2.5(-t)108
1573
408 S(imes)-2.5 E F0(list of times \(relati)133 420 Q .3 -.15(ve t)-.25
1574
H 2.5(ot).15 G(he referencec date\). T)-2.5 E(imes are gi)-.35 E -.15
1575
(ve)-.25 G 2.5(ni).15 G 2.5(nf)-2.5 G(ormat: HH.MM.)-2.5 E F2 2.5(-s)108
1576
436.8 S(tartdate)-2.5 E F0(starting date for the trajectory calculation\
1577
 \(in format YYYYMMD_HHMM\).)133 448.8 Q F2 2.5(-e)108 465.6 S(nddate)
1578
-2.5 E F0
1579
(end date for the trajectory calculation \(in format YYYYMMD_HHMM\).)133
1580
477.6 Q F2 2.5(-t)108 494.4 S(imerange)-2.5 E F0
1581
(time range \(in minutes\) of the trajectories.)133 506.4 Q F2 2.5(-l)
1582
108 523.2 S(ist)-2.5 E F0(list all trajectories.)133 523.2 Q F1
1583
(EXAMPLES)72 540 Q F2([1] trainf)108 552 Q 2.5<6f8c>-.25 G(le dim)-2.5 E
1584
F0(Gi)133 564 Q -.15(ve)-.25 G 3.129(nat).15 G .629(rajectory \214le wi\
1585
th name "\214le", write the three dimensions of the \214le to screen. T\
1586
he dimen-)-3.129 F(sions are: number of trajectories, number of times a\
1587
nd number of columns.)133 576 Q F1 -.548(AU)72 592.8 S(THOR).548 E F0
1588
(Written by Michael Sprenger and Heini W)108 604.8 Q
1589
(ernli \(January 2011\))-.8 E(1)535 768 Q 0 Cg EP
1590
%%Page: 1 26
1591
%%BeginPageSetup
1592
BP
1593
%%EndPageSetup
1594
/F0 10/Times-Roman@0 SF 389.98(wrfmap\(\) wrfmap\(\))72 48 R/F1 10.95
1595
/Times-Bold@0 SF -.219(NA)72 84 S(ME).219 E/F2 10/Times-Bold@0 SF
1596
(wrfmap -)108 96 Q F0
1597
(transform between geographical and WRF coordinates)2.5 E F1(SYNOPSIS)72
1598
112.8 Q F2(wrfmap)108 124.8 Q/F3 10/Times-Italic@0 SF 2.5([-)2.79 G(cr)
1599
-2.5 E(eate \214lename | -ll2xy coor)-.37 E(ds | -xy2ll coor)-.37 E 2.5
1600
(d])-.37 G F1(DESCRIPTION)72 141.6 Q F0 .032
1601
(Create a mapping netCDF \214le)108 153.6 R F3(wrfmap.nc)2.872 E F0 .032
1602
(which contains the forw)2.842 F .032(ard and the backw)-.1 F .033
1603
(ard coordinate transfor)-.1 F(-)-.2 E .529(mation. The forw)108 165.6 R
1604
.528(ard transformation is from WRF grid indices \(x/y\) to geographica\
1605
l longitude/latitude, the)-.1 F(backw)108 177.6 Q 1.522
1606
(ard transformation the other w)-.1 F 1.522(ay round.)-.1 F 1.522
1607
(If the mapping \214le)6.522 F F3(wrfmap.nc)4.362 E F0 -.15(ex)4.332 G
1608
1.522(ists, the tool allo).15 F 1.522(ws to)-.25 F
1609
(transform between the tw)108 189.6 Q 2.5(oc)-.1 G(oordinates.)-2.5 E F1
1610
-.81(PA)72 206.4 S(RAMETERS).81 E F3(-cr)108.01 218.4 Q(eate \214lename)
1611
-.37 E F0 1.133
1612
(create the mapping \214le wrfmap.nc based on the WRF output \214le)183
1613
218.4 R F3(\214lename)5.543 E(.)-.15 E F0 1.133(The WRF)6.133 F .477
1614
(output \214le must contain the \214elds XLA)183 230.4 R 2.977(Ta)-1.11
1615
G .477(nd XLONG which gi)-2.977 F .778 -.15(ve f)-.25 H .478
1616
(or each grid point of).15 F 1.686
1617
(the WRF domain the corresponding longitude and latitude. This forw)183
1618
242.4 R 1.686(ard transfoation)-.1 F 1.277
1619
(\214elds are written as LON and LA)183 254.4 R 3.778(Tt)-1.11 G 3.778
1620
(ow)-3.778 G 1.278(rfmap.nc. F)-3.778 F 1.278(or the re)-.15 F -.15(ve)
1621
-.25 G 1.278(rse transformation, tw).15 F(o)-.1 E .394
1622
(additional \214elds X and Y are calculated and written, which gi)183
1623
266.4 R .694 -.15(ve f)-.25 H .394(or each possible longi-).15 F
1624
(tude/latitude the corresponding grid point \(x/y\) in the WRF domain.)
1625
183 278.4 Q F3(-ll2xy [ lon lat | inp\214le out\214le ])108.01 295.2 Q
1626
F0 .322(write the x/y coordinates for either the single position lon/la\
1627
t to screen, or for all coordi-)183 307.2 R 1.209(nates on the input \
1628
\214le 'inp\214le' to 'out\214le'. The input and output \214le can be t\
1629
rajectory)183 319.2 R(\214les, in which case the column headers are acc\
1630
ordingly adapted to x/y or lon/lat.)183 331.2 Q F3
1631
(-xy2ll [ x y | inp\214le out\214le ])108.01 348 Q F0 .322(write the lo\
1632
n/lat coordinates for either the single position x/y to screen, or for \
1633
all coordi-)183 360 R 1.209(nates on the input \214le 'inp\214le' to 'o\
1634
ut\214le'. The input and output \214le can be trajectory)183 372 R(\214\
1635
les, in which case the column headers are accordingly adapted to x/y or\
1636
 lon/lat.)183 384 Q F1(EXAMPLES)72 400.8 Q F2([1] wrfmap -cr)108 412.8 Q
1637
(eate wrf)-.18 E(out_d01_2011-08-26_12:00:00)-.25 E F0 -.8(Ta)133 424.8
1638
S 2.371 -.1(ke t).8 H 2.172(he netCDF \214le 'wrfout_d01_2011-08-26_12:\
1639
00:00' from WRF and write the forw).1 F 2.172(ard and)-.1 F(backw)133
1640
436.8 Q(ard mapping \214elds to netCDF \214le 'wrfmap.nc')-.1 E F2
1641
([2] wrfmap -ll2xy -70 56)108 453.6 Q F0
1642
(Get the grid indices for \(lon=-70/lat=56\) and write them to screen.)
1643
133 465.6 Q F2([3] wrfmap -xy2ll traj.1 traj.1)108 482.4 Q F0 -.8(Ta)133
1644
494.4 S .206 -.1(ke t).8 H .006(he trajectory \214le ').1 F .006
1645
(traj.1' and con)-.18 F -.15(ve)-.4 G .005
1646
(rt all grid coordinates x/y into lon/lat coordinates. The output).15 F
1647
(will o)133 506.4 Q -.15(ve)-.15 G(rwrite the input \214le.).15 E F1
1648
-.548(AU)72 523.2 S(THOR).548 E F0
1649
(Written by Michael Sprenger and Heini W)108 535.2 Q
1650
(ernli \(January 2013\))-.8 E(1)535 768 Q 0 Cg EP
1651
%%Trailer
1652
end
1653
%%EOF