Changeset 6027 for anuga_work/production/geraldton/run_geraldton.py

Timestamp:

Nov 28, 2008, 11:48:45 AM (16 years ago)

Author:

kristy

Message:

general maintenance of scripts, added boundary to mesh in run_geraldton.py

File:

• anuga_work/production/geraldton/run_geraldton.py (modified) (4 diffs)

anuga_work/production/geraldton/run_geraldton.py

@@ -38 +38 @@
 from anuga.shallow_water.data_manager import export_grid, create_sts_boundary
 from anuga.pmesh.mesh_interface import create_mesh_from_regions
+from anuga.pmesh.mesh import importMeshFromFile
 from anuga.shallow_water.data_manager import start_screen_catcher, copy_code_files,store_parameters
 from anuga_parallel.parallel_abstraction import get_processor_name
@@ -44 +45 @@
 from anuga.fit_interpolate.benchmark_least_squares import mem_usage 
 from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon
-from anuga.geospatial_data.geospatial_data import find_optimal_smoothing_parameter
+from anuga.geospatial_data.geospatial_data import find_optimal_smoothing_parameter,Geospatial_data
 from polygon import Polygon_function
-    
+
+
 # Application specific imports
 import project  # Definition of file names and polygons
@@ -109 +111 @@
                          interior_regions=project.interior_regions,
                          filename=project.meshes_dir_name,
-                         use_cache=True,
+                         use_cache=False,
                          verbose=True)
+
+    m = importMeshFromFile(project.meshes_dir_name)
+    # load points from a points file
+    m.add_points_and_segments(Geospatial_data(project.barrier_dir_name))
+    m.add_points_and_segments(Geospatial_data(project.barrier_dir_name1))
+    m.generate_mesh()
+    m.export_mesh_file(project.meshes_dir_name)
    
     #-------------------------------------------------------------------------
@@ -224 +233 @@
     
     kwargs={}
+    kwargs['file_name']=project.dir_comment
     kwargs['finaltime']=project.finaltime
     kwargs['output_dir']=project.output_run_time_dir