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 /schedule.ccl | |
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 'schedule.ccl')
-rw-r--r-- | schedule.ccl | 31 |
1 files changed, 24 insertions, 7 deletions
diff --git a/schedule.ccl b/schedule.ccl index 284c773..d9e8552 100644 --- a/schedule.ccl +++ b/schedule.ccl @@ -68,10 +68,17 @@ STORAGE:GRHydro_reflevel STORAGE:densrhs STORAGE:taurhs STORAGE:srhs -STORAGE:Bvecrhs -if (clean_divergence) +if(CCTK_Equals(Bvec_evolution_method,"GRHydro")) { - STORAGE:psidcrhs + STORAGE:Bvecrhs + if (clean_divergence) + { + STORAGE:psidcrhs + } + if (track_divB) + { + STORAGE:divB + } } STORAGE:GRHydro_eos_scalars STORAGE:GRHydro_minima @@ -97,12 +104,22 @@ schedule group GRHydro_Initial IN HydroBase_Initial BEFORE SetTmunu { } "GRHydro initial data group" -if(CCTK_Equals(Bvec_evolution_method,"GRHydro") && clean_divergence) +if(CCTK_Equals(Bvec_evolution_method,"GRHydro")) { - schedule GRHydro_InitDivergenceClean IN HydroBase_Initial + if (clean_divergence) { - LANG: Fortran - } "Set psi for divergence cleaning initially to zero" + schedule GRHydro_InitDivergenceClean IN HydroBase_Initial + { + LANG: Fortran + } "Set psi for divergence cleaning initially to zero" + } + if (track_divB) + { + schedule GRHydro_DivBInit IN HydroBase_Initial + { + LANG: Fortran + } "Set divB initially to zero" + } } ################################################# |