Changeset 2155

Timestamp:

Dec 19, 2005, 4:34:15 PM (19 years ago)

Author:

ole

Message:

Simplified example

File:

• inundation/examples/bedslope.py (modified) (3 diffs)

inundation/examples/bedslope.py

@@ -19 +19 @@
 domain = Domain(points, vertices, boundary)
 domain.smooth = False
-domain.visualise = True
+domain.visualise = False
 domain.store = True
 domain.filename = 'bedslope'
@@ -25 +25 @@
 
 #######################
-#Bed-slope and friction
+# Initial conditions
 domain.set_quantity('elevation', lambda x,y: -x/3)
 domain.set_quantity('friction', 0.1)
+
+h = 0.05  # Constant depth
+domain.set_quantity('stage', expression = 'elevation + %f' %h)
 
 
@@ -39 +42 @@
                    f=lambda t: [(0.1*sin(t*2*pi)), 0.0, 0.0])
 
+
 domain.set_boundary({'left': Bd, 'right': Br, 'top': Br, 'bottom': Br})
+#domain.set_boundary({'left': Bw, 'right': Br, 'top': Br, 'bottom': Br})
 
 
 ######################
-#Initial condition
-h = 0.05
-#h = 0.0
-elevation = domain.quantities['elevation'].vertex_values
-domain.set_quantity('stage', elevation + h)
+#Evolution
 
 domain.check_integrity()
 
-  
-######################
-#Evolution
-for t in domain.evolve(yieldstep = 0.1, finaltime = 3.0):
+for t in domain.evolve(yieldstep = 0.1, finaltime = 4.0):
     domain.write_time()
 
-for t in domain.evolve(yieldstep = 0.01, finaltime = 4.0):
-    domain.write_time()    
 
+