User: rhaas
Date: 2014/04/15 02:49 PM
Modified:
/trunk/src/
GRHydro_PPMM.F90
Log:
GRHydro: fix reconstruct_Wv in F90 PPMM routine
was using velx in some places where vx should be used
File Changes:
Directory: /trunk/src/
======================
File [modified]: GRHydro_PPMM.F90
Delta lines: +7 -7
===================================================================
--- trunk/src/GRHydro_PPMM.F90 2014-04-15 19:49:37 UTC (rev 628)
+++ trunk/src/GRHydro_PPMM.F90 2014-04-15 19:49:40 UTC (rev 629)
@@ -137,9 +137,9 @@
dmx(i) = 0.d0 &&\
end if
STEEP(rho, drho, dmrho)
- STEEP(velx, dvelx, dmvelx)
- STEEP(vely, dvely, dmvely)
- STEEP(velz, dvelz, dmvelz)
+ STEEP(vx, dvelx, dmvelx)
+ STEEP(vy, dvely, dmvely)
+ STEEP(vz, dvelz, dmvelz)
STEEP(Bvcx, dBvcx, dmBvcx)
STEEP(Bvcy, dBvcy, dmBvcy)
STEEP(Bvcz, dBvcz, dmBvcz)
@@ -362,7 +362,7 @@
do i = 3, nx - 2
dpress2 = press(i+2) - press(i-2)
- dvel = velx(i+1) - velx(i-1)
+ dvel = vx(i+1) - vx(i-1)
if ( (abs(dpress(i)) > ppm_epsilon * min(press(i-1),press(i+1))) .and. &
(dvel < 0.d0) ) then
w = 1.d0
@@ -472,9 +472,9 @@
end if
MON(rhominus,rho,rhoplus)
- MON(velxminus,velx,velxplus)
- MON(velyminus,vely,velyplus)
- MON(velzminus,velz,velzplus)
+ MON(velxminus,vx,velxplus)
+ MON(velyminus,vy,velyplus)
+ MON(velzminus,vz,velzplus)
MON(Bvcxminus,Bvcx,Bvcxplus)
MON(Bvcyminus,Bvcy,Bvcyplus)
MON(Bvczminus,Bvcz,Bvczplus)