Changeset 3972 for anuga_work/production/hobart_2006/make_report_boundary_BOM.py

Timestamp:

Nov 13, 2006, 3:03:26 PM (16 years ago)

Author:

sexton

Message:

(i) incorporating new supply of interpolated data for Broome (ii) updating report to look at MOST versus ANUGA for Hobart scenario

File:

• anuga_work/production/hobart_2006/make_report_boundary_BOM.py (modified) (6 diffs)

anuga_work/production/hobart_2006/make_report_boundary_BOM.py

@@ -2 +2 @@
 import project
 from anuga.abstract_2d_finite_volumes.util import sww2timeseries, get_gauges_from_file
+
+def print_elev(elevfiles,production_dirs,figname):
+    from pylab import plot, xlabel, ylabel, hold, ion, legend, savefig, close
+    ion()
+    hold(True)
+    cstr = ['g', 'r', 'b', 'c', 'm', 'y', 'k']
+    j = -1
+    leg_label = []
+    for elevfile in elevfiles.keys():
+        x, elev = read_file(elevfile)
+        leg_label.append(production_dirs[elevfiles[elevfile]])
+        j += 1
+        plot(x, elev, '-', c = cstr[j])
+    xlabel('Easting')
+    ylabel('Water Depth')
+    legend((leg_label),loc='upper right')
+    savefig(figname)
+    close('all')
+    return
+
+def read_file(filename):
+    fid = open(filename)
+    lines = fid.readlines()
+    fid.close()
+    x = []
+    elev = []
+    i = -1
+    for line in lines[:]:
+        i += 1
+        fields = line.split(',')
+        x.append(float(fields[3].strip(' ')))
+        elev.append(float(fields[2].strip(' ')))
+    return x, elev
 
 # Derive scenario name
@@ -17 +50 @@
     
 # 
-production_dirs = {#'20061008_234702': 'Mw 8-7', # 2500 res, 750000 other res
-                   #'20061016_065743': 'Mw 8-7 refined', #new interior res
-                   #'20061017_004409': 'Mw 8-7 refined', #new interior res
-                   '20061022_224422': 'Mw 8-7', #new interior res 
+production_dirs = {'20061022_224422': 'Mw 8-7', # refined around paleo sites 
                    'MOST': 'MOST'} # MOST input for Mw 8-7
 
-#production_dirs = {'20061006_062319': 'Mw 8-5', # 2500 res, 750000 other res
-#                   'MOST': 'MOST'} # MOST input for Mw 8-5
-
-gauge_map = 'boundary_gauges_bom.jpg'
+gauge_map = 'gauges_map_bom.jpg'
 
 # Create sections and graphs for each designated production directory
@@ -44 +71 @@
                                       report = True,
                                       reportname = 'latexoutput_boundary',
-                                      #reportname = 'latexoutput_boundary_event2',
                                       plot_quantity = ['stage', 'speed'],
                                       surface = False,
@@ -54 +80 @@
 latex_output.append(texname)
 
+elevfiles = {}
+for label_id in production_dirs.keys():
+    
+    file_loc = project.outputdir + label_id + sep
+    elevfile = file_loc + 'gauges_maxmins' + '.csv'
+    if label_id == 'MOST':
+        elevfile = project.boundarydir + 'gauges_maxmins' + '.csv'
+    elevfiles[elevfile] = label_id
+
+figname = 'compare_elev.png'
+print_elev(elevfiles,production_dirs,report_dir+figname)
+
 # Start report generation
-# Future: generate_report(reportdir, scenario, report_title, 
-# project.gauge_filename, max_maps, damage_maps, production_dirs, latex_output)
 report_name = reportdir + 'boundary_report.tex'
-#report_name = reportdir + 'boundary_report_event2.tex'
 fid = open(report_name, 'w')
 
@@ -107 +142 @@
 
 \\begin{document}
-\\title{Comparison between ANUGA and MOST}
+\\title{Comparison between ANUGA and MOST - Hobart}
 \maketitle
 """
 fid.write(s)
 
-# Generate latex output for location points
-##s = '\\begin{table} \\begin{center} \n'
-##fid.write(s)
-##s = '\caption{Defined point locations for %s study area.}' %scenario_name
-##fid.write(s)
-##s = """
-##\label{table:locations}
-##\\begin{tabular}{|l|l|l|l|}\hline
-##\\bf{Point Name} & \\bf{Easting} & \\bf{Northing} & \\bf{Elevation}\\\\ \hline
-##"""
-##fid.write(s)
-##
-##gauges, locations, elevation = get_gauges_from_file(project.gauge_filename)
-##
-##for name, gauges, elev in zip(locations, gauges, elevation):
-##    east = gauges[0]
-##    north = gauges[1]
-##    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 = '\\begin{figure}[hbt] \n \centerline{ \includegraphics[width=\paperwidth]{../report_figures/%s}}' %gauge_map
+s = '\\begin{figure}[hbt] \n \centerline{ \includegraphics[scale=0.6]{../report_figures/%s}}' %gauge_map
 fid.write(s)
 
@@ -144 +156 @@
 """
 fid.write(s)
-    
+
+s = '\\begin{figure}[hbt] \n \centerline{ \includegraphics[width=\paperwidth]{%s}}' %figname
+fid.write(s)
+
+s  = """
+\caption{Elevation data for ANUGA and MOST.}  
+\label{fig:elevation}
+\end{figure}
+"""
+fid.write(s)
+
 s = '\input{%s} \n \clearpage \n \n' %latex_output[0]
 fid.write(s)