Adding Eric's standing plane wave solution.

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDLinearWaves/trunk@76 5c0f84ea-6048-4d6e-bfa4-55cd5f2e0dd7

5 files changed, 168 insertions, 2 deletions

diff --git a/param.ccl b/param.ccl
index ab334be..5793b63 100644
--- a/param.ccl
+++ b/param.ccl
@@ -4,8 +4,9 @@ shares: ADMBase
 
 EXTENDS KEYWORD initial_data 
 {
-"teukwaves"  :: "linear waves initial data- Teukolsky waves"
-"planewaves" :: "linear waves initial data- plane waves"
+"teukwaves"           :: "linear waves initial data- Teukolsky waves"
+"planewaves"          :: "linear waves initial data- plane waves"
+"standing_planewaves" :: "linear waves initial data- standing plane waves"
 }
 
 USES KEYWORD metric_type 
diff --git a/schedule.ccl b/schedule.ccl
index ecb50f6..2debdbe 100644
--- a/schedule.ccl
+++ b/schedule.ccl
@@ -13,6 +13,19 @@ if (CCTK_Equals(initial_data,"planewaves"))
   } "Construct linear planewave initial data"  
 }
 
+if (CCTK_Equals(initial_data,"standing_planewaves")) 
+{
+  schedule IDLinearWaves_ParamChecker at CCTK_PARAMCHECK
+  {
+    LANG: C
+  } "Check that the metric_type is recognised"
+
+  schedule IDLinearWaves_Standing_PlaneWaves in ADMBase_InitialData
+  {
+    LANG: Fortran
+  } "Construct linear planewave initial data"  
+}
+
 if (CCTK_Equals(initial_data,"teukwaves")) 
 {
   schedule IDLinearWaves_ParamChecker at CCTK_PARAMCHECK
diff --git a/src/make.code.defn b/src/make.code.defn
index 9969001..ff57d67 100644
--- a/src/make.code.defn
+++ b/src/make.code.defn
@@ -4,6 +4,7 @@
 # Source files in this directory
 SRCS = 	ParamCheck.c\
 	planewaves.F77\
+	standing_planewaves.F77\
 	teukwaves.F77
 
 # Subdirectories containing source files
diff --git a/src/standing_planewaves.F77 b/src/standing_planewaves.F77
new file mode 100644
index 0000000..c4af49b
--- /dev/null
+++ b/src/standing_planewaves.F77
@@ -0,0 +1,97 @@
+c     Written 2002-05-22 by Erik Schnetter <schnetter@uni-tuebingen.de>
+c     $Header$
+
+#include "cctk.h"
+#include "cctk_Arguments.h"
+#include "cctk_Parameters.h"
+      
+      subroutine IDLinearWaves_Standing_PlaneWaves (CCTK_ARGUMENTS)
+      
+      implicit none
+      
+      DECLARE_CCTK_ARGUMENTS
+      DECLARE_CCTK_FUNCTIONS
+      DECLARE_CCTK_PARAMETERS
+      
+      CCTK_REAL pi
+      parameter (pi = 3.14159265358979d0)
+      CCTK_REAL aa, kk
+      CCTK_REAL zz, tt
+      CCTK_REAL bb, bbdot
+      integer i,j,k
+      
+      aa = amplitude
+      kk = 2*pi / wavelength
+      
+      tt = cctk_time
+      
+      do k=1,cctk_lsh(3)
+         do j=1,cctk_lsh(2)
+            do i=1,cctk_lsh(1)
+               
+               zz = z(i,j,k)
+               
+               bb = aa * sin(kk*zz) * cos(kk*tt)
+               bbdot = -kk*aa * sin(kk*zz) * sin(kk*tt)
+               
+               gxx(i,j,k) = 1 + bb
+               gxy(i,j,k) = 0
+               gxz(i,j,k) = 0
+               gyy(i,j,k) = 1 - bb
+               gyz(i,j,k) = 0
+               gzz(i,j,k) = 1
+               
+               kxx(i,j,k) = bbdot / (-2)
+               kxy(i,j,k) = 0
+               kxz(i,j,k) = 0
+               kyy(i,j,k) = - bbdot / (-2)
+               kyz(i,j,k) = 0
+               kzz(i,j,k) = 0
+               
+            end do 
+         end do
+      end do
+      
+c     initialise the conformal factor
+      
+      if (CCTK_EQUALS(metric_type, "static conformal")) then 
+         
+         conformal_state = 1
+         
+         if (CCTK_EQUALS(conformal_storage,"factor+derivs")) then
+            
+            conformal_state = 2
+            
+         else if (CCTK_EQUALS(conformal_storage,"factor+derivs+2nd derivs")) then
+            conformal_state = 3
+            
+         end if
+         
+         do k=1,cctk_lsh(3)
+            do j=1,cctk_lsh(2)
+               do i=1,cctk_lsh(1)		 
+                  
+                  psi(i,j,k) = 1
+                  
+                  if (conformal_state .gt. 1) then
+                     psix(i,j,k) = 0
+                     psiy(i,j,k) = 0
+                     psiz(i,j,k) = 0
+                  end if
+                  
+                  if (conformal_state .gt. 2) then
+                     psixy(i,j,k) = 0
+                     psixz(i,j,k) = 0
+                     psiyz(i,j,k) = 0
+                     psixx(i,j,k) = 0
+                     psiyy(i,j,k) = 0
+                     psizz(i,j,k) = 0
+                  end if
+                  
+               end do
+            end do
+         end do
+         
+      end if
+      
+      end
diff --git a/test/test_spw_ADM_icn_par b/test/test_spw_ADM_icn_par
new file mode 100644
index 0000000..7aac2f6
--- /dev/null
+++ b/test/test_spw_ADM_icn_par
@@ -0,0 +1,54 @@
+#########################################################
+!DESC "Standing linearised plane wave along y-line, ADM, ICN"
+#########################################################
+
+# Required thorns
+ActiveThorns = "time ADMconstraints pugh cartgrid3d admbase staticconformal admcoupling admmacros boundary coordgauge pughslab spacemask ADM IDLinearWaves
+ ioascii ioutil iobasic pughreduce"
+
+# GENERAL 
+
+driver::global_nx  =  3
+driver::global_ny  =  3
+driver::global_nz  =  200
+time::dtfac    =  0.25   
+
+grid::type         =  "byspacing"
+grid::dxyz         =  0.01
+grid::domain       =  "full"
+grid::avoid_origin = no
+
+cactus::cctk_initial_time  =  0.     
+cactus::cctk_itlast        =  100      
+
+
+# Einstein 
+
+admbase::initial_data     = "standing_planewaves"
+admbase::evolution_method = "ADM"     
+
+adm::method        =  "icn" 
+adm::bound         =  "static"
+admbase::lapse_evolution_method  =  "geodesic"  
+
+# IDLinearwaves
+idlinearwaves::amplitude  = 0.001
+
+# Output
+
+IO::outdir         = "test_spw_ADM_icn"
+IOBasic::outScalar_every =  10      
+IOBasic::outScalar_vars  = "admbase::gzz admbase::kzz 
+	admconstraints::ADMconstraints 
+	admbase::alp"
+IOASCII::out1D_every =  10      
+IOASCII::out1D_vars  = "admbase::gxx admbase::gxy admbase::gxz
+                        admbase::gyy admbase::gyz admbase::gzz
+                        admbase::kxx admbase::kxy admbase::kxz
+                        admbase::kyy admbase::kyz admbase::kzz"
+
+
+IOBasic::outInfo_every = 5
+IOBasic::outInfo_vars = "admconstraints::ham admbase::gzz"
+
+##################################################################