Changeset 4314
- Timestamp:
- Mar 22, 2007, 5:25:44 PM (18 years ago)
- Location:
- anuga_work/production
- Files:
-
- 15 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/broome_2006/run_building_inundation.py
r4310 r4314 19 19 from os import sep 20 20 21 from damage.inundation_damage import add_depth_and_momentum2csv21 from anuga.damage_modelling.inundation_damage import add_depth_and_momentum2csv, inundation_damage 22 22 import project 23 23 … … 32 32 buildings_filename_out = file_loc + project.buildings_filename_out 33 33 34 add_depth_and_momentum2csv(swwfile, buildings_filename, buildings_filename_out)34 #add_depth_and_momentum2csv(swwfile, buildings_filename, buildings_filename_out) 35 35 inundation_damage(swwfile, buildings_filename, buildings_filename_out) 36 36 -
anuga_work/production/busselton_2006/run_building_inundation.py
r4310 r4314 19 19 from os import sep 20 20 21 from damage.inundation_damage import add_depth_and_momentum2csv21 from anuga.damage_modelling.inundation_damage import add_depth_and_momentum2csv, inundation_damage 22 22 import project 23 23 … … 32 32 buildings_filename_out = file_loc + project.buildings_filename_out 33 33 34 add_depth_and_momentum2csv(swwfile, buildings_filename, buildings_filename_out)34 #add_depth_and_momentum2csv(swwfile, buildings_filename, buildings_filename_out) 35 35 inundation_damage(swwfile, buildings_filename, buildings_filename_out) 36 36 -
anuga_work/production/dampier_2006/run_building_inundation.py
r4310 r4314 19 19 from os import sep 20 20 21 from damage.inundation_damage import add_depth_and_momentum2csv21 from anuga.damage_modelling.inundation_damage import add_depth_and_momentum2csv, inundation_damage 22 22 import project 23 23 … … 32 32 buildings_filename_out = file_loc + project.buildings_filename_out 33 33 34 add_depth_and_momentum2csv(swwfile, buildings_filename, buildings_filename_out)34 #add_depth_and_momentum2csv(swwfile, buildings_filename, buildings_filename_out) 35 35 inundation_damage(swwfile, buildings_filename, buildings_filename_out) 36 36 -
anuga_work/production/exmouth_2006/run_building_inundation.py
r4310 r4314 19 19 from os import sep 20 20 21 from damage.inundation_damage import add_depth_and_momentum2csv21 from anuga.damage_modelling.inundation_damage import add_depth_and_momentum2csv, inundation_damage 22 22 import project 23 23 … … 32 32 buildings_filename_out = file_loc + project.buildings_filename_out 33 33 34 add_depth_and_momentum2csv(swwfile, buildings_filename, buildings_filename_out)34 #add_depth_and_momentum2csv(swwfile, buildings_filename, buildings_filename_out) 35 35 inundation_damage(swwfile, buildings_filename, buildings_filename_out) 36 36 -
anuga_work/production/newcastle_2006/export_results.py
r4276 r4314 20 20 #time_dir = '20070219_060722' #bulli 21 21 #time_dir = '20070219_060808' #shovel 22 time_dir = '20070219_060854' #yacaaba 22 #time_dir = '20070219_060854' #yacaaba 23 24 #corrected yacaaba and birubi 25 time_dir = '20070322_004007' # yacaaba_c1 - potential 26 #time_dir = '20070322_003255' # yacaaba_c2 - historical 27 #time_dir = '20070322_032705' # birubi_c1 - historical 28 #time_dir = '20070322_032618' # birubi_c2 - potential 29 30 # new region 31 #time_dir = '20070322_033350' # birubi 32 #time_dir = '20070322_033545' # yacaaba 33 #time_dir = '20070322_033716' # shovel 34 #time_dir = '20070322_033608' # bulli 23 35 24 36 directory = project_slide.outputdir … … 39 51 40 52 if which_var == 2: # Depth 41 outname = name + '_depth-yacaaba -michelle'53 outname = name + '_depth-yacaaba' 42 54 quantityname = 'stage-elevation' 43 55 -
anuga_work/production/newcastle_2006/project_slide.py
r4276 r4314 79 79 # bounding polygon for study area 80 80 polyAll = read_polygon(polygondir+'extentC3.csv') 81 polyAll_c0 = read_polygon(polygondir+'extent_c0.csv') 81 82 82 83 print 'Area of bounding polygon', polygon_area(polyAll)/1000000.0 … … 93 94 print 'Area of CBD polygon', polygon_area(poly_newcastle)/1000000.0 94 95 96 # region c0 97 poly_local_c0 = read_polygon(polygondir+'local_c0.csv') 98 poly_entrance_c0 = read_polygon(polygondir+'entrance.csv') 99 95 100 #plot_polygons([polyAll,poly_local,poly_newcastle],'fig',verbose=True) 96 101 ################################################################### … … 103 108 northingmin_nsw = 6253970 104 109 northingmax_nsw = 6424600 110 111 # clipping for nsw grid for conversion 112 eastingminc0_nsw = 322440 113 eastingmaxc0_nsw = 436250 114 northingminc0_nsw = 6280000 115 northingmaxc0_nsw = 6410000 105 116 106 117 poly_surveyclip = read_polygon(polygondir+'surveyclip.csv') … … 125 136 slide_origin_bulli = [370110,6189489] 126 137 slide_origin_shovel = [354753,6187397] 127 slide_origin_yacaaba = [463012,6334735] 138 #slide_origin_yacaaba = [463012,6334735] original work was actually Birubi 139 slide_origin_yacaaba = [476465, 6362424] 140 slide_origin_birubi = [462267,6335100] 128 141 129 142 # potential slides: these correspond approximately to modelled depths 130 143 slide_origin_bulli_b = [405617,6249852] 131 slide_origin_bulli_c = [465283,6335210] 144 slide_origin_bulli_c0 = [438583,6308510] 145 slide_origin_bulli_c1 = [465283,6335210] 146 slide_origin_bulli_c2 = [483783,6358810] 132 147 slide_origin_shovel_b = [396262,6252144] 133 slide_origin_shovel_c = [458694,6340352] 134 slide_origin_yacaaba_a = [361753,6188786] 135 slide_origin_yacaaba_b = [401897,6250643] 148 slide_origin_shovel_c0 = [431683,6311110] 149 slide_origin_shovel_c1 = [458694,6340352] 150 slide_origin_shovel_c2 = [475589,6363030] 151 slide_origin_yacaaba_a = [359796,6190943] 152 slide_origin_yacaaba_b = [397330,6251857] 153 slide_origin_yacaaba_c0 = [431283,6312110] 154 slide_origin_yacaaba_c1 = [459983,6339110] 155 slide_origin_birubi_a = [362811,6191183] 156 slide_origin_birubi_b = [403573,6250268] 157 slide_origin_birubi_c0 = [437183,6309410] 158 slide_origin_birubi_c2 = [483392,6359044] 136 159 137 160 #bulli_depth = 2087.0 … … 154 177 155 178 #yacaaba_depth = 1119.0 156 yacaaba_depth = 1222.5157 yacaaba_length = 7050.0158 yacaaba_thickness = 144.0159 yacaaba_width = 3080.0179 yacaaba_depth = 938 180 yacaaba_length = 4189 181 yacaaba_thickness = 53 182 yacaaba_width = 2898 160 183 yacaaba_density = 1.48 161 yacaaba_slope = 3.0184 yacaaba_slope = 2.3 162 185 yacaaba_alpha = 133.0 - 90.0 186 187 birubi_depth = 1320. 188 birubi_length = 9903.0 189 birubi_width = 4150. 190 birubi_density = 1.48 191 birubi_slope = 3.7 192 birubi_alpha = 133.0 - 90.0 193 birubi_thickness = 140.0 -
anuga_work/production/newcastle_2006/run_newcastle_slide.py
r4276 r4314 143 143 from smf import slide_tsunami 144 144 145 tsunami_source = slide_tsunami(length=project_slide. yacaaba_length,146 width=project_slide. yacaaba_width,147 depth=project_slide. yacaaba_depth,148 slope=project_slide. yacaaba_slope,149 thickness=project_slide. yacaaba_thickness,150 x0=project_slide.slide_origin_ yacaaba[0],151 y0=project_slide.slide_origin_ yacaaba[1],152 alpha=project_slide. yacaaba_alpha,145 tsunami_source = slide_tsunami(length=project_slide.shovel_length, 146 width=project_slide.shovel_width, 147 depth=project_slide.shovel_depth, 148 slope=project_slide.shovel_slope, 149 thickness=project_slide.shovel_thickness, 150 x0=project_slide.slide_origin_shovel_c2[0], 151 y0=project_slide.slide_origin_shovel_c2[1], 152 alpha=project_slide.shovel_alpha, 153 153 domain=domain) 154 154 -
anuga_work/production/onslow_2006/run_building_inundation.py
r4310 r4314 19 19 from os import sep 20 20 21 from damage.inundation_damage import add_depth_and_momentum2csv21 from anuga.damage_modelling.inundation_damage import add_depth_and_momentum2csv, inundation_damage 22 22 import project 23 23 … … 32 32 buildings_filename_out = file_loc + project.buildings_filename_out 33 33 34 add_depth_and_momentum2csv(swwfile, buildings_filename, buildings_filename_out)34 #add_depth_and_momentum2csv(swwfile, buildings_filename, buildings_filename_out) 35 35 inundation_damage(swwfile, buildings_filename, buildings_filename_out) 36 36 -
anuga_work/production/perth_2006/run_building_inundation.py
r4310 r4314 19 19 from os import sep 20 20 21 from damage.inundation_damage import add_depth_and_momentum2csv21 from anuga.damage_modelling.inundation_damage import add_depth_and_momentum2csv, inundation_damage 22 22 import project 23 23 … … 32 32 buildings_filename_out = file_loc + project.buildings_filename_out 33 33 34 add_depth_and_momentum2csv(swwfile, buildings_filename, buildings_filename_out)34 #add_depth_and_momentum2csv(swwfile, buildings_filename, buildings_filename_out) 35 35 inundation_damage(swwfile, buildings_filename, buildings_filename_out) 36 36 -
anuga_work/production/pt_hedland_2006/run_building_inundation.py
r4310 r4314 19 19 from os import sep 20 20 21 from damage.inundation_damage import add_depth_and_momentum2csv21 from anuga.damage_modelling.inundation_damage import add_depth_and_momentum2csv, inundation_damage 22 22 import project 23 23 … … 32 32 buildings_filename_out = file_loc + project.buildings_filename_out 33 33 34 add_depth_and_momentum2csv(swwfile, buildings_filename, buildings_filename_out)34 #add_depth_and_momentum2csv(swwfile, buildings_filename, buildings_filename_out) 35 35 inundation_damage(swwfile, buildings_filename, buildings_filename_out) 36 36 -
anuga_work/production/sydney_2006/export_results.py
r4294 r4314 8 8 #time_dir = '20061212_012705' #shovel 9 9 #time_dir = '20061212_064807' #yacaaba 10 10 # correct depths 11 11 #time_dir = '20070121_213509' #shovel 12 12 #time_dir = '20070119_055222' #yacaaba 13 13 #time_dir = '20070121_213444' #bulli 14 #corrected yacaaba and birubi 15 time_dir = '20070322_035214' # yacaaba 16 #time_dir = '20070322_035135' # birubi 14 17 15 #directory = project_slide.outputdir16 #name = directory + time_dir + sep + project_slide.basename17 name = 'test500000.035000'18 directory = project_slide.outputdir 19 name = directory + time_dir + sep + project_slide.basename 20 #name = 'test500000.035000' 18 21 #name = 'test' 19 22 … … 32 35 33 36 if which_var == 2: # Depth 34 outname = name + '_depth- eq'37 outname = name + '_depth-yacaaba' 35 38 quantityname = 'stage-elevation' 36 39 -
anuga_work/production/sydney_2006/project_slide.py
r4294 r4314 82 82 # Domain definitions 83 83 ############################### 84 """ 84 85 85 # bounding polygon for study area 86 86 polyAll = read_polygon(polygondir+'extent3.csv') … … 133 133 northingmin = 6244800 134 134 northingmax = 6274550 135 """ 135 136 136 # exporting asc grid 137 137 eastingmin = 318922+10000 … … 152 152 slide_origin_bulli = [370110,6189489] 153 153 slide_origin_shovel = [354753,6187397] 154 slide_origin_yacaaba = [463012,6334735] 154 #slide_origin_yacaaba = [463012,6334735] original work was actually Birubi 155 slide_origin_yacaaba = [476465, 6362424] 156 slide_origin_birubi = [462267,6335100] 155 157 156 158 # potential slides: these correspond approximately to modelled depths 157 159 slide_origin_bulli_b = [405617,6249852] 158 slide_origin_bulli_c = [465283,6335210] 160 slide_origin_bulli_c0 = [438583,6308510] 161 slide_origin_bulli_c1 = [465283,6335210] 162 slide_origin_bulli_c2 = [483783,6358810] 159 163 slide_origin_shovel_b = [396262,6252144] 160 slide_origin_shovel_c = [458694,6340352] 161 slide_origin_yacaaba_a = [361753,6188786] 162 slide_origin_yacaaba_b = [401897,6250643] 164 slide_origin_shovel_c0 = [431683,6311110] 165 slide_origin_shovel_c1 = [458694,6340352] 166 slide_origin_shovel_c2 = [475589,6363030] 167 slide_origin_yacaaba_a = [359796,6190943] 168 slide_origin_yacaaba_b = [397330,6251857] 169 slide_origin_yacaaba_c0 = [431283,6312110] 170 slide_origin_yacaaba_c1 = [459983,6339110] 171 slide_origin_birubi_a = [362811,6191183] 172 slide_origin_birubi_b = [403573,6250268] 173 slide_origin_birubi_c0 = [437183,6309410] 174 slide_origin_birubi_c2 = [483392,6359044] 163 175 164 176 #bulli_depth = 2087.0 … … 181 193 182 194 #yacaaba_depth = 1119.0 183 yacaaba_depth = 1222.5184 yacaaba_length = 7050.0185 yacaaba_thickness = 144.0186 yacaaba_width = 3080.0195 yacaaba_depth = 938. 196 yacaaba_length = 4189. 197 yacaaba_thickness = 53. 198 yacaaba_width = 2898. 187 199 yacaaba_density = 1.48 188 yacaaba_slope = 3.0200 yacaaba_slope = 2.3 189 201 yacaaba_alpha = 133.0 - 90.0 202 203 birubi_depth = 1320. 204 birubi_length = 9903.0 205 birubi_width = 4150. 206 birubi_density = 1.48 207 birubi_slope = 3.7 208 birubi_alpha = 133.0 - 90.0 209 birubi_thickness = 140. -
anuga_work/production/wollongong_2006/export_results.py
r4276 r4314 21 21 #time_dir = '20070213_035110' #bulli 22 22 #time_dir = '20070213_035146' #shovel 23 time_dir = '20070213_035302' #yacaaba 23 #time_dir = '20070213_035302' #yacaaba 24 25 #corrected yacaaba and birubi 26 time_dir = '20070322_035513' # yacaaba 27 #time_dir = '20070322_035640' # birubi 24 28 25 29 directory = project_slide.outputdir … … 40 44 41 45 if which_var == 2: # Depth 42 outname = name + '_depth-yacaaba -amy'46 outname = name + '_depth-yacaaba' 43 47 quantityname = 'stage-elevation' 44 48 -
anuga_work/production/wollongong_2006/project_slide.py
r4276 r4314 130 130 slide_origin_bulli = [370110,6189489] 131 131 slide_origin_shovel = [354753,6187397] 132 slide_origin_yacaaba = [463012,6334735] 132 #slide_origin_yacaaba = [463012,6334735] original work was actually Birubi 133 slide_origin_yacaaba = [476465, 6362424] 134 slide_origin_birubi = [462267,6335100] 133 135 134 136 # potential slides: these correspond approximately to modelled depths 135 137 slide_origin_bulli_b = [405617,6249852] 136 slide_origin_bulli_c = [465283,6335210] 138 slide_origin_bulli_c0 = [438583,6308510] 139 slide_origin_bulli_c1 = [465283,6335210] 140 slide_origin_bulli_c2 = [483783,6358810] 137 141 slide_origin_shovel_b = [396262,6252144] 138 slide_origin_shovel_c = [458694,6340352] 139 slide_origin_yacaaba_a = [361753,6188786] 140 slide_origin_yacaaba_b = [401897,6250643] 142 slide_origin_shovel_c0 = [431683,6311110] 143 slide_origin_shovel_c1 = [458694,6340352] 144 slide_origin_shovel_c2 = [475589,6363030] 145 slide_origin_yacaaba_a = [359796,6190943] 146 slide_origin_yacaaba_b = [397330,6251857] 147 slide_origin_yacaaba_c0 = [431283,6312110] 148 slide_origin_yacaaba_c1 = [459983,6339110] 149 slide_origin_birubi_a = [362811,6191183] 150 slide_origin_birubi_b = [403573,6250268] 151 slide_origin_birubi_c0 = [437183,6309410] 152 slide_origin_birubi_c2 = [483392,6359044] 141 153 142 154 #bulli_depth = 2087.0 … … 159 171 160 172 #yacaaba_depth = 1119.0 161 yacaaba_depth = 1222.5162 yacaaba_length = 7050.0163 yacaaba_thickness = 144.0164 yacaaba_width = 3080.0173 yacaaba_depth = 938. 174 yacaaba_length = 4189. 175 yacaaba_thickness = 53. 176 yacaaba_width = 2898. 165 177 yacaaba_density = 1.48 166 yacaaba_slope = 3.0178 yacaaba_slope = 2.3 167 179 yacaaba_alpha = 133.0 - 90.0 168 180 181 birubi_depth = 1320. 182 birubi_length = 9903.0 183 birubi_width = 4150. 184 birubi_density = 1.48 185 birubi_slope = 3.7 186 birubi_alpha = 133.0 - 90.0 187 birubi_thickness = 140. -
anuga_work/production/wollongong_2006/run_gong_slide.py
r4276 r4314 147 147 from smf import slide_tsunami 148 148 149 tsunami_source = slide_tsunami(length=project_slide. yacaaba_length,150 width=project_slide. yacaaba_width,151 depth=project_slide. yacaaba_depth,152 slope=project_slide. yacaaba_slope,153 thickness=project_slide. yacaaba_thickness,154 x0=project_slide.slide_origin_ yacaaba_a[0],155 y0=project_slide.slide_origin_ yacaaba_a[1],156 alpha=project_slide. yacaaba_alpha,149 tsunami_source = slide_tsunami(length=project_slide.birubi_length, 150 width=project_slide.birubi_width, 151 depth=project_slide.birubi_depth, 152 slope=project_slide.birubi_slope, 153 thickness=project_slide.birubi_thickness, 154 x0=project_slide.slide_origin_birubi_a[0], 155 y0=project_slide.slide_origin_birubi_a[1], 156 alpha=project_slide.birubi_alpha, 157 157 domain=domain) 158 158
Note: See TracChangeset
for help on using the changeset viewer.