Changeset 8668

Timestamp:

Dec 8, 2012, 8:26:13 PM (12 years ago)

Author:

mungkasi

Message:

Adding radial_dam_break_dry and modifying the wet one.

Location:

trunk/anuga_core/validation_tests/Tests/Other_References

Files:

• radial_dam_break_dry (added)
• radial_dam_break_dry/Ver_numerical_2.000000.csv (added)
• radial_dam_break_dry/plot_results.py (added)
• radial_dam_break_dry/produce_results.py (added)
• radial_dam_break_dry/radial_dam_break.py (added)
• radial_dam_break_dry/results.tex (added)
• radial_dam_break_wet/radial_dam_break.py (modified) (3 diffs)

trunk/anuga_core/validation_tests/Tests/Other_References/radial_dam_break_wet/radial_dam_break.py

@@ -44 +44 @@
 
 #------------------------------------------------------------------------------
-# Setup Algorithm
+# Setup Algorithm, either using command line arguments
+# or override manually yourself
 #------------------------------------------------------------------------------
-domain.set_timestepping_method('rk2')
-domain.set_default_order(2)
-print domain.get_timestepping_method()
+from anuga.utilities.argparsing import parse_standard_args
+alg, cfl = parse_standard_args()
+domain.set_flow_algorithm(alg)
+domain.set_CFL(cfl)
 
 #------------------------------------------------------------------------------
@@ -79 +81 @@
 Br = anuga.Reflective_boundary(domain)      # Solid reflective wall
 
-
-
 # Associate boundary tags with boundary objects
 domain.set_boundary({'left': Br, 'right': Br, 'top': Br, 'bottom': Br})
-
 
 #===============================================================================
@@ -97 +96 @@
 # and print to file
 #---------------------------------------------
-Point = (0.0, 0.0)
-for n in range(len(domain.triangles)):
-    tri = domain.get_vertex_coordinates(n)
-    if is_inside_triangle(Point,tri):
-        #print 'Point is within triangle with vertices '+'%s'%tri
-        n_point = n
-        break
-print 'The triangle ID containing the point of origin is = ',n_point
+##Point = (0.0, 0.0)
+##for n in range(len(domain.triangles)):
+##    tri = domain.get_vertex_coordinates(n)
+##    if is_inside_triangle(Point,tri):
+##        #print 'Point is within triangle with vertices '+'%s'%tri
+##        n_point = n
+##        break
+##print 'The triangle ID containing the point of origin is = ',n_point