Changeset 3408

Timestamp:

Jul 24, 2006, 11:49:21 AM (19 years ago)

Author:

duncan

Message:

setting up a new simulation

Location:

development/dam_2006

Files:

• create_mesh.py (modified) (2 diffs)
• run_dam.py (modified) (2 diffs)

development/dam_2006/create_mesh.py

@@ -13 +13 @@
     #Basic geometry
     
-    xleft   = -5.4
-    xright  = .4
+    xleft   = -2.25
+    xright  = .75
     ybottom = 0
     ytop    = 0.4
@@ -56 +56 @@
         
     m.addVertsSegs(dict)
-
-    dam = m.add_region(0.2,0.2)
+    
+    dam = m.add_region(0.2,0.2) # this is the location of the region.
     dam.setTag("dam")
-    #m.generate_mesh(maximum_triangle_area=0.01)
-    m.generate_mesh(maximum_triangle_area=0.0001)
+    m.generate_mesh(maximum_triangle_area=0.01)
+    #m.generate_mesh(maximum_triangle_area=0.0001)
 
     m.export_mesh_file(project.mesh_filename)

development/dam_2006/run_dam.py

@@ -75 +75 @@
 
 tide = 0.0
+slope = -0.05
 
 def elevation_tilt(x, y):
-    return -x*0.15
+    return x*slope
         
 domain.set_quantity('stage', elevation_tilt)
@@ -85 +86 @@
 print 'Available boundary tags', domain.get_boundary_tags()
 
-domain.set_region(Set_region('dam','stage',0.4,location = 'unique vertices'))
+domain.set_region(Set_region('dam','stage',0.2,location = 'unique vertices'))
 
 Br = Reflective_boundary(domain)