Changeset 4219

Timestamp:

Feb 6, 2007, 12:00:06 PM (17 years ago)

Author:

ole

Message:

Diagnostic work on flux and slope limiters

Files:

• anuga_core/source/anuga/config.py (modified) (1 diff)
• anuga_core/source/anuga/examples/limiter_test.py (modified) (3 diffs)
• anuga_validation/okushiri_2005/run_okushiri.py (modified) (1 diff)
• anuga_work/development/dam_test_from_brad_2007 (added)
• anuga_work/development/dam_test_from_brad_2007/dam.csv (added)
• anuga_work/development/dam_test_from_brad_2007/dam_sample.py (added)

anuga_core/source/anuga/config.py

@@ -81 +81 @@
 beta_h      = 0.2
 
-alpha_balance = 2.0
+# Alpha_balance controls how limiters are balanced between deep and shallow.
+# A large value will favour the deep water limiters, allowing the a closer hug to the coastline.
+# This will minimise 'creep' but at the same time cause smaller time steps
+# Range:
+
+alpha_balance = 2.0 
+
 
 CFL = 1.0  #FIXME (ole): Is this in use yet??

anuga_core/source/anuga/examples/limiter_test.py

@@ -26 +26 @@
 domain.set_name('limiter_test')
 domain.set_store_vertices_uniquely(True)
+domain.set_default_order(2)
 
 #------------------------------------------------------------------------------
@@ -44 +45 @@
 
 domain.set_quantity('elevation', topography) # Use function for elevation
-domain.set_quantity('stage', 0.5)             # Constant negative initial stage
+domain.set_quantity('stage', 0.5)           # Constant initial stage
 
 
@@ -58 +59 @@
 #------------------------------------------------------------------------------
 
-for t in domain.evolve(yieldstep = 0.1, finaltime = 1.0):
+for t in domain.evolve(yieldstep = 0.1, finaltime = 50.0):
     domain.write_time()
+
+print 'done'    
 
     

anuga_validation/okushiri_2005/run_okushiri.py

@@ -80 +80 @@
 
 for t in domain.evolve(yieldstep = 0.05, finaltime = 22.5):
-    domain.write_time()
+    domain.write_time(track_speeds=True)
 
 print 'That took %.2f seconds' %(time.time()-t0)