Changeset 3188 for production/onslow_2006/MOST_timeseries.py

Timestamp:

Jun 20, 2006, 5:59:42 PM (18 years ago)

Author:

sexton

Message:

make gauges to compare MOST and ANUGA time series

File:

• production/onslow_2006/MOST_timeseries.py (modified) (6 diffs)

production/onslow_2006/MOST_timeseries.py

@@ -4 +4 @@
 Inputs:
 
-output dirs: dictionary of output directories with a
-             association to simulation, eg 50m boundary, 100m boundary
-                    
+production dirs: dictionary of output directories with a
+                 association to simulation, eg 50m boundary, 100m boundary
+
+gauge_map:       graphical output from make_gauges.py
+
 Outputs:
 
@@ -22 +24 @@
 
 # User defined inputs
-production_dirs = {'SU-AU_clip': 'MOST'}
+production_dirs = {'20060515_001733': '100m boundary',
+                   #'20060530_102753': '50m boundary'}#,
+                   'MOST': 'MOST'}
 
 gauge_map = 'boundarycomparison.png'
-plot_quantity = ['stage', 'speed']
 
 # Create sections and graphs for each designated production directory
@@ -32 +35 @@
 for label_id in production_dirs.keys():
     
-    file_loc = project.boundarydir
-    swwfile = file_loc + project.boundary_basename + '.sww'
+    file_loc = project.outputdir + label_id + sep
+    swwfile = file_loc + project.basename + '.sww'
+    if label_id == 'MOST':
+        swwfile = project.boundarydir + project.boundary_basename + '.sww'
     swwfiles[swwfile] = label_id
+
+print 'swwfiles', swwfiles
 
 texname, vec = sww2timeseries(swwfiles,
@@ -40 +47 @@
                               production_dirs,
                               report = True,
-                              reportname = 'latexoutput_MOST_input',
-                              plot_quantity = plot_quantity,
+                              reportname = 'latexoutput_MOST_ANUGA_comparison',
+                              plot_quantity = ['stage', 'speed'],
+                              surface = False,
                               time_min = None,
                               time_max = None,
@@ -49 +57 @@
 latex_output.append(texname)
 
+from shutil import copy, move
+copy ('report' + sep + texname + '.tex', project.comparereportdir + sep + texname + '.tex')
 
 # Start report generation
-reportdir = getcwd()+sep+'report'+sep
-if access(reportdir,F_OK) == 0:
-    mkdir (reportdir)
-input_name = reportdir + sep + 'MOST_input_onslow.tex'
+input_name = project.comparereportdir + sep + 'comparison_onslow.tex'
 fid = open(input_name, 'w')
+
+s = '\\begin{figure}[hbt] \n \centerline{ \includegraphics[width=150mm, height=100mm]{../report_figures/%s}}' %gauge_map
+fid.write(s)
+
+s  = """
+\caption{Point locations used to compare MOST and ANUGA for Onslow region.}  
+\label{fig:comparisonpoints}
+\end{figure}
+"""
+fid.write(s)
+    
+s = '\input{%s} \n \clearpage \n \n' %latex_output[0]
+fid.write(s)
+
 
 # Generate latex output for location points
@@ -81 +102 @@
 #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)