Changeset 4351 for anuga_work/development/friction_UA_flume_2006/create_mesh.py

Timestamp:

Apr 3, 2007, 4:22:40 PM (18 years ago)

Author:

duncan

Message:

checking in a scenario

File:

• anuga_work/development/friction_UA_flume_2006/create_mesh.py (modified) (5 diffs)

anuga_work/development/friction_UA_flume_2006/create_mesh.py

@@ -7 +7 @@
 
 
-def generate(mesh_filename, gate_position, is_coarse = False):
+def generate(mesh_filename, is_coarse = False):
     """
     Generate mesh for University of Aberbeen dam break flume.
@@ -16 +16 @@
     #Basic geometry
     
-    xright  = gate_position 
+    xright  = 19.0 #  web site says a 20 m long flume
+                   #(http://www.eng.abdn.ac.uk/envhrg/facilities/swashrig.hti)
+                   # Therefore 19 +1 (for reservoir)= 20
     ybottom = 0
     ytop    = 0.45
     xdam = 0.0
-    xleft = gate_position  - 13.0
+    xleft = -1.0
+    xslope = 4.0 # note, get in gen meshparrallel run dam as well!
 
     #Outline 
@@ -32 +35 @@
     point_dam_bottom = [xdam, ybottom]    
 
+    # slope seperation (middle)
+    point_slope_top = [xslope, ytop]
+    point_slope_bottom = [xslope, ybottom]    
+
     m = Mesh()
 
@@ -38 +45 @@
               point_nw,
               point_dam_top,
+              point_slope_top,
               point_ne, 
-              point_se, 
+              point_se,
+              point_slope_bottom,
               point_dam_bottom
               ]
     
     segments = [[0,1], [1,2], [2,3],
-                        [3,4 ],[4,5], [5,0],  #The outer border
-                        [2,5]]         #dam Separator
+                [3,4 ],[4,5], [5,6],[6,7],[7,0],  #The outer border
+                [2,7],    #dam Separator
+                [3,6]]        # slope separator
     
-    segment_tags = {'wall':[0,1,2,4,5],'edge':[3]} # '':[6]
+    segment_tags = {'wall':[0,1,2,3,5,6,7],'edge':[4]} # '':[6]
         
     m.add_points_and_segments(points, segments, segment_tags)
     
     dam = m.add_region(-0.0000001,(ytop - ybottom)/2)
-    # this is the location of the region.
+    # this is the location of the reservoir region.
     dam.setTag("dam")
+    
+    slope = m.add_region(xslope + 0.0000001,(ytop - ybottom)/2)
+    # this is the location of the slope region.
+    slope.setTag("slope")
     
     if is_coarse:
@@ -65 +79 @@
 #-------------------------------------------------------------
 if __name__ == "__main__":
-    generate("aa.tsh", 0.75, is_course = True)
+    generate("aa.tsh", is_coarse = True)