diff options
author | bmundim <bmundim@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2010-10-12 16:43:41 +0000 |
---|---|---|
committer | bmundim <bmundim@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2010-10-12 16:43:41 +0000 |
commit | dbf4be4d3522c08f81ed6d303f8f9f36e1d50cc7 (patch) | |
tree | a2ea60f96f040a709719daa89f0e1bc95a6dc827 /src/GRHydro_CalcUpdateM.F90 | |
parent | 1a80cb7aa17044e97ad435c3ac7f7744ba22504b (diff) |
RIT's divergence cleaning development.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@161 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src/GRHydro_CalcUpdateM.F90')
-rw-r--r-- | src/GRHydro_CalcUpdateM.F90 | 10 |
1 files changed, 8 insertions, 2 deletions
diff --git a/src/GRHydro_CalcUpdateM.F90 b/src/GRHydro_CalcUpdateM.F90 index fab2b55..ab5e55d 100644 --- a/src/GRHydro_CalcUpdateM.F90 +++ b/src/GRHydro_CalcUpdateM.F90 @@ -1,6 +1,6 @@ /*@@ @file GRHydro_CalcUpdateM.F90 - @date Aug 30, 2010 + @date Oct 10, 2010 @author Joshua Faber, Scott Noble, Bruno Mundim, Ian Hawke @desc Calculates the update terms given the fluxes. Moved to here so that @@ -15,7 +15,7 @@ /*@@ @routine UpdateCalculationM - @date Aug 30, 2010 + @date Oct 10, 2010 @author Joshua Faber, Scott Noble, Bruno Mundim, Ian Hawke @desc Calculates the update terms from the fluxes. @@ -89,6 +89,9 @@ subroutine UpdateCalculationM(CCTK_ARGUMENTS) Bvecrhs(i,j,k,3) = Bvecrhs(i,j,k,3) + & (alp_l * Bveczflux(i-xoffset,j-yoffset,k-zoffset) - & alp_r * Bveczflux(i,j,k)) * idx + psidcrhs(i,j,k) = psidcrhs(i,j,k) + & + (alp_l * psidcflux(i-xoffset,j-yoffset,k-zoffset) - & + alp_r * psidcflux(i,j,k)) * idx if (evolve_tracer .ne. 0) then @@ -234,6 +237,9 @@ subroutine UpdateCalculationM(CCTK_ARGUMENTS) Bvecrhs(i,j,k,3) = Bvecrhs(i,j,k,3) + & (Bveczflux(i-xoffset,j-yoffset,k-zoffset) - & Bveczflux(i,j,k)) * idx + psidcrhs(i,j,k) = psidcrhs(i,j,k) + & + (psidcflux(i-xoffset,j-yoffset,k-zoffset) - & + psidcflux(i,j,k)) * idx enddo enddo |