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
|
<HTML>
<HEAD>
<TITLE>grdvolume</TITLE>
</HEAD>
<BODY>
<H1>grdvolume</H1>
<HR>
<PRE>
<!-- Manpage converted by man2html 3.0.1 -->
grdvolume - Calculating volume under a surface within a
contour
</PRE>
<H2>SYNOPSIS</H2><PRE>
<B>grdvolume</B> <I>grdfile</I> [ <B>-C</B><I>cval</I> or <B>-C</B><I>low/high/delta</I> ] [ <B>-L</B><I>base</I>
] [ <B>-R</B><I>west/east/south/north</I>[<B>r</B>] ] [ <B>-S</B>[<B>k</B>] ] [ <B>-T</B> ] [ <B>-V</B>[<B>l</B>]
] [ <B>-Z</B><I>fact</I>[/<I>delta</I>] ]
</PRE>
<H2>DESCRIPTION</H2><PRE>
<B>grdvolume</B> reads a 2-D binary grd file and calculates the
volume contained between the surface and the plane speci
fied by the given contour (or zero if not given). Alter
natively, specify a range of contours to be tried and <B>grd</B>
<B>volume</B> will determine the volume and area inside the con
tour for all contour values. The contour that produced the
maximum mean height (volume/area) is reported as well.
This feature may be used with <B><A HREF="grdfilter.html">grdfilter</A></B> in designing an
Optimal Robust Separator [<I>Wessel</I>, 1998].
<I>grdfile</I>
The name of the input 2-D binary grd file.
</PRE>
<H2>OPTIONS</H2><PRE>
No space between the option flag and the associated argu
ments.
<B>-C</B> find area and volume inside the <I>cval</I> contour.
Alternatively, search using all contours from <I>low</I>
to <I>high</I> in steps of <I>delta</I>. [Default returns entire
area and volume of grid].
<B>-L</B> Also add in the volume from the level of the con
tour down to <I>base</I> [Default base is contour].
<B>-S</B> Convert degrees to meters, append <B>k</B> for km [Default
is Cartesian].
<B>-T</B> Use curvature minimum rather than maximum height to
find best contour value (when contour search is
selected with <B>-C</B>).
<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.
<B>-V</B> Selects verbose mode, which will send progress
reports to stderr [Default runs "silently"].
Append <B>l</B> to see all the results for each contour
level tested (when contour search has been
selected).
<I>fact</I>. [Default is no scaling]. (Numbers in <B>-C,</B> <B>-L</B>
refer to values after this scaling has occurred).
</PRE>
<H2>EXAMPLES</H2><PRE>
To determine the volume in km^3 under the <A HREF="surface.html">surface</A>
hawaii_topo.grd (height in km), try
grdvolume hawaii_topo.grd <B>-Sk</B>
To find the volume between the surface peaks.grd and the
contour z = 250, use
grdvolume peaks.grd <B>-Sk</B> <B>-C</B>250
To search for the contour, between 100 and 300 in steps of
10, that maximizes the ratio of volume to surface area for
the file peaks.grd, use
grdvolume peaks.grd <B>-Sk</B> <B>-C</B>100/300/10 > results.d
To see the areas and volumes for all the countours in the
previous example, try
grdvolume peaks.grd <B>-Sk</B> <B>-Vl</B> <B>-C</B>100/300/10 > results.d
</PRE>
<H2>SEE ALSO</H2><PRE>
<I>gmt</I>(l), <I><A HREF="grdfilter.html">grdfilter</A></I>(l)
</PRE>
<H2>REFERENCES</H2><PRE>
Wessel, P., 1998, An empirical method for optimal robust
regional-residual separation of geophysical data, <I>Math.</I>
<I>Geol.,</I> <I>30(4)</I>, 391-408.
</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">
|
