Initialize all active timelevels.

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinBase/ADMBase/trunk@23 d576a68a-b34a-40ae-82fc-004fa1a9d16f

1 files changed, 101 insertions, 0 deletions

diff --git a/src/Initialisation.c b/src/Initialisation.c
index 271ef87..ba0bba8 100644
--- a/src/Initialisation.c
+++ b/src/Initialisation.c
@@ -88,6 +88,70 @@ void ADMBase_CartesianMinkowski(CCTK_ARGUMENTS)
     kxz[i] = 0.0;
     kyz[i] = 0.0;
   }
+  
+  if (CCTK_ActiveTimeLevels(cctkGH, "ADMBase::metric") > 1)
+  {
+    for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
+    {
+
+      gxx_p[i] = 1.0;
+      gyy_p[i] = 1.0;
+      gzz_p[i] = 1.0;
+      
+      gxy_p[i] = 0.0;
+      gxz_p[i] = 0.0;
+      gyz_p[i] = 0.0;
+    }
+  }
+  
+  if (CCTK_ActiveTimeLevels(cctkGH, "ADMBase::metric") > 2)
+  {
+    for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
+    {
+
+      gxx_p_p[i] = 1.0;
+      gyy_p_p[i] = 1.0;
+      gzz_p_p[i] = 1.0;
+      
+      gxy_p_p[i] = 0.0;
+      gxz_p_p[i] = 0.0;
+      gyz_p_p[i] = 0.0;
+      
+    }
+    
+  }
+    
+  if (CCTK_ActiveTimeLevels(cctkGH, "ADMBase::curv") > 1)
+  {
+    for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
+    {
+
+      kxx_p[i] = 0.0;
+      kyy_p[i] = 0.0;
+      kzz_p[i] = 0.0;
+      
+      kxy_p[i] = 0.0;
+      kxz_p[i] = 0.0;
+      kyz_p[i] = 0.0;
+    }
+  }
+  
+  if (CCTK_ActiveTimeLevels(cctkGH, "ADMBase::curv") > 2)
+  {
+    for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
+    {
+      
+      kxx_p_p[i] = 0.0;
+      kyy_p_p[i] = 0.0;
+      kzz_p_p[i] = 0.0;
+      
+      kxy_p_p[i] = 0.0;
+      kxz_p_p[i] = 0.0;
+      kyz_p_p[i] = 0.0;
+      
+    }
+    
+  }
 
 }
 
@@ -116,6 +180,23 @@ void ADMBase_LapseOne(CCTK_ARGUMENTS)
   {
     alp[i] = 1.0;
   }
+  
+  if (CCTK_ActiveTimeLevels(cctkGH, "ADMBase::lapse") > 1)
+  {
+    for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
+    {
+      alp_p[i] = 1.0;
+    }
+  }
+  
+  if (CCTK_ActiveTimeLevels(cctkGH, "ADMBase::lapse") > 2)
+  {
+    for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
+    {
+      alp_p_p[i] = 1.0;
+    }
+  }
+  
 }
 
 /*@@
@@ -144,6 +225,26 @@ void ADMBase_ShiftZero(CCTK_ARGUMENTS)
     betay[i] = 0.0;
     betaz[i] = 0.0;
   }
+  
+  if (CCTK_ActiveTimeLevels(cctkGH, "ADMBase::shift") > 1)
+  {
+    for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
+    {
+      betax_p[i] = 0.0;
+      betay_p[i] = 0.0;
+      betaz_p[i] = 0.0;
+    }
+  }
+  
+  if (CCTK_ActiveTimeLevels(cctkGH, "ADMBase::shift") > 2)
+  {
+    for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
+    {
+      betax_p_p[i] = 0.0;
+      betay_p_p[i] = 0.0;
+      betaz_p_p[i] = 0.0;
+    }
+  }
 
 }