- Timestamp:
- Nov 13, 2006, 3:03:26 PM (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/hobart_2006/make_report_boundary_BOM.py
r3965 r3972 2 2 import project 3 3 from anuga.abstract_2d_finite_volumes.util import sww2timeseries, get_gauges_from_file 4 5 def print_elev(elevfiles,production_dirs,figname): 6 from pylab import plot, xlabel, ylabel, hold, ion, legend, savefig, close 7 ion() 8 hold(True) 9 cstr = ['g', 'r', 'b', 'c', 'm', 'y', 'k'] 10 j = -1 11 leg_label = [] 12 for elevfile in elevfiles.keys(): 13 x, elev = read_file(elevfile) 14 leg_label.append(production_dirs[elevfiles[elevfile]]) 15 j += 1 16 plot(x, elev, '-', c = cstr[j]) 17 xlabel('Easting') 18 ylabel('Water Depth') 19 legend((leg_label),loc='upper right') 20 savefig(figname) 21 close('all') 22 return 23 24 def read_file(filename): 25 fid = open(filename) 26 lines = fid.readlines() 27 fid.close() 28 x = [] 29 elev = [] 30 i = -1 31 for line in lines[:]: 32 i += 1 33 fields = line.split(',') 34 x.append(float(fields[3].strip(' '))) 35 elev.append(float(fields[2].strip(' '))) 36 return x, elev 4 37 5 38 # Derive scenario name … … 17 50 18 51 # 19 production_dirs = {#'20061008_234702': 'Mw 8-7', # 2500 res, 750000 other res 20 #'20061016_065743': 'Mw 8-7 refined', #new interior res 21 #'20061017_004409': 'Mw 8-7 refined', #new interior res 22 '20061022_224422': 'Mw 8-7', #new interior res 52 production_dirs = {'20061022_224422': 'Mw 8-7', # refined around paleo sites 23 53 'MOST': 'MOST'} # MOST input for Mw 8-7 24 54 25 #production_dirs = {'20061006_062319': 'Mw 8-5', # 2500 res, 750000 other res 26 # 'MOST': 'MOST'} # MOST input for Mw 8-5 27 28 gauge_map = 'boundary_gauges_bom.jpg' 55 gauge_map = 'gauges_map_bom.jpg' 29 56 30 57 # Create sections and graphs for each designated production directory … … 44 71 report = True, 45 72 reportname = 'latexoutput_boundary', 46 #reportname = 'latexoutput_boundary_event2',47 73 plot_quantity = ['stage', 'speed'], 48 74 surface = False, … … 54 80 latex_output.append(texname) 55 81 82 elevfiles = {} 83 for label_id in production_dirs.keys(): 84 85 file_loc = project.outputdir + label_id + sep 86 elevfile = file_loc + 'gauges_maxmins' + '.csv' 87 if label_id == 'MOST': 88 elevfile = project.boundarydir + 'gauges_maxmins' + '.csv' 89 elevfiles[elevfile] = label_id 90 91 figname = 'compare_elev.png' 92 print_elev(elevfiles,production_dirs,report_dir+figname) 93 56 94 # Start report generation 57 # Future: generate_report(reportdir, scenario, report_title,58 # project.gauge_filename, max_maps, damage_maps, production_dirs, latex_output)59 95 report_name = reportdir + 'boundary_report.tex' 60 #report_name = reportdir + 'boundary_report_event2.tex'61 96 fid = open(report_name, 'w') 62 97 … … 107 142 108 143 \\begin{document} 109 \\title{Comparison between ANUGA and MOST }144 \\title{Comparison between ANUGA and MOST - Hobart} 110 145 \maketitle 111 146 """ 112 147 fid.write(s) 113 148 114 # Generate latex output for location points 115 ##s = '\\begin{table} \\begin{center} \n' 116 ##fid.write(s) 117 ##s = '\caption{Defined point locations for %s study area.}' %scenario_name 118 ##fid.write(s) 119 ##s = """ 120 ##\label{table:locations} 121 ##\\begin{tabular}{|l|l|l|l|}\hline 122 ##\\bf{Point Name} & \\bf{Easting} & \\bf{Northing} & \\bf{Elevation}\\\\ \hline 123 ##""" 124 ##fid.write(s) 125 ## 126 ##gauges, locations, elevation = get_gauges_from_file(project.gauge_filename) 127 ## 128 ##for name, gauges, elev in zip(locations, gauges, elevation): 129 ## east = gauges[0] 130 ## north = gauges[1] 131 ## s = '%s & %.2f & %.2f & %.2f \\\\ \hline \n' %(name.replace('_',' '), east, north, elev) 132 ## fid.write(s) 133 ## 134 ##s = '\\end{tabular} \n \end{center} \n \end{table} \n \n' 135 ##fid.write(s) 136 137 s = '\\begin{figure}[hbt] \n \centerline{ \includegraphics[width=\paperwidth]{../report_figures/%s}}' %gauge_map 149 s = '\\begin{figure}[hbt] \n \centerline{ \includegraphics[scale=0.6]{../report_figures/%s}}' %gauge_map 138 150 fid.write(s) 139 151 … … 144 156 """ 145 157 fid.write(s) 146 158 159 s = '\\begin{figure}[hbt] \n \centerline{ \includegraphics[width=\paperwidth]{%s}}' %figname 160 fid.write(s) 161 162 s = """ 163 \caption{Elevation data for ANUGA and MOST.} 164 \label{fig:elevation} 165 \end{figure} 166 """ 167 fid.write(s) 168 147 169 s = '\input{%s} \n \clearpage \n \n' %latex_output[0] 148 170 fid.write(s)
Note: See TracChangeset
for help on using the changeset viewer.