Changeset 5785 for anuga_work/production/perth/run_perth_250m.py

Timestamp:

Sep 25, 2008, 1:34:14 PM (16 years ago)

Author:

kristy

Message:

Correlating with Janes updates

File:

• anuga_work/production/perth/run_perth_250m.py (modified) (8 diffs)

anuga_work/production/perth/run_perth_250m.py

@@ -1 @@
-"""Script for running tsunami inundation scenario for Perth, WA, Australia.
+"""Script for running a tsunami inundation scenario for Perth, WA, Australia.
 
 The scenario is defined by a triangular mesh created from project_250m.polygon,
-the elevation data is combiled into a pts file through build_perth.py
-and a simulated tsunami is generated through an sts file from build_boundary.py
-
-Input: sts file (build_boundary.py for resepective event)
+the elevation data is compiled into a pts file through build_perth.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_perth.py)
-       information from project_250m file
+       information from project file
 Outputs: sww file stored in project_250m.output_run_time_dir 
 The export_results_all.py and get_timeseries.py is reliant
@@ -48 +48 @@
     
 # Application specific imports
-import project_250m  # Definition of file names and polygons
+import project  # Definition of file names and polygons
 numprocs = 1
 myid = 0
@@ -90 +90 @@
     N = len(urs_bounding_polygon)-1
 
-    #boundary tags refer to project_250m.landward 4 points equals 5 segments start at N
+    # boundary tags refer to project_250m.landward 4 points equals 5 segments start at N
     boundary_tags={'back': [N+1,N+2,N+3], 'side': [N,N+4], 'ocean': range(N)}
 
-    
     #--------------------------------------------------------------------------
     # Create the triangular mesh based on overall clipping polygon with a tagged
@@ -100 +99 @@
     #--------------------------------------------------------------------------
 
-    #IMPORTANT don't cache create_mesh_from_region and Domain(mesh....) as it
+    # IMPORTANT don't cache create_mesh_from_region and Domain(mesh....) as it
     # causes problems with the ability to cache set quantity which takes alot of times
         
@@ -132 +131 @@
     print 'Setup initial conditions'
 
-    #following sets the stage/water to be offcoast only
+    # sets the initial stage in the offcoast region only
     IC = Polygon_function( [(project_250m.poly_mainland, 0)], default = kwargs['tide'],
                              geo_reference = domain.geo_reference)
@@ -161 +160 @@
     domain.set_name(kwargs['scenario_name'])
     domain.set_datadir(kwargs['output_dir'])
-    domain.set_default_order(2) # Apply second order scheme
-    domain.set_minimum_storable_height(0.01) # Don't store anything less than 1cm
+    domain.set_default_order(2)                 # Apply second order scheme
+    domain.set_minimum_storable_height(0.01)    # Don't store anything less than 1cm
     domain.set_store_vertices_uniquely(False)
     domain.set_quantities_to_be_stored(['stage', 'xmomentum', 'ymomentum'])
@@ -206 +205 @@
         domain.write_boundary_statistics(tags = 'ocean')
 
+    # these outputs should be checked with the resultant inundation map
     x, y = domain.get_maximum_inundation_location()
     q = domain.get_maximum_inundation_elevation()
-
     print 'Maximum runup observed at (%.2f, %.2f) with elevation %.2f' %(x,y,q)
 
-    print 'That took %.2f seconds' %(time.time()-t0)
+    print 'Simulation took %.2f seconds' %(time.time()-t0)
 
     #kwargs 'completed' must be added to write the final parameters to file
@@ -232 +231 @@
     kwargs['alpha'] = project_250m.alpha
     kwargs['friction']=project_250m.friction
-       
-
+     
     run_model(**kwargs)