Changeset 5149 for anuga_work/production/busselton/build_busselton.py

Timestamp:

Mar 10, 2008, 3:59:09 PM (16 years ago)

Author:

nick

Message:

update busselton

File:

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

anuga_work/production/busselton/build_busselton.py

@@ -27 +27 @@
 from anuga.shallow_water.data_manager import convert_dem_from_ascii2netcdf, dem2pts
 from anuga.geospatial_data.geospatial_data import *
-from anuga.abstract_2d_finite_volumes.util import start_screen_catcher, copy_code_files
+from anuga.shallow_water.data_manager import start_screen_catcher, copy_code_files,store_parameters
 
 # Application specific imports
@@ -39 +39 @@
 start_screen_catcher(project.output_build_time_dir)
 
-print 'time stamp: ',project.gtime
+print 'time stamp: ',project.time
 print 'USER: ', project.user
 
@@ -59 +59 @@
 onshore_in_dir_name = project.onshore_in_dir_name
 coast_in_dir_name = project.coast_in_dir_name
+coast_in_dir_name1 = project.coast_in_dir_name1
 #island_in_dir_name = project.island_in_dir_name
 offshore_in_dir_name = project.offshore_in_dir_name
@@ -64 +65 @@
 onshore_dir_name = project.onshore_dir_name
 coast_dir_name = project.coast_dir_name
+coast_dir_name1 = project.coast_dir_name1
 #island_dir_name = project.island_dir_name
 offshore_dir_name = project.offshore_dir_name
@@ -86 +88 @@
 
 print'create Geospatial data1 objects from topographies'
-G1 = Geospatial_data(file_name = onshore_dir_name + '.pts')
-G2 = Geospatial_data(file_name = coast_in_dir_name + '.xya')
-#G3 = Geospatial_data(file_name = island_dir_name + '.pts')
-G_off = Geospatial_data(file_name = offshore_in_dir_name + '.xya')
+G1 = Geospatial_data(file_name = onshore_dir_name + '.pts',verbose=True)
+G2 = Geospatial_data(file_name = coast_in_dir_name + '.txt',verbose=True)
+G3 = Geospatial_data(file_name = coast_in_dir_name1 + '.txt',verbose=True)
+#G3 = Geospatial_data(file_name = island_dir_name + '.pts',verbose=True)
+G_off = Geospatial_data(file_name = offshore_in_dir_name + '.txt',verbose=True)
 
 print'add all geospatial objects'
-G = G1 + G2 + G_off 
+G = G1 + G2 + G3 + G_off 
 
 print'clip combined geospatial object by bounding polygon'
@@ -102 +105 @@
 if access(project.topographies_dir,F_OK) == 0:
     mkdir (project.topographies_dir)
-G.export_points_file(project.combined_dir_name + '.xya')
+G.export_points_file(project.combined_dir_name + '.txt')
 #G_clipped.export_points_file(project.combined_dir_name + '.xya')
 
+print'project.combined_dir_name + .txt',project.combined_dir_name + '.txt'
+#G_all=Geospatial_data(file_name = project.combined_dir_name + '.txt')
+print'split'
+G_all_1, = G.split(.10)
+print'export 1'
+G_all_1.export_points_file(project.combined_dir_name+'_small' + '.txt')
+#print'export 2'
+#G_all_2.export_points_file(project.combined_dir_name+'_other1' + '.xya')
+
 '''
-print'project.combined_dir_name + 1.xya',project.combined_dir_name + '1.xya'
-G_all=Geospatial_data(file_name = project.combined_dir_name + '1.xya')
-print'split'
-G_all_1, G_all_2 = G_all.split(.10)
-print'export 1'
-G_all_1.export_points_file(project.combined_dir_name+'_small1' + '.xya')
-print'export 2'
-G_all_2.export_points_file(project.combined_dir_name+'_other1' + '.xya')
-
-
 #-------------------------------------------------------------------------
 # Convert URS to SWW file for boundary conditions