Changeset 4429 for anuga_work/production/broome_2006/build_broome.py

Timestamp:

May 15, 2007, 4:13:58 PM (17 years ago)

Author:

nick

Message:

update broome

File:

• anuga_work/production/broome_2006/build_broome.py (modified) (5 diffs)

anuga_work/production/broome_2006/build_broome.py

@@ -83 +83 @@
 dem2pts(offshore_dir_name1, use_cache=True, verbose=True)
 dem2pts(offshore_dir_name2, use_cache=True, verbose=True)
-'''
+
 print'create Geospatial data1 objects from topographies',onshore_dir_name + '.pts'
 G1 = Geospatial_data(file_name = onshore_dir_name + '.pts')
@@ -110 +110 @@
 print'export',project.combined_dir_name+'_unclipped' + '.txt'
 G.export_points_file(project.combined_dir_name+'_unclipped' + '.txt')
-'''
+
 print'split'
 G_small, G_other = G_clipped.split(.10,verbose=True)
@@ -120 +120 @@
 G_small.export_points_file(project.combined_small_dir_name + '.txt')
 #G_clipped.export_points_file(project.combined_dir_name + '.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')
@@ -132 +132 @@
 G_all_2.export_points_file(project.combined_dir_name+'_other1' + '.xya')
 
-
+'''
 #-------------------------------------------------------------------------
 # Convert URS to SWW file for boundary conditions 
 #-------------------------------------------------------------------------
-'''
+print 'starting to create boundary conditions'
+boundaries_in_dir_name = project.boundaries_in_dir_name
+
+from anuga.shallow_water.data_manager import urs2sww, urs_ungridded2sww
+
+print 'boundaries_in_dir_name',boundaries_in_dir_name
+print 'project.boundaries_dir_name',project.boundaries_dir_name
+
+urs_ungridded2sww(boundaries_in_dir_name, project.boundaries_dir_name,
+                  verbose=True, mint=4000, maxt=35000, zscale=1)
 
 
@@ -143 +152 @@
 
 
-