Index: /DVD_images/extra_files/Geraldton/analysis.html
===================================================================
--- /DVD_images/extra_files/Geraldton/analysis.html (revision 7619)
+++ /DVD_images/extra_files/Geraldton/analysis.html (revision 7620)
@@ -23,5 +23,5 @@
Analysis Scripts
- There are number of analysis scripts on this disk that you can use to extract information from the SWW files
+ There are number of analysis scripts on this disk that can be used to extract information from the ANUGA output (SWW) files
produced by a simulation run. These scripts are not meant to be the definitive answer on how to get information
from SWW files, but are offered as an example of one way to do it.
@@ -32,10 +32,10 @@
export_results_max.py
- The export_results_max.py script can be used to create an ASC file that contains the maximum
- of an ANUGA variable or expression for each point in a user-defined ArcGIS grid on the simulation region.
+ The export_results_max.py script can be used to create an ASCI (.asc) file that contains the maximum
+ 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.
- You must change certain values within the export_results_max.py file to get what you want: + Certain values within the export_results_max.py file must be changed to align with the correct output files:
time_dir1 = '20090505_150430_run_final_0.8_58292_None_kvanputt' time_dir2 = '20090505_150517_run_final_0_58292_None_kvanputt'@@ -43,7 +43,7 @@
Then modify the line:
time_dirs = [time_dir1, time_dir2]- to include each of the variables you defined above. + to include each of the variables defined above.
-
- Note we do not have to specify the first SWW file, only the extra file start times. If there are no extra SWW files, + 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, just do:
times = []
-
cellsize = 20 # ArcGIS grid cell size in metres+
cellsize = 20 # Grid cell size in metres
Note: That this value should not go beyond the most refined part of the mesh i.e. if mesh resolution is 500m2, cellsize = (500 x 2)1/2 = 32m
-
timestep = None # over all timesteps #timestep = 0
- Note that if you want the elevation use timestep = 0. Elevation does not change over time and you only need the - first SWW file. + Note that if extracting elevation use timestep = 0. Elevation does not change over time and therefore only the + first SWW file is needed.
-
area = ['geraldton', 'NW', 'South']
- Note that the name strings you put into the area list must match the names used in the
+ Note that the name strings in the area list must match the names used in the
project.py file. That is, if you use a name such as 'NW', then the
project.py file must have lines like:
@@ -91,9 +91,9 @@
area = ['All'] # no region of interest
-
var = ['depth', 'speed']- Note that the strings you supply must be defined in the var_equations dictionary - prior to your definition of var: + Note that the strings supplied must be defined in the var_equations dictionary + prior to definition of var:
var_equations = {'stage': 'stage', 'momentum': '(xmomentum**2 + ymomentum**2)**0.5', @@ -101,5 +101,5 @@ 'speed': '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-6)', 'elevation': 'elevation' }- This dictionary maps your var strings to a variable or expression. + This dictionary maps the var strings to a variable or expression.
- You must change certain values within the get_timeseries.py file to get what you want: + Certain values within the export_results_max.py file must be changed to align with the correct output files:
time_dir1 = '20090505_150430_run_final_0.8_58292_None_kvanputt' time_dir2 = '20090505_150517_run_final_0_58292_None_kvanputt'- as necessary. Then modify the line: + Then modify the line:
time_dirs = [time_dir1, time_dir2]- to include each of the variables you defined above. + to include each of the variables defined above.
easting,northing,name,elevation @@ -124,4 +124,10 @@ 532441.5164,5239639.686,Opossum,2
+