Changeset 4838 for anuga_validation/convergence_study/convergence_structured.py

Timestamp:

Nov 21, 2007, 9:40:40 AM (16 years ago)

Author:

steve

Message:

Changed exception in fit_interpolation from ToFewPointsError? to
TooFewPointsError?

File:

• anuga_validation/convergence_study/convergence_structured.py (modified) (4 diffs)

anuga_validation/convergence_study/convergence_structured.py

@@ -23 +23 @@
 # Setup computational domain
 #------------------------------------------------------------------------------
-dx = 200.
+dx = 100.
 dy = dx
 L = 100000.
-W = 3000.
+W = dx
 
 # structured mesh
@@ -35 +35 @@
 #points, vertices, boundary = rectangular_cross(1000, 100, 20, 3) # Basic mesh
 domain = Domain(points, vertices, boundary) 
-domain.set_timestepping_method('rk2')
-domain.set_default_order(2)
+
+
 ## # unstructured mesh
 ## poly_domain = [[0,-W],[0,W],[L,W],[L,-W]]
@@ -50 +50 @@
 ## # Create domain
 ## domain = Domain(meshname, use_cache=True, verbose = True)
-domain.set_name('myexample6')                
+
+domain.set_timestepping_method('euler')
+domain.set_default_order(2)
+domain.set_name('myexample9')                
 domain.set_datadir('.')                     # Use current directory for output
+
+domain.beta_w      = 100.0
+domain.beta_w_dry  = 0.2
+domain.beta_uh     = 100.0
+domain.beta_uh_dry = 0.2
+domain.beta_vh     = 100.0
+domain.beta_vh_dry = 0.2
+domain.beta_h      = 100.0
 
 #------------------------------------------------------------------------------
@@ -83 +94 @@
 domain.set_boundary({'left': Bw, 'right': Bt, 'top': Br, 'bottom': Br})
 
-#from anuga.visualiser import RealtimeVisualiser
-#vis = RealtimeVisualiser(domain)
-#vis.render_quantity_height("elevation", dynamic=False)
-#vis.render_quantity_height("stage", dynamic=True)
-#vis.colour_height_quantity('stage', (0.0, 0.0, 0.8))
-#vis.start()
+## from anuga.visualiser import RealtimeVisualiser
+## vis = RealtimeVisualiser(domain)
+## vis.render_quantity_height("elevation", dynamic=False)
+## vis.render_quantity_height("stage", dynamic=True, zScale=100.0)
+## #vis.colour_height_quantity('stage', (0.0, 0.0, 0.8))
+## vis.colour_height_quantity('stage', (lambda q:q['stage'], -1, 1))
+## vis.start()
 
 #------------------------------------------------------------------------------