Changeset 6298 for anuga_work/production/australia_ph2/sydney/run_model.py

Timestamp:

Feb 9, 2009, 4:26:06 PM (16 years ago)

Author:

jgriffin

Message:

scripts modified to include file_length function

File:

• anuga_work/production/australia_ph2/sydney/run_model.py (modified) (4 diffs)

anuga_work/production/australia_ph2/sydney/run_model.py

@@ -1 @@
-"""Run a tsunami inundation scenario for Busselton, WA, Australia.
+"""Run a tsunami inundation scenario for sydney region, NSW, Australia.
 
 The scenario is defined by a triangular mesh created from project.polygon, the
-elevation data is compiled into a pts file through build_busselton.py and a
+elevation data is compiled into a pts file through build_sydney.py and a
 simulated tsunami is generated through an sts file from build_boundary.py.
 
 Input: sts file (build_boundary.py for respective event)
-       pts file (build_busselton.py)
+       pts file (build_sydney.py)
        information from project file
 Outputs: sww file stored in project.output_run_time_dir 
@@ -32 +32 @@
 from anuga.interface import create_sts_boundary
 from anuga.interface import csv2building_polygons
+from file_length import file_length
 
 from anuga.shallow_water.data_manager import start_screen_catcher
@@ -70 +71 @@
 # Number of boundary segments
 N = len(event_sts)-1
+# Number of landward_boundary points
+M = file_length(project.landward_boundary)
+print 'M',M
 
 # Boundary tags refer to project.landward_boundary
 # 4 points equals 5 segments start at N
-boundary_tags={'back': [N+1,N+2,N+3,N+4, N+5],
-               'side': [N,N+6],
+boundary_tags={'back': range(N+1,N+M),
+               'side': [N,N+M],
                'ocean': range(N)}
 
@@ -80 +84 @@
 domain = create_domain_from_regions(bounding_polygon_sts,
                                     boundary_tags=boundary_tags,
-                                    maximum_triangle_area=project.res_poly_all,
+                                    maximum_triangle_area=project.bounding_maxarea,
                                     interior_regions=project.interior_regions,
                                     mesh_filename=project.meshes,