[Commits] [svn:einsteintoolkit] GRHydro/branches/divclean/src/ (Rev. 237)

[tanja.bode at physics.gatech.edu]

User: tbode
Date: 2011/04/27 06:05 PM

Modified:
 /branches/divclean/src/
  GRHydro_SourceM.F90

Log:
 Fix variable names.

File Changes:

Directory: /branches/divclean/src/
==================================

File [modified]: GRHydro_SourceM.F90
Delta lines: +7 -7
===================================================================
--- branches/divclean/src/GRHydro_SourceM.F90	2011-04-27 22:58:37 UTC (rev 236)
+++ branches/divclean/src/GRHydro_SourceM.F90	2011-04-27 23:05:18 UTC (rev 237)
@@ -444,7 +444,7 @@
         taurhs(i,j,k) = alp(i,j,k)*sqrtdet*tau_source
         if(clean_divergence.ne.0) then
            psidcrhs(i,j,k) = -1.d0 * ( ch_dc*ch_dc/cp_dc/cp_dc*alp(i,j,k) + &
-                dx_shiftx + dy_shifty + dz_shiftz ) * psidc(i,j,k) + &
+                dx_betax + dy_betay + dz_betaz ) * psidc(i,j,k) + &
                 Bconsx(i,j,k) * ( dx_alp - half*alp(i,j,k) * &
                  ( uxx*dx_gxx + uyy*dx_gyy + uzz*dx_gzz + 2.d0*uxy*dx_gxy + &
                    2.d0*uxz*dx_gxz + 2.d0*uyz*dx_gyz ) )/ sqrtdet + &
@@ -463,20 +463,20 @@
            dz_det_bydet = uxx*dz_gxx + uyy*dz_gyy + uzz*dz_gzz + &
                 2.d0*(uxy*dz_gxy+uxz*dz_gxz+uyz*dz_gyz)
 
-           bvcx_source = -1.d0 * ( Bconsx(i,j,k)*dx_shiftx + &
-                Bconsy(i,j,k)*dy_shiftx + Bconsz(i,j,k)*dz_shiftx ) + &
+           bvcx_source = -1.d0 * ( Bconsx(i,j,k)*dx_betax + &
+                Bconsy(i,j,k)*dy_betax + Bconsz(i,j,k)*dz_betax ) + &
                 psidc(i,j,k)*sqrtdet*( uxx*dx_alp+uxy*dy_alp+uxz*dz_alp ) + &
                 psidc(i,j,k)*3.d0*half*alp(i,j,k)*sqrtdet*( uxx*dx_det_bydet + &
                   uxy*dy_det_bydet + uxz*dz_det_bydet )
 
-           bvcy_source = -1.d0 * ( Bconsx(i,j,k)*dx_shifty + &
-                Bconsy(i,j,k)*dy_shifty + Bconsz(i,j,k)*dz_shifty ) + &
+           bvcy_source = -1.d0 * ( Bconsx(i,j,k)*dx_betay + &
+                Bconsy(i,j,k)*dy_betay + Bconsz(i,j,k)*dz_betay ) + &
                 psidc(i,j,k)*sqrtdet*( uxy*dx_alp+uyy*dy_alp+uyz*dz_alp ) + &
                 psidc(i,j,k)*3.d0*half*alp(i,j,k)*sqrtdet*( uxy*dx_det_bydet + &
                   uyy*dy_det_bydet + uyz*dz_det_bydet )
 
-           bvcz_source = -1.d0 * ( Bconsx(i,j,k)*dx_shiftz + &
-                Bconsy(i,j,k)*dy_shiftz + Bconsz(i,j,k)*dz_shiftz ) + &
+           bvcz_source = -1.d0 * ( Bconsx(i,j,k)*dx_betaz + &
+                Bconsy(i,j,k)*dy_betaz + Bconsz(i,j,k)*dz_betaz ) + &
                 psidc(i,j,k)*sqrtdet*( uxz*dx_alp+uyz*dy_alp+uzz*dz_alp ) + &
                 psidc(i,j,k)*3.d0*half*alp(i,j,k)*sqrtdet*( uxz*dx_det_bydet + &
                   uyz*dy_det_bydet + uzz*dz_det_bydet )