source: anuga_work/production/geraldton/geraldton_rerun/Geraldton4distribution.py @ 7118

Last change on this file since 7118 was 7118, checked in by kristy, 15 years ago
File size: 5.2 KB
Line 
1# ---------------------------------------------------------------------------
2# Carnarvon4distribution.py
3# Created on: Thu Apr 30 2009 03:51:03 PM
4#   (generated by ArcGIS/ModelBuilder)
5# ---------------------------------------------------------------------------
6
7# Import system modules
8import sys, string, os, arcgisscripting
9
10# Create the Geoprocessor object
11gp = arcgisscripting.create()
12
13# Check out any necessary licenses
14gp.CheckOutExtension("spatial")
15
16# Load required toolboxes...
17gp.AddToolbox("C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Spatial Analyst Tools.tbx")
18gp.AddToolbox("C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Data Management Tools.tbx")
19gp.overwriteoutput = 1
20
21community = 'Geraldton'
22output_loc = "\\\\nas2\\gemd\\georisk_models\\inundation\\data\\western_australia\\FESA_outputs_Dec_08\\ARCgrids\\April_2009\\"
23input_scenario = "\\\\nas2\\gemd\\georisk_models\\inundation\\data\\western_australia\\geraldton_tsunami_scenario\\"
24input_output = input_scenario + "\\anuga\\outputs\\"
25
26location = 'CBD'
27in_Tsunami_1_HAT = input_output + "20090324_140411_run_final_0.6_27283_Bt_kvanputt\\raster.gdb\\"
28in_Tsunami_1_MSL = input_output + "20090327_160022_run_final_0_27283_Bt_kvanputt\\raster.gdb\\"
29in_Tsunami_2_HAT = input_output + "20090324_134504_run_final_0.6_27255_Bt_kvanputt\\raster.gdb\\"
30in_Tsunami_2_MSL = input_output + "20090327_160043_run_final_0_27255_Bt_kvanputt\\raster.gdb\\"
31in_Tsunami_3_HAT = input_output + "20090324_140741_run_final_0.6_68693_Bt_kvanputt\\raster.gdb\\"
32in_Tsunami_3_MSL = input_output + "20090327_155230_run_final_0_68693_Bt_kvanputt\\raster.gdb\\"
33
34Tsunami_1_HAT = output_loc + "Tsunami1-Exmouth\\" + community + "_HAT.gdb\\"
35Tsunami_1_MSL = output_loc + "Tsunami1-Exmouth\\" + community + "_MSL.gdb\\"
36Tsunami_2_HAT = output_loc + "Tsunami2-Perth\\" + community + "_HAT.gdb\\"
37Tsunami_2_MSL = output_loc + "Tsunami2-Perth\\" + community + "_MSL.gdb\\"
38Tsunami_3_HAT = output_loc + "Tsunami3-1833\\" + community + "_HAT.gdb\\"
39Tsunami_3_MSL = output_loc + "Tsunami3-1833\\" + community + "_MSL.gdb\\"
40
41depth = 'CDB_depth_E'
42speed = 'CDB_speed_'
43hr = community + '_runup'
44
45# Local variables...
46# polygon clip
47if location == 'CBD':
48    extent = input_scenario + "map_work\\Geraldton.gdb\Internal_Polygons\CBD_500m"
49
50# Output locations
51
52
53
54hr_1_HAT = Tsunami_1_HAT +"highest_runup"
55hr_1_MSL = Tsunami_1_MSL +"highest_runup"
56hr_2_HAT = Tsunami_2_HAT +"highest_runup"
57hr_2_MSL = Tsunami_2_MSL +"highest_runup"
58hr_3_HAT = Tsunami_3_HAT +"highest_runup"
59hr_3_MSL = Tsunami_3_MSL +"highest_runup"
60
61depth_1_HAT = Tsunami_1_HAT + location + "_depth"
62depth_1_MSL = Tsunami_1_MSL + location + "_depth"
63depth_2_HAT = Tsunami_2_HAT + location + "_depth"
64depth_2_MSL = Tsunami_2_MSL + location + "_depth"
65depth_3_HAT = Tsunami_3_HAT + location + "_depth"
66depth_3_MSL = Tsunami_3_MSL + location + "_depth"
67
68speed_1_HAT = Tsunami_1_HAT + location + "_speed"
69speed_1_MSL = Tsunami_1_MSL + location + "_speed"
70speed_2_HAT = Tsunami_2_HAT + location + "_speed"
71speed_2_MSL = Tsunami_2_MSL + location + "_speed"
72speed_3_HAT = Tsunami_3_HAT + location + "_speed"
73speed_3_MSL = Tsunami_3_MSL + location + "_speed"
74
75
76# Input locations
77
78in_depth_1_HAT = in_Tsunami_1_HAT + depth
79in_depth_1_MSL = in_Tsunami_1_MSL + depth
80in_depth_2_HAT = in_Tsunami_2_HAT + depth
81in_depth_2_MSL = in_Tsunami_2_MSL + depth
82in_depth_3_HAT = in_Tsunami_3_HAT + depth
83in_depth_3_MSL = in_Tsunami_3_MSL + depth
84
85in_speed_1_HAT = in_Tsunami_1_HAT + speed
86in_speed_1_MSL = in_Tsunami_1_MSL + speed
87in_speed_2_HAT = in_Tsunami_2_HAT + speed
88in_speed_2_MSL = in_Tsunami_2_MSL + speed
89in_speed_3_HAT = in_Tsunami_3_HAT + speed
90in_speed_3_MSL = in_Tsunami_3_MSL + speed
91
92in_hr_1_HAT = in_Tsunami_1_HAT + hr
93in_hr_1_MSL = in_Tsunami_1_MSL + hr
94in_hr_2_HAT = in_Tsunami_2_HAT + hr
95in_hr_2_MSL = in_Tsunami_2_MSL + hr
96in_hr_3_HAT = in_Tsunami_3_HAT + hr
97in_hr_3_MSL = in_Tsunami_3_MSL + hr
98
99
100# Process: Extract by Mask
101
102print 'Extract by mask for depth'
103gp.ExtractByMask_sa(in_depth_1_HAT, extent, depth_1_HAT)
104gp.ExtractByMask_sa(in_depth_1_MSL, extent, depth_1_MSL)
105gp.ExtractByMask_sa(in_depth_2_HAT, extent, depth_2_HAT)
106gp.ExtractByMask_sa(in_depth_2_MSL, extent, depth_2_MSL)
107gp.ExtractByMask_sa(in_depth_3_HAT, extent, depth_3_HAT)
108gp.ExtractByMask_sa(in_depth_3_MSL, extent, depth_3_MSL)
109
110print 'Extract by mask for speed'
111gp.ExtractByMask_sa(in_speed_1_HAT, extent, speed_1_HAT)
112gp.ExtractByMask_sa(in_speed_1_MSL, extent, speed_1_MSL)
113gp.ExtractByMask_sa(in_speed_2_HAT, extent, speed_2_HAT)
114gp.ExtractByMask_sa(in_speed_2_MSL, extent, speed_2_MSL)
115gp.ExtractByMask_sa(in_speed_3_HAT, extent, speed_3_HAT)
116gp.ExtractByMask_sa(in_speed_3_MSL, extent, speed_3_MSL)
117
118
119# Process: Copy Features...
120
121print 'Copy Highest Runup'
122gp.CopyFeatures_management(in_hr_1_HAT, hr_1_HAT, "", "0", "0", "0")
123gp.CopyFeatures_management(in_hr_1_MSL, hr_1_MSL, "", "0", "0", "0")
124gp.CopyFeatures_management(in_hr_2_HAT, hr_2_HAT, "", "0", "0", "0")
125gp.CopyFeatures_management(in_hr_2_MSL, hr_2_MSL, "", "0", "0", "0")
126gp.CopyFeatures_management(in_hr_3_HAT, hr_3_HAT, "", "0", "0", "0")
127gp.CopyFeatures_management(in_hr_3_MSL, hr_3_MSL, "", "0", "0", "0")
128
Note: See TracBrowser for help on using the repository browser.