Changeset 7331


Ignore:
Timestamp:
Aug 7, 2009, 1:24:59 PM (15 years ago)
Author:
kristy
Message:

updated all scripts to reflect batemans bay changes

Location:
DVD_images/extra_files/Hobart
Files:
1 deleted
4 edited

Legend:

Unmodified
Added
Removed
  • DVD_images/extra_files/Hobart/analysis.html

    r7301 r7331  
    1616      </td>
    1717      <td>
    18         <H1>Tsunami Inundation Models for the HOBART region</H1>
     18        <H1>Tsunami Inundation Models for south east Tasmania</H1>
    1919      </td>
    2020    </tr>
     
    2727  from SWW files, but are offered as an example of one way to do it.
    2828  <p>
    29   This page describes the scripts and shows how to use them.  Note that due to the ad-hoc nature of the scripts, 
     29  This page describes the scripts and shows how to use them.  Note that due to the ad-hoc nature of the scripts,
    3030  you <b>must</b> modify these scripts before they can be used to extract information from any SWW files you
    3131  generate.
     
    3333  <a name="export_results_max.py scripts"><h3><b>export_results_max.py</b></h3></a>
    3434  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 raster on the simulation region.
     35  of an ANUGA variable or expression for each point in a user-defined ArcGIS grid on the simulation region.
    3636  <p>
    3737  You must change certain values within the <font color="red">export_results_max.py</font> file to get what you want:
    3838  <ul>
    39     <li> You must determine how many SWW files your simulation produced and create as many lines like this:
     39    <li> List the output folder names here, with one for each simulation you have 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>
    42          as necessary.  Then modify the line:
     42         For information on  <a href="modifications.html#output">output_folder_name</a>.
     43  <p>Then modify the line:
    4344<pre><font color="brown">    time_dirs = [time_dir1, time_dir2]</font></pre>
    4445         to include each of the variables you defined above.
    45          <p>   
    46     <li> Your SWW output will have a name like <font color="red">hobart.sww</font>. For a large model you 
    47          may have more than one output file with extra filenames like <font color="red">hobart_37860_0.sww</font>,
    48          where the <b>37860</b> tells us that the timestep at the start of this file is 37860 sec. The 
     46         <p>
     47    <li> Your SWW output will have a name like <font color="red">hobart.sww</font>. For a large model you
     48         may have more than one output file with extra filenames like <font color="red">hobart_16980_0.sww</font>,
     49         where the <b>16980</b> tells us that the timestep at the start of this file is 16980 sec. The
    4950         <font color="red">export_results_max.py</font> script needs to examine all SWW output files, so we specify
    5051         all the files produced by the model run to create a maximum over the entire time.
    51          We do this by specifying the start times of all the extra SWW files:
    52          <p>   
    53 <pre><font color="brown">    times = [37860]</font></pre>
     52         We do this by specifying the start times of all the extra SWW files:
    5453         <p>
    55          Note we do not have to specify the first SWW file, only the extra file start times. If there are no extra SWW files,
    56          just do:
    57          <p>   
     54<pre><font color="brown">    times = [16980, 33960, 50940] </font></pre>
     55         <p>
     56         <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         just do:</i>
     58         <p>
    5859<pre><font color="brown">    times = []</font></pre>
    5960         <p>
    60     <li> Modify the <b>cellsize</b> value to set the size of the raster you require.  For example:
    61 <pre><font color="brown">    cellsize = 20              # raster cell size in metres</font></pre>
    62          <p>   
    63     <li> You can get the maximum of a variable or expression over the entire model time, or for a single time in the simulation.
     61    <li> Modify the <b>cellsize</b> value to set the size of the ArcGIS grid you require.  For example:
     62<pre><font color="brown">    cellsize = 20              # ArcGIS grid cell size in metres</font></pre>
     63         <p>
     64         <i>Note: That this value should not go beyond the most refined part of the mesh
     65         i.e. if mesh resolution is 500m<sup>2</sup>, cellsize = (500 </i>x<i> 2)<sup>1/2</sup> = 32m</i>
     66         <p>
     67    <li> You can get the maximum of a variable or expression over the entire model time, or for a single time in the simulation.
    6468         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
    65          the maximum values in the raster over all timesteps:
     69         the maximum values over all timesteps:
    6670<pre><font color="brown">    timestep = None  # over all timesteps
    6771    #timestep = 0</font></pre>
    68          <p>         
    69          Note that if you want the elevation use <b>timestep = 0</b>. Elevation does not change over time and you only need the
    70          first SWW file.               
    7172         <p>
    72     <li> If you want to clip the raster to one or more small regions then put the names of the regions of interest
     73         <i>Note that if you want the elevation use <b>timestep = 0</b>. Elevation does not change over time and you only need the
     74         first SWW file.</i>
     75         <p>
     76    <li> If you want to clip the ArcGIS grid to one or more small regions then put the names of the regions of interest
    7377         into the <b>area</b> list:
    7478<pre><font color="brown">    area = ['Hobart', 'NW', 'South']</font></pre>
    75          Note that the name strings you put into the <b>area</b> list must match the names used in the
     79         <p>
     80         <i>Note that the name strings you put into the <b>area</b> list must match the names used in the
    7681         <font color="red">project.py</font> file.  That is, if you use a name such as 'NW', then the
    77          <font color="red">project.py</font> file must have lines like:
     82         <font color="red">project.py</font> file must have lines like:</i>
    7883<pre><font color="brown">    xminNW = 548000
    7984    xmaxNW = 561000
     
    8590<pre><font color="brown">    area = ['All']             # no region of interest</font></pre>
    8691         <p>
    87     <li> Finally you must decide which variable or expression values you want sampled on your raster.
     92    <li> Finally you must decide which variable or expression values you want sampled on your ArcGIS grid.
    8893         Define a list <b>var</b> that contains strings defining the required variable or expression:
    8994<pre><font color="brown">    var = ['depth', 'speed']</font></pre>
    90          Note that the strings you supply must be defined in the <b>var_equations</b> dictionary
    91          prior to your definition of <b>var</b>:
     95         <i>Note that the strings you supply must be defined in the <b>var_equations</b> dictionary
     96         prior to your definition of <b>var</b>:</i>
    9297<pre><font color="brown">    var_equations = {'stage':     'stage',
    9398                     'momentum':  '(xmomentum**2 + ymomentum**2)**0.5',
     
    100105  <a name="get_timeseries.py"><h3><b>get_timeseries.py</b></h3></a>
    101106  The <font color="red">get_timeseries.py</font> script is used to get timeseries data for a selection of variable
    102   data at one or more gauge points in a simulation.
     107  data at one or more timeseries gauge points in a simulation.
    103108  <p>
    104109  You must change certain values within the <font color="red">get_timeseries.py</font> file to get what you want:
    105110  <ul>
    106     <li> You must determine how many SWW files your simulation produced and create as many lines like this:
     111    <li> Depending on how many simulations you have run will determine the number of output folder names you insert here:
    107112<pre><font color="brown">    time_dir1 = '20090505_150430_run_final_0.8_58292_None_kvanputt'
    108113    time_dir2 = '20090505_150517_run_final_0_58292_None_kvanputt'</font></pre>
     
    111116         to include each of the variables you defined above.
    112117         <p>
    113     <li> Make sure that the <b>gauges</b> list in <font color="red">projects.py</font> contains one or more CSV files
     118    <li> Make sure that the <b>gauges</b> list specified in <font color="red">projects.py</font> contains one or more CSV files
    114119         defining the gauges within the simulation that you want the timeseries data for.
    115120         The gauge file must have this format:
  • DVD_images/extra_files/Hobart/index.html

    r7300 r7331  
    1616      </td>
    1717      <td>
    18         <H1>Tsunami Inundation Models for the HOBART region</H1>
     18        <H1>Tsunami Inundation Models for south east Tasmania</H1>
    1919      </td>
    2020    </tr>
     
    2323
    2424  <a name="introduction"><h3><b>Introduction</b></h3></a>
    25 
    26   The information on this DVD and the associated report is intended to assist emergency managers in developing preparation and
    27   response plans that can be used during a tsunami threat or emergency.
    28   <p>
    29   The Attorney General's Department (AGD) has supported Geoscience Australia (GA) in developing a range of products to support
    30   the understanding of tsunami hazard through the Australian Probabilistic Tsunami Warning System Project. The work reported here is intended
    31   to further build the capacity of the Jurisdictions in developing inundation models for prioritised locations.
     25  The information within this DVD is intended to assist emergency managers in tsunami planning and preparation
     26  activities. The Australian Government Attorney General’s Department (AGD) has supported Geoscience Australia (GA) in developing a
     27  range of products to support the understanding of tsunami hazard through the Australian Tsunami Warning System Project. The work
     28  reported here is intended to further build the capacity of the Tasmanian State Government in developing inundation models for
     29  prioritised locations.
    3230  <p>
    3331  While the associated report provides background, model results and interpretations, the DVD contains all data, scripts and
    34   software necessary to reproduce and potentially augment the models underpinning the report. This will allow the Jurisdiction
     32  instructions for installing software necessary to reproduce and potentially augment the models underpinning the report. This will allow Tasmanian State Government
    3533  to rerun the models with minor modifications or new elevation data as needed and also potentially replicate the methodology
    36   to other locations within the model extent, if desired. 
     34  for other locations within the model extent, if desired.
    3735  <p>
    3836  The data provided on this DVD is:
    3937  <ul>
    40     <li> The Professional Opinion Report 2009/XX entitled <blink><i><a href="documents/not_there.pdf">Capacity Building for Tsunami Planning and Preparation:
    41          Inundation Models for Four East Coast Australian Communities</a> <b><font color="red">UNFINISHED</font></b></i></blink>
    42     <li>  <a href="documents/Figures">Figures</a> of maximum inundation and maximum speed
    43     <li> Generated <a href="data/tasmania/hobart_tsunami_scenario_2009/anuga/outputs">Arc grids</a> of maximum inundation and maximum speed for the areas of interest
    44     <li> Generated <a href="data/tasmania/hobart_tsunami_scenario_2009/anuga/outputs">timeseries</a> from the models 
     38    <li> The Professional Opinion Report 2009/07 entitled <blink><i><a href="documents/PO_2009_09.pdf">Capacity building for tsunami planning and preparation:
     39         inundation models for south east Tasmania</a></i></blink>
     40    <li>  <a href="documents/Figures">Figures</a> of maximum inundation depth and maximum flow speed
     41    <li> Generated <a href="data/tasmania/hobart_tsunami_scenario_2009/anuga/outputs">ArcGIS grids</a> of maximum inundation depth and maximum flow speed for the areas of interest
     42    <li> Generated <a href="data/tasmania/hobart_tsunami_scenario_2009/anuga/outputs">timeseries</a> from the models
    4543    <li> Combined <a href="data/tasmania/hobart_tsunami_scenario_2009/anuga/topographies">elevation data</a> used by the simulation
    46     <li> The <a href="data/tasmania/hobart_tsunami_scenario_2009/anuga/outputs/elevation">Arc grid</a> of the resulting elevation generated by ANUGA
    47     <li> The Python scripts used to run the models
     44    <li> The <a href="data/tasmania/hobart_tsunami_scenario_2009/anuga/outputs/elevation">ArcGIS grid</a> of the resulting elevation
     45    generated by the model software; <a href="https://datamining.anu.edu.au/anuga">ANUGA</a>
     46    <li> The <a href="project">Python scripts</a> used to run the models
    4847  </ul>
    4948  <p>
    50   The tsunami wave data on this disk are all derived from events available in the Australian Tsunami
    51   Hazard Map. As this dataset is very large, this DVD only contains the events
     49  The input tsunami wave data on this disk are all derived from events available in the Probabilistic Tsunami Hazard Assessment (PTHA) for Australia. As this dataset is very large, this DVD only contains the events
    5250  described in the report.
    5351
  • DVD_images/extra_files/Hobart/installation.html

    r7301 r7331  
    1616      </td>
    1717      <td>
    18         <H1>Tsunami Inundation Models for the HOBART region</H1>
     18        <H1>Tsunami Inundation Models for south east Tasmania</H1>
    1919      </td>
    2020    </tr>
     
    3131
    3232  The tsunami inundation simulations are based on the Open Source software package called ANUGA.
    33   For more information about ANUGA visit <a href="https://datamining.anu.edu.au/anuga">https://datamining.anu.edu.au/anuga</a>. 
     33  For more information about ANUGA visit <a href="https://datamining.anu.edu.au/anuga">https://datamining.anu.edu.au/anuga</a>.
    3434  <p>
    35   ANUGA has been developed and tested under the Windows XP and Linux (Ubuntu, Mint and Red Hat) operating systems. 
     35  ANUGA has been developed and tested under the Windows XP and Linux (Ubuntu, Mint and Red Hat) operating systems.
    3636  ANUGA may also be installed on a Windows Vista system and Debian Linux, but these have not been extensively tested.
    3737  <p>
     
    4040
    4141  <a name="installing"><h3><b>Installing ANUGA</b></h3></a>
    42 
     42  <a href="http://sourceforge.net/projects/anuga/">Download Software</a><i>: please note Internet Explorer has the habit of renaming the
     43  .tgz file to .gz - the remedy is to rename them back or use another browser such as Firefox.</i>
     44  <p>
    4345  You must install ANUGA according to the <a href="documents/anuga_installation_guide.pdf">Installation Guide</a>.
    4446  As the installation procedure comprises a number of steps we suggest that it is carried out by someone with a
    45   basic understanding of how to install software packages on the operating system of choice. 
     47  basic understanding of how to install software packages on the operating system of choice.
    4648
    4749  <a name="simulating"><h3><b>How to run a simulation</b></h3></a>
     
    5254  can be used to get help.
    5355  <p>
    54  
     56
    5557  <p>
    5658  <hr align="left" width="100">
     
    6163  <hr align="left" width="100">
    6264  <p>
    63   Next, you must copy the entire contents of the DVD to a place in your filesystem. Let's assume
    64   that you want to put the data into <b>C:\ANUGA</b>. You would do this in a Command Prompt window:
    65   <table class="code"><tr><td>
    66     C:<br>
    67     cd \<br>
    68     mkdir ANUGA<br>
     65  Next, you must copy the entire contents of the DVD to a place in your file system. Let's assume
     66  that you want to put the data into <b>C:\ANUGA</b>. You could drag all files from the DVD into this folder with Explorer.
     67  Or through the Command Prompt window:
     68  <pre><font color="brown">
     69    C:
     70    cd
     71    mkdir ANUGA
    6972    xcopy /e X:\ ANUGA
    70   </td></tr></table>
    71   <p>
    72   Of course, you could just drag all files on the DVD (X:\ is the assumed drive letter) to the desired target directory with Explorer.
     73  </font></pre>
     74  X:\ is the assumed to be the CD drive letter.
    7375  <p>
    7476  <hr align="left" width="100">
    7577  <p>
    7678  Now you must create an environment variable <b>ANUGADATA</b> that points to the directory you just created:
    77   <table class="code"><tr><td>
     79  <pre><font color="brown">
    7880    set ANUGADATA=C:\ANUGA
    79   </td></tr></table>
    80   <p>
     81  </font></pre>
    8182  Doing the above in a Command Prompt window is temporary - if you open another
    8283  window you will not have the ANUGADATA environment variable defined. To make
     
    8485  tool (<b>Advanced</b> tab).  See the <a href="documents/anuga_installation_guide.pdf">Installation Guide</a>
    8586  if you don't know how to do this.
    86   <p>
    8787  <hr align="left" width="100">
    88   <p>
    8988  Next, you prepare and execute the model:
    90   <table class="code"><tr><td>
    91     cd C:\ANUGA<br>
    92     cd project<br>
    93     python setup_model.py<br >
     89  <pre><font color="brown">
     90    cd C:\ANUGA
     91    cd project
     92    python setup_model.py
    9493    python run_model.py
    95   </td></tr></table>
    96   <p>
     94  </font></pre>
    9795  <hr align="left" width="100">
    98   <p>
    9996  This should all run to completion though it may take several days, depending on
    10097  the hardware used.
     
    105102  drilling down from that point until you find a directory called <b>outputs</b>.
    106103  The latest directory in <b>outputs</b> should contain a file called <font color="red">screen_error.txt</font>
    107   with an explanation of the error at the end. 
     104  with an explanation of the error at the end.
    108105  <p>
    109106  Once the simulation has completed you may view the results using the ANUGA
    110   viewer as described in the <a href="documents/anuga_installation_guide.pdf">Installation Guide</a> 
     107  viewer as described in the <a href="documents/anuga_installation_guide.pdf">Installation Guide</a>
    111108  and <a href="documents/anuga_user_manual.pdf">User Manual</a>.  You may also
    112109  generate rasters suitable for GIS mapping or write dedicated scripts extracting timeseries at selected
     
    120117  After you have run a simulation you may wish to extract information from the generated SWW file.
    121118  On this disk there are a selection of example analysis scripts that you can use or modify.
    122   The <a href="analysis.html">analysis scripts</a> page explains how to use the scripts.
     119  The <a href="analysis.html">analysis</a> page explains how to use the scripts.
    123120
    124121  <a name="modifications"><h3><b>How to modify a simulation</b></h3></a>
     
    138135  <a name="FAQ"><h3><b>Frequently Asked Questions</b></h3></a>
    139136
    140   There is a <b>Frequently Asked Questions</b> (FAQ) page at
    141   <a href="https://datamining.anu.edu.au/anuga/wiki/FrequentlyAskedQuestions">https://datamining.anu.edu.au/anuga/wiki/FrequentlyAskedQuestions</a>.
    142   Look in there and search the mailing list archives before asking questions on the ANUGA user's mailing list.
    143   Someone may already have asked your question!
     137  Look in <b><a href="https://datamining.anu.edu.au/anuga/wiki/FrequentlyAskedQuestions">Frequently Asked Questions</a></b> (FAQ) and search the mailing list archives before asking questions on the ANUGA user's mailing list.
     138  Someone may already have asked and answered your question!
    144139  <p>
    145140  Note that you can edit the FAQ pages, so it is helpful if you can update the FAQ with your question and its answer if you feel
  • DVD_images/extra_files/Hobart/modifications.html

    r7301 r7331  
    1616      </td>
    1717      <td>
    18         <H1>Tsunami Inundation Models for the HOBART region</H1>
     18        <H1>Tsunami Inundation Models for south east Tasmania</H1>
    1919      </td>
    2020    </tr>
     
    2424  <a name="modifications"><h2><b>Modifying a simulation</b></h2></a>
    2525
    26   Here we talk about how you change a simulation.  Why would you need to change a simulation?  Well,
     26  Here we talk about how to change a simulation.  Why would you need to change a simulation?  Well,
    2727  you might have some better elevation data, you want to run the simulation on a finer mesh in certain
    2828  areas, or maybe you want to see the result for different events.
    2929  <p>
    30   First we describe the uses of the script files you might change, then we walk through a few examples
     30  First we describe the usage of the script files you might change, then we walk through a few examples
    3131  of changes you might make.
    3232  <p>
     
    3737  <table class="code">
    3838    <tr><td><font color="red">project.py</font></td><td>Defines the input data used, where to place output, etc.</td></tr>
    39     <tr><td><font color="red">build_elevation.py</font></td><td>Builds a PTS file of all elevation data specified by <font color="red">project.py</font></td></tr>
     39    <tr><td><font color="red">build_elevation.py</font></td><td>Combines the elevation data specified by <font color="red">project.py</font> into one file; with the extension .PTS</td></tr>
    4040    <tr><td><font color="red">setup_model.py</font></td><td>Prepares the simulation before actually running it</td></tr>
    4141    <tr><td><font color="red">run_model.py</font></td><td>Runs the simulation</td></tr>
     
    4848    For further details on changing parameters see <a href="#change">Making changes to a simulation</a>.
    4949
    50 <!--  This file is the heart of the simulation.  Here we specify many things, including the tide level, the simulation
    51   start and stop times, input data files, and so on.  Because we are using python as the programming language in
    52   these files, we can use the power of the language to speed up our development.  We use this by <i>parameterising</i>
    53   data as much as possible.
    54   <p>
    55   Parameters that you might change
    56   As an example, let's look at the <b>scenario_name</b> variable we define in <font color="red">project.py</font>:
    57   <table class="code">
    58     <tr><td>scenario_name = 'hobart'&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;# name of the scenario</td></tr>
    59   </table>
    60   <p>
    61   We define this variable solely to quickly change other input filenames.  Where we want a different input files for
    62   a different scenario,
    63   depending on whether we want a 'small' or 'large' model, we build the filenames using the <b>model</b> variable.
    64   Some examples are:
    65   <p>
    66   <table class="code">
    67     <tr height="25"><td>bounding_polygon_filename = 'bounding_polygon_'+ model + '.csv'</td></tr>
    68     <tr height="25"><td>land_initial_conditions_filename = [['initial_conditions_' + model +'.csv', 0]]</td></tr>
    69     <tr height="25"><td>landward_boundary_filename = 'landward_boundary_'+ model +'.csv'</td></tr>
    70   </table>
    71   <p>
    72   We must, of course, ensure that there are input files <font color="red">bounding_polygon_small.csv</font>
    73   and <font color="red">bounding_polygon_large.csv</font>, and so on.
    74   <p>
    75   Now, if we need to change the simulation model from 'small' to 'large', we need only change the
    76   value of the <b>model</b> variable, and all the input files that depend on the model will change.
    77   This speeds up development and make errors less likely.  If you want to add another model to the two
    78   we show above, just set <b>model</b> to 'medium' and create appropriate input files with the string
    79   'medium' in the filename. -->
    80 
    8150  <a name="build_elevation.py"><h4><b>build_elevation.py</b></h4></a>
    8251
    83   This script combines all input elevation files into one PTS file.
     52  This script combines all input elevation files into a single elevation point file (PTS file).
    8453  We have provided you with the PTS file used to create the outputs on this DVD.
    85   If you would like to change the elevation see <a href="#change">Making changes to a simulation</a>.
     54  If you would like to change the elevation, see <a href="#change">Making changes to a simulation</a>.
    8655  <p>
    8756
     
    9564
    9665  This script runs a tsunami inundation scenario.  It relies on the parameters set in <font color="red">project.py</font>
    97   as well as the PTS and STS files.
    98   An STS file has been generated for each event listed in the boundaries directory. For further details on events see
     66  as well as the elevation and event input files (PTS and STS files respectively).
     67  An STS file has been generated for each event listed in the <u>boundaries</u> directory. For further details on events see
    9968  <a href="#change">Making changes to a simulation</a>.
    10069  <p>
     
    10675  <p>
    10776
     77  <a name="output"><h4><b>Output Folder Name</b></h4></a>
     78
     79  The <b>output folder name</b> should be unique between different runs on different data.
     80  The list of items below will be used to create the folder in your <u>output</u> directory.
     81  Your user name and time+date will be automatically added.  For example,
     82  <pre><font color="brown">
     83  output_comments = [setup, tide, event_number]</font></pre>
     84  will result in a folder name like
     85  <pre><font color="brown">
     86  20090212_091046_run_final_0_27283_rwilson</font>
     87Where you (<u>rwilson</u>) ran a <u>run</u> script at <u>9:10.46</u> in the morning on the <u>2/12/09</u>, <b>setup</b> = <u>final</u>,
     88<b>tide</b> = <u>0</u>, <b>event_number</b> = <u>27283</u>  - refer below for more information on these parameters </pre>
     89  <p>
     90  You can also add strings to this list
     91  <pre><font color="brown">
     92  output_comments = [setup, tide, event_number, 'large']</font></pre>
     93  will result in a folder name like
     94  <pre><font color="brown">
     95  20090212_091046_run_final_0_27283_large_rwilson</font></pre>
     96  <p>
     97  <p>
     98    <a name="setup"><h4><b>Setup</b></h4></a>
     99
     100  The <b>setup</b> parameter determines the type of run. This can be one of three values:
     101  <pre><font color="brown">
     102  'trial' <font color="black">- coarsest mesh, fast </font>
     103  'basic' <font color="black">- coarse mesh</font>
     104  'final' <font color="black">- fine mesh, slowest</font>
     105  </pre></font>
     106  Note: <b>'final'</b> must be used if determining the best estimate of inundation for your area of interest.
     107  <p>
    108108  <a name="tide"><h4><b>Tide</b></h4></a>
    109109
    110   The <b>tide</b> parameter is used to change the mean stage of the simulation.  When <b>tide</b> is set to 0
     110  The <b>tide</b> parameter is used to change the mean inital water level of the simulation.  When <b>tide</b> is set to 0
    111111  the initial water level will be at Mean Sea Level.  If you increase the <b>tide</b> value the water level will become deeper.
    112   However land is masked out using initial conditions.  Within ANUGA <b>tide</b> is modelled as a constant.
     112  This setting will also increase non tidal lakes and rivers inside the model.  To compensate a mask is used on land called
     113  <b>initial conditions</b> which brings the internal water bodies back to 0.  Within ANUGA <b>tide</b> is modelled as a constant.
    113114  <p>
    114115
     
    116117
    117118  The <b>event_number</b> variable contains the event number that initiates the tsunami we are modelling.
    118   You can change <b>event_number</b> to any event number in the anuga/boundaries directory.
    119   As we said above, a new STS file must be generated for <font color="red">run_model.py</font> to work.
    120   <p>
    121   The event numbers correspond to a quake ID from the probabilistic tsunami hazard map of Australia.
    122   <table class="code">
    123     <tr><td><pre>event_number = 51436    # 1 in 10000 yr event from New Hebrides</pre></td></tr>
    124   </table>
     119  You can change <b>event_number</b> to any event number in the <u>boundaries</u> directory.
     120  <a href="data/tasmania/hobart_tsunami_scenario_2009/anuga/boundaries/events_chosen.xls">Click here</a> for more information on the events.
     121  An STS file has been generated for all events listed in the <u>boundaries</u> directory.
     122  <p>
     123  The event numbers correspond to a quake ID from the Probabilistic Tsunami Hazard Map Assessment of Australia.
     124
     125  <pre><font color="brown">
     126  event_number = 58260    # 1 in 10000 yr event from Puysegur Trench</font></pre>
     127
    125128  <p>
    126129
    127130  <a name="elevation"><h4><b>Elevation</b></h4></a>
    128131
    129   Elevation data can be changed in the <font color="red">project.py</font> elevation script.
     132  Elevation data can be changed in the <font color="red">project.py</font> script under ELEVATION DATA.
    130133  Elevation data can be read as either a point file, comma delimited, or as an ASCII grid file
    131   (ASC) with an accompanying projection file (PRJ). All elevation input must be projected in the correct UTM zone.
     134  (ASC) with an accompanying projection file (PRJ). All elevation input should sit in <u>topographies</u> and must be projected in the correct UTM zone.
    132135  <p>
    133136
    134137  A header for a CSV file has the format:
    135   <table class="code">
    136     <tr><td><pre>x,y,elevation</pre></td></tr>
    137   </table>
     138
     139    <pre><font color="brown">
     140  x,y,elevation</font></pre>
     141
    138142  <p>
    139143
    140144  An ASC file header has the format:
    141   <table class="code">
    142     <tr><td><pre>ncols         868
    143 nrows         856
    144 xllcorner     418933.86055096
    145 yllcorner     5151810.6668096
    146 cellsize      250
    147 NODATA_value  -9999</pre></td></tr>
    148   </table>
     145
     146    <pre><font color="brown">
     147  ncols         868
     148  nrows         856
     149  xllcorner     418933.86055096
     150  yllcorner     5151810.6668096
     151  cellsize      250
     152  NODATA_value  -9999</font></pre>
     153
    149154  <p>
    150155
    151156  The header of a PRJ file has the format:
    152   <table class="code">
    153     <tr><td><pre>Projection    UTM
    154 Zone          55
    155 Datum         D_GDA_1994
    156 Zunits        NO
    157 Units         METERS
    158 Spheroid      GRS_1980
    159 Xshift        500000
    160 Yshift        10000000
    161 Parameters</pre></td></tr>
    162   </table>
    163   <p>
    164 
    165   The elevation filenames must be listed in either <b>point_filenames</b> or <b>ascii_grid_filenames</b>
     157
     158    <pre><font color="brown">
     159  Projection    UTM
     160  Zone          55
     161  Datum         D_GDA_1994
     162  Zunits        NO
     163  Units         METERS
     164  Spheroid      GRS_1980
     165  Xshift        500000
     166  Yshift        10000000
     167  Parameters</font></pre>
     168
     169  <p>
     170
     171  The elevation filenames in <font color="red">project.py</font> must be listed in either <b>point_filenames</b> or <b>ascii_grid_filenames</b>
    166172  depending on their format. Point files need to have their extension shown however the ascii grid files have the .asc extension assumed:
    167   <table class="code">
    168     <tr><td><pre>point_filenames = ['SD100031996_jgriffin_clip.csv',
    169                    'tomaga_offshore_AHD_MGA_1997.csv',
    170                    'hobart_BBHD_MGA_1995.csv',
    171                    'hobart_AHD_MGA_2000.csv']
    172 
    173 ascii_grid_filenames = ['sd100031996_p',
    174                         'sd100031996_p2',
    175                         'sd100031996_p3',
    176                         'sd100031996_p4']</pre></td></tr>
    177   </table>
     173
     174    <pre><font color="brown">
     175  point_filenames = ['point1.csv',
     176                     'point2.csv',
     177                     'point3.csv']
     178
     179  ascii_grid_filenames = ['grid1',
     180                          'grid2',
     181                          'grid3']</font></pre>
     182
    178183  <p>
    179184
     
    183188  <a name="interior_regions"><h4><b>Interior regions</b></h4></a>
    184189
    185   The <b>interior_regions</b> parameter allows you to change the mesh of the model.
     190  The user can specify a number of internal polygons within each of which the resolution of the mesh can be specified.
     191  Mesh resolution is the maximum allowable area specified for each region, defining the largest area an indivdual
     192  triangular element of the mesh can take (and therefore the minimum mesh resolution).
     193  These polygons need to be nested within each other with no overlapping edges.
     194  <p>
     195  The <b>interior regions</b> can be changed in the <font color="red">project.py</font> script under INTERIOR REGIONS.
     196  Interior regions can be read as either seperate CSV files for each polygon displayed as a listed paired with its
     197  resolution and/or one CSV file for all polygons, where its resolution is defined within the csv under 'id'.
     198  All file inputs should sit in <u>polygons</u> and must be projected in the correct UTM zone.
     199  <p>
     200  The format for a CSV file with ONE polygon has the format:
     201  <pre><font color="brown">
     202  easting,northing  </font><font color="black"> Note: NO Header  </font></pre>
     203
     204
     205  The header for a CSV file with MANY polygons has the format:
     206  <pre><font color="brown">
     207  easting,northing,id,value</font></pre>
     208  <p>
     209  Where id = polygon number and value = maximum allowable area.
     210  <p>
     211  The <font color="red">project.py</font> script for this section looks like this:
     212
     213
     214  <pre><font color="brown">
     215  interior_regions_list = [['aos1.csv', 1500],
     216                           ['aos2.csv', 1500],
     217                           ['sw.csv', 30000]]
     218  interior_regions_multiple_csv = 'PriorityAreas.csv'</font></pre>
     219
     220  <p>
     221
     222  For further information on ANUGA file formats please see the ANUGA User Manual, section 5.1.
    186223  <p>
    187224</body>
    188225</html>
     226
Note: See TracChangeset for help on using the changeset viewer.