Changeset 4715

Timestamp:

Sep 10, 2007, 5:17:33 PM (17 years ago)

Author:

ole

Message:

Work on Will's convergence problem

Location:

anuga_validation/convergence_study

Files:

• animatesww2d_alt.py (modified) (1 diff)
• convergence.py (modified) (3 diffs)

anuga_validation/convergence_study/animatesww2d_alt.py

@@ -124 +124 @@
     close('all')
 
-swwfile = join(getenv('INUNDATIONHOME'),'data','validation','convergence_study','5000','myexample2.sww')
-animatesww2d(swwfile)
+if __name__ == '__main__':
+    import sys
+    swwfile = sys.argv[1]
+    
+    animatesww2d(swwfile)
+    
+#swwfile = join(getenv('INUNDATIONHOME'),'data','validation','convergence_study','5000','myexample2.sww')
+#animatesww2d(swwfile)
 
 # once png files have been created, they can then be dragged into

anuga_validation/convergence_study/convergence.py

@@ -44 +44 @@
                          boundary_tags={'left': [0], 'top': [1],
                                         'right': [2], 'bottom': [3]},
-                         maximum_triangle_area = 1000,
+                         maximum_triangle_area = 100000,
                          filename=meshname)
 
@@ -50 +50 @@
 domain = Domain(meshname, use_cache=True, verbose = True)
 domain.set_name('myexample2')                
+domain.set_default_order(2) # Use second order spatial scheme
 domain.set_datadir('.')                     # Use current directory for output
 
@@ -56 +57 @@
 #------------------------------------------------------------------------------
 #domain.set_quantity('elevation', topography) # Use function for elevation
-domain.set_quantity('elevation',-100)
+domain.set_quantity('elevation', -100)
 domain.set_quantity('friction', 0.00)
 domain.set_quantity('stage', 0.0)