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

[tanja.bode at physics.gatech.edu]

User: tbode
Date: 2011/04/27 07:47 PM

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

Log:
 Fix Bcons sources.

File Changes:

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

File [modified]: GRHydro_SourceM.F90
Delta lines: +26 -6
===================================================================
--- branches/divclean/src/GRHydro_SourceM.F90	2011-04-28 00:37:56 UTC (rev 238)
+++ branches/divclean/src/GRHydro_SourceM.F90	2011-04-28 00:47:43 UTC (rev 239)
@@ -85,6 +85,7 @@
   !!CCTK_REAL :: dx_psidc, dy_psidc, dz_psidc
   CCTK_REAL :: bvcx_source, bvcy_source, bvcz_source
   CCTK_REAL :: dx_det_bydet, dy_det_bydet, dz_det_bydet
+  CCTK_REAL :: gdg_x, gdg_y, gdg_z !! g^{ik} d_k g_{ij}
 
   CCTK_REAL :: Bvecxlow,Bvecylow,Bveczlow,bdotv,b2,dum,bxlow,bylow,bzlow
   CCTK_REAL :: bt,bx,by,bz,rhohstarW2,pstar
@@ -463,23 +464,42 @@
            dz_det_bydet = uxx*dz_gxx + uyy*dz_gyy + uzz*dz_gzz + &
                 2.d0*(uxy*dz_gxy+uxz*dz_gxz+uyz*dz_gyz)
 
+           !!$ g^{ik} d_k g_{li}
+           gdg_x = uxx*dx_gxx + uxy*dy_gxx + uxz*dz_gxx + &
+                   uxy*dx_gxy + uyy*dy_gxy + uyz*dz_gxy + &
+                   uxz*dx_gxz + uyz*dy_gxz + uzz*dz_gxz
+
+           gdg_y = uxx*dx_gxy + uxy*dy_gxy + uxz*dz_gxy + &
+                   uxy*dx_gyy + uyy*dy_gyy + uyz*dz_gyy + &
+                   uxz*dx_gyz + uyz*dy_gyz + uzz*dz_gyz
+
+           gdg_z = uxx*dx_gxz + uxy*dy_gxz + uxz*dz_gxz + &
+                   uxy*dx_gyz + uyy*dy_gyz + uyz*dz_gyz + &
+                   uxz*dx_gzz + uyz*dy_gzz + uzz*dz_gzz
+
            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 )
+                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 )
 
            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 )
+                psidc(i,j,k)*alp(i,j,k)*sqrtdet*( 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 )
 
            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 )
+                psidc(i,j,k)*alp(i,j,k)*sqrtdet*( 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 )
 
            Bconsrhsx(i,j,k) = bvcx_source
            Bconsrhsy(i,j,k) = bvcy_source