diff options
author | bmundim <bmundim@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2011-04-13 03:29:13 +0000 |
---|---|---|
committer | bmundim <bmundim@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2011-04-13 03:29:13 +0000 |
commit | 4adfcb611a12689096a604e2a510b7ff12c145f1 (patch) | |
tree | 37b8c6c6627e625b3fe7d7194518e983350884d1 /src/GRHydro_FluxM.F90 | |
parent | 7e9323a4e1f5ccedc9d11f4a425b9b03d2325ed7 (diff) |
RIT MHD development:
1) Fix a bug in the divergence cleaning implementation:
a psidc term in the induction equation was implemented
as a source term where it was supposed to be coded as
part of the flux calculation.
2) Introduce divB as a diagnostic grid function.
3) Remove an old file: GRHydro_CalcUpdateM.F90
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@228 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src/GRHydro_FluxM.F90')
-rw-r--r-- | src/GRHydro_FluxM.F90 | 13 |
1 files changed, 1 insertions, 12 deletions
diff --git a/src/GRHydro_FluxM.F90 b/src/GRHydro_FluxM.F90 index 77bee7a..fadbfe6 100644 --- a/src/GRHydro_FluxM.F90 +++ b/src/GRHydro_FluxM.F90 @@ -47,21 +47,10 @@ subroutine num_x_fluxM(dens,sx,sy,sz,tau,Bx,By,Bz,& !!$ [psi^6 tau] vtilde^i +p* v^i - alp b^0 B^i/w tauf = tau*vxt + psipstar*velm - ab0*psiBx/w -!!$ [psi^6 (B^k vtilde^i - B^i vtilde^k] +!!$ [psi^6 (B^k vtilde^i - B^i vtilde^k)] bxf = 0.0 byf = psiBy * vxt - psiBx*vyt bzf = psiBz * vxt - psiBx*vzt end subroutine num_x_fluxM - - - - - - - - - - - |