GRHydro: Introduce variable atmopshere level / tolerance as function of radius. Off by default.

From: Christian Reisswig <reisswig@scriwalker.(none)>

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

1 files changed, 20 insertions, 0 deletions

diff --git a/param.ccl b/param.ccl
index 0d17914..47ca995 100644
--- a/param.ccl
+++ b/param.ccl
@@ -455,6 +455,26 @@ REAL GRHydro_atmo_tolerance "A point is set to atmosphere in the Con2Prim's if i
    0.0: :: "Zero or larger. A useful value could be 0.0001"
 }  0.0
 
+REAL atmo_falloff_radius "The radius for which the atmosphere starts to falloff as (atmo_falloff_radius/r)**atmo_falloff_power"
+{
+  0:* :: "Anything positive"
+} 50.0
+
+REAL atmo_falloff_power "The power at which the atmosphere level falls off as (atmo_falloof_radius/r)**atmo_falloff_power"
+{
+  0:* :: "Anything positive"
+} 0.0
+
+REAL atmo_tolerance_radius "The radius for which the atmosphere tolerance starts to increase as (r/atmo_tolerance_radius)**atmo_tolerance_power"
+{
+  0:* :: "Anything positive"
+} 50.0
+
+REAL atmo_tolerance_power "The power at which the atmosphere tolerance increases as (r/atmo_tolerance_radius)**atmo_tolerance_power"
+{
+  0:* :: "Anything positive"
+} 0.0
+
 BOOLEAN wk_atmosphere "Use some of Wolfgang Kastauns atmosphere tricks"
 {
 } "no"