diff options
| author | bmundim <bmundim@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2011-09-26 15:50:04 +0000 |
|---|---|---|
| committer | bmundim <bmundim@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2011-09-26 15:50:04 +0000 |
| commit | 56b2e300fb3c5e7ff100d5f2813393492679f99e (patch) | |
| tree | d6e06b7a49f70106c047754ecc22170948c64f05 /src/GRHydro_SourceM.F90 | |
| parent | e76b9040dbc377e03def86e744503523a41e0ed6 (diff) | |
RIT dev: fix source term in divergence cleaning. Thanks Tanja
for catching this up too.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@280 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src/GRHydro_SourceM.F90')
| -rw-r--r-- | src/GRHydro_SourceM.F90 | 23 |
1 files changed, 10 insertions, 13 deletions
diff --git a/src/GRHydro_SourceM.F90 b/src/GRHydro_SourceM.F90 index ed0963e..918b308 100644 --- a/src/GRHydro_SourceM.F90 +++ b/src/GRHydro_SourceM.F90 @@ -457,27 +457,24 @@ subroutine SourceTermsM(CCTK_ARGUMENTS) 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 |