Implemented enhanced PPM scheme by Colella & Sekora 2008, McCorquodale & Colella 2011.

Can be actived by setting use_enhanced_ppm = yes.

Also fixed a bug in GRHydro_Reconstruct.F90 when using MP.



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

1 files changed, 21 insertions, 0 deletions

diff --git a/param.ccl b/param.ccl
index a4f037c..537c2a8 100644
--- a/param.ccl
+++ b/param.ccl
@@ -193,6 +193,27 @@ int ppm_mppm_debug_eigenvalues "write eigenvalues into debug grid variables"
   0:1 :: "0 (off, default) or 1 (on)" 
 } 0
 
+
+boolean use_enhanced_ppm "Use the enhanced ppm reconstruction method proposed by Colella & Sekora 2008 and McCorquodale & Colella 2011" STEERABLE=RECOVER
+{
+} no
+
+real enhanced_ppm_C2 "Parameter for enhancecd ppm limiter. Default from McCorquodale & Colella 2011"
+{
+   *:* :: "must be greater than 1. According to Colella&Sekora 2008, enhanced ppm is insensitive to C in [1.25,5]"
+} 1.25
+
+real enhanced_ppm_C3 "Parameter for enhancecd ppm limiter. Default from McCorquodale & Colella 2011"
+{
+   0:* :: "must be greater than 0."
+} 0.1
+
+real GRHydro_ppm_atmo_tolerance "A relative value for rho below which we switch to standard ppm. This may be necessary for very sharp transitions in the density to atmosphere values." STEERABLE=RECOVER
+{
+   0:* :: "anything positive"
+} 1e4
+
+
 int eno_order "The order of accuracy of the ENO reconstruction"
 {
   1:*		:: "Anything above 1, but above 5 is pointless"