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