Changeset 4058 for anuga_work/production/sydney_2006
- Timestamp:
- Dec 5, 2006, 5:30:41 PM (18 years ago)
- Location:
- anuga_work/production/sydney_2006
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/sydney_2006/project_slide.py
r4056 r4058 31 31 nsw100 = 'dem_utm' 32 32 33 # AHO33 # survey data 34 34 offshore_name1 = 'surveyAreaA' 35 35 offshore_name2 = 'surveyAreaB' 36 36 offshore_name3 = 'surveyAreaC' 37 38 # AHO data 39 offshore_name4 = '1000003611export' 40 offshore_name5 = '1000003613export' 41 offshore_name6 = '1000003614export' 42 offshore_name7 = '1000003627export' 43 offshore_name8 = '1000003628export' 44 offshore_name9 = 'AHDexport' 37 45 38 46 #swollen/ all data output … … 63 71 offshore_dem_name2 = datadir + offshore_name2 64 72 offshore_dem_name3 = datadir + offshore_name3 73 offshore_dem_name4 = datadir + offshore_name4 74 offshore_dem_name5 = datadir + offshore_name5 75 offshore_dem_name6 = datadir + offshore_name6 76 offshore_dem_name7 = datadir + offshore_name7 77 offshore_dem_name8 = datadir + offshore_name8 78 offshore_dem_name9 = datadir + offshore_name9 65 79 combined_dem_name = datadir + 'sydney_combined_elevation' 66 80 … … 134 148 bulli_density = 1.46 135 149 bulli_slope = 4.0 136 bulli_alpha = 126.0 +90.0150 bulli_alpha = 126.0 - 90.0 137 151 138 152 shovel_depth = 968.0 … … 142 156 shovel_density = 1.49 143 157 shove_slope = 4.0 144 shovel_alpha = 118.0 +90.0158 shovel_alpha = 118.0 - 90.0 145 159 146 160 yacaaba_depth = 1119.0 … … 150 164 yacaaba_density = 1.48 151 165 yacaaba_slope = 3.0 152 yacaaba_alpha = 133.0 +90.0166 yacaaba_alpha = 133.0 - 90.0 -
anuga_work/production/sydney_2006/run_sydney_slide.py
r4056 r4058 83 83 84 84 print 'create offshore' 85 G1 = Geospatial_data(file_name = project_slide.offshore_dem_name1 + '.xya')+\ 86 Geospatial_data(file_name = project_slide.offshore_dem_name2 + '.xya')+\ 87 Geospatial_data(file_name = project_slide.offshore_dem_name3 + '.xya')+\ 85 G11 = Geospatial_data(file_name = project_slide.offshore_dem_name1 + '.xya')+\ 86 Geospatial_data(file_name = project_slide.offshore_dem_name2 + '.xya')+\ 87 Geospatial_data(file_name = project_slide.offshore_dem_name3 + '.xya') 88 G12 = Geospatial_data(file_name = project_slide.offshore_dem_name4 + '.xya')+\ 89 Geospatial_data(file_name = project_slide.offshore_dem_name5 + '.xya')+\ 90 Geospatial_data(file_name = project_slide.offshore_dem_name6 + '.xya')+\ 91 Geospatial_data(file_name = project_slide.offshore_dem_name7 + '.xya')+\ 92 Geospatial_data(file_name = project_slide.offshore_dem_name8 + '.xya')+\ 93 Geospatial_data(file_name = project_slide.offshore_dem_name9 + '.xya') 88 94 print 'create onshore' 89 95 G2 = Geospatial_data(file_name = project_slide.on_offshore25_dem_name + '.pts') … … 91 97 G4 = Geospatial_data(file_name = project_slide.nsw_dem_name + '.pts') 92 98 print 'add' 93 G = G1.clip(Geospatial(project_slide.poly_surveyclip)) +\ 99 G = G11.clip(Geospatial_data(project_slide.poly_surveyclip)) +\ 100 G12.clip(Geospatial_data(project_slide.polyAll)) 94 101 G2.clip(Geospatial_data(project_slide.poly_25mclip)) +\ 95 102 G3.clip(Geospatial_data(project_slide.poly_origsyd)) +\ 96 (G4.clip(Geospatial_data(project_slide.polyAll))).clip_outside(Geospatial_data(project_slide.poly_surveyclip)) 103 (G4.clip(Geospatial_data(project_slide.polyAll))).clip_outside(Geospatial_data(project_slide.poly_surveyclip)).clip_outside(Geospatial_data(project_slide.poly_origsyd)) 97 104 print 'export points' 98 105 G.export_points_file(project_slide.combined_dem_name + '.pts')
Note: See TracChangeset
for help on using the changeset viewer.