Add noise initial data

darcs-hash:20040507205108-07bb3-666d03915673b2dd38382818cc4dfe13f29c78c1.gz

2 files changed, 18 insertions, 2 deletions

diff --git a/CarpetExtra/IDScalarWaveFO/param.ccl b/CarpetExtra/IDScalarWaveFO/param.ccl
index 58c4c07ca..e6bc7b617 100644
--- a/CarpetExtra/IDScalarWaveFO/param.ccl
+++ b/CarpetExtra/IDScalarWaveFO/param.ccl
@@ -1,11 +1,12 @@
 # Parameter definitions for thorn IDScalarWaveFO
-# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/IDScalarWaveFO/param.ccl,v 1.2 2004/01/11 12:20:39 schnetter Exp $
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/IDScalarWaveFO/param.ccl,v 1.3 2004/05/07 22:51:08 schnetter Exp $
 
 RESTRICTED:
 
 KEYWORD initial_data "Type of initial data"
 {
   "plane" :: "Plane wave"
+  "noise" :: "Noise"
 } "plane"
 
 PRIVATE:
diff --git a/CarpetExtra/IDScalarWaveFO/src/initialdata.F77 b/CarpetExtra/IDScalarWaveFO/src/initialdata.F77
index 9e193be6b..e2808d54f 100644
--- a/CarpetExtra/IDScalarWaveFO/src/initialdata.F77
+++ b/CarpetExtra/IDScalarWaveFO/src/initialdata.F77
@@ -1,4 +1,4 @@
-c     $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/IDScalarWaveFO/src/initialdata.F77,v 1.4 2004/01/11 12:20:39 schnetter Exp $
+c     $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/IDScalarWaveFO/src/initialdata.F77,v 1.5 2004/05/07 22:51:08 schnetter Exp $
 
 #include "cctk.h"
 #include "cctk_Arguments.h"
@@ -14,6 +14,7 @@ c     $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/IDScalarWaveFO/src/ini
       parameter (pi = 3.141592653589793238462643383279502884197169399375105820974944592307816406286208998628034825342117068d0)
       CCTK_REAL omega
       integer i, j, k
+      
       if (CCTK_EQUALS(initial_data, "plane")) then
          omega = sqrt(wave_number(1)**2 + wave_number(2)**2 + wave_number(3)**2)
          do k=1,cctk_lsh(3)
@@ -42,5 +43,19 @@ c     $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/IDScalarWaveFO/src/ini
                end do
             end do
          end do
+         
+      else if (CCTK_EQUALS(initial_data, "noise")) then
+         
+         do k=1,cctk_lsh(3)
+            do j=1,cctk_lsh(2)
+               do i=1,cctk_lsh(1)
+                  call random_number (phi(i,j,k))
+                  call random_number (psix(i,j,k))
+                  call random_number (psiy(i,j,k))
+                  call random_number (psiz(i,j,k))
+               end do
+            end do
+         end do
+         
       end if
       end