Ignore:
Timestamp:
Feb 11, 2010, 4:43:41 PM (15 years ago)
Author:
griffin
Message:

Edits to html pages

File:
1 edited

Legend:

Unmodified
Added
Removed

  • DVD_images/extra_files/Geraldton/analysis.html

    r7514 r7620  
    2323
    2424  <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
     25  There are number of analysis scripts on this disk that can be used to extract information from the ANUGA output (SWW) files
    2626  produced by a simulation run.  These scripts are not meant to be the definitive answer on how to get information
    2727  from SWW files, but are offered as an example of one way to do it.
    … …  
    3232
    3333  <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.
     34  The <font color="red">export_results_max.py</font> script can be used to create an ASCI (.asc) file that contains the maximum
     35  of an ANUGA variable or expression for each point in a user-defined grid on the simulation region. ASCI files can then be imported into a GIS environment.
    3636  <p>
    37   You must change certain values within the <font color="red">export_results_max.py</font> file to get what you want:
     37  Certain values within the <font color="red">export_results_max.py</font> file must be changed to align with the correct output files:
    3838  <ul>
    39     <li> List the output folder names here, with one for each simulation you have run:
     39    <li> List the output folder names here, with one for each simulation run:
    4040<pre><font color="brown">    time_dir1 = '20090505_150430_run_final_0.8_58292_None_kvanputt'
    4141    time_dir2 = '20090505_150517_run_final_0_58292_None_kvanputt'</font></pre>
    … …  
    4343  <p>Then modify the line:
    4444<pre><font color="brown">    time_dirs = [time_dir1, time_dir2]</font></pre>
    45          to include each of the variables you defined above.
     45         to include each of the variables defined above.
    4646         <p>
    47     <li> Your SWW output will have a name like <font color="red">geraldton.sww</font>. For a large model you
     47    <li> The SWW output will have a name like <font color="red">geraldton.sww</font>. For a large model you
    4848         may have more than one output file with extra filenames like <font color="red">geraldton_37860_0.sww</font>,
    4949         where the <b>37860</b> tells us that the timestep at the start of this file is 37860 sec. The
    … …  
    5555    times = [20000, 40000, 60000]</font></pre>
    5656         <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,
     57         <i>Note that it is not necessary to specify the first SWW file, only the extra file start times. If there are no extra SWW files,
    5858         just do:</i>
    5959         <p>
    6060<pre><font color="brown">    times = []</font></pre>
    6161         <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>
     62    <li> Modify the <b>cellsize</b> value to set the resolution of the output grid required.  For example:
     63<pre><font color="brown">    cellsize = 20              # Grid cell size in metres</font></pre>
    6464         <p>
    6565         <i>Note: That this value should not go beyond the most refined part of the mesh
    6666         i.e. if mesh resolution is 500m<sup>2</sup>, cellsize = (500 </i>x<i> 2)<sup>1/2</sup> = 32m</i>
    6767         <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
     68    <li> Maximum values of a variable or expression may be extracted over the entire model time or for a single time in the simulation.
     69         Set <b>timestep</b> to the required time to extract the maximum for, or set it to <b>None</b> to extract
    7070         the maximum values over all timesteps:
    7171<pre><font color="brown">    timestep = None  # over all timesteps
    7272    #timestep = 0</font></pre>
    7373         <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>
     74         <i>Note that if extracting elevation use <b>timestep = 0</b>. Elevation does not change over time and therefore only the
     75         first SWW file is needed.</i>
    7676         <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
     77    <li> To clip the ASCI grid to one or more small regions then put the names of the regions of interest
    7878         into the <b>area</b> list:
    7979<pre><font color="brown">    area = ['geraldton', 'NW', 'South']</font></pre>
    8080         <p>
    81          <i>Note that the name strings you put into the <b>area</b> list must match the names used in the
     81         <i>Note that the name strings in the <b>area</b> list must match the names used in the
    8282         <font color="red">project.py</font> file.  That is, if you use a name such as 'NW', then the
    8383         <font color="red">project.py</font> file must have lines like:</i>
    … …  
    9191<pre><font color="brown">    area = ['All']             # no region of interest</font></pre>
    9292         <p>
    93     <li> Finally you must decide which variable or expression values you want sampled on your ArcGIS grid.
     93    <li> Finally the variable or expression values to be sampled on your ASCI grid must be specified.
    9494         Define a list <b>var</b> that contains strings defining the required variable or expression:
    9595<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>
     96         <i>Note that the strings supplied must be defined in the <b>var_equations</b> dictionary
     97         prior to definition of <b>var</b>:</i>
    9898<pre><font color="brown">    var_equations = {'stage':     'stage',
    9999                     'momentum':  '(xmomentum**2 + ymomentum**2)**0.5',
    … …  
    101101                     'speed':     '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-6)',
    102102                     'elevation': 'elevation' }</font></pre>
    103          This dictionary maps your <b>var</b> strings to a variable or expression.
     103         This dictionary maps the <b>var</b> strings to a variable or expression.
    104104  </ul>
    105105
    … …  
    108108  data at one or more timeseries gauge points in a simulation.
    109109  <p>
    110   You must change certain values within the <font color="red">get_timeseries.py</font> file to get what you want:
     110  Certain values within the <font color="red">export_results_max.py</font> file must be changed to align with the correct output files:
    111111  <ul>
    112     <li> Depending on how many simulations you have run will determine the number of output folder names you insert here:
     112    <li> List the output folder names here, with one for each simulation run:
    113113<pre><font color="brown">    time_dir1 = '20090505_150430_run_final_0.8_58292_None_kvanputt'
    114114    time_dir2 = '20090505_150517_run_final_0_58292_None_kvanputt'</font></pre>
    115          as necessary. Then modify the line:
     115          Then modify the line:
    116116<pre><font color="brown">    time_dirs = [time_dir1, time_dir2]</font></pre>
    117          to include each of the variables you defined above.
     117         to include each of the variables defined above.
    118118         <p>
    119119    <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.
     120         defining the synthetic gauges within the simulation that the timeseries data will be extracted for.
    121121         The gauge file must have this format:
    122122<pre><font color="brown">    easting,northing,name,elevation
    … …  
    124124    532441.5164,5239639.686,Opossum,2</font></pre>
    125125  </ul>
     126    <p>
     127  <hr>
     128   <a rel="license" href="http://creativecommons.org/licenses/by-nc/3.0/"><img alt="Creative Commons License" style="border-width:0" src="http://i.creativecommons.org/l/by-nc/3.0/88x31.png" /></a><br />
     129  © Commonwealth of Australia (Geoscience Australia) 2010.
     130  This material is licensed under a <a rel="license" href="http://creativecommons.org/licenses/by-nc/3.0/">Creative Commons Attribution-Noncommercial
     131  3.0 Australia License</a>. 
    126132</body>
    127133</html>
Note: See TracChangeset for help on using the changeset viewer.