Ignore:
Timestamp:
Aug 14, 2009, 9:03:48 AM (16 years ago)
Author:
Leharne
Message:

changes to gold coast html files

File:
1 edited

Legend:

Unmodified
Added
Removed

  • DVD_images/extra_files/GoldCoast/analysis.html

    r7301 r7363  
    1616      </td>
    1717      <td>
    18         <H1>Tsunami Inundation Models for the GOLD COAST region</H1>
     18        <H1>Tsunami Inundation Models for the Gold Coast region</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:
    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          as necessary.  Then modify the line:
     39    <li> List the output folder names here, with one for each simulation you have run:
     40<pre><font color="brown">    time_dir1 = '20090522_165600_run_final_1.1_51423_lfountai'
     41    time_dir2 = '20090522_164640_run_final_1.1_51392_lfountai'</font></pre>
     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">gold_coast.sww</font>. For a large model you 
    47          may have more than one output file with extra filenames like <font color="red">gold_coast_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">gold_coast.sww</font>. For a large model you
     48         may have more than one output file with extra filenames like <font color="red">gold_coast_28860_0.sww</font>,
     49         where the <b>28860</b> tells us that the timestep at the start of this file is 28860 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 = [28860, 57720] </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:
    74 <pre><font color="brown">    area = ['Gold_Coast', '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
    76          <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:
    78 <pre><font color="brown">    xminNW = 548000
    79     xmaxNW = 561000
    80     yminNW = 5250000
    81     ymaxNW = 5258000</font></pre>
     78<pre><font color="brown">    area = ['PalmBeach', 'MainBeach', 'CBD']</font></pre>
     79         <p>
     80         <i>Note that the name strings you put into the <b>area</b> list must match the names used in the
     81         <font color="red">project.py</font> file.  That is, if you use a name such as 'PalmBeach', then the
     82         <font color="red">project.py</font> file must have lines like:</i>
     83<pre><font color="brown">    xminPalmBeach = 548000
     84    xmaxPalmBeach = 561000
     85    yminPalmBeach = 5250000
     86    ymaxPalmBeach = 5258000</font></pre>
    8287         which set the maximum and minimum eastings and northings that define a rectangular region.
    8388         <p>
    … …  
    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:
    107 <pre><font color="brown">    time_dir1 = '20090505_150430_run_final_0.8_58292_None_kvanputt'
    108     time_dir2 = '20090505_150517_run_final_0_58292_None_kvanputt'</font></pre>
     111    <li> Depending on how many simulations you have run will determine the number of output folder names you insert here:
     112<pre><font color="brown">    time_dir1 = '20090522_165600_run_final_1.1_51423_lfountai'
     113    time_dir2 = '20090522_164640_run_final_1.1_51392_lfountai'</font></pre>
    109114         as necessary.  Then modify the line:
    110115<pre><font color="brown">    time_dirs = [time_dir1, time_dir2]</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:
    116121<pre><font color="brown">    easting,northing,name,elevation
    117     559134.8963,5251176.875,Connelly,2
    118     532441.5164,5239639.686,Opossum,2</font></pre>
     122    550154.557,6885684.489,Coolangatta,-5
     123    546917.144,6889565.504,PalmBeach,-5</font></pre>
    119124  </ul>
    120125</body>
Note: See TracChangeset for help on using the changeset viewer.