Initial revision

darcs-hash:20010301114010-f6438-12fb8a9ffcc80e86c0a97e37b5b0dae0dbc59b79.gz

8 files changed, 191 insertions, 0 deletions

diff --git a/CarpetExtra/WaveToyExpl/src/finalboundaries.F77 b/CarpetExtra/WaveToyExpl/src/finalboundaries.F77
new file mode 100644
index 000000000..c00960fc3
--- /dev/null
+++ b/CarpetExtra/WaveToyExpl/src/finalboundaries.F77
@@ -0,0 +1,25 @@
+c     $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyExpl/src/finalboundaries.F77,v 1.2 2003/10/27 15:31:40 schnetter Exp $
+
+#include "cctk.h"
+#include "cctk_Arguments.h"
+#include "cctk_Functions.h"
+#include "cctk_Parameters.h"
+
+      subroutine WaveToyExpl_FinalBoundaries (CCTK_ARGUMENTS)
+      implicit none
+      DECLARE_CCTK_ARGUMENTS
+      DECLARE_CCTK_FUNCTIONS
+      DECLARE_CCTK_PARAMETERS
+      CCTK_INT bndwidth
+      parameter (bndwidth = 1)
+      character fbound*100
+      CCTK_INT  fboundlen
+      integer ierr
+      call CCTK_FortranString (fboundlen, bound, fbound)
+      if (fboundlen.lt.0) call CCTK_WARN (0, "internal error")
+      ierr = Boundary_SelectGroupForBC (cctkGH, CCTK_ALL_FACES, bndwidth, -1,
+     $     "WaveToyExpl::scalarevolve", fbound);
+      if (ierr.ne.0) call CCTK_WARN (0, "internal error")
+      call CartSymGN (ierr, cctkGH, "WaveToyExpl::scalarevolve")
+      if (ierr.ne.0) call CCTK_WARN (0, "internal error")
+      end
diff --git a/CarpetExtra/WaveToyExpl/src/finalstep.F77 b/CarpetExtra/WaveToyExpl/src/finalstep.F77
new file mode 100644
index 000000000..b4304d60e
--- /dev/null
+++ b/CarpetExtra/WaveToyExpl/src/finalstep.F77
@@ -0,0 +1,36 @@
+c     $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyExpl/src/finalstep.F77,v 1.2 2003/10/27 15:31:40 schnetter Exp $
+
+#include "cctk.h"
+#include "cctk_Arguments.h"
+#include "cctk_Functions.h"
+#include "cctk_Parameters.h"
+
+      subroutine WaveToyExpl_FinalStep (CCTK_ARGUMENTS)
+      implicit none
+      DECLARE_CCTK_ARGUMENTS
+      DECLARE_CCTK_FUNCTIONS
+      DECLARE_CCTK_PARAMETERS
+      CCTK_REAL dx(3), dt
+      integer i, j, k
+      integer d
+      do d=1,3
+         dx(d) = CCTK_DELTA_SPACE(d)
+      end do
+      dt = CCTK_DELTA_TIME
+      do k=1+cctk_nghostzones(3),cctk_lsh(3)-cctk_nghostzones(3)
+         do j=1+cctk_nghostzones(2),cctk_lsh(2)-cctk_nghostzones(2)
+            do i=1+cctk_nghostzones(1),cctk_lsh(1)-cctk_nghostzones(1)
+c     d/dt phi = psi
+               phi(i,j,k) = phi_p(i,j,k) + dt/2 * (psi_p(i,j,k) + psiint(i,j,k))
+c     d/dt psi = d^2/dx^2 phi
+               psi(i,j,k) = psi_p(i,j,k) + dt/2 *
+     $              ( (phi_p(i-1,j,k) - 2*phi_p(i,j,k) + phi_p(i+1,j,k)) / dx(1)**2
+     $              + (phi_p(i,j-1,k) - 2*phi_p(i,j,k) + phi_p(i,j+1,k)) / dx(2)**2
+     $              + (phi_p(i,j,k-1) - 2*phi_p(i,j,k) + phi_p(i,j,k+1)) / dx(3)**2
+     $              + (phiint(i-1,j,k) - 2*phiint(i,j,k) + phiint(i+1,j,k)) / dx(1)**2
+     $              + (phiint(i,j-1,k) - 2*phiint(i,j,k) + phiint(i,j+1,k)) / dx(2)**2
+     $              + (phiint(i,j,k-1) - 2*phiint(i,j,k) + phiint(i,j,k+1)) / dx(3)**2)
+            end do
+         end do
+      end do
+      end
diff --git a/CarpetExtra/WaveToyExpl/src/initsymbound.F77 b/CarpetExtra/WaveToyExpl/src/initsymbound.F77
new file mode 100644
index 000000000..2f4f182fc
--- /dev/null
+++ b/CarpetExtra/WaveToyExpl/src/initsymbound.F77
@@ -0,0 +1,23 @@
+c     $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyExpl/src/initsymbound.F77,v 1.2 2003/10/27 15:31:40 schnetter Exp $
+
+#include "cctk.h"
+#include "cctk_Arguments.h"
+#include "cctk_Functions.h"
+#include "cctk_Parameters.h"
+
+      subroutine WaveToyExpl_InitSymBound (CCTK_ARGUMENTS)
+      implicit none
+      DECLARE_CCTK_ARGUMENTS
+      DECLARE_CCTK_FUNCTIONS
+      DECLARE_CCTK_PARAMETERS
+      integer sym(3)
+      integer d
+      integer ierr
+      do d=1,3
+         sym(d) = +1
+      end do
+      call SetCartSymGN (ierr, cctkGH, sym, "WaveToyExpl::scalarevolve")
+      if (ierr.ne.0) call CCTK_WARN (0, "internal error")
+      call SetCartSymGN (ierr, cctkGH, sym, "WaveToyExpl::scalarevolveint")
+      if (ierr.ne.0) call CCTK_WARN (0, "internal error")
+      end
diff --git a/CarpetExtra/WaveToyExpl/src/intboundaries.F77 b/CarpetExtra/WaveToyExpl/src/intboundaries.F77
new file mode 100644
index 000000000..838729f73
--- /dev/null
+++ b/CarpetExtra/WaveToyExpl/src/intboundaries.F77
@@ -0,0 +1,25 @@
+c     $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyExpl/src/intboundaries.F77,v 1.3 2003/11/05 16:18:40 schnetter Exp $
+
+#include "cctk.h"
+#include "cctk_Arguments.h"
+#include "cctk_Functions.h"
+#include "cctk_Parameters.h"
+
+      subroutine WaveToyExpl_IntBoundaries (CCTK_ARGUMENTS)
+      implicit none
+      DECLARE_CCTK_ARGUMENTS
+      DECLARE_CCTK_FUNCTIONS
+      DECLARE_CCTK_PARAMETERS
+      CCTK_INT bndwidth
+      parameter (bndwidth = 1)
+      character fbound*100
+      CCTK_INT  fboundlen
+      integer ierr
+      call CCTK_FortranString (fboundlen, bound, fbound)
+      if (fboundlen.lt.0) call CCTK_WARN (0, "internal error")
+      ierr = Boundary_SelectGroupForBC (cctkGH, CCTK_ALL_FACES, bndwidth, -1,
+     $     "WaveToyExpl::scalarevolveint", fbound);
+      if (ierr.ne.0) call CCTK_WARN (0, "internal error")
+      call CartSymGN (ierr, cctkGH, "WaveToyExpl::scalarevolveint")
+      if (ierr.ne.0) call CCTK_WARN (0, "internal error")
+      end
diff --git a/CarpetExtra/WaveToyExpl/src/intinit.F77 b/CarpetExtra/WaveToyExpl/src/intinit.F77
new file mode 100644
index 000000000..b2ebbaea2
--- /dev/null
+++ b/CarpetExtra/WaveToyExpl/src/intinit.F77
@@ -0,0 +1,26 @@
+c     $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyExpl/src/intinit.F77,v 1.2 2003/10/27 15:31:40 schnetter Exp $
+
+#include "cctk.h"
+#include "cctk_Arguments.h"
+#include "cctk_Functions.h"
+#include "cctk_Parameters.h"
+
+      subroutine WaveToyExpl_IntInit (CCTK_ARGUMENTS)
+      implicit none
+      DECLARE_CCTK_ARGUMENTS
+      DECLARE_CCTK_FUNCTIONS
+      DECLARE_CCTK_PARAMETERS
+      integer i, j, k
+      do k=1,cctk_lsh(3)
+         do j=1,cctk_lsh(2)
+            do i=1,cctk_lsh(1)
+               phiint(i,j,k) = phi(i,j,k)
+               phiint_p(i,j,k) = phi_p(i,j,k)
+               phiint_p_p(i,j,k) = phi_p_p(i,j,k)
+               psiint(i,j,k) = psi(i,j,k)
+               psiint_p(i,j,k) = psi_p(i,j,k)
+               psiint_p_p(i,j,k) = psi_p_p(i,j,k)
+            end do
+         end do
+      end do
+      end
diff --git a/CarpetExtra/WaveToyExpl/src/intstep.F77 b/CarpetExtra/WaveToyExpl/src/intstep.F77
new file mode 100644
index 000000000..909842c2a
--- /dev/null
+++ b/CarpetExtra/WaveToyExpl/src/intstep.F77
@@ -0,0 +1,33 @@
+c     $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyExpl/src/intstep.F77,v 1.2 2003/10/27 15:31:40 schnetter Exp $
+
+#include "cctk.h"
+#include "cctk_Arguments.h"
+#include "cctk_Functions.h"
+#include "cctk_Parameters.h"
+
+      subroutine WaveToyExpl_IntStep (CCTK_ARGUMENTS)
+      implicit none
+      DECLARE_CCTK_ARGUMENTS
+      DECLARE_CCTK_FUNCTIONS
+      DECLARE_CCTK_PARAMETERS
+      CCTK_REAL dx(3), dt
+      integer i, j, k
+      integer d
+      do d=1,3
+         dx(d) = CCTK_DELTA_SPACE(d)
+      end do
+      dt = CCTK_DELTA_TIME
+      do k=1+cctk_nghostzones(3),cctk_lsh(3)-cctk_nghostzones(3)
+         do j=1+cctk_nghostzones(2),cctk_lsh(2)-cctk_nghostzones(2)
+            do i=1+cctk_nghostzones(1),cctk_lsh(1)-cctk_nghostzones(1)
+c     d/dt phi = psi
+               phiint(i,j,k) = phi_p(i,j,k) + dt * psi_p(i,j,k)
+c     d/dt psi = d^2/dx^2 phi
+               psiint(i,j,k) = psi_p(i,j,k) + dt *
+     $              ( (phi_p(i-1,j,k) - 2*phi_p(i,j,k) + phi_p(i+1,j,k)) / dx(1)**2
+     $              + (phi_p(i,j-1,k) - 2*phi_p(i,j,k) + phi_p(i,j+1,k)) / dx(2)**2
+     $              + (phi_p(i,j,k-1) - 2*phi_p(i,j,k) + phi_p(i,j,k+1)) / dx(3)**2)
+            end do
+         end do
+      end do
+      end
diff --git a/CarpetExtra/WaveToyExpl/src/make.code.defn b/CarpetExtra/WaveToyExpl/src/make.code.defn
new file mode 100644
index 000000000..d828a9051
--- /dev/null
+++ b/CarpetExtra/WaveToyExpl/src/make.code.defn
@@ -0,0 +1,9 @@
+# Main make.code.defn file for thorn WaveToyExpl
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyExpl/src/make.code.defn,v 1.1 2003/06/18 18:24:30 schnetter Exp $
+
+# Source files in this directory
+SRCS = finalboundaries.F77 finalstep.F77 initsymbound.F77 intboundaries.F77 intinit.F77 intstep.F77 startup.F77
+
+# Subdirectories containing source files
+SUBDIRS = 
+
diff --git a/CarpetExtra/WaveToyExpl/src/startup.F77 b/CarpetExtra/WaveToyExpl/src/startup.F77
new file mode 100644
index 000000000..ce33d9a1c
--- /dev/null
+++ b/CarpetExtra/WaveToyExpl/src/startup.F77
@@ -0,0 +1,14 @@
+c     $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyExpl/src/startup.F77,v 1.2 2003/10/27 15:31:40 schnetter Exp $
+
+#include "cctk.h"
+#include "cctk_Functions.h"
+#include "cctk_Parameters.h"
+
+      subroutine WaveToyExpl_Startup
+      implicit none
+      DECLARE_CCTK_FUNCTIONS
+      DECLARE_CCTK_PARAMETERS
+      integer ierr
+      call CCTK_RegisterBanner (ierr, "WaveToyExpl: Evolving a scalar field with explicit intermediate integration steps")
+      if (ierr.ne.0) call CCTK_WARN (0, "internal error")
+      end