Fix a severe bug: forgot to divide by the grid spacing.

Remove spurious factor of 1/16.

Add second order dissipation, which is first order correct.

Checked stability interval.  It should now be as described.


git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/Dissipation/trunk@4 850bcc8b-0e4f-0410-8c26-8d28fbf1eda9

1 files changed, 7 insertions, 1 deletions

diff --git a/param.ccl b/param.ccl
index b50b7da..5640ea1 100644
--- a/param.ccl
+++ b/param.ccl
@@ -5,9 +5,15 @@ BOOLEAN verbose "produce log output" STEERABLE=always
 {
 } "no"
 
+INT order "Dissipation order" STEERABLE=always
+{
+  2 :: "second order dissipation (first order accurate)"
+  4 :: "fourth order dissipation (third order accurate)"
+} 4
+
 REAL epsdis "Dissipation strength" STEERABLE=always
 {
-  *:* :: "0 for no dissipation.  Unstable for epsdis<0 and epsdis>1/"
+  *:* :: "0 for no dissipation.  Unstable for epsdis<0 and epsdis>1/3"
 } 0.2
 
 STRING vars "List of evolved grid functions that should have dissipation added" STEERABLE=always