1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
|
<HTML>
<HEAD>
<TITLE>grdcontour</TITLE>
</HEAD>
<BODY>
<H1>grdcontour</H1>
<HR>
<PRE>
<!-- Manpage converted by man2html 3.0.1 -->
grdcontour - contouring of 2-D gridded data sets
</PRE>
<H2>SYNOPSIS</H2><PRE>
<B>grdcontour</B> <I>grdfile</I> <B>-C</B><I>cont</I><B>_</B><I>int</I> <B>-J</B><I>parameters</I> [
<B>-A</B>[<B>-</B>][<I>anot</I><B>_</B><I>int</I>][<B>f</B><I>font</I><B>_</B><I>size</I>][<B>a</B><I>angle</I>][<I>/r/g/b</I>][<B>o|t</B>]] ] [
<B>-B</B><I>tickinfo</I> ] [ <B>-D</B><I>dumpfile</I> ] [ <B>-E</B><I>azimuth/elevation</I> ] [
<B>-G</B><I>gap/width</I> ] [ <B>-K</B> ] [ <B>-L</B><I>low/high</I> ] [ <B>-M</B>[<I>flag</I>] ] [
<B>-N</B>[[-]<I>unit</I>] ] [ <B>-O</B> ] [ <B>-P</B> ] [ <B>-Q</B><I>cut</I> ] [
<B>-R</B><I>west/east/south/north</I>[<B>r</B>] ] [ <B>-S</B><I>smoothfactor</I> ] [
<B>-T</B>[<B>+|-</B>][<I>gap/length</I>][<B>:</B><I>LH</I>] ] [ <B>-U</B>[<I>/dx/dy/</I>][<I>label</I>] ] [ <B>-V</B> ]
[ <B>-W</B>[<B>+</B>][<I>type</I>]<I>pen</I> ] [ <B>-X</B><I>x-shift</I> ] [ <B>-Y</B><I>y-shift</I> ] [ <B>-Z</B>[<I>fac</I>
<I>tor</I>[/<I>shift</I>]][<B>p</B>] ] [ <B>-c</B><I>copies</I> ] [ <B>-bo</B>[<B>s</B>] ]
</PRE>
<H2>DESCRIPTION</H2><PRE>
<B>grdcontour</B> reads a 2-D gridded file and produces a contour
map by tracing each contour through the grid. As an
option, the x/y/z positions of the contour lines may be
dumped to a single multisegment file or many separate
files. <I>PostScript</I> code is generated and sent to standard
output. Various options that affect the plotting are
available.
No space between the option flag and the associ
ated arguments. Use upper case for the option flags and
lower case for modifiers.
<I>grdfile</I>
2-D gridded data set to be contoured
<B>-C</B> The contours to be drawn may be specified in one of
three possible ways:
If <I>cont</I><B>_</B><I>int</I> has the suffix ".cpt" and can be opened
as a file, it is assumed to be a color palette
table. The color boundaries are then used as con
tour levels. If the cpt-file has anotation flags
in the last column then those contours will be ano
tated. By default all contours are labeled; use <B>-A-</B>
to disable all anotations.
If <I>cont</I><B>_</B><I>int</I> is a file but not a cpt-file, it is
expected to contain contour levels in column 1 and
a C(ontour) OR A(nnotate) in col 2. The levels
marked C (or c) are contoured, the levels marked A
(or a) are contoured and annotated. Optionally, a
third column may be present and contain the fixed
anotation angle for this contour level.
If no file is found, then <I>cont</I><B>_</B><I>int</I> is interpreted
as a constant contour interval. If <B>-A</B> is set and <B>-C</B>
is not, then the contour interval is set equal to
the specified anotation interval.
If a file is given and <B>-T</B> is set, then only con
tours marked with upper case C or A will have tick
marks.
UNIT is cm, inch, or m, depending on the MEA
SURE_UNIT setting in .<A HREF="gmtdefaults.html">gmtdefaults</A>, but this can be
overridden on the command line by appending the c,
i, or m to the scale/width value.
<B>CYLINDRICAL</B> <B>PROJECTIONS:</B>
<B>-Jc</B><I>lon0/lat0/scale</I> (Cassini)
<B>-Jj</B><I>lon0/scale</I> (Miller)
<B>-Jm</B><I>scale</I> (Mercator - Greenwich and Equator as ori
gin)
<B>-Jm</B><I>lon0/lat0/scale</I> (Mercator - Give meridian and
standard parallel)
<B>-Joa</B><I>lon0/lat0/azimuth/scale</I> (Oblique Mercator -
point and azimuth)
<B>-Job</B><I>lon0/lat0/lon1/lat1/scale</I> (Oblique Mercator -
two points)
<B>-Joc</B><I>lon0/lat0/lonp/latp/scale</I> (Oblique Mercator -
point and pole)
<B>-Jq</B><I>lon0/scale</I> (Equidistant Cylindrical Projection
(Plate Carree))
<B>-Jt</B><I>lon0/scale</I> (TM - Transverse Mercator, with Equa
tor as y = 0)
<B>-Jt</B><I>lon0/lat0/scale</I> (TM - Transverse Mercator, set
origin)
<B>-Ju</B><I>zone/scale</I> (UTM - Universal Transverse Mercator)
<B>-Jy</B><I>lon0/lats/scale</I> (Basic Cylindrical Projection)
<B>AZIMUTHAL</B> <B>PROJECTIONS:</B>
<B>-Ja</B><I>lon0/lat0/scale</I> (Lambert).
<B>-Je</B><I>lon0/lat0/scale</I> (Equidistant).
<B>-Jf</B><I>lon0/lat0/horizon/scale</I> (Gnomonic).
<B>-Jg</B><I>lon0/lat0/scale</I> (Orthographic).
<B>-Js</B><I>lon0/lat0/</I>[<I>slat/</I>]<I>scale</I> (General Stereographic)
<B>CONIC</B> <B>PROJECTIONS:</B>
<B>-Jb</B><I>lon0/lat0/lat1/lat2/scale</I> (Albers)
<B>-Jd</B><I>lon0/lat0/lat1/lat2/scale</I> (Equidistant)
<B>-Jl</B><I>lon0/lat0/lat1/lat2/scale</I> (Lambert)
<B>MISCELLANEOUS</B> <B>PROJECTIONS:</B>
<B>-Jh</B><I>lon0/scale</I> (Hammer)
<B>-Ji</B><I>lon0/scale</I> (Sinusoidal)
<B>-Jk</B>[<B>f|s</B>]<I>lon0/scale</I> (Eckert IV (f) and VI (s))
<B>-Jn</B><I>lon0/scale</I> (Robinson)
<B>-Jr</B><I>lon0/scale</I> (Winkel Tripel)
<B>-Jv</B><I>lon0/scale</I> (Van der Grinten)
<B>-Jw</B><I>lon0/scale</I> (Mollweide)
optional <B>a</B> for azimuths and offset theta [0])
<B>-Jx</B><I>x-scale</I>[<B>l|p</B><I>pow</I>][<I>/y-scale</I>[<B>l|p</B><I>pow</I>]] (Linear, log,
and power scaling)
More details can be found in the <B><A HREF="psbasemap.html">psbasemap</A></B> man
pages.
</PRE>
<H2>OPTIONS</H2><PRE>
No space between the option flag and the associated argu
ments.
<B>-A</B> <I>anot</I><B>_</B><I>int</I> is annotation interval in data units.
Ignored if contour levels are given in a file.
[Default is no annotations]. Several options can be
set to modify the form of the annotation. Give - to
disable all anotations. Append f<I>font</I><B>_</B><I>size</I> to
change font size [9], append <I>/r/g/b</I> to change color
of text fill box [PAGE_COLOR], append a<I>angle</I> to fix
annotation angle [Default follows contour], append
<B>o</B> to draw the outline of the surrounding text box
[Default is no outline], or append <B>t</B> to use a
transparent label box [Default is opaque].
<B>-B</B> Sets map boundary tickmark intervals. See <B><A HREF="psbasemap.html">psbasemap</A></B>
for details.
<B>-D</B> Dump the (x,y,z) coordinates of each contour to
separate files, one for each contour segment. The
files will be named <I>dumpfile</I><B>_</B><I>cont</I><B>_</B><I>segment[</I><B>_</B><I>i]</I>.xyz
(or <I>.b</I> is <B>-b</B> is selected), where <I>cont</I> is the con
tour value and <I>segment</I> is a running segment number
for each contour interval (for closed contours we
append _i.) If the prefix is given as '-' the file
names are instead <I>C#</I><B>_</B><I>i</I> (interior) or <I>C#</I><B>_</B><I>e</I> (exter
nal) plus extension, and # is just a running num
ber. This allows us to make short file names that
will work with GNU utilities under DOS. However,
when <B>-M</B> is used in conjunction with <B>-D</B> a single
multisegment file is created instead.
<B>-E</B> Sets the viewpoint's azimuth and elevation for per
spective view [180/90]'
<B>-G</B> <I>gap</I> is distance between each annotation along the
same contour. <I>width</I> is number of points over which
to estimate the best fitting slope for contour
labels [Default is 10<B>c</B>/10 or 4<B>i</B>/10].
<B>-K</B> More <I>PostScript</I> code will be appended later
[Default terminates the plot system].
<B>-L</B> Limit range: Do not draw contours for data values
below <I>low</I> or above <I>high</I>.
ment file is created, and each contour section is
preceeded by a header record whose first column is
<I>flag</I> followed by the contour level.
<B>-N</B> Appends <I>unit</I> to all contour labels (assuming <B>-A</B> is
set). If no name is specified, the unit given in
the <I>grdfile</I> will be used. If <I>unit</I> starts with a
leading - then there will be no space between con
tour value and the unit. [Default is no unit].
<B>-O</B> Selects Overlay plot mode [Default initializes a
new plot system].
<B>-P</B> Selects Portrait plotting mode [<B><A HREF="GMT.html">GMT</A></B> Default is
Landscape, see gmtdefaults to change this].
<B>-Q</B> Do not draw contours with less than <I>cut</I> number of
points [Draw all contours].
<B>-R</B> <I>west,</I> <I>east,</I> <I>south,</I> and <I>north</I> specify the Region of
interest. To specify boundaries in degrees and min
utes [and seconds], use the dd:mm[:ss] format.
Append <B>r</B> if lower left and upper right map coordi
nates are given instead of wesn. [Default is
region defined in the grd file].
<B>-S</B> Used to resample the contour lines at roughly every
(gridbox_size/<I>smoothfactor</I>) interval.
<B>-T</B> Will draw tickmarks pointing in the downward direc
tion every <I>gap</I> along the innermost closed contours.
Append <I>gap</I> and tickmark length or use defaults
[0.5<B>c</B>/0.1c or 0.2<B>i</B>/0.04i]. User may choose to tick
only local highs or local lows by specifying <B>-T+</B> or
<B>-T-</B>, respectively. Appending :LH will plot the
characters L and H at the center of closed inner
most contours (local lows and highs). L and H can
be any single character (e.g., LH, -+, etc.) If a
file is given by <B>-C</B> and <B>-T</B> is set, then only con
tours marked with upper case C or A will have tick
marks [and anotation].
<B>-U</B> Draw Unix System time stamp on plot. User may spec
ify where the lower left corner of the stamp should
fall on the page relative to lower left corner of
plot. Optionally, append a label, or <B>c</B> (which will
plot the command string.). The GMT parameters
UNIX_TIME and UNIX_TIME_POS can affect the appear
ance; see the <B><A HREF="gmtdefaults.html">gmtdefaults</A></B> man page for details.
<B>-V</B> Selects verbose mode, which will send progress
reports to stderr [Default runs "silently"].
or 'c' for regular contours [Default]. <I>pen</I> sets
the attributes for the particular line. Default
values for annotated contours: width = 3, color =
black, texture = solid. Regular contours have
default width = 1. If the <B>+</B> flag is specified then
the color of the contour lines are taken from the
cpt file (see <B>-C</B>).
<B>-X</B> <B>-Y</B> Shift origin of plot by (<I>x-shift,y-shift</I>). Prepend
<B>a</B> for absolute coordinates; the default (<B>r</B>) will
reset plot origin.
<B>-Z</B> Use to subtract <I>shift</I> from the data and multiply
the results by <I>factor</I> before contouring starts
[1/0]. (Numbers in <B>0</B> <B>-L</B> refer to values after this
scaling has occurred.) Append <B>p</B> to indicate that
this grid file contains z-values that are periodic
in 360 degrees (e.g., phase data, angular distribu
tions) and that special precautions must be taken
when determining 0-contours.
<B>-c</B> Specifies the number of plot copies. [Default is 1]
<B>-bo</B> Selects binary output. Append <B>s</B> for single preci
sion [Default is double].
</PRE>
<H2>EXAMPLES</H2><PRE>
To contour the file hawaii_grav.grd every 25 mGal on a
Mercator map at 0.5 inch/degree, annotate every 50 mGal
(using fontsize = 10), using 1 degree tickmarks, and draw
30 minute gridlines, try
grdcontour hawaii_grav.grd <B>-Jm</B>0.5<B>i</B> <B>-C</B>25 <B>-A</B>50<B>f</B>10 <B>-B</B>1<B>g</B>30<B>m</B> >
hawaii_grav.ps
To contour the file image.grd using the levels in the file
cont.d on a linear projection at 0.1 cm/x-unit and 50
cm/y-unit, using 20 (x) and 0.1 (y) tickmarks, smooth the
contours a bit, use "RMS Misfit" as plot-title, use a
thick red pen for annotated contours, and a thin, dashed,
blue pen for the rest, try
grdcontour image.grd <B>-Jx</B>0.1<B>c</B>/50.0<B>c</B> <B>-C</B>cont.d <B>-S</B>4
<B>-B</B>20/0.1:."RMS Misfit": <B>-Wa</B>5/255/0/0 <B>-Wc</B>1/0/0/255<B>ta</B> | lp
</PRE>
<H2>BUGS</H2><PRE>
The text-outline box has not been implemented for perspec
tive view plots. Anotations will plot correctly, but con
tours will go right through the text.
The labeling of local highs and lows may plot outside the
innermost contour since only the mean value of the contour
coordinates is used to position the label.
ANOT_OFFSET, respectively. You must thus change this value
in .gmtdefaults to affect the box size.
There are too many options.
</PRE>
<H2>SEE ALSO</H2><PRE>
<I>gmt</I>(l), <I><A HREF="gmtdefaults.html">gmtdefaults</A></I>(l), <I><A HREF="psbasemap.html">psbasemap</A></I>(l), <I><A HREF="grdimage.html">grdimage</A></I>(l), <I>grd</I>
<I>view</I>(l), <I><A HREF="pscontour.html">pscontour</A></I>(l)
</PRE>
<HR>
<ADDRESS>
Man(1) output converted with
<a href="http://www.oac.uci.edu/indiv/ehood/man2html.html">man2html</a>
</ADDRESS>
</BODY>
</HTML>
<body bgcolor="#ffffff">
|
