Add Frank Herrmann's extra dissipation switches.

Two new parameters, Dissipation::extra_dissipation_in_horizons and
Dissipation::extra_dissipation_at_outerbound, increase the disspation
coefficient inside horizons and near the outer boundary.


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

1 files changed, 48 insertions, 0 deletions

diff --git a/param.ccl b/param.ccl
index ca4cbcb..3b4725f 100644
--- a/param.ccl
+++ b/param.ccl
@@ -20,3 +20,51 @@ STRING vars "List of evolved grid functions that should have dissipation added"
 {
   .* :: "must be a valid list of grid functions"
 } ""
+
+
+
+# Extra dissipation in horizons
+
+BOOLEAN extra_dissipation_in_horizons "extra dissipation in horizons" STEERABLE=always
+{
+} "no"
+
+INT update_ah_every "how often to update the AH information for dissipation" STEERABLE=always
+{
+  0:* :: "positive iteration number"
+} 1
+
+REAL ah_slope "Slope (epsdis / distance) inside AH" STEERABLE=always
+{
+  *:* :: "0 for no dissipation.  Unstable when total epsdis<0 or total epsdis>1/3"
+} 0.2
+
+INT surface_number[100] "horizon number for extra dissipation in horizons -- SphericalSurface number" STEERABLE=always
+{
+  -1  :: "do not use a spherical surface"
+  0:* :: "surface number (from SphericalSurface); starts from 0"
+} -1
+
+INT horizon_number[100] "horizon number for extra dissipation in horizons -- AHFinderDirect number" STEERABLE=always
+{
+  -1  :: "do not use a horizon"
+  1:* :: "horizon number (from AHFinderDirect); starts from 1"
+} -1
+
+
+
+# Extra dissipation at outer boundary
+
+BOOLEAN extra_dissipation_at_outerbound "increase dissipation at outer boundary" STEERABLE=always
+{
+} "no"
+
+INT outer_bound_npoints "number of points in which dissipation should be increased" STEERABLE=always
+{
+  0:* :: "positive number"
+} 3
+
+REAL ob_slope "slope (epsdis / distance) at outer boundary" STEERABLE=always
+{
+  0:* :: "increase dissipation"
+} 5