Changeset 718 for inundation/ga/storm_surge/examples/run_tsh.py

Timestamp:

Dec 22, 2004, 2:58:16 PM (20 years ago)

Author:

duncan

Message:

adding regional modifiers

File:

• inundation/ga/storm_surge/examples/run_tsh.py (modified) (3 diffs)

inundation/ga/storm_surge/examples/run_tsh.py

@@ -18 +18 @@
 
 from shallow_water import Domain, Reflective_boundary, Dirichlet_boundary,\
-     Transmissive_boundary, Time_boundary, Constant_height, Weir
-
+     Transmissive_boundary, Time_boundary, Set_region, Add_value_to_region
 from mesh_factory import rectangular
 from pmesh2domain import pmesh_to_domain, pmesh_to_domain_instance
@@ -80 +79 @@
     #
     print 'Boundaries'
-    Br = Reflective_boundary(domain)
+    reflective = Reflective_boundary(domain)
     Bt = Transmissive_boundary(domain)
 
     #Constant inflow
-    Bd = Dirichlet_boundary(array([inflow_stage, 0.0, 0.0]))
+    Bd = Dirichlet_boundary(array([7, 0.0, 0.0]))
 
     #Time dependent inflow
@@ -105 +104 @@
     tags['internal'] = None
     tags['levee'] = None
-    tags['0'] = Br
-    tags['wall'] = Br 
-    tags['external'] = Br 
-    tags['open'] = Bd 
+    tags['0'] = reflective
+    tags['wall'] = reflective 
+    tags['external'] = reflective 
+    tags['open'] = Bd  
+    tags['opening'] = None 
 
     domain.set_boundary(tags)
 
-
+    # region tags
+    #domain.set_region(Set_region('mound', 'elevation', 100, location='unique vertices'))
+    domain.set_region(Add_value_to_region('mound', 'elevation', 100, location='unique vertices', initial_quantity='elevation'))
+    domain.set_region(Add_value_to_region('reservoir', 'level', 100.0,initial_quantity='elevation'))
+    
     #print domain.quantities['elevation'].vertex_values
     #print domain.quantities['level'].vertex_values