Changeset 5646

Timestamp:

Aug 13, 2008, 8:22:25 AM (17 years ago)

Author:

kristy

Message:

updated with more data within the 500m2 mesh area.

File:

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

anuga_work/production/busselton/build_busselton.py

@@ -62 +62 @@
 offshore_in_dir_name2 = project.offshore_in_dir_name2
 offshore_in_dir_name3 = project.offshore_in_dir_name3
+offshore_in_dir_name4 = project.offshore_in_dir_name4
+offshore_in_dir_name5 = project.offshore_in_dir_name5
 
 onshore_dir_name = project.onshore_dir_name
@@ -70 +72 @@
 offshore_dir_name2 = project.offshore_dir_name2
 offshore_dir_name3 = project.offshore_dir_name3
+offshore_dir_name4 = project.offshore_dir_name4
+offshore_dir_name5 = project.offshore_dir_name5
 
 # creates DEM from asc data
@@ -96 +100 @@
 G_off2 = Geospatial_data(file_name = offshore_in_dir_name2 + '.txt',verbose=True)
 G_off3 = Geospatial_data(file_name = offshore_in_dir_name3 + '.txt',verbose=True)
-                         
+G_off4 = Geospatial_data(file_name = offshore_in_dir_name4 + '.txt',verbose=True)
+G_off5 = Geospatial_data(file_name = offshore_in_dir_name5 + '.txt',verbose=True)
 print'add all geospatial objects'
-G = G1 + G2 + G3 + G_off + G_off1 + G_off2 + G_off3 
+G = G1 + G2 + G3 + G_off + G_off1 + G_off2 + G_off3 + G_off4 + G_off5 
 
 print'clip combined geospatial object by bounding polygon'
@@ -110 +115 @@
 print'project.combined_dir_name + .txt',project.combined_dir_name + '.txt'
 
-#-------------------------------------------------------------------------
-# Convert URS to SWW file for boundary conditions 
-#-------------------------------------------------------------------------
-print 'starting to create boundary conditions'
-from anuga.shallow_water.data_manager import urs2sww, urs_ungridded2sww
-
-boundaries_in_dir_name = project.boundaries_in_dir_name
-print 'boundaries_in_dir_name',project.boundaries_in_dir_name
-
-
-#import sys; sys.exit()
-
-urs_ungridded2sww(project.boundaries_in_dir_name, project.boundaries_in_dir_name,
-                  verbose=True, mint=4000, maxt=80000, zscale=1)
-
+###-------------------------------------------------------------------------
+### Convert URS to SWW file for boundary conditions 
+###-------------------------------------------------------------------------
+##print 'starting to create boundary conditions'
+##from anuga.shallow_water.data_manager import urs2sww, urs_ungridded2sww
+##
+##boundaries_in_dir_name = project.boundaries_in_dir_name
+##print 'boundaries_in_dir_name',project.boundaries_in_dir_name
+##
+##
+###import sys; sys.exit()
+##
+##urs_ungridded2sww(project.boundaries_in_dir_name, project.boundaries_in_dir_name,
+##                  verbose=True, mint=4000, maxt=80000, zscale=1)
+##