GRHydro: Add basic vector potential support

   Basic cell-centered, algebraic gauge vector potential
   method with place-holders for lorenz gauge.
   Initial Avec constrained to poloidal at the moment.

From: Tanja Bode <tanja.bode@physics.gatech.edu>

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

1 files changed, 16 insertions, 3 deletions

diff --git a/param.ccl b/param.ccl
index 409cad3..6eb77df 100644
--- a/param.ccl
+++ b/param.ccl
@@ -32,7 +32,8 @@ EXTENDS KEYWORD evolution_method ""
 }
 EXTENDS KEYWORD Bvec_evolution_method ""
 {
-  "GRHydro" :: "Use GRHydro to evolve the MHD variables"
+  "GRHydro"      :: "Use GRHydro to evolve the MHD variables, evolving the Magnetic Field"
+  "GRHydro_Avec" :: "Use GRHydro to evolve the MHD variables, evolving the Vector Potential"
 }
 EXTENDS KEYWORD Y_e_evolution_method ""
 {
@@ -46,6 +47,10 @@ EXTENDS KEYWORD entropy_evolution_method ""
 {
   "GRHydro" :: "Use GRHydro to evolve entropy"
 }
+EXTENDS KEYWORD initial_Bvec ""
+{
+  "Bvec_from_Avec" :: "Calculate B^i for an initially specified A_i"
+}
 
 #########################################
 ###  PARAMETRES SHARED FROM MoL       ###
@@ -94,13 +99,13 @@ CCTK_INT GRHydro_hydro_excision "Turns excision automatically on in HydroBase" A
 CCTK_INT GRHydro_MaxNumEvolvedVars "The maximum number of evolved variables used by GRHydro" ACCUMULATOR-BASE=MethodofLines::MoL_Num_Evolved_Vars
 {
   0             :: "when using multirate"
-  5:10		:: "dens scon[3] tau Bvec[3] psidc ye"
+  5:12		:: "dens scon[3] tau (B/A)vec[3] psidc ye entropy Aphi"
 } 5
 
 CCTK_INT GRHydro_MaxNumEvolvedVarsSlow "The maximum number of evolved variables used by GRHydro" ACCUMULATOR-BASE=MethodofLines::MoL_Num_Evolved_Vars_Slow
 {
   0             :: "do not use multirate"
-  5:10		:: "dens scon[3] tau Bvec[3] psidc ye"
+  5:12		:: "dens scon[3] tau (B/A)vec[3] psidc ye entropy Aphi"
 } 0
 
 #  7  primitives (rho,press,eps,w_lorentz,vel)
@@ -657,6 +662,14 @@ boolean calculate_bcom  "Calculate the comoving contravariant magnetic 4-vector
 {
 } "no"
 
+############### Parameters for Avec evolution ###########################################
+
+KEYWORD Avec_gauge "Which gauge condition to use when evolving the vector potential" 
+{
+  "algebraic"  :: "Algebraic gauge"
+  "lorenz"     :: "Lorenz gauge"
+} "lorenz"
+
 ############### Parameters for Tmunu damping (useful in atmosphere) #####################
 
 CCTK_REAL Tmunu_damping_radius_min "damping radius at which we start to damp with a tanh function"