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)