Changeset 7657 for anuga_work/production/bunbury_storm_surge_2009/run_model.py

Timestamp:

Mar 9, 2010, 1:54:17 PM (13 years ago)

Author:

fountain

Message:

updates to bunbury model to refine elevation around leschenault inlet and estuary and include storm gates

File:

• anuga_work/production/bunbury_storm_surge_2009/run_model.py (modified) (2 diffs)

anuga_work/production/bunbury_storm_surge_2009/run_model.py

@@ -113 +113 @@
 print domain.statistics()
 
+domain.set_starttime(project.starttime)
 domain.set_name(project.scenario_name)
 domain.set_datadir(project.output_run) 
@@ -137 +138 @@
                     verbose=True,
                     alpha=project.alpha)
+
+#--------------------------
+# Add storm surge gate
+#--------------------------
+storm_gate_polygon, storm_gate_height = csv2building_polygons(project.storm_gate, floor_height=4)
+
+gate = []
+for key in storm_gate_polygon:
+    poly = storm_gate_polygon[key]
+    elev = storm_gate_height[key]
+    gate.append((poly, elev))
+
+domain.add_quantity('elevation', 
+                    Polygon_function(gate, geo_reference=domain.geo_reference),
+                    use_cache=True,
+                    verbose=True)
 
 #-------------------------------------------------------------------------------