Changeset 6560

Timestamp:

Mar 20, 2009, 4:02:42 PM (16 years ago)

Author:

kristy

Message:

added datasets for bunbury

File:

• anuga_work/production/busselton/build_busselton.py (modified) (4 diffs)

anuga_work/production/busselton/build_busselton.py

@@ -62 +62 @@
 offshore_in_dir_name4 = project.offshore_in_dir_name4
 offshore_in_dir_name5 = project.offshore_in_dir_name5
+offshore_in_dir_name6 = project.offshore_in_dir_name6
+offshore_in_dir_name7 = project.offshore_in_dir_name7
 
 # output elevation directory filenames
@@ -73 +75 @@
 offshore_dir_name4 = project.offshore_dir_name4
 offshore_dir_name5 = project.offshore_dir_name5
+offshore_dir_name6 = project.offshore_dir_name6
+offshore_dir_name7 = project.offshore_dir_name7
+
 # creates DEM from asc data
 print "creates DEMs from asc data"
 convert_dem_from_ascii2netcdf(onshore_in_dir_name, basename_out=onshore_dir_name, use_cache=True, verbose=True)
-convert_dem_from_ascii2netcdf(offshore_in_dir_name5, basename_out=offshore_dir_name5, use_cache=True, verbose=True)
+convert_dem_from_ascii2netcdf(offshore_in_dir_name7, basename_out=offshore_dir_name7, use_cache=True, verbose=True)
 
 # creates pts file for onshore DEM
 print "creates pts file for onshore DEM"
 dem2pts(onshore_dir_name ,use_cache=True,verbose=True)
-dem2pts(offshore_dir_name5 ,use_cache=True,verbose=True)
+dem2pts(offshore_dir_name7 ,use_cache=True,verbose=True)
 
 # create onshore pts files
@@ -103 +108 @@
 G_off4 = Geospatial_data(file_name = offshore_in_dir_name4)
 print'create Geospatial data9 objects from topographies'
-G_off5 = Geospatial_data(file_name = offshore_dir_name5 + '.pts')
+G_off5 = Geospatial_data(file_name = offshore_in_dir_name5)
+print'create Geospatial data9 objects from topographies'
+G_off6 = Geospatial_data(file_name = offshore_in_dir_name6)
+print'create Geospatial data9 objects from topographies'
+G_off7 = Geospatial_data(file_name = offshore_dir_name7 + '.pts')
+
+
 
 
@@ -111 +122 @@
 
 print'add all geospatial objects'
-G = G1 + G2 + G3 + G_off + G_off1 + G_off2 + G_off3 + G_off4
+G = G1 + G2 + G3 + G_off + G_off1 + G_off2 + G_off3 + G_off4 
 
 print'clip combined geospatial object by bounding polygon'
 G_clip = G.clip_outside(project.poly_aoi1)
-G_all = G_clip + G_off5
+G_all = G_clip + G_off5 + G_off6 + G_off7
 G_clipped = G_all.clip(project.poly_all)