Changeset 8696

Timestamp:

Feb 14, 2013, 9:17:51 PM (12 years ago)

Author:

davies

Message:

Adjusting the tsunami algorithm defaults to always store vertices uniquely

File:

• trunk/anuga_core/source/anuga/shallow_water/shallow_water_domain.py (modified) (3 diffs)

trunk/anuga_core/source/anuga/shallow_water/shallow_water_domain.py

@@ -361 +361 @@
 
         # We demand that vertex values are stored uniquely
-        #self.set_store_vertices_smoothly(False)
+        # NOTE: Without this, wet-dry velocity spikes tend to occur in output files.
+        # The tsunami algorithm is probably particularly sensitive to this.
+        # Hence, this seems like a reasonable default.
+        self.set_store_vertices_smoothly(False)
 
         self.maximum_allowed_speed=0.0
@@ -371 +374 @@
 #        print '##########################################################################'
 #        print '#'
-#        print '# Using anuga_tsunami solver in anuga_work/development/gareth/experimental/anuga_tsunami/'
 #        print "#"
-#        print "# This solver is experimental. Here are some tips on using it"
+#        print "# Here are some tips on using the 'tsunami' solver"
 #        print "#"
 #        print "# 1) When plotting outputs, I strongly suggest you examine centroid values, not vertex values"
@@ -399 +401 @@
 #        print "# You will probably be able to tell this is causing you problems by convergence testing"
 #        print "#"
-#        print '# 5) Note that many options in config.py have been overridden by the solver -- I have '
+#        print '# 5) Note that many options in config.py have been overridden by the solver -- we have '
 #        print '# deliberately attempted to get the solver to perform well with consistent values of '
-#        print '# these parameters -- so I would advise against changing them unless you at least check that '
+#        print '# these parameters -- so I advise against changing them unless you at least check that '
 #        print '# it really does improve things'
 #        print '#'