Changeset 6282

Timestamp:

Feb 5, 2009, 3:27:30 PM (16 years ago)

Author:

jgriffin

Message:

 

File:

• anuga_work/production/busselton/standardised_version/project.py (modified) (6 diffs)

anuga_work/production/busselton/standardised_version/project.py

@@ -40 +40 @@
 # Model specific parameters. One or all can be changed each time the
 # run_scenario script is executed
-tide = 0              #0.6
-#event_number = 27255 # Java 9.3 worst case for Perth
-#event_number = 68693 # Sumatra 9.2
-event_number = 27283  # Java 9.3 original
+tide = 0                #0.6
+#event_number = 27255   # Java 9.3 worst case for Perth
+#event_number = 68693   # Sumatra 9.2
+event_number = 27283    # Java 9.3 original
 alpha = 0.1             # smoothing parameter for mesh
 friction=0.01           # manning's friction coefficient
@@ -137 +137 @@
 # Check the run_busselton.py for boundary_tags
 
-# thinned ordering file from Hazard Map (geographic)
+# Thinned ordering file from Hazard Map (geographic)
 # Format is index,latitude,longitude (with header)
 urs_order_filename = 'thinned_boundary_ordering_extend.csv'
 
-#landward bounding points
+# Landward bounding points
 # Format easting,northing (no header)
 landward_boundary_filename = 'landward_boundary_extend.csv'
@@ -209 +209 @@
 output_run_time = join(output_run, scenario_name) 
 
-# The absolute pathname for the gauges file, used for get_timeseries.py
+# The absolute pathname for the gauges file
+# Used for get_timeseries.py
 gauges = join(gauges_folder, gauges_filename)       
 
-# The absolute pathname for the building file,
+# The absolute pathname for the building file
 # Used for run_building_inundation.py
 building_exposure = join(gauges_folder, building_exposure_filename)
@@ -220 +221 @@
 #------------------------------------------------------------------------------
 
+# Create list of land polygons with initial conditions
 land_initial_conditions = []
 for filename, MSL in land_initial_conditions_filename:
@@ -225 +227 @@
     land_initial_conditions.append([filename, MSL])
 
+# Create list of interior polygons with scaling factor
+interior_regions = []
+for filename, maxarea in interior_regions_data:
+    polygon = read_polygon(join(polygons_folder, filename))
+    interior_regions.append([polygon, maxarea*scale_factor])
+
 # Initial bounding polygon for data clipping 
 bounding_polygon = read_polygon(join(polygons_folder,
@@ -230 +238 @@
 bounding_maxarea = 100000*scale_factor
 
-interior_regions = []
-for filename, maxarea in interior_regions_data:
-    polygon = read_polygon(join(polygons_folder, filename))
-    interior_regions.append([polygon, maxarea*scale_factor])
-                     
+# Estimate the number of triangles                     
 trigs_min = number_mesh_triangles(interior_regions,
                                   bounding_polygon, bounding_maxarea)