Ignore:
Timestamp:
Nov 13, 2006, 3:03:26 PM (17 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:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/hobart_2006/make_report_boundary_BOM.py

    r3965 r3972  
    22import project
    33from anuga.abstract_2d_finite_volumes.util import sww2timeseries, get_gauges_from_file
     4
     5def 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
     24def 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
    437
    538# Derive scenario name
     
    1750   
    1851#
    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
     52production_dirs = {'20061022_224422': 'Mw 8-7', # refined around paleo sites
    2353                   'MOST': 'MOST'} # MOST input for Mw 8-7
    2454
    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'
     55gauge_map = 'gauges_map_bom.jpg'
    2956
    3057# Create sections and graphs for each designated production directory
     
    4471                                      report = True,
    4572                                      reportname = 'latexoutput_boundary',
    46                                       #reportname = 'latexoutput_boundary_event2',
    4773                                      plot_quantity = ['stage', 'speed'],
    4874                                      surface = False,
     
    5480latex_output.append(texname)
    5581
     82elevfiles = {}
     83for 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
     91figname = 'compare_elev.png'
     92print_elev(elevfiles,production_dirs,report_dir+figname)
     93
    5694# Start report generation
    57 # Future: generate_report(reportdir, scenario, report_title,
    58 # project.gauge_filename, max_maps, damage_maps, production_dirs, latex_output)
    5995report_name = reportdir + 'boundary_report.tex'
    60 #report_name = reportdir + 'boundary_report_event2.tex'
    6196fid = open(report_name, 'w')
    6297
     
    107142
    108143\\begin{document}
    109 \\title{Comparison between ANUGA and MOST}
     144\\title{Comparison between ANUGA and MOST - Hobart}
    110145\maketitle
    111146"""
    112147fid.write(s)
    113148
    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
     149s = '\\begin{figure}[hbt] \n \centerline{ \includegraphics[scale=0.6]{../report_figures/%s}}' %gauge_map
    138150fid.write(s)
    139151
     
    144156"""
    145157fid.write(s)
    146    
     158
     159s = '\\begin{figure}[hbt] \n \centerline{ \includegraphics[width=\paperwidth]{%s}}' %figname
     160fid.write(s)
     161
     162s  = """
     163\caption{Elevation data for ANUGA and MOST.} 
     164\label{fig:elevation}
     165\end{figure}
     166"""
     167fid.write(s)
     168
    147169s = '\input{%s} \n \clearpage \n \n' %latex_output[0]
    148170fid.write(s)
Note: See TracChangeset for help on using the changeset viewer.