Changeset 9292 for trunk/anuga_work/development/gareth/tests/dam_break

Timestamp:

Aug 12, 2014, 12:10:51 PM (11 years ago)

Author:

davies

Message:

Adding local extrapolation and flux updating

Location:

trunk/anuga_work/development/gareth/tests/dam_break

Files:

• dam_break.py (modified) (3 diffs)
• plotme.py (modified) (1 diff)

trunk/anuga_work/development/gareth/tests/dam_break/dam_break.py

@@ -14 +14 @@
 from numpy import zeros, float
 from time import localtime, strftime, gmtime
-from bal_and import *
+#from bal_and import *
 #from anuga_tsunami import *
 
@@ -42 +42 @@
 points, vertices, boundary = anuga.rectangular_cross(int(L/dx), int(W/dy), L, W, (0.0, -W/2))
 
-#domain = anuga.Domain(points, vertices, boundary) 
-domain = Domain(points, vertices, boundary) 
+domain = anuga.Domain(points, vertices, boundary) 
+#domain = Domain(points, vertices, boundary) 
 
 domain.set_name(output_file)                
@@ -56 +56 @@
 print domain.get_timestepping_method()
 
+domain.set_flow_algorithm('DE0')
+#domain.set_local_extrapolation_and_flux_updating(nlevels=8)
 #domain.use_edge_limiter = True
 #domain.use_edge_limiter = False

trunk/anuga_work/development/gareth/tests/dam_break/plotme.py

@@ -4 +4 @@
 """
 from anuga.utilities import plot_utils as util
-from bal_and import plot_utils as util2
 from matplotlib import pyplot as pyplot
 
-p_st = util.get_output('dam_break_20120305_124724/dam_break.sww')
+p_st = util.get_output('dam_break_20140811_130209/dam_break.sww')
 p2_st=util.get_centroids(p_st)
 
 
-p_dev = util2.get_output('dam_break_20131204_145613/dam_break.sww', 0.001)
-p2_dev=util2.get_centroids(p_dev, velocity_extrapolation=True)
+p_dev = util.get_output('dam_break_20140811_130049/dam_break.sww', 0.001)
+p2_dev=util.get_centroids(p_dev, velocity_extrapolation=True)
 
 pyplot.clf()