Making the initial data consistent between the initial data thorns -

in particular adding "none" to the C and C++ versions.

Tom


git-svn-id: http://svn.cactuscode.org/arrangements/CactusWave/IDScalarWaveCXX/trunk@30 7ec00dc9-1e7a-42ad-bf73-a299ee168e72

2 files changed, 22 insertions, 0 deletions

diff --git a/param.ccl b/param.ccl
index fe6f8e7..3d3c0cb 100644
--- a/param.ccl
+++ b/param.ccl
@@ -12,10 +12,13 @@ KEYWORD initial_data "Type of initial data"
   "plane"      :: "Plane wave"
   "gaussian"   :: "Gaussian wave"
   "box"        :: "Box wave"
+  "none"       :: "No initial data, zero phi"
 } "gaussian"
 
 private:
 
+## Parameter for initial wavepulses
+
 REAL radius "The radius of the gaussian wave"
 {
  0:* :: "Positive"
@@ -30,10 +33,12 @@ REAL kx "The wave number in the x-direction"
 {
  *:* :: "No restriction"
 } 4.0
+
 REAL ky "The wave number in the y-direction"
 {
  *:* :: "No restriction"
 } 0.0
+
 REAL kz "The wave number in the z-direction"
 {
  *:* :: "No restriction"
diff --git a/src/InitialData.cc b/src/InitialData.cc
index 1f659c6..5e457d8 100644
--- a/src/InitialData.cc
+++ b/src/InitialData.cc
@@ -115,6 +115,23 @@ extern "C" void IDScalarWaveCXX_InitialData(CCTK_ARGUMENTS)
       }
     }
   }
+  else if (CCTK_Equals(initial_data, "none"))
+  {
+    for(int k=0; k<cctk_lsh[2]; k++)
+    {
+      for(int j=0; j<cctk_lsh[1]; j++)
+      {
+        for(int i=0; i<cctk_lsh[0]; i++)
+        {
+          int index =  CCTK_GFINDEX3D(cctkGH,i,j,k);
+
+          phi[index] = 0.0;
+
+          phi_p[index] = 0.0;
+        }
+      }
+    }
+  }
 
   return;
 }