add option to smoothly dampen Tmunu to vacuum outside of a given radius

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@266 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45

1 files changed, 16 insertions, 0 deletions

diff --git a/param.ccl b/param.ccl
index e042819..b93b451 100644
--- a/param.ccl
+++ b/param.ccl
@@ -532,6 +532,22 @@ boolean track_divB "Track the magnetic field constraint violations"
 {
 } "no"
 
+############### Parameters for Tmunu damping (useful in atmosphere) #####################
+
+CCTK_REAL Tmunu_damping_radius_min "damping radius at which we start to damp with a tanh function"
+{
+   -1  :: "damping switched off"
+   0:* :: "damping radius at which we start to damp"
+} -1
+
+
+CCTK_REAL Tmunu_damping_radius_max "damping radius at which Tmunu becomes 0"
+{
+   -1  :: "damping switched off"
+   0:* :: "greater than minimum radius above"
+} -1
+
+
 ############### temporary parameters to be removed once connected issues are fixed.
 
 boolean con2prim_oct_hack "Disregard c2p failures in oct/rotsym90 boundary cells with xyz < 0" STEERABLE=ALWAYS