Changeset 6617 for anuga_work/production/australia_ph2/darwin/project.py

Timestamp:

Mar 25, 2009, 4:26:41 PM (16 years ago)

Author:

myall

Message:

putting in the option for a internal polygon with same triangle area (to test discretisation issues) in project scripts.
adding required scripts for broome

File:

• anuga_work/production/australia_ph2/darwin/project.py (modified) (4 diffs)

anuga_work/production/australia_ph2/darwin/project.py

@@ -27 +27 @@
 # One or all can be changed each time the run_model script is executed
 tide = 0                # difference between MSL and HAT
-zone = 52               # specify zone of model
+central_meridian = 129.0 # for projection (necessary for nonstandard)
 event_number = 70906    # the event number or the mux file name
 # contributing events: 17898 (Flores), 64977 (Seram), 70906 (Sumba_normal),
@@ -48 +48 @@
 # for 1860 - point to north of Darwin, large wave from event off Timor
 # event_number 71065    #(timor - mag 7.4, prob 1e-4)
+
+internal_polygon = False
 
 #-------------------------------------------------------------------------------
@@ -60 +62 @@
 #-------------------------------------------------------------------------------
 
-output_comment = [setup, tide, event_number, index, wave]
+if internal_polygon:
+    internal_poly_comment = 'internal'
+else:
+    internal_poly_comment = None
+
+output_comment = [setup, tide, event_number, index, wave, internal_poly_comment]
 
 #-------------------------------------------------------------------------------
@@ -91 +98 @@
 #interior_regions_data = []
 interior_regions_data = [['poly_island.csv',200000000]]  # large triangles on island
+
+# add an internal polygon to force different mesh generation
+# used to test for discretisation error when building elevation
+# make sure file is in same folder as interior regions and bouding polygon;
+# format is same (2 column .csv; easting, northing; no header)
+if internal_polygon:
+    interior_regions_data.append(['internal_polygon.csv',
+                                  bounding_polygon_maxarea])
 
 # LAND - used to set the initial stage/water to be offcoast only