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

diff --git a/src/paramcheck.c b/src/paramcheck.c
new file mode 100644
index 0000000..d3629fe
--- /dev/null
+++ b/src/paramcheck.c
@@ -0,0 +1,29 @@
+/* $Header$ */
+
+#include "cctk.h"
+#include "cctk_Arguments.h"
+#include "cctk_Parameters.h"
+
+void
+dissipation_paramcheck (CCTK_ARGUMENTS)
+{
+  DECLARE_CCTK_ARGUMENTS;
+  DECLARE_CCTK_PARAMETERS;
+  
+  int i;
+  int want_horizon;
+  
+  if (extra_dissipation_in_horizons)
+  {
+    want_horizon = 0;
+    for (i=0; i<100; ++i)
+    {
+      want_horizon = want_horizon | (horizon_number[i] >= 0);
+    }
+    
+    if (want_horizon && ! CCTK_IsFunctionAliased ("HorizonRadiusInDirection"))
+    {
+      CCTK_PARAMWARN ("The aliased function \"HorizonRadiusInDirection\" must be defined when the parameter \"extra_dissipation_in_horizons\" is set and one of the sources is AHFinderDirect");
+    }
+  }
+}