GRHydro: fixes for GRMHD

All by Philipp Moesta.

1) Fix parity of psidc and divb
2) Fix a wrong index in the source terms of scon
3) Fix wrong indices of derivatives of space-time metric in the source of the
   divergence cleaning scalar.
4) Calculate divergence of B in MoL PseudoEvolution and set its
   Prolongation="Restrict".
5) Correct the source terms and fluxes for the Bfield and the divergence
   cleaning field when having a non-flat space-time.
6) Make sure alpha factors match between UpdateCalculation and fluxes
   definition.
7) Include 1/sqrt(detg) factor in calculation of \epsilon^{\muijk} in the cross
   product to obtain the Bfield form the vector potential.


git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@330 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45

1 files changed, 19 insertions, 0 deletions

diff --git a/schedule.ccl b/schedule.ccl
index 392086a..7e9b765 100644
--- a/schedule.ccl
+++ b/schedule.ccl
@@ -1162,3 +1162,22 @@ schedule GROUP SetTmunu AT POSTPOSTINITIAL AFTER Con2Prim BEFORE ADMConstraintsG
 {
 } "Calculate the stress-energy tensor"
 
+# Calculate divergence in PseudoEvolution
+
+if(track_divB)
+{
+  # bit of a misnomer. Currently only contains divB
+  schedule group GRHydroAnalysis IN MoL_PseudoEvolution
+  {
+  } "Calculate analysis quantities"
+
+  schedule GRHydro_Analysis_Init IN GRHydroAnalysis
+  {
+    LANG: Fortran
+  } "Initialize variables"
+
+  schedule GRHydro_CalcDivB IN GRHydroAnalysis AFTER GRHydro_Analysis_Init
+  {
+    LANG: Fortran
+  } "Calculate divB"
+}