[Commits] [svn:einsteintoolkit] GRHydro_InitData/trunk/src/ (Rev. 128)

[roland.haas at physics.gatech.edu]

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