source: DVD_images/extra_files/BatemansBay/analysis.html @ 7329

Last change on this file since 7329 was 7329, checked in by kristy, 15 years ago

updated htmls to reflect editing
File size: 7.3 KB
Line 
1<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01//EN" "http://www.w3.org/TR/html4/strict.dtd">
2<html>
3<head>
4  <title>Geoscience Australia</title>
5  <link rel="stylesheet" href="browser_files/style.css">
6  <link REL="SHORTCUT ICON" HREF="browser_files/favicon.ico">
7</head>
8<body>
9  <table cellspacing="10" border="0">
10    <tr>
11      <td>
12        <img src="browser_files/logo.jpg" alt="Australian Government, Geoscience Australia" width="327" height="80" border="0"/>
13      </td>
14      <td width="50">
15        &nbsp;
16      </td>
17      <td>
18        <H1>Tsunami Inundation Models for the BATEMANS BAY region</H1>
19      </td>
20    </tr>
21  </table>
22  <hr>
23
24  <a name="analysis scripts"><h2><b>Analysis Scripts</b></h2></a>
25  There are number of analysis scripts on this disk that you can use to extract information from the SWW files
26  produced by a simulation run.  These scripts are not meant to be the definitive answer on how to get information
27  from SWW files, but are offered as an example of one way to do it.
28  <p>
29  This page describes the scripts and shows how to use them.  Note that due to the ad-hoc nature of the scripts,
30  you <b>must</b> modify these scripts before they can be used to extract information from any SWW files you
31  generate.
32
33  <a name="export_results_max.py scripts"><h3><b>export_results_max.py</b></h3></a>
34  The <font color="red">export_results_max.py</font> script can be used to create an ASC file that contains the maximum
35  of an ANUGA variable or expression for each point in a user-defined ArcGIS grid on the simulation region.
36  <p>
37  You must change certain values within the <font color="red">export_results_max.py</font> file to get what you want:
38  <ul>
39    <li> List the output folder names here, with one for each simulation you have run:
40<pre><font color="brown">    time_dir1 = '20090505_150430_run_final_0.8_58292_None_kvanputt'
41    time_dir2 = '20090505_150517_run_final_0_58292_None_kvanputt'</font></pre>
42         For information on  <a href="modifications.html#output">output_folder_name</a>.
43  <p>Then modify the line:
44<pre><font color="brown">    time_dirs = [time_dir1, time_dir2]</font></pre>
45         to include each of the variables you defined above.
46         <p>
47    <li> Your SWW output will have a name like <font color="red">batemans_bay.sww</font>. For a large model you
48         may have more than one output file with extra filenames like <font color="red">batemans_bay_37860_0.sww</font>,
49         where the <b>37860</b> tells us that the timestep at the start of this file is 37860 sec. The
50         <font color="red">export_results_max.py</font> script needs to examine all SWW output files, so we specify
51         all the files produced by the model run to create a maximum over the entire time.
52         We do this by specifying the start times of all the extra SWW files:
53         <p>
54<pre><font color="brown">    times = [37860] <font color="black"> or </font>
55    times = [20000, 40000, 60000]</font></pre>
56         <p>
57         <i>Note we do not have to specify the first SWW file, only the extra file start times. If there are no extra SWW files,
58         just do:</i>
59         <p>
60<pre><font color="brown">    times = []</font></pre>
61         <p>
62    <li> Modify the <b>cellsize</b> value to set the size of the ArcGIS grid you require.  For example:
63<pre><font color="brown">    cellsize = 20              # ArcGIS grid cell size in metres</font></pre>
64         <p>
65         <i>Note: That this value should not go beyond the most refined part of the mesh
66         i.e. if mesh resolution is 500m<sup>2</sup>, cellsize = (500 </i>x<i> 2)<sup>1/2</sup> = 32m</i>
67         <p>
68    <li> You can get the maximum of a variable or expression over the entire model time, or for a single time in the simulation.
69         Set <b>timestep</b> to the required time you want the maximum for, or set it to <b>None</b> to indicate that you want
70         the maximum values over all timesteps:
71<pre><font color="brown">    timestep = None  # over all timesteps
72    #timestep = 0</font></pre>
73         <p>
74         <i>Note that if you want the elevation use <b>timestep = 0</b>. Elevation does not change over time and you only need the
75         first SWW file.</i>
76         <p>
77    <li> If you want to clip the ArcGIS grid to one or more small regions then put the names of the regions of interest
78         into the <b>area</b> list:
79<pre><font color="brown">    area = ['Batemans_Bay', 'NW', 'South']</font></pre>
80         <p>
81         <i>Note that the name strings you put into the <b>area</b> list must match the names used in the
82         <font color="red">project.py</font> file.  That is, if you use a name such as 'NW', then the
83         <font color="red">project.py</font> file must have lines like:</i>
84<pre><font color="brown">    xminNW = 548000
85    xmaxNW = 561000
86    yminNW = 5250000
87    ymaxNW = 5258000</font></pre>
88         which set the maximum and minimum eastings and northings that define a rectangular region.
89         <p>
90         If you don't want to clip to a region of interest, then do this:
91<pre><font color="brown">    area = ['All']             # no region of interest</font></pre>
92         <p>
93    <li> Finally you must decide which variable or expression values you want sampled on your ArcGIS grid.
94         Define a list <b>var</b> that contains strings defining the required variable or expression:
95<pre><font color="brown">    var = ['depth', 'speed']</font></pre>
96         <i>Note that the strings you supply must be defined in the <b>var_equations</b> dictionary
97         prior to your definition of <b>var</b>:</i>
98<pre><font color="brown">    var_equations = {'stage':     'stage',
99                     'momentum':  '(xmomentum**2 + ymomentum**2)**0.5',
100                     'depth':     'stage-elevation',
101                     'speed':     '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-6)',
102                     'elevation': 'elevation' }</font></pre>
103         This dictionary maps your <b>var</b> strings to a variable or expression.
104  </ul>
105
106  <a name="get_timeseries.py"><h3><b>get_timeseries.py</b></h3></a>
107  The <font color="red">get_timeseries.py</font> script is used to get timeseries data for a selection of variable
108  data at one or more timeseries gauge points in a simulation.
109  <p>
110  You must change certain values within the <font color="red">get_timeseries.py</font> file to get what you want:
111  <ul>
112    <li> Depending on how many simulations you have run will determine the number of output folder names you insert here:
113<pre><font color="brown">    time_dir1 = '20090505_150430_run_final_0.8_58292_None_kvanputt'
114    time_dir2 = '20090505_150517_run_final_0_58292_None_kvanputt'</font></pre>
115         as necessary.  Then modify the line:
116<pre><font color="brown">    time_dirs = [time_dir1, time_dir2]</font></pre>
117         to include each of the variables you defined above.
118         <p>
119    <li> Make sure that the <b>gauges</b> list specified in <font color="red">projects.py</font> contains one or more CSV files
120         defining the gauges within the simulation that you want the timeseries data for.
121         The gauge file must have this format:
122<pre><font color="brown">    easting,northing,name,elevation
123    559134.8963,5251176.875,Connelly,2
124    532441.5164,5239639.686,Opossum,2</font></pre>
125  </ul>
126</body>
127</html>
Note: See TracBrowser for help on using the repository browser.