Changeset 4294

Timestamp:

Mar 5, 2007, 4:13:51 PM (18 years ago)

Author:

sexton

Message:

slide report updates and gold coast hypothetical

Location:

anuga_work/production

Files:

• gold_coast_2007 (added)
• gold_coast_2007/export_results.py (added)
• gold_coast_2007/project.py (added)
• gold_coast_2007/project.pyc (added)
• gold_coast_2007/run_goldcoast.py (added)
• sydney_2006/export_results.py (modified) (4 diffs)
• sydney_2006/project_slide.py (modified) (2 diffs)
• sydney_2006/report/anuga_setup.tex (modified) (1 diff)
• sydney_2006/report/bulli_A_map.tex (modified) (1 diff)
• sydney_2006/report/bulli_B_map.tex (modified) (1 diff)
• sydney_2006/report/bulli_C_map.tex (modified) (1 diff)
• sydney_2006/report/execsum.tex (modified) (1 diff)
• sydney_2006/report/interpretation.tex (modified) (4 diffs)
• sydney_2006/report/shovel_A_map.tex (modified) (1 diff)
• sydney_2006/report/shovel_B_map.tex (modified) (1 diff)
• sydney_2006/report/shovel_C_map.tex (modified) (1 diff)
• sydney_2006/report/summary.tex (modified) (1 diff)
• sydney_2006/report/yacaaba_A_map.tex (modified) (1 diff)
• sydney_2006/report/yacaaba_B_map.tex (modified) (1 diff)
• sydney_2006/report/yacaaba_C_map.tex (modified) (1 diff)
• sydney_2006/smf_approx.py (modified) (2 diffs)

anuga_work/production/sydney_2006/export_results.py

@@ -11 +11 @@
 #time_dir = '20070121_213509' #shovel
 #time_dir = '20070119_055222' #yacaaba
-time_dir = '20070121_213444' #bulli
+#time_dir = '20070121_213444' #bulli
 
-#correct depths
-time_dir = '20070213_031322' #  shovel
-#time_dir = '20070213_033012' # bulli
-#time_dir = '20070213_033058' # yacaaba
-
-directory = project_slide.outputdir
-name = directory + time_dir + sep + project_slide.basename
+#directory = project_slide.outputdir
+#name = directory + time_dir + sep + project_slide.basename
+name = 'test500000.035000'
+#name = 'test'
 
 is_parallel = False
@@ -25 +22 @@
 print 'output dir:', name
 
-which_var = 2
+which_var = 0
 if which_var == 0:  # Stage
-    outname = name + '_stage'
+    outname = name + '_stage-eq-thurs-test'
     quantityname = 'stage'
 
@@ -35 +32 @@
 
 if which_var == 2:  # Depth
-    outname = name + '_depth-shovel'
+    outname = name + '_depth-eq'
     quantityname = 'stage-elevation'  
 
@@ -67 +64 @@
     sww2dem(name, basename_out = outname,
                 quantity = quantityname,
-                cellsize = 20,      
+                timestep = 0,
+                cellsize = 50,      
                 easting_min = project_slide.eastingmin,
                 easting_max = project_slide.eastingmax,

anuga_work/production/sydney_2006/project_slide.py

@@ -82 +82 @@
 # Domain definitions
 ###############################
-
+"""
 # bounding polygon for study area
 polyAll = read_polygon(polygondir+'extent3.csv')
@@ -133 +133 @@
 northingmin = 6244800
 northingmax = 6274550
+"""
+# exporting asc grid
+eastingmin = 318922+10000
+eastingmax = 318922+45000
+northingmin = 6231749+7500
+northingmax = 6231749+48000
+
+#eastingmin = -40
+#eastingmax = 40
+#northingmin = -20
+#northingmax = 20
 
 ###################################################################

anuga_work/production/sydney_2006/report/anuga_setup.tex

@@ -120 +120 @@
 $\backslash {\rm inundation \backslash data \backslash new\_south\_wales
 \backslash }$ \\ \hline
-A & Bulli &  $ {\rm wollongong \_tsunami\_scenario\_2006 \backslash anuga
-\backslash outputs \backslash 20061211\_060105}$ \\ \hline
-A & Shovel & ${\rm wollongong \_tsunami\_scenario\_2006 \backslash anuga
+A & Bulli - Historical &  $ {\rm wollongong \_tsunami\_scenario\_2006 \backslash anuga
+\backslash outputs \backslash 20070213\_035110}$ \\ \hline
+A & Shovel - Historical & ${\rm wollongong \_tsunami\_scenario\_2006 \backslash anuga
 \backslash outputs \backslash
-20061212\_012715}$ \\ \hline
+20070213\_035146}$ \\ \hline
 A & Yacaaba & ${\rm wollongong \_tsunami\_scenario\_2006 \backslash anuga
 \backslash outputs \backslash
-20061212\_064735}$ \\ \hline
+20070213\_035302}$ \\ \hline
 B & Bulli &  ${\rm sydney \_tsunami\_scenario\_2006 \backslash anuga
 \backslash outputs \backslash
-20061211\_071516}$ \\ \hline
+20070213\_033012}$ \\ \hline
 B & Shovel & ${\rm sydney \_tsunami\_scenario\_2006 \backslash anuga
 \backslash outputs \backslash
-20061212\_012705}$ \\ \hline
+20070213\_031322}$ \\ \hline
 B & Yacaaba & $\rm{ sydney \_tsunami\_scenario\_2006 \backslash anuga
 \backslash outputs \backslash
-20061212\_064807}$ \\ \hline
+20070213\_033058}$ \\ \hline
 C & Bulli &  ${\rm newcastle \_tsunami\_scenario\_2006 \backslash anuga
 \backslash outputs \backslash
-20061211\_073709}$ \\ \hline
+20070219\_060722}$ \\ \hline
 C & Shovel &  ${\rm newcastle \_tsunami\_scenario\_2006 \backslash anuga
 \backslash outputs \backslash
-20061212\_012802}$ \\ \hline
-C & Yacaaba & ${\rm newcastle \_tsunami\_scenario\_2006 \backslash anuga
+20070219\_060808}$ \\ \hline
+C & Yacaaba - Historical & ${\rm newcastle \_tsunami\_scenario\_2006 \backslash anuga
 \backslash outputs \backslash
-20061212\_064757}$ \\ \hline
+20070219\_060854}$ \\ \hline
 \end{tabular}
 

anuga_work/production/sydney_2006/report/bulli_A_map.tex

@@ -1 +1 @@
 %\begin{sidewaysfigure}
 
+%\begin{figure}
+%\centering
+%\begin{tabular}{cc}
+%\includegraphics[scale=0.4]{../report_figures/bulliA.jpg} & 
+%\includegraphics[scale=0.4]{../report_figures/bulliAimagery.jpg} \\
+%\end{tabular}
 \begin{figure}
 \centering
-\begin{tabular}{cc}
-\includegraphics[scale=0.4]{../report_figures/bulliA.jpg} & 
-\includegraphics[scale=0.4]{../report_figures/bulliAimagery.jpg} \\
-\end{tabular}
+\includegraphics[scale=0.5]{..../report_figures/bulliAimageryreviseddepth.jpg}
 \caption{Maximum inundation map for the Bulli slide in region A. Data: NSW Department of Lands and AHO.}  
 \label{fig:bulliA}

anuga_work/production/sydney_2006/report/bulli_B_map.tex

@@ -3 +3 @@
 \begin{figure}
 \centering
-\begin{tabular}{cc}
-\includegraphics[scale=0.4]{../report_figures/bulliB.jpg} & 
-\includegraphics[scale=0.4]{../report_figures/bulliBimagery.jpg} \\
-\end{tabular}
+\includegraphics[scale=0.4]{../report_figures/bulliBimageryreviseddepth.jpg}
 \caption{Maximum inundation map for the Bulli slide in region B. Data: NSW Department of Lands and AHO.}  
 \label{fig:bulliB}

anuga_work/production/sydney_2006/report/bulli_C_map.tex

@@ -3 +3 @@
 \begin{figure}
 \centering
-\begin{tabular}{cc}
-\includegraphics[scale=0.4]{../report_figures/bulliC.jpg} &
-\includegraphics[scale=0.4]{../report_figures/bulliCimagery.jpg} \\
-\end{tabular}
+\includegraphics[scale=0.4]{../report_figures/bulliCimageryreviseddepth.jpg} 
 \caption{Maximum inundation map for the Bulli slide in region C. Data: NSW Department of Lands and AHO.}  
 \label{fig:bulliC}

anuga_work/production/sydney_2006/report/execsum.tex

@@ -10 +10 @@
 for submarine mass failures based on results from the October 2006 Continental
 Slope Survey. In particular, three slides are modelled in three regions; 
-Wollongong, Sydney and Newcastle. One of these slides represents a historical event
-and the remainder represent potential slide scenarios.
+Wollongong, Sydney and Newcastle. Three of the nine slides modelled
+represent historical events and the remainder represent potential slide scenarios.
 
 

anuga_work/production/sydney_2006/report/interpretation.tex

@@ -1 +1 @@
 
 The modelling has considered three separate regions based on the location
-of the potential slide failures. The regions A, B and C refer
+of the historical and potential slide failures. The regions A, B and C refer
 to Wollongong, Sydney and Newcastle. Figure \ref{fig:overallmodel} shows
 each model with respect to each other and highlights the location
@@ -8 +8 @@
 
 \begin{figure}[h]
-\centerline{\includegraphics[scale=0.4]{../report_figures/overallmodel.jpg}}
+\centerline{\includegraphics[scale=0.4]{../report_figures/overallmodelv2.jpg}}
 \caption{NSW region showing the location of each study region and the location
-of the origin of the slide.}
+of the origin of the historical and potential slide.}
 \label{fig:overallmodel}
 \end{figure}
 
 The purpose of this section is show the maximum inundation map for each
-slide scenario for region A, B and C. The Bulli slide is considered to
-have occurred with the Shovel and Yacaaba considered as potential slides.
+slide scenario for region A, B and C. The Bulli and Shovel slide occurring
+in region A and the Yacaaba slide in region C are considered to have occurred.
+The potential slides are based on the volumes of these slides located in 
+appropriate locations in the other regions.
+The outputs in region B should thus be interpretted as potential tsunami inundation
+for slides of varied volume. 
 The modelling effectively describes the output as a function of slide
-volume with Bulli the largest. The outputs should thus be interpretted
-as potential outputs with the Bulli volume occurring in region A as the 
-only to have occurred. This event may be used in paleo-tsunami 
-investigations due to occur in early 2007. The remaining events may
-then be used to described the potential tsunami inundation for slides 
-of varied volume.
+volume with Bulli the largest. This event may be used in paleo-tsunami 
+investigations due to occur in early 2007. 
 
 The time to impact after the initiation of the slide is approximately 25 mins 
@@ -40 +40 @@
 Thickness (m) & 424  & 165  & 144  \\ \hline
 Density g m$^{-3}$ & 1.46  & 1.49  & 1.48  \\ \hline
-Water depth to centre of mass (m) & 2087  & 968  & 1119  \\ \hline
+Water depth to centre of mass (m) & 1470  & 877  & 1223  \\ \hline
 Bathymetric slope (\degree) & 4  & 4  & 3  \\ \hline
 Angular orientation (\degree) & 126 & 118  & 133  \\ \hline
@@ -48 +48 @@
 \end{table}
 
-
-
 Figures \ref{fig:bulliA}, \ref{fig:bulliB} and \ref{fig:bulliC}
 describe the inundation extent to regions A, B and C respectively
-for the Bulli volume. Similarly, Figures \ref{fig:shovelA}, \ref{fig:shovelB}
+for the Bulli volume, noting that \ref{fig:bulliA} is the historical event.
+Similarly, Figures \ref{fig:shovelA}, \ref{fig:shovelB}
 and \ref{fig:shovelC} describe the inundation extent for the Shovel
-volume with Figures \ref{fig:yacaabaA}, \ref{fig:yacaabaB} and 
-\ref{fig:yacaabaC} describing the inundation extent for the Yacaaba volume.
-It is immediately apparent that there is little difference in the
-inundation extents between the Bulli and Shovel slide volumes.  The reason
-for this may possibly be explained in Figure \ref{fig:depthvsamp}. Here, the 
+volume, noting that \ref{fig:shovelA} is the historical event. 
+Figures \ref{fig:yacaabaA}, \ref{fig:yacaabaB} and 
+\ref{fig:yacaabaC} describe the inundation extent for the Yacaaba volume, noting
+that \ref{fig:yacaabaC} is the historical event.
+It appears that the inundation extent grows with the volume of the slide. This may
+be explained through Figure \ref{fig:depthvsamp} where the 
 characteristic three dimensional 
-amplitude is determined as a function of depth for each
-slide volume (see Appendix \ref{sec:smfmodel} for details) 
-with the location of each specific slide marked on its
-corresponding volume. 
-As Figure \ref{fig:depthvsamp} shows, even though the Bulli and Shovel
-slides are considered to occur in different depths of water, the
-resultant amplitude are approximately equal.
+amplitude is shown as a function of depth for each
+slide volume (see Appendix \ref{sec:smfmodel} for details). The depth
+for each specific slide is marked on the curve corresponding to its volume. From this we
+can see that the initial amplitude is greatest for the Bulli slide, followed by the Shove
+and Yacaaba.
 
 \begin{figure}[h]

anuga_work/production/sydney_2006/report/shovel_A_map.tex

@@ -3 +3 @@
 \begin{figure}
 \centering
-\begin{tabular}{cc}
-\includegraphics[scale=0.4]{../report_figures/shovelA.jpg} &
-\includegraphics[scale=0.4]{../report_figures/shovelAimagery.jpg} \\
-\end{tabular}
+\includegraphics[scale=0.4]{../report_figures/shovelAimageryreviseddepth.jpg}
 \caption{Maximum inundation map for the Shovel slide in region A. Data: NSW Department of Lands and AHO.}  
 \label{fig:shovelA}

anuga_work/production/sydney_2006/report/shovel_B_map.tex

@@ -3 +3 @@
 \begin{figure}
 \centering
-\begin{tabular}{cc}
-\includegraphics[scale=0.4]{../report_figures/shovelB.jpg} &
-\includegraphics[scale=0.4]{../report_figures/shovelBimagery.jpg} \\
-\end{tabular}
+\includegraphics[scale=0.4]{../report_figures/shovelBimageryreviseddepth.jpg}
 \caption{Maximum inundation map for the Shovel slide in region B. Data: NSW Department of Lands and AHO.}  
 \label{fig:shovelB}

anuga_work/production/sydney_2006/report/shovel_C_map.tex

@@ -3 +3 @@
 \begin{figure}
 \centering
-\begin{tabular}{cc}
-\includegraphics[scale=0.4]{../report_figures/shovelC.jpg} &
-\includegraphics[scale=0.4]{../report_figures/shovelCimagery.jpg} \\
-\end{tabular}
+\includegraphics[scale=0.4]{../report_figures/shovelCimageryreviseddepth.jpg}
 \caption{Maximum inundation map for the Shovel slide in region C. Data: NSW Department of Lands and AHO.}  
 \label{fig:shovelC}

anuga_work/production/sydney_2006/report/summary.tex

@@ -1 +1 @@
 This report has described the tsunami inundation occurring at Wollongong,
 Sydney and Newcastle for a range of different slides 
-occurring at Mean Sea Level. 
+occurring at Mean Sea Level. Three of the nine modelled slides are 
+historical; Bulli, Shovel at Wollongong and Yacaaba at Newcastle. The 
+remaining slides are considered potential scenarios.
 The results described must be tempered with the modelling assumptions 
 and the accuracy and coverage of the underlying elevation data.

anuga_work/production/sydney_2006/report/yacaaba_A_map.tex

@@ -3 +3 @@
 \begin{figure}
 \centering
-\begin{tabular}{cc}
-\includegraphics[scale=0.4]{../report_figures/yacaabaA.jpg} &
-\includegraphics[scale=0.4]{../report_figures/yacaabaAimagery.jpg} \\
-\end{tabular}
+\includegraphics[scale=0.4]{../report_figures/yacaabaAimageryreviseddepth.jpg} 
 \caption{Maximum inundation map for the Yacaaba slide in region A. Data: NSW Department of Lands and AHO.}  
 \label{fig:yacaabaA}

anuga_work/production/sydney_2006/report/yacaaba_B_map.tex

@@ -3 +3 @@
 \begin{figure}
 \centering
-\begin{tabular}{cc}
-\includegraphics[scale=0.4]{../report_figures/yacaabaB.jpg} &
-\includegraphics[scale=0.4]{../report_figures/yacaabaBimagery.jpg} \\
-\end{tabular}
+\includegraphics[scale=0.4]{../report_figures/yacaabaBimageryreviseddepth.jpg} 
 \caption{Maximum inundation map for the Yacaaba slide in region B. Data: NSW Department of Lands and AHO.}  
 \label{fig:yacaabaB}

anuga_work/production/sydney_2006/report/yacaaba_C_map.tex

@@ -3 +3 @@
 \begin{figure}
 \centering
-\begin{tabular}{cc}
-\includegraphics[scale=0.4]{../report_figures/yacaabaC.jpg} &
-\includegraphics[scale=0.4]{../report_figures/yacaabaCimagery.jpg} \\
-\end{tabular}
+\includegraphics[scale=0.4]{../report_figures/yacaabaCimageryreviseddepth.jpg} 
 \caption{Maximum inundation map for the Yacaaba slide in region C. Data: NSW Department of Lands and AHO.}  
 \label{fig:yacaabaC}

anuga_work/production/sydney_2006/smf_approx.py

@@ -88 +88 @@
 #par = 'depth'
 #par = 'density'
-par = 'slope'
+par = 'depth'
 
 depths = range(750,2500,50)
@@ -100 +100 @@
 gammas = [1.46,1.49,1.48]
 smftype = ['Bulli','Shovel','Yacaaba']
-smfdepths = [2087,968,1119]
+#smfdepths = [2087,968,1119]
+smfdepths = [1470,877,1222.5]
 
 if par == 'depth':