Changeset 3015 for production/pt_hedland_2006/make_report.py

Timestamp:

May 30, 2006, 12:00:56 PM (18 years ago)

Author:

sexton

Message:

including references and more updates to report words

File:

• production/pt_hedland_2006/make_report.py (modified) (9 diffs)

production/pt_hedland_2006/make_report.py

@@ -34 +34 @@
 * an introduction must be written in introduction.tex; a basic outline and
   some of the core inputs are already in place
-* an results.tex file needs to be written for the particular scenario
-* the tsunami-genic event should be discussed in tsunami_scenario.tex
+  * the tsunami-genic event should be discussed in tsunami_scenario.tex
+* a computational_setup.tex file needs to be written for the particular scenario
+* the interpretation of the results needs to be written to interpretation.tex
+* maximum inundation maps need to be included in HAT_map and LAT_map etc.
+* damage modelling maps need to be included in HAT_damage and LAT_damage etc.
 * a summary must be written into summary.tex
 * metadata for the scenario data to be included in metadata.tex
@@ -44 +47 @@
 from os import getcwd, sep, altsep, mkdir, access, F_OK
 import project
-from pyvolution.util import sww2timeseries
+from pyvolution.util import sww2timeseries, get_gauges_from_file
 
 # Derive scenario name
@@ -54 +57 @@
     scenario_name = '%s %s' %(test[0], test[1])
 
+# Create report directory
+reportdir = getcwd()+sep+'report'+sep
+if access(reportdir,F_OK) == 0:
+    mkdir (reportdir)
+    
 # User defined inputs
 report_title = 'Tsunami impact modelling for the North West shelf: %s' %scenario_name.title()
@@ -60 +68 @@
                    '': 'Lowest Astronomical Tide'}
 
+max_maps = {'Highest Astronomical Tide': 'HAT_map',
+            'Lowest Astronomical Tide': 'LAT_map'}
+
+damage_maps = {'Highest Astronomical Tide': 'HAT_damage',
+               'Lowest Astronomical Tide': 'LAT_damage'}
+  
 # Create sections and graphs for each designated production directory
 latex_output = []
@@ -71 +85 @@
                              label_id,
                              report = True,
-                             plot_quantity = ['stage', 'velocity', 'bearing'],
+                             plot_quantity = ['stage', 'velocity'],
                              time_min = None,
                              time_max = None,
-                             title_on = False,   
+                             title_on = False,
                              verbose = True)
 
@@ -80 +94 @@
 
 
+
 # Start report generation
-reportdir = getcwd()+sep+'report'+sep
-if access(reportdir,F_OK) == 0:
-    mkdir (reportdir)
 report_name = reportdir + scenario + '_report.tex'
 fid = open(report_name, 'w')
@@ -94 +106 @@
 % * an introduction must be written in introduction.tex; a basic outline and
 %   some of the core inputs are already in place
-% * an results.tex file needs to be written for the particular scenario
 % * the tsunami-genic event should be discussed in tsunami_scenario.tex
+% * an computational_setup.tex file needs to be written for the particular scenario
+% * the interpretation of the results needs to be written to interpretation.tex
+% * maximum inundation maps need to be included in HAT_map.tex and LAT_map.tex etc.
+% * damage modelling maps need to be included in HAT_damage and LAT_damage etc.
 % * a summary must be written into summary.tex
 % * metadata for the scenario data to be included in metadata.tex
@@ -151 +166 @@
   \section{Inundation modelling results}
      \label{sec:results}
-     \input{results}
+     \input{computational_setup}
      
 """
+fid.write(s)
+
+# Generate latex output for gauges
+s = '\\begin{table} \label{table:gaugelocations} \n'
+fid.write(s)
+s = '\caption{Defined gauge locations for %s study area.}' %scenario_name.title()
+fid.write(s)
+s = """
+\\begin{center}
+\\begin{tabular}{|l|l|l|l|}\hline
+\\bf{Gauge Name} & \\bf{Easting} & \\bf{Northing} & \\bf{Elevation}\\\\ \hline
+"""
+fid.write(s)
+
+gauges, locations = get_gauges_from_file(project.gauge_filename)
+
+for name, gauges in zip(locations, gauges):
+    east = gauges[0]
+    north = gauges[1]
+    elev = 0.0
+    #elev = gauges[2]
+    s = '%s & %.2f & %.2f & %.2f \\\\ \hline \n' %(name.replace('_',' '), east, north, elev)
+    fid.write(s)
+
+s = '\\end{tabular} \n  \end{center} \n \end{table} \n \n'
+fid.write(s)
+
+s = '\input{interpretation} \n'
 fid.write(s)
 
@@ -162 +205 @@
     s = '\subsection{%s} \n \n' %production_dirs[name]      
     fid.write(s)
-    
+
+    s = '\input{%s} \n \clearpage \n \n' %max_maps[production_dirs[name]]
+    fid.write(s)
+
     s = '\input{%s} \n \clearpage \n \n' %latex_output[i]
     fid.write(s)
 
-
-# Closing    
+# Closing
+
 s = """
    \section{Damage modelling}
      \input{damage}
-
-    \section{Summary}
+"""
+fid.write(s)
+
+for i, name in enumerate(production_dirs.keys()):
+
+    s = '\subsection{%s} \n \n' %production_dirs[name]      
+    fid.write(s)
+
+    s = '\input{%s} \n \clearpage \n \n' %damage_maps[production_dirs[name]]
+    fid.write(s)
+
+
+s = """
+   \section{Summary}
      \input{summary}
      
    \section{References}
-
+    \input{references}
+    
    \section{Metadata}
      \label{sec:metadata}