Changeset 4151 for anuga_work/production/onslow_2006

Timestamp:

Jan 9, 2007, 3:26:53 PM (18 years ago)

Author:

sexton

Message:

(1) updates to make_report script to handle multiple events for multiple return periods (2) take Onslow changes through to Port Hedland and Dampier

Location:

anuga_work/production/onslow_2006

Files:

• make_report.py (modified) (7 diffs)
• project.py (modified) (1 diff)
• report/anuga_setup.tex (modified) (3 diffs)
• report/metadata.tex (modified) (1 diff)
• report/modelling_methodology.tex (modified) (1 diff)

anuga_work/production/onslow_2006/make_report.py

@@ -25 +25 @@
 * Tsunami scenario
 * Data sources
-
 * Inundation modelling results
 * Impact modelling
@@ -81 +80 @@
 
 # WA DLI data
-production_dirs = {'20060704_063234': 'MSL',
-                   '20060704_063005': 'HAT',
-                   '20060706_235246': 'LAT'}
-
-max_maps = {'MSL': 'MSL_map',
+production_dirs_test1 = {'20060704_063234': 'MSL',
+                         '20060704_063005': 'HAT',
+                         '20060706_235246': 'LAT'}
+
+production_dirs_test2 = production_dirs_test1
+
+return_periods = [1000, 1000000]
+
+all_prod_dirs = [production_dirs_test1, production_dirs_test2]
+#all_prod_dirs = [production_dirs_test1]
+
+is_parallel = False
+if is_parallel == True:
+    nodes = 8
+    
+max_maps1 = {'MSL': 'MSL_map',
             'HAT': 'HAT_map',
             'LAT': 'LAT_map'}
 
+max_maps2 = max_maps1
+
+all_maps = [max_maps1, max_maps2]
+
 gauge_map = 'onslow_dli_gauge.jpg'
 
@@ -95 +109 @@
 # Create sections and graphs for each designated production directory
 latex_output = []
+report_name = 'latexoutput'
 swwfiles = {}
-for label_id in production_dirs.keys():
-    
-    file_loc = project.outputdir + label_id + sep
-    swwfile = file_loc + project.basename + '.sww'
-    swwfiles[swwfile] = label_id
-
-texname, elev_output = sww2timeseries(swwfiles,
-                                      project.gauge_filename,
-                                      #project.gauge_filename_bindi,
-                                      production_dirs,
-                                      report = True,
-                                      reportname = 'latexoutput',
-                                      plot_quantity = ['stage', 'speed'],
-                                      surface = False,
-                                      time_min = None,
-                                      time_max = None,
-                                      title_on = False,
-                                      verbose = True)
-
-latex_output.append(texname)
-
+if is_parallel == True:
+
+    for j, production_dirs in enumerate(all_prod_dirs):
+        
+        for i in range(nodes):
+            print 'Sending node %d of %d' %(i,nodes)
+            swwfiles = {}
+            reportname = report_name + '%s' + 'rp%s' %(i,j)
+            
+            for label_id in production_dirs.keys():
+                file_loc = project.output_dir + label_id + sep
+                sww_extra = '_P%s_%s' %(i,nodes)
+                swwfile = file_loc + project.basename + sww_extra + '.sww'
+                #swwfile = file_loc + project.scenario_name + sww_extra + '.sww'
+                swwfiles[swwfile] = label_id
+
+                texname, elev_output = sww2timeseries(swwfiles,
+                                                      project.gauge_filename,
+                                                      production_dirs,
+                                                      report = True,
+                                                      reportname = reportname,
+                                                      plot_quantity = ['stage', 'momentum'],
+                                                      surface = False,
+                                                      time_min = None,
+                                                      time_max = None,
+                                                      title_on = False,
+                                                      verbose = True)
+                
+                latex_output.append(texname)
+    
+else:
+
+    for i, production_dirs in enumerate(all_prod_dirs):
+
+        swwfiles = {}
+        reportname = report_name + 'rp%s' %(i)
+        
+        for label_id in production_dirs.keys():
+
+            file_loc = project.outputdir + label_id + sep
+            swwfile = file_loc + project.basename + '.sww'
+            #swwfile = file_loc + project.scenario_name + sww_extra + '.sww'
+            swwfiles[swwfile] = label_id
+            
+        texname, elev_output = sww2timeseries(swwfiles,
+                                              project.gauge_filename,
+                                              production_dirs,
+                                              report = True,
+                                              reportname = reportname,
+                                              plot_quantity = ['stage', 'speed'],
+                                              surface = False,
+                                              time_min = None,
+                                              time_max = None,
+                                              title_on = False,
+                                              verbose = True)
+        
+        latex_output.append(texname)
+    
 # Start report generation
 # Future: generate_report(reportdir, scenario, report_title, 
@@ -152 +204 @@
 \usepackage{setspace}
 \usepackage{rotating}
+\usepackage{pdfpages}
 \include{appendix}
 \setstretch{1.25}
@@ -164 +217 @@
 fid.write(s)
 
-#s = '\\title{%s} \n' %report_title
-#fid.write(s)
-
 s = """
 \date{\\today}
@@ -218 +268 @@
 # Assign titles to each production section 
 # Must specify one name per section
-for i, name in enumerate(production_dirs.keys()):
-#
-#    s = '\subsection{%s} \n \n' %production_dirs[name]      
-#    fid.write(s)
-
-    s = '\input{%s} \n \clearpage \n \n' %max_maps[production_dirs[name]]
+for i, max_maps in enumerate(all_maps):
+
+    s = '\subsection{Return Period: %s years} \n \n' %return_periods[i]      
     fid.write(s)
+
+    production_dirs = all_prod_dirs[i]
+    for i, name in enumerate(production_dirs.keys()):
+
+        s = '\input{%s} \n \clearpage \n \n' %max_maps[production_dirs[name]]
+        fid.write(s)
 
 # Generate latex output for location points
@@ -308 +361 @@
 fid.write(s)
 
-s = '\input{%s} \n \clearpage \n \n' %latex_output[0]
-fid.write(s)
+for i in range(len(latex_output)):
+    if latex_output[i] <> '':
+        s = '\subsection{Return Period: %s years} \n' %return_periods[i]
+        fid.write(s)
+        s = '\input{%s} \n \clearpage \n \n' %latex_output[i]   
+        fid.write(s)
 
 s="""

anuga_work/production/onslow_2006/project.py

@@ -17 +17 @@
 state = 'western_australia'
 scenario_dir_name = 'onslow_tsunami_scenario_2006'
+scenario_name = 'onslow'
 
 # 250m data to be provided

anuga_work/production/onslow_2006/report/anuga_setup.tex

@@ -12 +12 @@
  
 The initial conditions used for this scenario are MSL, HAT and LAT which were
-defined in Section \ref{sec:data}.
+defined in Section \ref{sec:data}. Figure \ref{fig:IC} shows the MSL, HAT and 
+LAT contours to illustrate the water level for each of these inial conditions.
 The dynamics of
 tidal effects (that is, the changes in water height over time for
@@ -23 +24 @@
 in the scenario. The 
 results are therefore likely to be over estimates.
+
+\begin{figure}[h]
+
+  \centerline{\includegraphics[width=\paperwidth]{../report_figures/onslow_dli_contour.jpg}}
+
+  \caption{Onslow region showing the initial conditions used for the study; 
+  -1.5m AHD (LAT), 0m AHD (MSL) and 1.5m AHD (HAT) contour lines.}
+  \label{fig:IC}
+\end{figure}
+
 
 To set up a model for the tsunami scenario, a study area is first 
@@ -69 +80 @@
 The final item to be addressed to complete the model setup is the 
 definition of the boundary condition. As 
-discussed in Section \ref{sec:methodology}, a range of events with
-different return periods will modelled. It must be 
-remembered that tsunamis are typically made up of a series
+discussed in Section \ref{sec:methodology}, a series of events corresponding
+to a range of return periods are selected as the tsunami sources.
+The resultant tsunami wave is made up of a series
 of waves with different amplitudes which is affected by the energy
 and style of the event as well as the bathymetry whilst it travels
-from its source to Onslow. The amplitude and velocity of each of these
+from its source to Dampier. The amplitude and velocity of each of these
 waves are then provided to ANUGA as boundary conditions and propagated
 inshore.

anuga_work/production/onslow_2006/report/metadata.tex

@@ -1 @@
-
-See Enclosure.
+%\includepdf[pages={1-6}]{MetadataforATWSOnslowScenario}

anuga_work/production/onslow_2006/report/modelling_methodology.tex

@@ -1 @@
-
 Geoscience Australia aims to define the economic and social threat posed to urban communities 
 by a range of rapid onset natural hazards. Through the integration of natural hazard research, defining national exposure and