[Commits] [svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 280)

[bcmsma at astro.rit.edu]

User: bmundim
Date: 2011/09/26 10:50 AM

Modified:
 /trunk/src/
  GRHydro_SourceM.F90

Log:
 RIT dev: fix source term in divergence cleaning. Thanks Tanja
 for catching this up too.

File Changes:

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

File [modified]: GRHydro_SourceM.F90
Delta lines: +10 -13
===================================================================
--- trunk/src/GRHydro_SourceM.F90	2011-09-25 22:59:03 UTC (rev 279)
+++ trunk/src/GRHydro_SourceM.F90	2011-09-26 15:50:04 UTC (rev 280)
@@ -457,27 +457,24 @@
 
            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)*alp(i,j,k)*sqrtdet*( uxx*dx_det_bydet + &
-                  uxy*dy_det_bydet + uxz*dz_det_bydet ) - &
-                psidc(i,j,k)*alp(i,j,k)*sqrtdet*( uxx*gdg_x + uxy*gdg_y + &
-                uxz*gdg_z )
+                psidc(i,j,k)*sqrtdet*(( uxx*dx_alp+uxy*dy_alp+uxz*dz_alp ) + &
+                alp(i,j,k)*(0.5d0*( uxx*dx_det_bydet + &
+                  uxy*dy_det_bydet + uxz*dz_det_bydet) - &
+                ( uxx*gdg_x + uxy*gdg_y + uxz*gdg_z )))
 
            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)*alp(i,j,k)*sqrtdet*( uxy*dx_det_bydet + &
+                psidc(i,j,k)*sqrtdet*(( uxy*dx_alp+uyy*dy_alp+uyz*dz_alp ) + &
+                alp(i,j,k)*(0.5d0*( uxy*dx_det_bydet + &
                   uyy*dy_det_bydet + uyz*dz_det_bydet ) - &
-                psidc(i,j,k)*alp(i,j,k)*sqrtdet*( uxy*gdg_x + uyy*gdg_y + &
-                uyz*gdg_z )
+                ( uxy*gdg_x + uyy*gdg_y + uyz*gdg_z )))
 
            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)*alp(i,j,k)*sqrtdet*( uxz*dx_det_bydet + &
+                psidc(i,j,k)*sqrtdet*(( uxz*dx_alp+uyz*dy_alp+uzz*dz_alp ) + &
+                alp(i,j,k)*(0.5d0*( uxz*dx_det_bydet + &
                   uyz*dy_det_bydet + uzz*dz_det_bydet ) - &
-                psidc(i,j,k)*alp(i,j,k)*sqrtdet*( uxz*gdg_x + uyz*gdg_y + &
-                uzz*gdg_z )
+                ( uxz*gdg_x + uyz*gdg_y + uzz*gdg_z )))
 
            Bconsrhsx(i,j,k) = bvcx_source
            Bconsrhsy(i,j,k) = bvcy_source