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@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"