From 985a7b06d2069e7cb281b7c76315a916d01354e4 Mon Sep 17 00:00:00 2001
From: goodale <goodale@7ec00dc9-1e7a-42ad-bf73-a299ee168e72>
Date: Wed, 16 Jan 2002 03:13:38 +0000
Subject: 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
---
 param.ccl          |  5 +++++
 src/InitialData.cc | 17 +++++++++++++++++
 2 files changed, 22 insertions(+)

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;
 }
-- 
cgit v1.2.3