Changeset 4212 for anuga_work/production/dampier_2006/build_dampier.py

Timestamp:

Feb 5, 2007, 11:28:49 AM (18 years ago)

Author:

nick

Message:

update to dampier

File:

• anuga_work/production/dampier_2006/build_dampier.py (modified) (6 diffs)

anuga_work/production/dampier_2006/build_dampier.py

@@ -61 +61 @@
 coast_in_dir_name = project.coast_in_dir_name
 offshore_in_dir_name = project.offshore_in_dir_name
-#offshore1_in_dir_name = project.offshore1_in_dir_name
-#offshore2_in_dir_name = project.offshore2_in_dir_name
+offshore1_in_dir_name = project.offshore1_in_dir_name
+offshore2_in_dir_name = project.offshore2_in_dir_name
 
 onshore_dir_name = project.onshore_dir_name
 coast_dir_name = project.coast_dir_name
 offshore_dir_name = project.offshore_dir_name
-#offshore1_dir_name = project.offshore1_dir_name
-#offshore2_dir_name = project.offshore2_dir_name
-
+offshore1_dir_name = project.offshore1_dir_name
+offshore2_dir_name = project.offshore2_dir_name
+'''
 # creates DEM from asc data
 print "creates DEMs from asc data"
@@ -75 +75 @@
                               basename_out=onshore_dir_name,
                               use_cache=True, verbose=True)
-'''
+
 convert_dem_from_ascii2netcdf(offshore1_in_dir_name,
                               basename_out=offshore1_dir_name,
@@ -82 +82 @@
                               basename_out=offshore2_dir_name,
                               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)
-'''
+
 #creates pts file for clipped DEMs
 dem2pts(offshore1_dir_name,
@@ -94 +94 @@
         basename_out=offshore2_dir_name,
         use_cache=True, verbose=True)
-'''
+
 print'create Geospatial onshore objects from topographies'
 #print'create Geospatial coastal objects from topographies'
@@ -101 +101 @@
 #print'create Geospatial offshore objects from topographies'
 G = Geospatial_data(file_name = onshore_dir_name + '.pts') +\
-     Geospatial_data(file_name = coast_in_dir_name + '.txt') +\
-     Geospatial_data(file_name = offshore_in_dir_name + '.txt')
+    Geospatial_data(file_name = coast_in_dir_name + '.txt') +\
+    Geospatial_data(file_name = offshore_in_dir_name + '.txt') +\
+    Geospatial_data(file_name = offshore1_dir_name + '.pts') +\
+    Geospatial_data(file_name = offshore2_dir_name + '.pts')
+     
 
 print'add all geospatial objects'
@@ -124 +127 @@
 G_small.export_points_file(project.combined_dir_name + '_small' + '.pts')
 G_other.export_points_file(project.combined_dir_name + '_other' + '.pts')
+'''
+print 'start reading:',project.combined_dir_name + '.pts'
+G = Geospatial_data(file_name = project.combined_dir_name + '.pts')
+print 'start export',project.combined_dir_name + '.txt'
+G.export_points_file(project.combined_dir_name + '.txt')
 
 
+
+
+