Changeset 4145
- Timestamp:
- Jan 8, 2007, 3:36:39 PM (17 years ago)
- Files:
-
- 39 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_core/documentation/planning/Inundation Scenarios - 06 07 Plan.Gnt
r4071 r4145 1 1 <?xml version="1.0" encoding="utf-8"?> 2 2 <!--****WARNING: DO NOT MODIFY THIS FILE DIRECTLY******--> 3 <Gantt Version="0.902" Checksum=" +R8DXMfKYgs=">3 <Gantt Version="0.902" Checksum="FkpD2Gk1zVI="> 4 4 <Display Width="0" Height="0" Split="178" /> 5 <Miscellaneous Saved="2006-12-1 1 18:35:09Z" Created="Monday, 11 December 2006 7:35:09" />5 <Miscellaneous Saved="2006-12-19 10:48:35Z" Created="Monday, 18 December 2006 11:48:35 " /> 6 6 <Globalization> 7 7 <Culture LCID="3081" ISO="ENA" DN="English (Australia)" Cal="System.Globalization.GregorianCalendar" Parent="English">English (Australia)</Culture> … … 10 10 <FirstDay Date="2006-08-28" /> 11 11 <LastDay Date="2007-06-30" /> 12 <Tasks Size=" 38">12 <Tasks Size="40"> 13 13 <!--S=Start D=Duration P=Predecessor H=Height L=Lag VA=Valign R=Remarks F=Flag I=Importance--> 14 14 <!--B=BarType BC=BarColor BTx=BarText--> … … 17 17 <Task ID="3" Type="1" S="2006-09-20" SH="0" P="2" In="0" B="2" BC="-16776961" D="17" H="35" L="1" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Hobart scenario</Task> 18 18 <Task ID="4" Type="1" S="2006-09-29" SH="12" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of Karratha/Dampier data to NM&I</Task> 19 <Task ID="5" Type="1" S="2006-10-12" SH="12" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of K/D data </Task>19 <Task ID="5" Type="1" S="2006-10-12" SH="12" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of K/D data for CIPMA scenario</Task> 20 20 <Task ID="6" Type="1" S="2006-09-15" SH="12" P="0" In="0" B="2" BC="-16776961" D="50.5" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">K/D scenario</Task> 21 21 <Task ID="7" Type="1" S="2006-10-06" SH="12" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of Broome offshore data to NM&I</Task> 22 <Task ID="8" Type="1" S="2006-1 0-26" SH="12" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of Broome data</Task>23 <Task ID="9" Type="1" S="200 6-10-27" SH="12" P="8" In="0" B="2" BC="-16776961" D="15.0" H="35" L="1" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Broome setup</Task>22 <Task ID="8" Type="1" S="2006-12-15" SH="12" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of Broome data</Task> 23 <Task ID="9" Type="1" S="2007-01-08" SH="0" P="8" In="0" B="2" BC="-16776961" D="3.0" H="35" L="15.5" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Broome setup</Task> 24 24 <Task ID="10" Type="1" S="2006-11-24" SH="12" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of Perth offshore data to NM&I</Task> 25 25 <Task ID="11" Type="1" S="2006-12-07" SH="12" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of Perth data</Task> 26 <Task ID="12" Type="1" S="2006-12-08" SH="0" P="11" In="0" B="2" BC="-16776961" D=" 3" H="35" L="0.5" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Perth setup</Task>26 <Task ID="12" Type="1" S="2006-12-08" SH="0" P="11" In="0" B="2" BC="-16776961" D="7" H="35" L="0.5" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Perth setup</Task> 27 27 <Task ID="13" Type="1" S="2006-11-24" SH="0" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of Busselton offshore data to NM&I</Task> 28 <Task ID="14" Type="1" S="2006-12-1 4" SH="12" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of Busselton data</Task>29 <Task ID="15" Type="1" S="2006-12-1 8" SH="0" P="14" In="0" B="2" BC="-16776961" D="2.5" H="35" L="1.5" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Busselton setup</Task>28 <Task ID="14" Type="1" S="2006-12-18" SH="0" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of Busselton data</Task> 29 <Task ID="15" Type="1" S="2006-12-19" SH="12" P="14" In="0" B="2" BC="-16776961" D="3.5" H="35" L="1.5" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Busselton setup</Task> 30 30 <Task ID="16" Type="1" S="2006-11-27" SH="0" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of Exmouth offshore data to NM&I</Task> 31 <Task ID="17" Type="1" S="2006-12-11" SH="0" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of Exmouth data</Task> 32 <Task ID="18" Type="1" S="2006-12-12" SH="12" P="17" In="0" B="2" BC="-16776961" D="3.5" H="35" L="1.5" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Exmouth setup</Task> 33 <Task ID="19" Type="1" S="2007-01-12" SH="12" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of Shark Bay offshore data to NM&I</Task> 34 <Task ID="20" Type="1" S="2007-02-02" SH="0" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of Shark Bay data</Task> 35 <Task ID="21" Type="1" S="2007-01-24" SH="12" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of BC for Shark Bay</Task> 36 <Task ID="22" Type="1" S="2007-02-05" SH="0" P="20" In="0" B="2" BC="-16776961" D="15" H="35" L="1" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Shark Bay validation exercise</Task> 37 <Task ID="23" Type="1" S="2007-01-17" SH="12" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Decision of number and type of events</Task> 38 <Task ID="24" Type="1" S="2007-01-26" SH="0" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of events</Task> 39 <Task ID="25" Type="1" S="2007-01-29" SH="12" P="0" In="0" B="2" BC="-16776961" D="15" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">K/D run</Task> 40 <Task ID="26" Type="1" S="2007-02-05" SH="0" P="0" In="0" B="2" BC="-16776961" D="15" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Broome run</Task> 41 <Task ID="27" Type="1" S="2007-02-12" SH="0" P="0" In="0" B="2" BC="-16776961" D="15" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Perth run</Task> 42 <Task ID="28" Type="1" S="2007-02-19" SH="0" P="0" In="0" B="2" BC="-16776961" D="15" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Busselton run</Task> 43 <Task ID="29" Type="1" S="2007-02-26" SH="0" P="0" In="0" B="2" BC="-16776961" D="14.5" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Exmouth run</Task> 44 <Task ID="30" Type="1" S="2007-03-05" SH="0" P="0" In="0" B="2" BC="-16776961" D="14.5" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Onslow run</Task> 45 <Task ID="31" Type="1" S="2007-03-12" SH="0" P="0" In="0" B="2" BC="-16776961" D="15" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Pt Hedland run</Task> 46 <Task ID="32" Type="1" S="2006-12-05" SH="12" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of offshore data for NSW slump scenarios</Task> 47 <Task ID="33" Type="1" S="2006-12-07" SH="0" P="0" In="0" B="2" BC="-16776961" D="10" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Wollongong slump scenario</Task> 48 <Task ID="34" Type="1" S="2006-12-07" SH="0" P="0" In="0" B="2" BC="-16776961" D="10" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Sydney slump scenario</Task> 49 <Task ID="35" Type="1" S="2006-12-07" SH="0" P="0" In="0" B="2" BC="-16776961" D="10" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Newcastle slump scenario</Task> 50 <Task ID="36" Type="1" S="2006-12-15" SH="0" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">ANUGA release</Task> 51 <Task ID="37" Type="1" S="2007-01-19" SH="0" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Beowulf operational</Task> 52 <Task ID="38" Type="1" S="2006-12-22" SH="0" P="0" In="0" B="2" BC="-16776961" D="11" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Christmas Break</Task> 31 <Task ID="17" Type="1" S="2006-12-21" SH="0" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of Exmouth data</Task> 32 <Task ID="18" Type="1" S="2007-01-11" SH="0" P="17" In="0" B="2" BC="-16776961" D="3" H="35" L="15" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Exmouth setup</Task> 33 <Task ID="19" Type="1" S="2006-12-22" SH="0" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Deliver of updated K/D data</Task> 34 <Task ID="20" Type="1" S="2007-01-03" SH="12" P="19" In="0" B="2" BC="-16776961" D="2.5" H="35" L="8.5" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">K/D resetup</Task> 35 <Task ID="21" Type="1" S="2007-01-12" SH="12" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of Shark Bay offshore data to NM&I</Task> 36 <Task ID="22" Type="1" S="2007-02-02" SH="0" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of Shark Bay data</Task> 37 <Task ID="23" Type="1" S="2007-01-24" SH="12" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of BC for Shark Bay</Task> 38 <Task ID="24" Type="1" S="2007-02-05" SH="0" P="22" In="0" B="2" BC="-16776961" D="15" H="35" L="1" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Shark Bay validation exercise</Task> 39 <Task ID="25" Type="1" S="2007-01-17" SH="12" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Decision of number and type of events</Task> 40 <Task ID="26" Type="1" S="2007-01-26" SH="0" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of events</Task> 41 <Task ID="27" Type="1" S="2007-01-29" SH="12" P="0" In="0" B="2" BC="-16776961" D="15" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">K/D run</Task> 42 <Task ID="28" Type="1" S="2007-02-05" SH="0" P="0" In="0" B="2" BC="-16776961" D="15" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Broome run</Task> 43 <Task ID="29" Type="1" S="2007-02-12" SH="0" P="0" In="0" B="2" BC="-16776961" D="15" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Perth run</Task> 44 <Task ID="30" Type="1" S="2007-02-19" SH="0" P="0" In="0" B="2" BC="-16776961" D="15" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Busselton run</Task> 45 <Task ID="31" Type="1" S="2007-02-26" SH="0" P="0" In="0" B="2" BC="-16776961" D="14.5" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Exmouth run</Task> 46 <Task ID="32" Type="1" S="2007-03-05" SH="0" P="0" In="0" B="2" BC="-16776961" D="14.5" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Onslow run</Task> 47 <Task ID="33" Type="1" S="2007-03-12" SH="0" P="0" In="0" B="2" BC="-16776961" D="15" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Pt Hedland run</Task> 48 <Task ID="34" Type="1" S="2006-12-05" SH="12" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Delivery of offshore data for NSW slump scenarios</Task> 49 <Task ID="35" Type="1" S="2006-12-07" SH="0" P="0" In="0" B="2" BC="-16776961" D="10" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Wollongong slump scenario</Task> 50 <Task ID="36" Type="1" S="2006-12-07" SH="0" P="0" In="0" B="2" BC="-16776961" D="10" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Sydney slump scenario</Task> 51 <Task ID="37" Type="1" S="2006-12-07" SH="0" P="0" In="0" B="2" BC="-16776961" D="10" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Newcastle slump scenario</Task> 52 <Task ID="38" Type="1" S="2006-12-22" SH="0" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">ANUGA release</Task> 53 <Task ID="39" Type="1" S="2007-01-19" SH="0" P="0" In="0" B="0" BC="-16776961" D="0" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Beowulf operational</Task> 54 <Task ID="40" Type="1" S="2006-12-22" SH="0" P="0" In="0" B="2" BC="-16776961" D="11" H="35" L="0" I="F" F="F" VA="1" R="" BTL="0" BTT="0" BTR="0" BTB="0">Christmas Break</Task> 53 55 </Tasks> 54 56 <TextStyles Size="8"> -
anuga_core/source/anuga/abstract_2d_finite_volumes/util.py
r4131 r4145 833 833 834 834 import sys 835 if not sys.version.startswith('2.4'):835 #if not sys.version.startswith('2.4'): 836 836 # Backwards compatibility 837 838 839 840 if sys.platform == 'win32':841 842 843 844 845 846 837 # from sets import Set as set 838 839 from sets import Set as set 840 841 for i,j in enumerate(quantity): 842 quantity[i] = quantity[i].lower() 843 p = list(set(quantity).difference(set(all_quantity))) 844 if len(p) <> 0: 845 msg = 'Quantities %s do not exist - please try again' %p 846 raise Exception, msg 847 847 848 848 return … … 928 928 max_momentums = [] 929 929 max_speeds = [] 930 c = 0930 max_depths = [] 931 931 model_time_plot3d = zeros((n0,m), Float) 932 932 stages_plot3d = zeros((n0,m), Float) … … 944 944 max_stage = 0 945 945 max_momentum = 0 946 max_speed = 0 946 max_speed = 0 947 max_depth = 0 947 948 gaugeloc = locations[k] 948 949 thisfile = file_loc[j]+sep+'gauges_time_series'+'_'+gaugeloc+'.csv' … … 983 984 if m > max_momentum: max_momentum = m 984 985 if vel > max_speed: max_speed = vel 986 if z > 0 and depth > max_depth: max_depth = depth 985 987 986 988 … … 990 992 min_stages.append(min_stage) 991 993 max_momentums.append(max_momentum) 992 max_speeds.append(max_speed) 994 max_speeds.append(max_speed) 995 max_depths.append(max_depth) 993 996 #### finished generating quantities for each swwfile ##### 994 997 … … 1026 1029 profilefig = 'solution_xprofile' 1027 1030 savefig('profilefig') 1028 1029 #stage_axis = axis([time_min/60.0, time_max/60.0, min(min_stages), max(max_stages)*1.1]) 1030 #stage_axis = axis([time_min/60.0, time_max/60.0, -3.0, 3.0]) 1031 #vel_axis = axis([time_min/60.0, time_max/60.0, min(max_speeds), max(max_speeds)*1.1]) 1032 #mom_axis = axis([time_min/60.0, time_max/60.0, min(max_momentums), max(max_momentums)*1.1]) 1031 1032 depth_axis = axis([time_min/60.0, time_max/60.0, 0, max(max_depths)*1.1]) 1033 stage_axis = axis([time_min/60.0, time_max/60.0, min(min_stages), max(max_stages)*1.1]) 1034 stage_axis = axis([time_min/60.0, time_max/60.0, -3.0, 3.0]) 1035 vel_axis = axis([time_min/60.0, time_max/60.0, min(max_speeds), max(max_speeds)*1.1]) 1036 mom_axis = axis([time_min/60.0, time_max/60.0, min(max_momentums), max(max_momentums)*1.1]) 1033 1037 1034 1038 cstr = ['g', 'r', 'b', 'c', 'm', 'y', 'k'] … … 1037 1041 elev_output = [] 1038 1042 if len(label_id) > 1: graphname_report = [] 1043 pp = 1 1044 div = 11. 1039 1045 for k in gauge_index: 1040 1046 g = gauges[k] 1041 1047 count1 = 0 1042 1048 if report == True and len(label_id) > 1: 1043 s = '\\begin{figure}[h bt] \n \\centering \n \\begin{tabular}{cc} \n'1049 s = '\\begin{figure}[ht] \n \\centering \n \\begin{tabular}{cc} \n' 1044 1050 fid.write(s) 1045 1051 if len(label_id) > 1: graphname_report = [] … … 1063 1069 plot(model_time[0:n[j]-1,k,j], depths[0:n[j]-1,k,j], '-', c = cstr[j]) 1064 1070 units = 'm' 1071 axis(depth_axis) 1065 1072 if which_quantity == 'stage': 1066 if elevations[0 :n[j]-1,k,j] < 0:1073 if elevations[0,k,j] < 0: 1067 1074 plot(model_time[0:n[j]-1,k,j], stages[0:n[j]-1,k,j], '-', c = cstr[j]) 1075 axis(stage_axis) 1068 1076 else: 1069 1077 plot(model_time[0:n[j]-1,k,j], depths[0:n[j]-1,k,j], '-', c = cstr[j]) 1070 #axis(stage_axis)1078 axis(depth_axis) 1071 1079 units = 'm' 1072 1080 if which_quantity == 'momentum': 1073 1081 plot(model_time[0:n[j]-1,k,j], momenta[0:n[j]-1,k,j], '-', c = cstr[j]) 1074 #axis(mom_axis)1082 axis(mom_axis) 1075 1083 units = 'm^2 / sec' 1076 1084 if which_quantity == 'xmomentum': 1077 1085 plot(model_time[0:n[j]-1,k,j], xmom[0:n[j]-1,k,j], '-', c = cstr[j]) 1078 #axis(mom_axis)1086 axis(mom_axis) 1079 1087 units = 'm^2 / sec' 1080 1088 if which_quantity == 'ymomentum': 1081 1089 plot(model_time[0:n[j]-1,k,j], ymom[0:n[j]-1,k,j], '-', c = cstr[j]) 1082 #axis(mom_axis)1090 axis(mom_axis) 1083 1091 units = 'm^2 / sec' 1084 1092 if which_quantity == 'speed': 1085 1093 plot(model_time[0:n[j]-1,k,j], speed[0:n[j]-1,k,j], '-', c = cstr[j]) 1086 #axis(vel_axis)1094 axis(vel_axis) 1087 1095 units = 'm / sec' 1088 1096 if which_quantity == 'bearing': … … 1170 1178 s = '\end{tabular} \n \\caption{%s} \n \label{fig:%s} \n \end{figure} \n \n' %(caption, label) 1171 1179 fid.write(s) 1172 c += 11173 if c % 6 == 0: fid.write('\\clearpage \n')1180 cc += 1 1181 if cc % 6 == 0: fid.write('\\clearpage \n') 1174 1182 savefig(graphname_latex) 1175 1183 … … 1186 1194 else: 1187 1195 word_quantity += plot_quantity[i] 1188 print 'hello', elevations[0,k,j]1189 1196 where1 = 0 1190 1197 count1 += 1 … … 1213 1220 s = '\end{tabular} \n \\caption{%s} \n \label{fig:%s} \n \end{figure} \n \n' %(caption, label) 1214 1221 fid.write(s) 1215 c += 1 1216 if c % 6 == 0: fid.write('\\clearpage \n') 1222 if float((k+1)/div - pp) == 0.: 1223 fid.write('\\clearpage \n') 1224 pp += 1 1217 1225 1218 1226 #### finished generating figures ### -
anuga_validation/solitary_waves/solitary_wave_runup.py
r4025 r4145 54 54 55 55 # Structured mesh 56 dx = 30# Resolution: Length of subdivisions on x axis (length)57 dy = 30# Resolution: Length of subdivisions on y axis (width)56 dx = 5 # Resolution: Length of subdivisions on x axis (length) 57 dy = 5 # Resolution: Length of subdivisions on y axis (width) 58 58 59 59 length = east-west -
anuga_work/production/onslow_2006/make_report.py
r4134 r4145 187 187 \input{execsum} 188 188 189 \clearpage 190 189 191 \\tableofcontents 190 192 … … 210 212 """ 211 213 fid.write(s) 212 213 # Generate latex output for location points214 s = '\\begin{table} \\begin{center} \n'215 fid.write(s)216 s = '\caption{Defined point locations for %s study area.}' %scenario_name217 fid.write(s)218 s = """219 \label{table:locations}220 \\begin{tabular}{|l|l|l|l|}\hline221 \\bf{Point Name} & \\bf{Easting} & \\bf{Northing} & \\bf{Elevation (m)}\\\\ \hline222 """223 fid.write(s)224 225 gauges, locations, elevation = get_gauges_from_file(project.gauge_filename)226 227 for name, gauges, elev in zip(locations, gauges, elevation):228 east = gauges[0]229 north = gauges[1]230 s = '%s & %.2f & %.2f & %.2f \\\\ \hline \n' %(name.replace('_',' '), east, north, elev)231 fid.write(s)232 233 s = '\\end{tabular} \n \end{center} \n \end{table} \n \n'234 fid.write(s)235 236 s = '\\begin{figure}[hbt] \n \centerline{ \includegraphics[width=\paperwidth]{../report_figures/%s}}' %gauge_map237 fid.write(s)238 239 s = """240 \caption{Point locations used for Onslow study.}241 \label{fig:points}242 \end{figure}243 """244 fid.write(s)245 214 246 215 s = '\input{interpretation} \n' … … 257 226 fid.write(s) 258 227 228 # Generate latex output for location points 229 s = '\\begin{table} \\begin{center} \n' 230 fid.write(s) 231 s = '\caption{Defined point locations for %s study area.}' %scenario_name 232 fid.write(s) 233 s = """ 234 \label{table:locations} 235 \\begin{tabular}{|l|l|l|l|}\hline 236 \\bf{Point Name} & \\bf{Easting} & \\bf{Northing} & \\bf{Elevation (m)}\\\\ \hline 237 """ 238 fid.write(s) 239 240 gauges, locations, elevation = get_gauges_from_file(project.gauge_filename) 241 242 for name, gauges, elev in zip(locations, gauges, elevation): 243 east = gauges[0] 244 north = gauges[1] 245 s = '%s & %.2f & %.2f & %.2f \\\\ \hline \n' %(name.replace('_',' '), east, north, elev) 246 fid.write(s) 247 248 s = '\\end{tabular} \n \end{center} \n \end{table} \n \n' 249 fid.write(s) 250 251 #s = '\\begin{figure}[h] \n \centerline{ \includegraphics[width=\paperwidth]{../report_figures/%s}}' %gauge_map 252 s = '\\begin{figure}[h] \n \centerline{ \includegraphics[scale=0.7]{../report_figures/%s}}' %gauge_map 253 fid.write(s) 254 255 s = """ 256 \caption{Point locations used for Onslow study.} 257 \label{fig:points} 258 \end{figure} 259 260 \clearpage 261 """ 262 fid.write(s) 263 259 264 # Closing 260 265 … … 265 270 """ 266 271 fid.write(s) 267 268 #for i, name in enumerate(production_dirs.keys()):269 270 # s = '\input{%s} \n \clearpage \n \n' %damage_maps[production_dirs[name]]271 # fid.write(s)272 272 273 273 s = """ … … 285 285 \input{references} 286 286 287 \clearpage 288 287 289 \\appendix 288 290 … … 302 304 \label{sec:timeseries} 303 305 \input{timeseriesdiscussion} 306 \clearpage 304 307 """ 305 308 fid.write(s) … … 309 312 310 313 s=""" 311 312 \pagebreak313 314 314 315 \section{Damage modelling inputs} -
anuga_work/production/onslow_2006/report/anuga_setup.tex
r4134 r4145 55 55 \begin{figure}[hbt] 56 56 57 \centerline{ \includegraphics[scale=0.1 5]{../report_figures/onslow_resolution_zones.jpg}}57 \centerline{ \includegraphics[scale=0.125]{../report_figures/onslow_resolution_zones.jpg}} 58 58 59 59 \caption{Study area for Onslow scenario highlighting four regions of increased refinement. … … 87 87 \caption{Parameters used in ANUGA for the Onslow scenario.} 88 88 \begin{tabular}{|l|l|l|}\hline 89 Mesh & & \ hline90 & resolution in Region 1 & 500 m$^2$ \ hline91 & resolution in Region 2 & 2500 m$^2$ \ hline92 & resolution in Region 3 & 20000 m$^2$ \ hline93 & remaining resolution & 100 000 m$^2$ \ hline94 Model parameters & & \ hline95 & friction & 0.017 96 & minimum stored height & 0.1 m \ hline89 Mesh & & \\ \hline 90 & resolution in Region 1 & 500 m$^2$ \\ \hline 91 & resolution in Region 2 & 2500 m$^2$ \\ \hline 92 & resolution in Region 3 & 20000 m$^2$ \\ \hline 93 & remaining resolution & 100 000 m$^2$ \\ \hline 94 Model parameters & & \\ \hline 95 & friction & 0.017 \\ \hline 96 & minimum stored height & 0.1 m \\ \hline 97 97 \end{tabular} 98 98 \end{center} -
anuga_work/production/onslow_2006/report/damage_inputs.tex
r4134 r4145 1 \begin{table}[ p]1 \begin{table}[h] 2 2 \begin{center} 3 3 \caption{Framed residential building collapse probability. $h$ is the … … 16 16 \end{table} 17 17 18 \begin{table} 18 \begin{table}[h] 19 19 \begin{center} 20 20 \caption{Mortality and injury state probability} … … 34 34 \end{table} 35 35 36 \begin{table} 36 \begin{table}[h] 37 37 \begin{center} 38 38 \caption{Injury level classificationse. Floor height is assumed to be 30cm} -
anuga_work/production/onslow_2006/report/introduction.tex
r3407 r4145 22 22 detailing critical infrastructure as well as damage modelling estimates. 23 23 24 This report is the first in a series of tsunami assessments 25 of the North West Shelf. The scenario used for this study has 26 an unknown return period, but is considered a plausible event (see 27 Section \ref{sec:methodology}). 28 Subsequent assessments will use refined hazard models with 29 associate return periods. A suite of assessments will be 30 made for Onslow and other localities, as advised by FESA. 24 This report details the impact assessments for a range of tsunami events. 25 These events are based on the probabilistic hazard assessment conducted 26 for the Western Australian coastline. A number of events are selected for 27 return periods of 500, 1000 and 2000 years, see Section \ref{sec:tsunamiscenario}. 28 %This report is the first in a series of tsunami assessments 29 %of the North West Shelf. The scenario used for this study has 30 %an unknown return period, but is considered a plausible event (see 31 %Section \ref{sec:methodology}). 32 %Subsequent assessments will use refined hazard models with 33 %associate return periods. A suite of assessments will be 34 %made for Onslow and other localities, as advised by FESA. 31 35 32 36 Onslow has a population of around 800 and … … 38 42 39 43 The modelling technique to simulate the 40 impact ashore is discussed in Section \ref{sec:methodology} and data inputs 41 discussed in Section \ref{sec:data}. 44 impact ashore is discussed in Section \ref{sec:methodology} with 45 event and data inputs 46 discussed in Sections \ref{sec:tsunamiscenario} and \ref{sec:data} respectively. 42 47 The inundation results are presented and discussed in Section \ref{sec:results} 43 48 and the impact modelling results outlined in Section \ref{sec:impact}. -
anuga_work/production/onslow_2006/report/latexoutput.tex
r4134 r4145 1 \begin{figure}[hbt] 1 \begin{figure}[ht] 2 \centering 3 \begin{tabular}{cc} 4 \includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBeadonPointLoadingBerthstage.png}& 5 \includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBeadonPointLoadingBerthspeed.png}\\ 6 \end{tabular} 7 \caption{Time series for stage and speed at Beadon Point Loading Berth location (elevation -8.70m)} 8 \label{fig:gaugeBeadonPointLoadingBerth} 9 \end{figure} 10 11 \begin{figure}[ht] 12 \centering 13 \begin{tabular}{cc} 14 \includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeHospitalstage.png}& 15 \includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeHospitalspeed.png}\\ 16 \end{tabular} 17 \caption{Time series for depth and speed at Hospital location (elevation 6.02m)} 18 \label{fig:gaugeHospital} 19 \end{figure} 20 21 \begin{figure}[ht] 2 22 \centering 3 23 \begin{tabular}{cc} … … 9 29 \end{figure} 10 30 31 \begin{figure}[ht] 32 \centering 33 \begin{tabular}{cc} 34 \includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugePowerStationstage.png}& 35 \includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugePowerStationspeed.png}\\ 36 \end{tabular} 37 \caption{Time series for depth and speed at Power Station location (elevation 5.81m)} 38 \label{fig:gaugePowerStation} 39 \end{figure} 40 41 \begin{figure}[ht] 42 \centering 43 \begin{tabular}{cc} 44 \includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeAirportRunwaystage.png}& 45 \includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeAirportRunwayspeed.png}\\ 46 \end{tabular} 47 \caption{Time series for depth and speed at Airport Runway location (elevation 5.49m)} 48 \label{fig:gaugeAirportRunway} 49 \end{figure} 50 51 \begin{figure}[ht] 52 \centering 53 \begin{tabular}{cc} 54 \includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBeadonCreekDocksstage.png}& 55 \includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBeadonCreekDocksspeed.png}\\ 56 \end{tabular} 57 \caption{Time series for depth and speed at Beadon Creek Docks location (elevation 1.65m)} 58 \label{fig:gaugeBeadonCreekDocks} 59 \end{figure} 60 61 \begin{figure}[ht] 62 \centering 63 \begin{tabular}{cc} 64 \includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeWestofGroynestage.png}& 65 \includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeWestofGroynespeed.png}\\ 66 \end{tabular} 67 \caption{Time series for stage and speed at West of Groyne location (elevation -2.11m)} 68 \label{fig:gaugeWestofGroyne} 69 \end{figure} 70 71 \begin{figure}[ht] 72 \centering 73 \begin{tabular}{cc} 74 \includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBeadonCreekmouthstage.png}& 75 \includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBeadonCreekmouthspeed.png}\\ 76 \end{tabular} 77 \caption{Time series for stage and speed at Beadon Creek mouth location (elevation -2.90m)} 78 \label{fig:gaugeBeadonCreekmouth} 79 \end{figure} 80 81 \begin{figure}[ht] 82 \centering 83 \begin{tabular}{cc} 84 \includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBeadonCreeksouthofdockstage.png}& 85 \includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBeadonCreeksouthofdockspeed.png}\\ 86 \end{tabular} 87 \caption{Time series for stage and speed at Beadon Creek south of dock location (elevation -1.81m)} 88 \label{fig:gaugeBeadonCreeksouthofdock} 89 \end{figure} 90 91 \begin{figure}[ht] 92 \centering 93 \begin{tabular}{cc} 94 \includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeCentredamwallstage.png}& 95 \includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeCentredamwallspeed.png}\\ 96 \end{tabular} 97 \caption{Time series for depth and speed at Centre dam wall location (elevation 2.09m)} 98 \label{fig:gaugeCentredamwall} 99 \end{figure} 100 101 \begin{figure}[ht] 102 \centering 103 \begin{tabular}{cc} 104 \includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeDamoverflowstage.png}& 105 \includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeDamoverflowspeed.png}\\ 106 \end{tabular} 107 \caption{Time series for depth and speed at Dam overflow location (elevation 0.88m)} 108 \label{fig:gaugeDamoverflow} 109 \end{figure} 110 111 \clearpage 112 \begin{figure}[ht] 113 \centering 114 \begin{tabular}{cc} 115 \includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeLightTowerstage.png}& 116 \includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeLightTowerspeed.png}\\ 117 \end{tabular} 118 \caption{Time series for depth and speed at Light Tower location (elevation 3.86m)} 119 \label{fig:gaugeLightTower} 120 \end{figure} 121 122 \begin{figure}[ht] 123 \centering 124 \begin{tabular}{cc} 125 \includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBeadonBayweststage.png}& 126 \includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBeadonBaywestspeed.png}\\ 127 \end{tabular} 128 \caption{Time series for stage and speed at Beadon Bay west location (elevation -4.62m)} 129 \label{fig:gaugeBeadonBaywest} 130 \end{figure} 131 132 \begin{figure}[ht] 133 \centering 134 \begin{tabular}{cc} 135 \includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBeadonBayeaststage.png}& 136 \includegraphics[width=0.49\linewidth, height=50mm]{../report_figures/gaugeBeadonBayeastspeed.png}\\ 137 \end{tabular} 138 \caption{Time series for stage and speed at Beadon Bay east location (elevation -3.56m)} 139 \label{fig:gaugeBeadonBayeast} 140 \end{figure} 141 -
anuga_work/production/onslow_2006/report/modelling_methodology.tex
r4134 r4145 29 29 evenly spaced along the Sunda Arc subduction zone. 30 30 The assessment used the Method of Splitting Tsunamis (MOST) 31 \cite{VT:MOST} model. A detailed probabilistic hazard map has now been completed 32 for the WA coastline \cite{prob:fesa} which is based on the paradigm used by the URS corporation's 33 Probabilistic Tsunami Hazard Analysis \cite{somerville:urs}. 31 \cite{VT:MOST} model. This preliminary hazard map had no probability attached to the event which is 32 required to conduct a tsunami risk assessment. Using 33 the Probabilistic Tsunami Hazard Analysis \cite{somerville:urs} paradigm 34 used by the URS corporation, a detailed probabilistic hazard map has now been completed 35 for the WA coastline \cite{prob:fesa}. 34 36 35 37 While MOST and URS are suitable for generating and propagating the tsunami wave from its source, … … 57 59 Further details 58 60 regarding the inundation modelling requirements for this study can be found in 59 Appendix \ref{ anugasetup}.61 Appendix \ref{sec:anugasetup}. 60 62 61 63 The risk of a given tsunami scenario can only be determined when the likelihood … … 77 79 for a range of probabilities (or return periods). As Figure \ref{fig:probonslow} 78 80 shows, for a given probability, a number of events are possible. The resulting 79 impact to Onslow would then vary depending on the source of the event. 81 impact to Onslow would then vary depending on the source of the event. Further detail 82 on the tsunami scenarios are outlined in Section \ref{sec:tsunamiscenario}. 80 83 81 84 % used for the 2005 report when looking at one event -
anuga_work/production/onslow_2006/report/onslow_2006_report.tex
r4134 r4145 53 53 \input{execsum} 54 54 55 \clearpage 56 55 57 \tableofcontents 56 58 … … 74 76 \label{sec:results} 75 77 78 \input{interpretation} 79 \input{HAT_map} 80 \clearpage 81 82 \input{MSL_map} 83 \clearpage 84 85 \input{LAT_map} 86 \clearpage 87 76 88 \begin{table} \begin{center} 77 89 \caption{Defined point locations for onslow study area.} … … 97 109 \end{table} 98 110 99 \begin{figure}[h bt]100 \centerline{ \includegraphics[ width=\paperwidth]{../report_figures/onslow_dli_gauge.jpg}}111 \begin{figure}[h] 112 \centerline{ \includegraphics[scale=0.7]{../report_figures/onslow_dli_gauge.jpg}} 101 113 \caption{Point locations used for Onslow study.} 102 114 \label{fig:points} 103 115 \end{figure} 104 \input{interpretation} 105 \input{HAT_map} 106 \clearpage 107 108 \input{MSL_map} 109 \clearpage 110 111 \input{LAT_map} 112 \clearpage 113 116 117 \clearpage 114 118 115 119 \section{Impact modelling} … … 130 134 \input{references} 131 135 136 \clearpage 137 132 138 \appendix 133 139 … … 147 153 \label{sec:timeseries} 148 154 \input{timeseriesdiscussion} 155 \clearpage 149 156 \input{latexoutput} 150 157 \clearpage 151 158 152 153 154 \pagebreak155 159 156 160 \section{Damage modelling inputs}
Note: See TracChangeset
for help on using the changeset viewer.