Ignore:
Timestamp:
Aug 7, 2009, 2:39:09 PM (16 years ago)
Author:
jgriffin
Message:

updated html pages to be consistent with other models

File:
1 edited

Legend:

Unmodified
Added
Removed

  • DVD_images/extra_files/Gosford/analysis.html

    r7301 r7336  
    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">gosford.sww</font>. For a large model you 
     46         <p>
     47    <li> Your SWW output will have a name like <font color="red">gosford.sww</font>. For a large model you
    4748         may have more than one output file with extra filenames like <font color="red">gosford_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 
     49         where the <b>37860</b> tells us that the timestep at the start of this file is 37860 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 = [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>
    5860<pre><font color="brown">    times = []</font></pre>
    5961         <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.
     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.
    6469         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:
     70         the maximum values over all timesteps:
    6671<pre><font color="brown">    timestep = None  # over all timesteps
    6772    #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.               
    7173         <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
     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
    7378         into the <b>area</b> list:
    74 <pre><font color="brown">    area = ['Batemans_Bay', '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<pre><font color="brown">    area = ['gosford', '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
    7682         <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:
     83         <font color="red">project.py</font> file must have lines like:</i>
    7884<pre><font color="brown">    xminNW = 548000
    7985    xmaxNW = 561000
    … …  
    8591<pre><font color="brown">    area = ['All']             # no region of interest</font></pre>
    8692         <p>
    87     <li> Finally you must decide which variable or expression values you want sampled on your raster.
     93    <li> Finally you must decide which variable or expression values you want sampled on your ArcGIS grid.
    8894         Define a list <b>var</b> that contains strings defining the required variable or expression:
    8995<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>:
     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>
    9298<pre><font color="brown">    var_equations = {'stage':     'stage',
    9399                     'momentum':  '(xmomentum**2 + ymomentum**2)**0.5',
    … …  
    100106  <a name="get_timeseries.py"><h3><b>get_timeseries.py</b></h3></a>
    101107  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.
     108  data at one or more timeseries gauge points in a simulation.
    103109  <p>
    104110  You must change certain values within the <font color="red">get_timeseries.py</font> file to get what you want:
    105111  <ul>
    106     <li> You must determine how many SWW files your simulation produced and create as many lines like this:
     112    <li> Depending on how many simulations you have run will determine the number of output folder names you insert here:
    107113<pre><font color="brown">    time_dir1 = '20090505_150430_run_final_0.8_58292_None_kvanputt'
    108114    time_dir2 = '20090505_150517_run_final_0_58292_None_kvanputt'</font></pre>
    … …  
    111117         to include each of the variables you defined above.
    112118         <p>
    113     <li> Make sure that the <b>gauges</b> list in <font color="red">projects.py</font> contains one or more CSV files
     119    <li> Make sure that the <b>gauges</b> list specified in <font color="red">projects.py</font> contains one or more CSV files
    114120         defining the gauges within the simulation that you want the timeseries data for.
    115121         The gauge file must have this format:
Note: See TracChangeset for help on using the changeset viewer.