Changeset 4570 for anuga_work/production/shark_bay_2007/build_shark_bay.py

Timestamp:

Jun 28, 2007, 5:54:00 PM (18 years ago)

Author:

ole

Message:

Setup shark bay initial runs

File:

• anuga_work/production/shark_bay_2007/build_shark_bay.py (modified) (5 diffs)

anuga_work/production/shark_bay_2007/build_shark_bay.py

@@ -28 +28 @@
 from anuga.shallow_water import File_boundary
 from anuga.shallow_water import Reflective_boundary
-from anuga.shallow_water.data_manager import convert_dem_from_ascii2netcdf, dem2pts
+from anuga.shallow_water.data_manager import convert_dem_from_ascii2netcdf
+from anuga.shallow_water.data_manager import dem2pts
 from anuga.pmesh.mesh_interface import create_mesh_from_regions
 from anuga.geospatial_data.geospatial_data import Geospatial_data
 from anuga.shallow_water.data_manager import start_screen_catcher
 from anuga.shallow_water.data_manager import copy_code_files
-
+from anuga.shallow_water.data_manager import urs_ungridded2sww
 from anuga_parallel.parallel_abstraction import get_processor_name
 
@@ -46 +47 @@
 
 copy_code_files(project.output_build_time_dir,__file__, 
-               dirname(project.__file__)+sep+ project.__name__+'.py' )
+               dirname(project.__file__)+sep+ project.__name__+'.py' )#
 
 start_screen_catcher(project.output_build_time_dir)
@@ -67 +68 @@
 geospatial_data = None
 # create DEMs from asc data
+
+# FIXME (Ole): Clip data by interior regions if possible
+
 print 'creating geospatial data objects from asc data (via dem and pts formats)'
 for filename in project.ascii_grid_filenames:
@@ -74 +78 @@
     dem2pts(filename, use_cache=True, verbose=True)
 
-    geospatial_data += Geospatial_data(file_name = filename + '.pts', verbose=True)
+    geospatial_data += Geospatial_data(file_name=filename + '.pts',
+                                       verbose=True)
 
 
 print 'creating geospatial data objects from txt data'
 for filename in project.point_filenames:
-    geospatial_data += Geospatial_data(file_name = filename + '.txt', verbose=True)
+    geospatial_data += Geospatial_data(file_name=filename + '.txt',
+                                       verbose=True)
 
 
 print 'clip combined geospatial object by bounding polygon'
-G = geospatial_data.clip(project.poly_all)
+G = geospatial_data.clip(project.bounding_polygon)
 
 
@@ -104 +110 @@
 #-------------------------------------------------------------------------
 print 'converting boundary conditions to sww format'
-#boundaries_in_dir_name = project.boundaries_in_dir_name
 
-#print 'minlat=project.north_boundary, maxlat=project.south_boundary',project.north_boundary, project.south_boundary
-#print 'minlon= project.west_boundary, maxlon=project.east_boundary',project.west_boundary, project.east_boundary
-       
-from anuga.shallow_water.data_manager import urs_ungridded2sww
+print 'boundary_dir', project.boundary_dir
+print 'project.boundary_name', project.boundary_name
 
-print 'boundaries_dir', project.boundaries_dir
-print 'project.boundaries_name', project.boundaries_name
-
-urs_ungridded2sww(project.boundaries_name,
-                  verbose=True, mint=5000, maxt=35000, zscale=1)
+urs_ungridded2sww(project.boundary_name, verbose=True);
+                  #verbose=True, mint=5000, maxt=35000, zscale=1)