-
Notifications
You must be signed in to change notification settings - Fork 7
/
fgmax_5_2_0.html
282 lines (270 loc) · 13.3 KB
/
fgmax_5_2_0.html
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
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN"
"http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<html xmlns="http://www.w3.org/1999/xhtml">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
<title>Fixed grid monitoring — Clawpack 5.1.0 documentation</title>
<link rel="stylesheet" href="_static/flasky.css" type="text/css" />
<link rel="stylesheet" href="_static/pygments.css" type="text/css" />
<link rel="stylesheet" href="_static/base.css" type="text/css" />
<link rel="stylesheet" href="_static/layout.css" type="text/css" />
<script type="text/javascript">
var DOCUMENTATION_OPTIONS = {
URL_ROOT: '',
VERSION: '5.1.0',
COLLAPSE_INDEX: false,
FILE_SUFFIX: '.html',
HAS_SOURCE: true
};
</script>
<script type="text/javascript" src="_static/jquery.js"></script>
<script type="text/javascript" src="_static/underscore.js"></script>
<script type="text/javascript" src="_static/doctools.js"></script>
<script type="text/javascript" src="http://cdn.mathjax.org/mathjax/latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML"></script>
<link rel="shortcut icon" href="_static/clawicon.ico"/>
<link rel="author" title="About these documents" href="about.html" />
<link rel="top" title="Clawpack 5.1.0 documentation" href="index.html" />
</head>
<body>
<div id="main-wrapper" class="sphinx">
<div id="header-wrapper">
<section id="header">
<!-- <h1><a href="http://clawpack.org/">Clawpack</a></h1> -->
<h1><a href="http://clawpack.org/">Clawpack-5</a></h1>
<nav>
<ul>
<li>
<a href="toc_condensed.html">Docs</a>
</li>
<li>
<a href="installing.html">Install</a>
</li>
<li>
<a class="" href="galleries.html">Gallery</a>
</li>
<li>
<a class="active" href="https://groups.google.com/forum/#!forum/claw-users">Support</a>
</li>
<li>
<a class="active" href="http://github.com/organizations/clawpack">Source</a>
</li>
<li>
<a class="" href="developers.html">Develop</a>
</li>
</ul>
</nav>
</section>
<div class="decoration"></div>
</div>
<div class="related">
<h3>Navigation</h3>
<ul>
<li class="right" style="margin-right: 10px">
<a href="genindex.html" title="General Index"
accesskey="I">index</a></li>
<li class="right" >
<a href="py-modindex.html" title="Python Module Index"
>modules</a> |</li>
<li><a href="contents.html">Clawpack 5.1.0 documentation</a> »</li>
</ul>
</div>
<div class="document">
<div class="documentwrapper">
<div class="bodywrapper">
<div class="body">
<div class="section" id="fixed-grid-monitoring">
<span id="fgmax-5-2-0"></span><h1>Fixed grid monitoring<a class="headerlink" href="#fixed-grid-monitoring" title="Permalink to this headline">¶</a></h1>
<div class="admonition warning">
<p class="first admonition-title">Warning</p>
<p class="last">This feature has been modified and this documentation describes
the version introduced in 5.2.0.
The documentation is also still incomplete.</p>
</div>
<p>GeoClaw has the capability to monitor certain quantities on a specified
“fixed grid” by interpolating from the AMR grids active at each time step,
or at specified time increments.
This is useful in particular to record the maximum flow depth observed at
each point over the course of a computation, or the maximum flow velocity,
momentum, or momentum flux. These quantities are often of interest in
hazard modeling.</p>
<p>It is also possible to record the <em>arrival time</em> of a flow or wave at each
point on the grid.</p>
<p>The “grids” do not have to be rectangular grids aligned with the
coordinate directions, but can consist of an arbitrary list of points
that could also be points along a one-dimensional transect or points
following a coastline, for example.</p>
<p>Each grid is specified by an input file in a specified form described below.
The list of file names for desired grids is specified in the <cite>setrun</cite>
function, see <a class="reference internal" href="setrun_geoclaw.html#setrun-fgmax"><em>Fixed grid maximum monitoring / arrival times</em></a>.</p>
<p>This is an improved version of the algorithms used in earlier versions of
GeoClaw, and now
correctly interpolates when a grid point lies near the junction of two
grid patches, which was not always handled properly before.
The earlier version can still be used for outputing results at intermediate
times on a fixed grid (see <a class="reference internal" href="fgout.html#fgout"><em>Fixed grid output</em></a>), but is not recommended for the purpose
of monitoring maxima or arrival times.</p>
<div class="section" id="input-file-specification">
<span id="fgmax-input"></span><h2>Input file specification<a class="headerlink" href="#input-file-specification" title="Permalink to this headline">¶</a></h2>
<p>(changed in Clawpack 5.2.0.)</p>
<p>The input file describing a grid of points has the following form:</p>
<div class="highlight-python"><div class="highlight"><pre><span class="n">tstart_max</span>
<span class="n">tend_max</span>
<span class="n">dt_check</span>
<span class="n">min_level_check</span>
<span class="n">arrival_tol</span>
<span class="n">point_style</span>
</pre></div>
</div>
<p>followed by additional lines that depend on the value of <cite>point_style</cite>.</p>
<p>If <cite>point_style == 0</cite>, an arbitrary collection of <cite>(x,y)</cite> points is allowed
and all must be listed, preceeded by the number of points:</p>
<div class="highlight-python"><pre>npts # number of points
x1 y1 # first point
x2 y2 # second point
... # etc.</pre>
</div>
<p>These points need not lie on a regular grid and can be specified in any order.</p>
<p>If <cite>point_style == 1</cite>, a 1-dimensional transect of points is specified by
three lines of the form:</p>
<div class="highlight-python"><div class="highlight"><pre><span class="n">npts</span> <span class="c"># number of points to generate</span>
<span class="n">x1</span><span class="p">,</span> <span class="n">y1</span> <span class="c"># first point</span>
<span class="n">x2</span><span class="p">,</span> <span class="n">y2</span> <span class="c"># last point</span>
</pre></div>
</div>
<p>If <cite>point_style == 2</cite>, a 2-dimensional cartesian of points is specified by
three lines of the form:</p>
<div class="highlight-python"><div class="highlight"><pre><span class="n">nx</span><span class="p">,</span> <span class="n">ny</span> <span class="c"># number of points in x and y (nx by ny grid)</span>
<span class="n">x1</span><span class="p">,</span> <span class="n">y1</span> <span class="c"># lower left corner of cartesian grid</span>
<span class="n">x2</span><span class="p">,</span> <span class="n">y2</span> <span class="c"># upper right corner of cartesian grid</span>
</pre></div>
</div>
<p>The output files will list values for the points in the same order as in the
input file. See <cite>fgmax_processing</cite> for some hints on processing and
plotting the results.</p>
<p>The other paramters in the input file are:</p>
<ul>
<li><p class="first"><cite>tstart_max</cite> : float</p>
<p>starting time to monitor maximum</p>
</li>
<li><p class="first"><cite>tend_max</cite> : float</p>
<p>ending time to monitor maximum</p>
</li>
<li><p class="first"><cite>dt_check</cite> : float</p>
<p>time increment for monitoring maximum and arrivals.
Interpolate to fixed grid and
update values only if the time since the last updating exceeds this time
increment. Set to 0 to monitor every time step.</p>
</li>
<li><p class="first"><cite>min_level_check</cite> : integer</p>
<p>Minimum AMR level to check for updating the maximum value observed and
the arrival time.
Care must be taken in selecting this value since the maximum observed
when interpolating to a point from a coarse AMR level may be much larger
than the value that would be seen on a fine grid that better resolves the
topography at this point. Often AMR “regions” are used to specify that a
fine grid at some level <cite>L</cite> should always be used in the region of
interest over the time period from <cite>start_max</cite> to <cite>tend_max</cite>, and then
it is natural to set <cite>min_level_check</cite> to <cite>L</cite>.</p>
</li>
<li><p class="first"><cite>arrival_tol</cite> : float</p>
<p>The time reported as the “arrival time” is the first time the value
of the surface elevation is greater than <cite>sea_level</cite> + <cite>arrival_tol</cite>.</p>
</li>
</ul>
</div>
<div class="section" id="values-to-monitor">
<span id="fgmax-values"></span><h2>Values to monitor<a class="headerlink" href="#values-to-monitor" title="Permalink to this headline">¶</a></h2>
<p>The values to be monitored are specified by the subroutine <cite>fgmax_values</cite>.
The default subroutine found in the library
<cite>$CLAW/geoclaw/src/2d/shallow/fgmax_values.f90</cite>
is now set up to monitor the
depth <cite>h</cite> (rather than the value <cite>eta_tilde</cite> used in Version 5.1)
and optionally will also monitor the speed <span class="math">\(s = \sqrt{u^2 + v^2}\)</span>
and three other quantities (the momentum <span class="math">\(hs\)</span>,
the momentum flux <span class="math">\(hs^2\)</span>, and <span class="math">\(-h\)</span>, which is useful to monitor
the minimum depth at each point).</p>
<p>The values monitored by the default routine described above is determined
by the value of the <cite>fgmax_module</cite> variable <cite>FG_NUM_VAL</cite>, which can be set
to 1, 2, or 5. This value is now read in from the data file <cite>fgmax.data</cite>
and can be set by specifying the value of
<cite>rundata.fgmax_data.num_fgmax_val</cite> in <cite>setrun.py</cite>.</p>
</div>
<div class="section" id="choice-of-interpolation-procedure">
<h2>Choice of interpolation procedure<a class="headerlink" href="#choice-of-interpolation-procedure" title="Permalink to this headline">¶</a></h2>
<p>The library routine <cite>geoclaw/src/2d/shallow/fgmax_interpolate.f90</cite> has
been improved in 5.2.0 to fix some bugs. This routine does bilinear
interpolation the finite volume grid centers to the fixed grid in
order to update the maximum of values such as depth or velocity.</p>
<p>An alternative version of this routine has been added in 5.2.0
that does piecewise constant interpolation instead. This simply uses the
value in the finite volume grid cell that contains the fixed grid
point (0 order extrapolation) and avoids problems sometimes seen when
doing linear interpolation near the margins of the flow.</p>
<p>This routine is in <cite>fgmax_interpolate0.f90</cite> and is now recommended.
To use this routine, modify the <cite>Makefile</cite> in an application directory to
replace the line</p>
<div class="highlight-python"><pre>$(GEOLIB)/fgmax_interpolate.f90 \</pre>
</div>
<p>by</p>
<div class="highlight-python"><pre>$(GEOLIB)/fgmax_interpolate0.f90 \</pre>
</div>
</div>
<div class="section" id="processing-and-plotting-fgmax-output">
<span id="fgmax-processing"></span><h2>Processing and plotting fgmax output<a class="headerlink" href="#processing-and-plotting-fgmax-output" title="Permalink to this headline">¶</a></h2>
<p>For an example see <cite>apps/tsunami/chile2010_fgmax</cite> in the <a class="reference internal" href="apps.html#apps"><em>Clawpack Applications repository</em></a>.</p>
<p><strong>Describe further.</strong></p>
</div>
</div>
</div>
</div>
</div>
<div class="sphinxsidebar">
<div class="sphinxsidebarwrapper">
<p><a href="http://clawpack.org/">
<img class="logo" src= "_static/clawlogo.jpg" alt="Logo"/>
</a>
<h2>Version 5.1.0</h2>
</p>
<h3><a href="contents.html">Table Of Contents</a></h3>
<ul>
<li><a class="reference internal" href="#">Fixed grid monitoring</a><ul>
<li><a class="reference internal" href="#input-file-specification">Input file specification</a></li>
<li><a class="reference internal" href="#values-to-monitor">Values to monitor</a></li>
<li><a class="reference internal" href="#choice-of-interpolation-procedure">Choice of interpolation procedure</a></li>
<li><a class="reference internal" href="#processing-and-plotting-fgmax-output">Processing and plotting fgmax output</a></li>
</ul>
</li>
</ul>
<h3>Related Topics</h3>
<ul>
<li><a href="contents.html">Documentation overview</a><ul>
</ul></li>
</ul>
<h3>This Page</h3>
<ul class="this-page-menu">
<li><a href="_sources/fgmax_5_2_0.txt"
rel="nofollow">Show Source</a></li>
</ul>
<div id="searchbox" style="display: none">
<h3>Quick search</h3>
<form class="search" action="search.html" method="get">
<input type="text" name="q" />
<input type="submit" value="Go" />
<input type="hidden" name="check_keywords" value="yes" />
<input type="hidden" name="area" value="default" />
</form>
<p class="searchtip" style="font-size: 90%">
Enter search terms or a module, class or function name.
</p>
</div>
<script type="text/javascript">$('#searchbox').show(0);</script>
</div>
</div>
<div class="clearer"></div>
</div>
<div class="footer">
© Copyright 2013, The Clawpack Development Team.
Created using <a href="http://sphinx.pocoo.org/">Sphinx</a>.
</div>
</body>
</html>