diff options
author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2012-05-11 22:17:18 +0000 |
---|---|---|
committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2012-05-11 22:17:18 +0000 |
commit | af6dfaa5b36603dc81d04e3c5a8e054f64c37293 (patch) | |
tree | 6d05161e5ec7780f44c88587f3c6d4f8a0486db4 /src | |
parent | d96d80fc26def7ffd86ab0f14cae0e48441aab94 (diff) |
GRHydro: correct reconstruction of velocities in ePPM
This is a bug fix to the ePPM scheme.
The velocity components were reconstructed from the plus face values!
This is clearly incorrect.
Patch by Christian Reisswig.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@318 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src')
-rw-r--r-- | src/GRHydro_PPM.F90 | 12 |
1 files changed, 6 insertions, 6 deletions
diff --git a/src/GRHydro_PPM.F90 b/src/GRHydro_PPM.F90 index 832cd38..6bc5880 100644 --- a/src/GRHydro_PPM.F90 +++ b/src/GRHydro_PPM.F90 @@ -229,15 +229,15 @@ trivial_rp = .true. LIMIT(rho, rhoplus(i), enhanced_ppm_C2, rhoplus(i)) rhominus(i+1) = rhoplus(i) - APPROX_AT_CELL_INTERFACE(velxplus, velxplus(i)) + APPROX_AT_CELL_INTERFACE(velx, velxplus(i)) LIMIT(velx, velxplus(i), enhanced_ppm_C2, velxplus(i)) velxminus(i+1) = velxplus(i) - APPROX_AT_CELL_INTERFACE(velyplus, velyplus(i)) + APPROX_AT_CELL_INTERFACE(vely, velyplus(i)) LIMIT(vely, velyplus(i), enhanced_ppm_C2, velyplus(i)) velyminus(i+1) = velyplus(i) - APPROX_AT_CELL_INTERFACE(velzplus, velzplus(i)) + APPROX_AT_CELL_INTERFACE(velz, velzplus(i)) LIMIT(velz, velzplus(i), enhanced_ppm_C2, velzplus(i)) velzminus(i+1) = velzplus(i) end do @@ -258,15 +258,15 @@ trivial_rp = .true. LIMIT(rho, rhoplus(i), enhanced_ppm_C2, rhoplus(i)) rhominus(i+1) = rhoplus(i) - APPROX_AT_CELL_INTERFACE_STENCIL4(velxplus, velxplus(i)) + APPROX_AT_CELL_INTERFACE_STENCIL4(velx, velxplus(i)) LIMIT(velx, velxplus(i), enhanced_ppm_C2, velxplus(i)) velxminus(i+1) = velxplus(i) - APPROX_AT_CELL_INTERFACE_STENCIL4(velyplus, velyplus(i)) + APPROX_AT_CELL_INTERFACE_STENCIL4(vely, velyplus(i)) LIMIT(vely, velyplus(i), enhanced_ppm_C2, velyplus(i)) velyminus(i+1) = velyplus(i) - APPROX_AT_CELL_INTERFACE_STENCIL4(velzplus, velzplus(i)) + APPROX_AT_CELL_INTERFACE_STENCIL4(velz, velzplus(i)) LIMIT(velz, velzplus(i), enhanced_ppm_C2, velzplus(i)) velzminus(i+1) = velzplus(i) end do |