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, 22 insertions, 0 deletions

diff --git a/interface.ccl b/interface.ccl
index 6e26e97..1ff5f1f 100644
--- a/interface.ccl
+++ b/interface.ccl
@@ -3,5 +3,27 @@
 
 IMPLEMENTS: Dissipation
 
+INHERITS: grid SphericalSurface
+
 CCTK_INT FUNCTION MoLQueryEvolvedRHS (CCTK_INT IN EvolvedIndex)
 REQUIRES FUNCTION MoLQueryEvolvedRHS
+
+CCTK_INT FUNCTION \
+    SymmetryTableHandleForGrid (CCTK_POINTER_TO_CONST IN cctkGH)
+USES FUNCTION SymmetryTableHandleForGrid
+
+CCTK_INT FUNCTION HorizonRadiusInDirection                              \
+   (CCTK_INT IN horizon_number,                                         \
+    CCTK_INT IN N_points,                                               \
+    CCTK_REAL IN ARRAY x, CCTK_REAL IN ARRAY y, CCTK_REAL IN ARRAY z,   \
+    CCTK_REAL OUT ARRAY radius)
+USES FUNCTION HorizonRadiusInDirection
+
+
+
+PRIVATE:
+
+REAL epsdisA_group TYPE=gf
+{
+  epsdisA
+} "dissipation array for spatially different dissipation"