Add gridfunctions and boilerplate for vector potential.

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinBase/HydroBase/trunk@59 57fe0bb3-ccba-405f-9b23-de0201f165b7

1 files changed, 40 insertions, 0 deletions

diff --git a/schedule.ccl b/schedule.ccl
index 2e8207e..3389936 100755
--- a/schedule.ccl
+++ b/schedule.ccl
@@ -15,6 +15,14 @@ if (timelevels == 1)
   {
     STORAGE: Bvec[1]
   }
+  if (!CCTK_EQUALS(initial_Avec, "none"))
+  {
+    STORAGE: Avec[1]
+  }
+  if (!CCTK_EQUALS(initial_Aphi, "none"))
+  {
+    STORAGE: Aphi[1]
+  }
   if (!CCTK_EQUALS(initial_temperature, "none"))
   {
     STORAGE: temperature[1]
@@ -39,6 +47,14 @@ else if (timelevels == 2)
   {
     STORAGE: Bvec[2]
   }
+  if (!CCTK_EQUALS(initial_Avec, "none"))
+  {
+    STORAGE: Avec[2]
+  }
+  if (!CCTK_EQUALS(initial_Aphi, "none"))
+  {
+    STORAGE: Aphi[2]
+  }
   if (!CCTK_EQUALS(initial_temperature, "none"))
   {
     STORAGE: temperature[2]
@@ -63,6 +79,14 @@ else if (timelevels == 3)
   {
     STORAGE: Bvec[3]
   }
+  if (!CCTK_EQUALS(initial_Avec, "none"))
+  {
+    STORAGE: Avec[3]
+  }
+  if (!CCTK_EQUALS(initial_Aphi, "none"))
+  {
+    STORAGE: Aphi[3]
+  }
   if (!CCTK_EQUALS(initial_temperature, "none"))
   {
     STORAGE: temperature[3]
@@ -166,6 +190,22 @@ if (CCTK_Equals(initial_Bvec, "zero"))
   } "Set magnetic field to 0"
 }
 
+if (CCTK_Equals(initial_Avec, "zero"))
+{
+  SCHEDULE HydroBase_Avec_zero in HydroBase_Initial
+  {
+    LANG: C
+  } "Set vector potential to 0"
+}
+
+if (CCTK_Equals(initial_Aphi, "zero"))
+{
+  SCHEDULE HydroBase_Aphi_zero in HydroBase_Initial
+  {
+    LANG: C
+  } "Set vector potential Phi to 0"
+}
+
 if (hydro_excision)
 {
   schedule GROUP HydroBase_ExcisionMaskSetup in HydroBase_Initial