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

Timestamp:

Jun 20, 2007, 6:34:21 PM (18 years ago)

Author:

ole

Message:

Got the shark bay scenario to build

File:

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

anuga_work/production/shark_bay_2007/build_shark_bay.py

@@ -3 +3 @@
 Source data such as elevation and boundary data is assumed to be available in
 directories specified by project_urs.py
-The output sww file is stored in project_urs.output_time_dir 
+The output sww file is stored in project.output_time_dir 
 
 The scenario is defined by a triangular mesh created from project_urs.polygon,
@@ -31 +31 @@
 from anuga.pmesh.mesh_interface import create_mesh_from_regions
 from anuga.geospatial_data.geospatial_data import Geospatial_data
-from anuga.abstract_2d_finite_volumes.util import start_screen_catcher, copy_code_files
+from anuga.shallow_water.data_manager import start_screen_catcher
+from anuga.shallow_water.data_manager import copy_code_files
+
 from anuga_parallel.parallel_abstraction import get_processor_name
 
 # Application specific imports
-import project_urs   # Definition of file names and polygons
+import project   # Definition of file names and polygons
+
 
 #------------------------------------------------------------------------------
@@ -42 +45 @@
 #------------------------------------------------------------------------------
 
-copy_code_files(project_urs.output_build_time_dir,__file__, 
-               dirname(project_urs.__file__)+sep+ project_urs.__name__+'.py' )
+copy_code_files(project.output_build_time_dir,__file__, 
+               dirname(project.__file__)+sep+ project.__name__+'.py' )
 
-start_screen_catcher(project_urs.output_build_time_dir)
+start_screen_catcher(project.output_build_time_dir)
 
-print "Processor Name:",get_processor_name()
-print 'USER: ', project_urs.user
+print 'Processor Name:', get_processor_name()
+print 'User: ', project.user
+
 
 #-------------------------------------------------------------------------------
@@ -57 +61 @@
 # Fine pts file to be clipped to area of interest
 #-------------------------------------------------------------------------------
-print"project_urs.combined_dir_name",project_urs.combined_dir_name
-'''
-# topography directory filenames
-onshore_in_dir_name = project_urs.onshore_in_dir_name
-coast_in_dir_name = project_urs.coast_in_dir_name
-offshore_in_dir_name = project_urs.offshore_in_dir_name
-offshore1_in_dir_name = project_urs.offshore1_in_dir_name
-offshore2_in_dir_name = project_urs.offshore2_in_dir_name
 
-onshore_dir_name = project_urs.onshore_dir_name
-coast_dir_name = project_urs.coast_dir_name
-offshore_dir_name = project_urs.offshore_dir_name
-offshore1_dir_name = project_urs.offshore1_dir_name
-offshore2_dir_name = project_urs.offshore2_dir_name
+print 'project.combined_dir_name', project.combined_dir_name
 
-# creates DEM from asc data
-print "creates DEMs from asc data"
-convert_dem_from_ascii2netcdf(onshore_in_dir_name,
-                              basename_out=onshore_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)
 
-print'create Geospatial onshore 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 = offshore1_dir_name + '.pts') +\
-#    Geospatial_data(file_name = offshore2_dir_name + '.pts')
-     
+geospatial_data = None
+# create DEMs from asc data
+print 'creating geospatial data objects from asc data (via dem and pts formats)'
+for filename in project.ascii_grid_filenames:
+    convert_dem_from_ascii2netcdf(filename,
+                                  basename_out=filename,
+                                  use_cache=True, verbose=True)
+    dem2pts(filename, use_cache=True, verbose=True)
 
-print'clip combined geospatial object by bounding polygon'
-G_clipped = G.clip(project_urs.poly_all)
+    geospatial_data += Geospatial_data(file_name = filename + '.pts', verbose=True)
 
-print'export combined DEM file'
-if access(project_urs.topographies_dir,F_OK) == 0:
-    mkdir (project_urs.topographies_dir)
-G.export_points_file(project_urs.combined_dir_name + '.txt')
 
-print'split combined data set'
-G_small, G_other = G_clipped.split(0.1, 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)
 
-print'export split DEM file'
-G_small.export_points_file(project_urs.combined_dir_name + '_small' + '.txt')
-#G_other.export_points_file(project_urs.combined_dir_name + '_other' + '.pts')
 
-'''
-print 'start reading:',project_urs.combined_dir_name + '.txt'
-G = Geospatial_data(file_name = (project_urs.combined_dir_name + '.txt'))
-G_clipped = G.clip(project_urs.poly_topo_clip)
+print 'clip combined geospatial object by bounding polygon'
+G = geospatial_data.clip(project.poly_all)
 
-print 'start split'
-#G_smallest, G_other = G.split(0.1,True)
 
-print 'start export',project_urs.combined_smallest_dir_name + '.txt'
-#G.export_points_file(project_urs.combined_smaller_dir_name + '.txt')
-G_clipped.export_points_file(project_urs.combined_smallest_dir_name + '.txt')
+print 'export combined geospatial data'
+if access(project.topographies_dir, F_OK) == 0:
+    mkdir (project.topographies_dir)
+G.export_points_file(project.combined_dir_name + '.txt')
 
+
+#print 'split combined data set'
+#G_small, _ = G.split(0.1, True)
+#
+#print 'export sub sampled DEM file'
+#G_small.export_points_file(project.combined_dir_name + '_small' + '.txt')
+
+
+ 
 #-------------------------------------------------------------------------
 # Convert URS to SWW file for boundary conditions 
 #-------------------------------------------------------------------------
-print 'starting to create boundary conditions'
-#boundaries_in_dir_name = project_urs.boundaries_in_dir_name
+print 'converting boundary conditions to sww format'
+#boundaries_in_dir_name = project.boundaries_in_dir_name
 
-#print 'minlat=project_urs.north_boundary, maxlat=project_urs.south_boundary',project_urs.north_boundary, project_urs.south_boundary
-#print 'minlon= project_urs.west_boundary, maxlon=project_urs.east_boundary',project_urs.west_boundary, project_urs.east_boundary
-'''       
+#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 'boundaries_in_dir_name',project_urs.boundaries_in_dir_name
-print 'project.boundaries_dir_name',project_urs.boundaries_dir_name
+print 'boundaries_dir', project.boundaries_dir
+print 'project.boundaries_name', project.boundaries_name
 
-urs_ungridded2sww(project_urs.boundaries_in_dir_name, project_urs.boundaries_dir_name,
+urs_ungridded2sww(project.boundaries_name,
                   verbose=True, mint=5000, maxt=35000, zscale=1)
-'''
 
 
+print 'Finished building the %s scenario' %project.scenario_name