source: production/onslow_2006/MOST_timeseries.py @ 3181

"""
Generate latex output comparing time series of given sww files

Inputs:

output dirs: dictionary of output directories with a
             association to simulation, eg 50m boundary, 100m boundary
                    
Outputs:

* latex output to generate figures for report 
* figures used for report stored in the report_figure directory
NOTE, this directory will need to be committed, as well as
the latex files.

June 2006                    
"""

from os import sep, getcwd, F_OK, access, mkdir
import project
from pyvolution.util import sww2timeseries, get_gauges_from_file

# User defined inputs
production_dirs = {'SU-AU_clip': 'MOST'}

gauge_map = 'boundarycomparison.png'
plot_quantity = ['stage', 'speed']

# Create sections and graphs for each designated production directory
latex_output = []
swwfiles = {}
for label_id in production_dirs.keys():
    
    file_loc = project.boundarydir
    swwfile = file_loc + project.boundary_basename + '.sww'
    swwfiles[swwfile] = label_id

texname, vec = sww2timeseries(swwfiles,
                              project.gauge_comparison,
                              production_dirs,
                              report = True,
                              reportname = 'latexoutput_MOST_input',
                              plot_quantity = plot_quantity,
                              time_min = None,
                              time_max = None,
                              title_on = False,
                              verbose = True)

latex_output.append(texname)


# Start report generation
reportdir = getcwd()+sep+'report'+sep
if access(reportdir,F_OK) == 0:
    mkdir (reportdir)
input_name = reportdir + sep + 'MOST_input_onslow.tex'
fid = open(input_name, 'w')

# Generate latex output for location points
#s = '\\begin{table} \label{table:MOSTlocationsonslow} \n'
#fid.write(s)
#s = '\caption{Defined point locations for comparison study for Onslow region.}' 
#fid.write(s)
#s = """
#\\begin{center}
#\\begin{tabular}{|l|l|l|l|}\hline
#\\bf{Point Name} & \\bf{Easting} & \\bf{Northing} & \\bf{Elevation (m)}\\\\ \hline
#"""
#fid.write(s)

#for i, thisvec in enumerate(vec):
#    name = thisvec[0]
#    east = thisvec[1]
#    north = thisvec[2]
#    elev = thisvec[3]
#    s = '%s & %.2f & %.2f & %.2f \\\\ \hline \n' %(name, east, north, elev)
#    fid.write(s)
#    for which_one in plot_quantity:
#        move ('report_figures' + sep + 'gauge'+name.replace(' ','')+which_one+'.png', project.comparereportfigdir)

#s = '\\end{tabular} \n  \end{center} \n \end{table} \n \n'
#fid.write(s)

s = '\\begin{figure}[hbt] \n \centerline{ \includegraphics[width=150mm, height=100mm]{../report_figures/%s}}' %gauge_map
fid.write(s)

s  = """
\caption{Point locations used for MOST input for Onslow region.}  
\label{fig:boundarypoints}
\end{figure}
"""
fid.write(s)
    
s = '\input{%s} \n \clearpage \n \n' %latex_output[0]
fid.write(s)