[Commits] [svn:einsteintoolkit] GRHydro/trunk/ (Rev. 398)

[rhaas at tapir.caltech.edu]

User: rhaas
Date: 2012/07/17 12:08 PM

Modified:
 /trunk/
  param.ccl
 /trunk/src/
  GRHydro_PPM.F90

Log:
 GRHydro: ePPM: Allow for tolerance in vel^2 limiting.
 
 From: Christian Reisswig <reisswig at scriwalker.(none)>

File Changes:

Directory: /trunk/src/
======================

File [modified]: GRHydro_PPM.F90
Delta lines: +1 -1
===================================================================
--- trunk/src/GRHydro_PPM.F90	2012-07-17 17:08:26 UTC (rev 397)
+++ trunk/src/GRHydro_PPM.F90	2012-07-17 17:08:28 UTC (rev 398)
@@ -791,7 +791,7 @@
 
 !! Check if velocity is below speed of light. If not, reduce reconstruction to first order!
 #define VELCHECK(vxminus,vx,vxplus,vyminus,vy,vyplus,vzminus,vz,vzplus) \
-   if (VEL2M(vxminus(i),vyminus(i),vzminus(i)) .ge. 1.0d0 .or. VEL2P(vxplus(i),vyplus(i),vzplus(i)) .ge. 1.0d0) then &&\
+   if (VEL2M(vxminus(i),vyminus(i),vzminus(i)) .ge. (1.0d0-enhanced_ppm_C4) .or. VEL2P(vxplus(i),vyplus(i),vzplus(i)) .ge. (1.0d0-enhanced_ppm_C4)) then &&\
       vxminus(i) = vx(i)  &&\
       vxplus(i) = vx(i)   &&\
       vyminus(i) = vy(i)  &&\

Directory: /trunk/
==================

File [modified]: param.ccl
Delta lines: +5 -0
===================================================================
--- trunk/param.ccl	2012-07-17 17:08:26 UTC (rev 397)
+++ trunk/param.ccl	2012-07-17 17:08:28 UTC (rev 398)
@@ -218,7 +218,12 @@
    0:* :: "must be greater than 0."
 } 0.1
 
+real enhanced_ppm_C4 "Parameter for enhancecd ppm limiter. This is for the additional velocity limiter" STEERABLE=ALWAYS
+{
+   0:1 :: "must be greater than or equal to 0, and not larger than 1"
+} 0.0
 
+
 int eno_order "The order of accuracy of the ENO reconstruction"
 {
   1:*		:: "Anything above 1, but above 5 is pointless"