User: rhaas Date: 2011/04/29 07:22 PM

Modified: /trunk/src/ GRHydro_ShockTubeM.F90

Log: GRHydro_InitialData: apply boundary conditions to Bcons for diagonal shocks

File Changes:

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

File [modified]: GRHydro_ShockTubeM.F90 Delta lines: +9 -9 =================================================================== --- trunk/src/GRHydro_ShockTubeM.F90 2011-04-28 01:14:57 UTC (rev 127) +++ trunk/src/GRHydro_ShockTubeM.F90 2011-04-30 00:22:45 UTC (rev 128) @@ -708,9 +708,9 @@ sy(i,j,k) = sy(inew,jnew,knew) sz(i,j,k) = sz(inew,jnew,knew) tau(i,j,k) = tau(inew,jnew,knew) - Bvecx(i,j,k)=Bvecx(inew,jnew,knew) - Bvecy(i,j,k)=Bvecy(inew,jnew,knew) - Bvecz(i,j,k)=Bvecz(inew,jnew,knew) + Bconsx(i,j,k)=Bconsx(inew,jnew,knew) + Bconsy(i,j,k)=Bconsy(inew,jnew,knew) + Bconsz(i,j,k)=Bconsz(inew,jnew,knew) if(clean_divergence.ne.0) then psidc(i,j,k)=psidc(inew,jnew,knew) endif @@ -789,9 +789,9 @@ sy(i,j,k) = sy(inew,jnew,k) sz(i,j,k) = sz(inew,jnew,k) tau(i,j,k) = tau(inew,jnew,k) - Bvecx(i,j,k)=Bvecx(inew,jnew,k) - Bvecy(i,j,k)=Bvecy(inew,jnew,k) - Bvecz(i,j,k)=Bvecz(inew,jnew,k) + Bconsx(i,j,k)=Bconsx(inew,jnew,k) + Bconsy(i,j,k)=Bconsy(inew,jnew,k) + Bconsz(i,j,k)=Bconsz(inew,jnew,k) if(clean_divergence.ne.0) then psidc(i,j,k)=psidc(inew,jnew,k) endif @@ -805,9 +805,9 @@ sy(i,j,k) = sy(i,j,kc) sz(i,j,k) = sz(i,j,kc) tau(i,j,k) = tau(i,j,kc) - Bvecx(i,j,k)=Bvecx(i,j,kc) - Bvecy(i,j,k)=Bvecy(i,j,kc) - Bvecz(i,j,k)=Bvecz(i,j,kc) + Bconsx(i,j,k)=Bconsx(i,j,kc) + Bconsy(i,j,k)=Bconsy(i,j,kc) + Bconsz(i,j,k)=Bconsz(i,j,kc) if(clean_divergence.ne.0) then psidc(i,j,k)=psidc(i,j,kc) endif