Initial revision

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

22 files changed, 1318 insertions, 0 deletions

diff --git a/CarpetExtra/WaveToyMoL/README b/CarpetExtra/WaveToyMoL/README
new file mode 100644
index 000000000..1cee06d63
--- /dev/null
+++ b/CarpetExtra/WaveToyMoL/README
@@ -0,0 +1,10 @@
+CVS info   : $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyMoL/README,v 1.1 2003/06/18 18:24:30 schnetter Exp $
+
+Cactus Code Thorn WaveToyMoL
+Thorn Author(s)     : Erik Schnetter <schnetter@uni-tuebingen.de>
+Thorn Maintainer(s) : Erik Schnetter <schnetter@uni-tuebingen.de>
+--------------------------------------------------------------------------
+
+Purpose of the thorn:
+
+Implement WaveToy using MoL as time integrator.
diff --git a/CarpetExtra/WaveToyMoL/configuration.ccl b/CarpetExtra/WaveToyMoL/configuration.ccl
new file mode 100644
index 000000000..f3c7db05c
--- /dev/null
+++ b/CarpetExtra/WaveToyMoL/configuration.ccl
@@ -0,0 +1,4 @@
+# Configuration definition for thorn Carpet/WaveToyMoL
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyMoL/configuration.ccl,v 1.2 2004/04/18 14:10:29 schnetter Exp $
+
+REQUIRES CartGrid3D
diff --git a/CarpetExtra/WaveToyMoL/doc/documentation.tex b/CarpetExtra/WaveToyMoL/doc/documentation.tex
new file mode 100644
index 000000000..0d275c478
--- /dev/null
+++ b/CarpetExtra/WaveToyMoL/doc/documentation.tex
@@ -0,0 +1,144 @@
+% *======================================================================*
+%  Cactus Thorn template for ThornGuide documentation
+%  Author: Ian Kelley
+%  Date: Sun Jun 02, 2002
+%  $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyMoL/doc/documentation.tex,v 1.1 2003/06/18 18:24:30 schnetter Exp $                                                             
+%
+%  Thorn documentation in the latex file doc/documentation.tex 
+%  will be included in ThornGuides built with the Cactus make system.
+%  The scripts employed by the make system automatically include 
+%  pages about variables, parameters and scheduling parsed from the 
+%  relevent thorn CCL files.
+%  
+%  This template contains guidelines which help to assure that your     
+%  documentation will be correctly added to ThornGuides. More 
+%  information is available in the Cactus UsersGuide.
+%                                                    
+%  Guidelines:
+%   - Do not change anything before the line
+%       % START CACTUS THORNGUIDE",
+%     except for filling in the title, author, date etc. fields.
+%        - Each of these fields should only be on ONE line.
+%        - Author names should be sparated with a \\ or a comma
+%   - You can define your own macros, but they must appear after
+%     the START CACTUS THORNGUIDE line, and must not redefine standard 
+%     latex commands.
+%   - To avoid name clashes with other thorns, 'labels', 'citations', 
+%     'references', and 'image' names should conform to the following 
+%     convention:          
+%       ARRANGEMENT_THORN_LABEL
+%     For example, an image wave.eps in the arrangement CactusWave and 
+%     thorn WaveToyC should be renamed to CactusWave_WaveToyC_wave.eps
+%   - Graphics should only be included using the graphix package. 
+%     More specifically, with the "includegraphics" command. Do
+%     not specify any graphic file extensions in your .tex file. This 
+%     will allow us (later) to create a PDF version of the ThornGuide
+%     via pdflatex. |
+%   - References should be included with the latex "bibitem" command. 
+%   - Use \begin{abstract}...\end{abstract} instead of \abstract{...}
+%   - Do not use \appendix, instead include any appendices you need as 
+%     standard sections. 
+%   - For the benefit of our Perl scripts, and for future extensions, 
+%     please use simple latex.     
+%
+% *======================================================================* 
+% 
+% Example of including a graphic image:
+%    \begin{figure}[ht]
+% 	\begin{center}
+%    	   \includegraphics[width=6cm]{MyArrangement_MyThorn_MyFigure}
+% 	\end{center}
+% 	\caption{Illustration of this and that}
+% 	\label{MyArrangement_MyThorn_MyLabel}
+%    \end{figure}
+%
+% Example of using a label:
+%   \label{MyArrangement_MyThorn_MyLabel}
+%
+% Example of a citation:
+%    \cite{MyArrangement_MyThorn_Author99}
+%
+% Example of including a reference
+%   \bibitem{MyArrangement_MyThorn_Author99}
+%   {J. Author, {\em The Title of the Book, Journal, or periodical}, 1 (1999), 
+%   1--16. {\tt http://www.nowhere.com/}}
+%
+% *======================================================================* 
+
+% If you are using CVS use this line to give version information
+% $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyMoL/doc/documentation.tex,v 1.1 2003/06/18 18:24:30 schnetter Exp $
+
+\documentclass{article}
+
+% Use the Cactus ThornGuide style file
+% (Automatically used from Cactus distribution, if you have a 
+%  thorn without the Cactus Flesh download this from the Cactus
+%  homepage at www.cactuscode.org)
+\usepackage{../../../../doc/latex/cactus}
+
+\begin{document}
+
+% The author of the documentation
+\author{Erik Schnetter \textless schnetter@uni-tuebingen.de\textgreater}
+
+% The title of the document (not necessarily the name of the Thorn)
+\title{}
+
+% the date your document was last changed, if your document is in CVS, 
+% please use:
+%    \date{$ $Date: 2003/06/18 18:24:30 $ $}
+\date{June 07 2003}
+
+\maketitle
+
+% Do not delete next line
+% START CACTUS THORNGUIDE
+
+% Add all definitions used in this documentation here 
+%   \def\mydef etc
+
+% Add an abstract for this thorn's documentation
+\begin{abstract}
+
+\end{abstract}
+
+% The following sections are suggestive only.
+% Remove them or add your own.
+
+\section{Introduction}
+
+\section{Physical System}
+
+\section{Numerical Implementation}
+
+\section{Using This Thorn}
+
+\subsection{Obtaining This Thorn}
+
+\subsection{Basic Usage}
+
+\subsection{Special Behaviour}
+
+\subsection{Interaction With Other Thorns}
+
+\subsection{Examples}
+
+\subsection{Support and Feedback}
+
+\section{History}
+
+\subsection{Thorn Source Code}
+
+\subsection{Thorn Documentation}
+
+\subsection{Acknowledgements}
+
+
+\begin{thebibliography}{9}
+
+\end{thebibliography}
+
+% Do not delete next line
+% END CACTUS THORNGUIDE
+
+\end{document}
diff --git a/CarpetExtra/WaveToyMoL/interface.ccl b/CarpetExtra/WaveToyMoL/interface.ccl
new file mode 100644
index 000000000..65e4cb5a0
--- /dev/null
+++ b/CarpetExtra/WaveToyMoL/interface.ccl
@@ -0,0 +1,42 @@
+# Interface definition for thorn WaveToyMoL
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyMoL/interface.ccl,v 1.4 2004/05/21 17:43:57 schnetter Exp $
+
+IMPLEMENTS: WaveToyMoL
+
+INHERITS: grid
+
+
+
+CCTK_INT FUNCTION Boundary_SelectGroupForBC(CCTK_POINTER_TO_CONST IN GH, \
+  CCTK_INT IN faces, CCTK_INT IN boundary_width, CCTK_INT IN table_handle, \
+  CCTK_STRING IN group_name, CCTK_STRING IN bc_name)
+
+REQUIRES FUNCTION Boundary_SelectGroupForBC
+
+
+
+CCTK_INT FUNCTION MoLRegisterEvolvedGroup(CCTK_INT IN EvolvedIndex, \
+                                          CCTK_INT IN RHSIndex)
+
+REQUIRES FUNCTION MoLRegisterEvolvedGroup
+
+
+
+PUBLIC:
+
+CCTK_REAL scalarevolve TYPE=gf TIMELEVELS=3
+{
+  phi psi
+} "Scalar field"
+
+PRIVATE:
+
+CCTK_REAL scalarevolvedot TYPE=gf
+{
+  phidot psidot
+} "Time derivative of the scalar field"
+
+CCTK_REAL scalarevolveerrorestimate TYPE=gf
+{
+  error
+} "Error estimate"
diff --git a/CarpetExtra/WaveToyMoL/par/sw1d.gnu b/CarpetExtra/WaveToyMoL/par/sw1d.gnu
new file mode 100644
index 000000000..834756543
--- /dev/null
+++ b/CarpetExtra/WaveToyMoL/par/sw1d.gnu
@@ -0,0 +1,55 @@
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyMoL/par/sw1d.gnu,v 1.2 2003/06/30 17:28:51 schnetter Exp $
+
+set grid
+
+dt=0.05
+f(x)=cos(2*pi*(x+t*dt))
+
+t=200
+
+
+p [0:1] "sw1d_1l_0020/phi.xl" i t u 10:13 w l, "sw1d_1l_0040/phi.xl" i t u 10:13 w l, "sw1d_1l_0080/phi.xl" i t u 10:13 w l, f(x)
+
+p [0:1] "sw1d_1l_0020/phi.xl" i t u 10:($13-f($10)) w l, "sw1d_1l_0040/phi.xl" i t u 10:(4*($13-f($10))) w l, "sw1d_1l_0080/phi.xl" i t u 10:(16*($13-f($10))) w l
+
+
+
+p [0:1] "sw1d_2l_0020/phi.xl" i 2*t u 10:13 w l, "sw1d_2l_0040/phi.xl" i 2*t u 10:13 w l, "sw1d_2l_0080/phi.xl" i 2*t u 10:13 w l, "sw1d_2l_0160/phi.xl" i 2*t u 10:13 w l, "sw1d_2l_0320/phi.xl" i 2*t u 10:13 w l, f(x)
+
+p [0:1] "sw1d_2l_0020/phi.xl" i 2*t u 10:($13-f($10)) w l, "sw1d_2l_0040/phi.xl" i 2*t u 10:(4*($13-f($10))) w l, "sw1d_2l_0080/phi.xl" i 2*t u 10:(16*($13-f($10))) w l, "sw1d_2l_0160/phi.xl" i 2*t u 10:(64*($13-f($10))) w l, "sw1d_2l_0320/phi.xl" i 2*t u 10:(256*($13-f($10))) w l
+
+
+
+
+
+set grid
+
+
+dt=0.1
+A=1.0/16.0
+
+f(x)=A*t*dt*sin(2*pi*x)*cos(2*pi*dt*t)
+
+
+
+t=400
+
+p [0:1] "sw1d_1l_0020/phierror.xl" i t u 10:13 w lp, "sw1d_1l_0040/phierror.xl" i 2*t u 10:(4*$13) w lp, "sw1d_1l_0080/phierror.xl" i 4*t u 10:(16*$13) w lp
+
+
+
+p [0:1] "sw1d_2l_0020/phi.xl" i 3*t ev :::2::2 u 10:13 w l, "sw1d_2l_0040/phi.xl" i 3*2*t ev :::2::2 u 10:13 w l, "sw1d_2l_0080/phi.xl" i 3*4*t ev :::2::2 u 10:13 w l, "sw1d_2l_0160/phi.xl" i 3*8*t ev :::2::2 u 10:13 w l, "sw1d_2l_0320/phi.xl" i 3*16*t ev :::2::2 u 10:13 w l, cos(2*pi*x)
+p [0:1] "sw1d_2l_20/phierror.xl" i 3*t u 10:13 w l, "sw1d_2l_40/phierror.xl" i 3*2*t u 10:(4*$13) w l, "sw1d_2l_80/phierror.xl" i 3*4*t u 10:(16*$13) w l, "sw1d_2l_160/phierror.xl" i 3*8*t u 10:(64*$13) w l, "sw1d_2l_320/phierror.xl" i 3*16*t u 10:(256*$13) w l
+
+p [0:1] "sw1d_2l_0020/phi.xl" i 2*t u 10:13 w l, "sw1d_2l_0040/phi.xl" i 2*t u 10:13 w l, "sw1d_2l_0080/phi.xl" i 2*t u 10:13 w l, "sw1d_2l_0160/phi.xl" i 2*t u 10:13 w l, "sw1d_2l_0320/phi.xl" i 2*t u 10:13 w l, cos(2*pi*x)
+p [0:1] "sw1d_2l_0020/phierror.xl" i 2*t u 10:13 w l, "sw1d_2l_0040/phierror.xl" i 2*t u 10:(4*$13) w l, "sw1d_2l_0080/phierror.xl" i 2*t u 10:(16*$13) w l, "sw1d_2l_0160/phierror.xl" i 2*t u 10:(64*$13) w l, "sw1d_2l_0320/phierror.xl" i 2*t u 10:(256*$13) w l
+p [0:1] "sw1d_2l_0020/phierror.xl" i 2*t u 10:($13-f($10)/20) w l, "sw1d_2l_0040/phierror.xl" i 2*t u 10:(4*($13-f($10)/40)) w l, "sw1d_2l_0080/phierror.xl" i 2*t u 10:(16*($13-f($10)/80)) w l, "sw1d_2l_0160/phierror.xl" i 2*t u 10:(64*($13-f($10)/160)) w l, "sw1d_2l_0320/phierror.xl" i 2*t u 10:(256*($13-f($10)/320)) w l
+
+
+
+p [:] "< tail +6 sw1d_1l_0020/phierror_norm2.xg" t "unigrid error, dx=1/20" w l, "< tail +6 sw1d_1l_0040/phierror_norm2.xg" u 1:(4*$2) t "4 * unigrid error, dx=1/40" w l, "< tail +6 sw1d_2l_0020/phierror_norm2.xg" t "FMR error, dx=1/20" w l, "< tail +6 sw1d_2l_0040/phierror_norm2.xg" u 1:(4*$2) t "4 * FMR error, dx=1/40" w l, "< tail +6 sw1d_2l_0080/phierror_norm2.xg" u 1:(16*$2) t "16 * FMR error, dx=1/80" w l, "< tail +6 sw1d_2l_0160/phierror_norm2.xg" u 1:(64*$2) t "64 * FMR error, dx=1/160" w l, "< tail +6 sw1d_2l_0320/phierror_norm2.xg" u 1:(256*$2) t "256 * FMR error, dx=1/320" w l
+, "< tail +6 sw1d_1l_0080/phierror_norm2.xg" u 1:(16*$2) t "16 * unigrid error, dx=1/80" w l
+
+p [:5] "< tail +6 sw1d_1l_20/phierror_norm2.xg" t "unigrid error, dx=1/20" w l, "< tail +6 sw1d_1l_40/phierror_norm2.xg" u 1:(4*$2) t "4 * unigrid error, dx=1/40" w l, "< tail +6 sw1d_1l_80/phierror_norm2.xg" u 1:(16*$2) t "16 * unigrid error, dx=1/80" w l, "< tail +6 sw1d_2l_20/phierror_norm2.xg" ev 2 t "FMR error, dx=1/20" w l, "< tail +6 sw1d_2l_40/phierror_norm2.xg" ev 2 u 1:(4*$2) t "4 * FMR error, dx=1/40" w l, "< tail +6 sw1d_2l_80/phierror_norm2.xg" ev 2 u 1:(16*$2) t "16 * FMR error, dx=1/80" w l, "< tail +6 sw1d_2l_160/phierror_norm2.xg" ev 2 u 1:(64*$2) t "64 * FMR error, dx=1/160" w l, "< tail +6 sw1d_2l_320/phierror_norm2.xg" ev 2 u 1:(256*$2) t "256 * FMR error, dx=1/320" w l
+
+p [:5] "< tail +6 sw1d_1l_20/psierror_norm2.xg" w l, "< tail +6 sw1d_1l_40/psierror_norm2.xg" u 1:(4*$2) w l, "< tail +6 sw1d_1l_80/psierror_norm2.xg" u 1:(16*$2) w l, "< tail +6 sw1d_2l_20/psierror_norm2.xg" ev 2 w l, "< tail +6 sw1d_2l_40/psierror_norm2.xg" ev 2 u 1:(4*$2) w l, "< tail +6 sw1d_2l_80/psierror_norm2.xg" ev 2 u 1:(16*$2) w l, "< tail +6 sw1d_2l_160/psierror_norm2.xg" ev 2 u 1:(64*$2) w l
diff --git a/CarpetExtra/WaveToyMoL/par/sw1d_1l_0020.par b/CarpetExtra/WaveToyMoL/par/sw1d_1l_0020.par
new file mode 100644
index 000000000..7cd2fb167
--- /dev/null
+++ b/CarpetExtra/WaveToyMoL/par/sw1d_1l_0020.par
@@ -0,0 +1,99 @@
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyMoL/par/sw1d_1l_0020.par,v 1.4 2004/07/07 16:15:31 tradke Exp $
+
+Cactus::cctk_initial_time = 0
+Cactus::cctk_itlast = 3200
+
+
+
+ActiveThorns = "IOUtil Carpet CarpetLib CarpetRegrid CarpetInterp CarpetReduce CarpetSlab"
+
+driver::global_nx = 24
+driver::global_ny = 5
+driver::global_nz = 5
+driver::ghost_size = 2
+
+#Carpet::verbose = yes
+#Carpet::veryverbose = yes
+Carpet::max_refinement_levels = 3
+Carpet::prolongation_order_space = 3
+Carpet::prolongation_order_time = 2
+#Carpet::init_3_timelevels = yes
+Carpet::init_each_timelevel = yes
+
+#CarpetLib::output_bboxes = yes
+
+CarpetRegrid::refinement_levels = 1
+CarpetRegrid::refined_regions = manual-gridpoint-list
+CarpetRegrid::gridpoints = "[[ ([28,4,4]:[68,12,12]:[2,2,2]) ]]"
+CarpetRegrid::outerbounds = "[[ [[0,0],[1,1],[1,1]] ]]"
+
+
+
+ActiveThorns = "LocalInterp NaNChecker Slab"
+
+NaNChecker::check_every = 16
+NaNChecker::check_vars = "all"
+NaNChecker::action_if_found = terminate
+
+
+
+ActiveThorns = "Boundary CoordBase SymBase CartGrid3D Periodic"
+
+grid::domain = full
+grid::type = byrange
+grid::xmin = -0.1
+grid::ymin = -0.15
+grid::zmin = -0.15
+grid::xmax = +1.05
+grid::ymax = +0.15
+grid::zmax = +0.15
+
+Periodic::periodic = yes
+
+
+
+ActiveThorns = "Time MoL"
+
+Time::dtfac = 0.25
+
+MoL::ODE_Method = ICN
+
+#MoL::ODE_Method = RK2
+#MoL::MoL_Intermediate_Steps = 2
+
+
+
+ActiveThorns = "IDScalarWaveMoL WaveToyMoL"
+
+IDScalarWaveMoL::wave_number[0] = 1.0
+
+
+
+#ActiveThorns = "Dissipation"
+#
+#Dissipation::epsdis = 0.1
+#Dissipation::timelevel = "n+1"
+#Dissipation::vars = "WaveToyMoL::scalarevolve"
+
+
+
+ActiveThorns = "IOBasic CarpetIOASCII"
+
+IO::out_dir = "sw1d_1l_0020"
+
+IOBasic::outInfo_every = 16
+IOBasic::outInfo_vars = "WaveToyMoL::phi IDScalarWaveMoL::phierror"
+ 
+IOBasic::outScalar_every = 16
+IOBasic::outScalar_reductions = "minimum maximum norm1 norm2 norm_inf"
+IOBasic::outScalar_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
+
+#IOASCII::output_all_timelevels = yes
+
+IOASCII::out1D_every = 16
+IO::out_yline_x = 0.5
+IO::out_zline_x = 0.5
+IOASCII::out1D_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
+
+#IOASCII::out3D_every = 16
+#IOASCII::out3D_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
diff --git a/CarpetExtra/WaveToyMoL/par/sw1d_1l_0040.par b/CarpetExtra/WaveToyMoL/par/sw1d_1l_0040.par
new file mode 100644
index 000000000..9c9f93248
--- /dev/null
+++ b/CarpetExtra/WaveToyMoL/par/sw1d_1l_0040.par
@@ -0,0 +1,99 @@
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyMoL/par/sw1d_1l_0040.par,v 1.4 2004/07/07 16:15:31 tradke Exp $
+
+Cactus::cctk_initial_time = 0
+Cactus::cctk_itlast = 6400
+
+
+
+ActiveThorns = "IOUtil Carpet CarpetLib CarpetRegrid CarpetInterp CarpetReduce CarpetSlab"
+
+driver::global_nx = 44
+driver::global_ny = 5
+driver::global_nz = 5
+driver::ghost_size = 2
+
+#Carpet::verbose = yes
+#Carpet::veryverbose = yes
+Carpet::max_refinement_levels = 3
+Carpet::prolongation_order_space = 3
+Carpet::prolongation_order_time = 2
+#Carpet::init_3_timelevels = yes
+Carpet::init_each_timelevel = yes
+
+#CarpetLib::output_bboxes = yes
+
+CarpetRegrid::refinement_levels = 1
+CarpetRegrid::refined_regions = manual-gridpoint-list
+CarpetRegrid::gridpoints = "[[ ([48,4,4]:[128,12,12]:[2,2,2]) ]]"
+CarpetRegrid::outerbounds = "[[ [[0,0],[1,1],[1,1]] ]]"
+
+
+
+ActiveThorns = "LocalInterp NaNChecker Slab"
+
+NaNChecker::check_every = 32
+NaNChecker::check_vars = "all"
+NaNChecker::action_if_found = terminate
+
+
+
+ActiveThorns = "Boundary CoordBase SymBase CartGrid3D Periodic"
+
+grid::domain = full
+grid::type = byrange
+grid::xmin = -0.05
+grid::ymin = -0.05
+grid::zmin = -0.05
+grid::xmax = +1.025
+grid::ymax = +0.05
+grid::zmax = +0.05
+
+Periodic::periodic = yes
+
+
+
+ActiveThorns = "Time MoL"
+
+Time::dtfac = 0.25
+
+MoL::ODE_Method = ICN
+
+#MoL::ODE_Method = RK2
+#MoL::MoL_Intermediate_Steps = 2
+
+
+
+ActiveThorns = "IDScalarWaveMoL WaveToyMoL"
+
+IDScalarWaveMoL::wave_number[0] = 1.0
+
+
+
+#ActiveThorns = "Dissipation"
+#
+#Dissipation::epsdis = 0.1
+#Dissipation::timelevel = "n+1"
+#Dissipation::vars = "WaveToyMoL::scalarevolve"
+
+
+
+ActiveThorns = "IOBasic CarpetIOASCII"
+
+IO::out_dir = "sw1d_1l_0040"
+
+IOBasic::outInfo_every = 32
+IOBasic::outInfo_vars = "WaveToyMoL::phi IDScalarWaveMoL::phierror"
+ 
+IOBasic::outScalar_every = 32
+IOBasic::outScalar_reductions = "minimum maximum norm1 norm2 norm_inf"
+IOBasic::outScalar_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
+
+#IOASCII::output_all_timelevels = yes
+
+IOASCII::out1D_every = 32
+IO::out_yline_x = 0.5
+IO::out_zline_x = 0.5
+IOASCII::out1D_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
+
+#IOASCII::out3D_every = 32
+#IOASCII::out3D_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
diff --git a/CarpetExtra/WaveToyMoL/par/sw1d_1l_0080.par b/CarpetExtra/WaveToyMoL/par/sw1d_1l_0080.par
new file mode 100644
index 000000000..c8a040c7c
--- /dev/null
+++ b/CarpetExtra/WaveToyMoL/par/sw1d_1l_0080.par
@@ -0,0 +1,99 @@
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyMoL/par/sw1d_1l_0080.par,v 1.4 2004/07/07 16:15:31 tradke Exp $
+
+Cactus::cctk_initial_time = 0
+Cactus::cctk_itlast = 12800
+
+
+
+ActiveThorns = "IOUtil Carpet CarpetLib CarpetRegrid CarpetInterp CarpetReduce CarpetSlab"
+
+driver::global_nx = 84
+driver::global_ny = 5
+driver::global_nz = 5
+driver::ghost_size = 2
+
+#Carpet::verbose = yes
+#Carpet::veryverbose = yes
+Carpet::max_refinement_levels = 3
+Carpet::prolongation_order_space = 3
+Carpet::prolongation_order_time = 2
+#Carpet::init_3_timelevels = yes
+Carpet::init_each_timelevel = yes
+
+#CarpetLib::output_bboxes = yes
+
+CarpetRegrid::refinement_levels = 1
+CarpetRegrid::refined_regions = manual-gridpoint-list
+CarpetRegrid::gridpoints = "[[ ([88,4,4]:[248,12,12]:[2,2,2]) ]]"
+CarpetRegrid::outerbounds = "[[ [[0,0],[1,1],[1,1]] ]]"
+
+
+
+ActiveThorns = "LocalInterp NaNChecker Slab"
+
+NaNChecker::check_every = 64
+NaNChecker::check_vars = "all"
+NaNChecker::action_if_found = terminate
+
+
+
+ActiveThorns = "Boundary CoordBase SymBase CartGrid3D Periodic"
+
+grid::domain = full
+grid::type = byrange
+grid::xmin = -0.025
+grid::ymin = -0.025
+grid::zmin = -0.025
+grid::xmax = +1.0125
+grid::ymax = +0.025
+grid::zmax = +0.025
+
+Periodic::periodic = yes
+
+
+
+ActiveThorns = "Time MoL"
+
+Time::dtfac = 0.25
+
+MoL::ODE_Method = ICN
+
+#MoL::ODE_Method = RK2
+#MoL::MoL_Intermediate_Steps = 2
+
+
+
+ActiveThorns = "IDScalarWaveMoL WaveToyMoL"
+
+IDScalarWaveMoL::wave_number[0] = 1.0
+
+
+
+#ActiveThorns = "Dissipation"
+#
+#Dissipation::epsdis = 0.1
+#Dissipation::timelevel = "n+1"
+#Dissipation::vars = "WaveToyMoL::scalarevolve"
+
+
+
+ActiveThorns = "IOBasic CarpetIOASCII"
+
+IO::out_dir = "sw1d_1l_0080"
+
+IOBasic::outInfo_every = 64
+IOBasic::outInfo_vars = "WaveToyMoL::phi IDScalarWaveMoL::phierror"
+ 
+IOBasic::outScalar_every = 64
+IOBasic::outScalar_reductions = "minimum maximum norm1 norm2 norm_inf"
+IOBasic::outScalar_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
+
+#IOASCII::output_all_timelevels = yes
+
+IOASCII::out1D_every = 64
+IO::out_yline_x = 0.5
+IO::out_zline_x = 0.5
+IOASCII::out1D_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
+
+#IOASCII::out3D_every = 64
+#IOASCII::out3D_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
diff --git a/CarpetExtra/WaveToyMoL/par/sw1d_2l_0020.par b/CarpetExtra/WaveToyMoL/par/sw1d_2l_0020.par
new file mode 100644
index 000000000..13f2de3ee
--- /dev/null
+++ b/CarpetExtra/WaveToyMoL/par/sw1d_2l_0020.par
@@ -0,0 +1,99 @@
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyMoL/par/sw1d_2l_0020.par,v 1.5 2004/07/07 16:15:31 tradke Exp $
+
+Cactus::cctk_initial_time = 0
+Cactus::cctk_itlast = 3200
+
+
+
+ActiveThorns = "IOUtil Carpet CarpetLib CarpetRegrid CarpetInterp CarpetReduce CarpetSlab"
+
+driver::global_nx = 24
+driver::global_ny = 5
+driver::global_nz = 5
+driver::ghost_size = 2
+
+#Carpet::verbose = yes
+#Carpet::veryverbose = yes
+Carpet::max_refinement_levels = 3
+Carpet::prolongation_order_space = 3
+Carpet::prolongation_order_time = 2
+#Carpet::init_3_timelevels = yes
+Carpet::init_each_timelevel = yes
+
+#CarpetLib::output_bboxes = yes
+
+CarpetRegrid::refinement_levels = 2
+CarpetRegrid::refined_regions = manual-gridpoint-list
+CarpetRegrid::gridpoints = "[[ ([30,4,4]:[66,12,12]:[2,2,2]) ]]"
+CarpetRegrid::outerbounds = "[[ [[0,0],[1,1],[1,1]] ]]"
+
+
+
+ActiveThorns = "LocalInterp NaNChecker Slab"
+
+NaNChecker::check_every = 16
+NaNChecker::check_vars = "all"
+NaNChecker::action_if_found = terminate
+
+
+
+ActiveThorns = "Boundary CoordBase SymBase CartGrid3D Periodic"
+
+grid::domain = full
+grid::type = byrange
+grid::xmin = -0.1
+grid::ymin = -0.1
+grid::zmin = -0.1
+grid::xmax = +1.05
+grid::ymax = +0.1
+grid::zmax = +0.1
+
+Periodic::periodic = yes
+
+
+
+ActiveThorns = "Time MoL"
+
+Time::dtfac = 0.25
+
+MoL::ODE_Method = ICN
+
+#MoL::ODE_Method = RK2
+#MoL::MoL_Intermediate_Steps = 2
+
+
+
+ActiveThorns = "IDScalarWaveMoL WaveToyMoL"
+
+IDScalarWaveMoL::wave_number[0] = 1.0
+
+
+
+#ActiveThorns = "Dissipation"
+#
+#Dissipation::epsdis = 0.1
+#Dissipation::timelevel = "n+1"
+#Dissipation::vars = "WaveToyMoL::scalarevolve"
+
+
+
+ActiveThorns = "IOBasic CarpetIOASCII"
+
+IO::out_dir = "sw1d_2l_0020"
+
+IOBasic::outInfo_every = 16
+IOBasic::outInfo_vars = "WaveToyMoL::phi IDScalarWaveMoL::phierror"
+ 
+IOBasic::outScalar_every = 16
+IOBasic::outScalar_reductions = "minimum maximum norm1 norm2 norm_inf"
+IOBasic::outScalar_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
+
+#IOASCII::output_all_timelevels = yes
+
+IOASCII::out1D_every = 16
+IO::out_yline_x = 0.5
+IO::out_zline_x = 0.5
+IOASCII::out1D_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
+
+#IOASCII::out3D_every = 16
+#IOASCII::out3D_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
diff --git a/CarpetExtra/WaveToyMoL/par/sw1d_2l_0040.par b/CarpetExtra/WaveToyMoL/par/sw1d_2l_0040.par
new file mode 100644
index 000000000..5f215164d
--- /dev/null
+++ b/CarpetExtra/WaveToyMoL/par/sw1d_2l_0040.par
@@ -0,0 +1,99 @@
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyMoL/par/sw1d_2l_0040.par,v 1.6 2004/07/07 16:15:31 tradke Exp $
+
+Cactus::cctk_initial_time = 0
+Cactus::cctk_itlast = 6400
+
+
+
+ActiveThorns = "IOUtil Carpet CarpetLib CarpetRegrid CarpetInterp CarpetReduce CarpetSlab"
+
+driver::global_nx = 44
+driver::global_ny = 5
+driver::global_nz = 5
+driver::ghost_size = 2
+
+#Carpet::verbose = yes
+#Carpet::veryverbose = yes
+Carpet::max_refinement_levels = 3
+Carpet::prolongation_order_space = 3
+Carpet::prolongation_order_time = 2
+#Carpet::init_3_timelevels = yes
+Carpet::init_each_timelevel = yes
+
+#CarpetLib::output_bboxes = yes
+
+CarpetRegrid::refinement_levels = 2
+CarpetRegrid::refined_regions = manual-gridpoint-list
+CarpetRegrid::gridpoints = "[[ ([50,4,4]:[126,12,12]:[2,2,2]) ]]"
+CarpetRegrid::outerbounds = "[[ [[0,0],[1,1],[1,1]] ]]"
+
+
+
+ActiveThorns = "CoordBase LocalInterp Slab SymBase"
+
+ActiveThorns = "NaNChecker"
+
+NaNChecker::check_every = 32
+NaNChecker::check_vars = "all"
+NaNChecker::action_if_found = terminate
+
+
+
+ActiveThorns = "Boundary CartGrid3D Periodic"
+
+grid::domain = full
+grid::type = byrange
+grid::xmin = -0.05
+grid::ymin = -0.05
+grid::zmin = -0.05
+grid::xmax = +1.025
+grid::ymax = +0.05
+grid::zmax = +0.05
+
+Periodic::periodic = yes
+
+
+
+ActiveThorns = "Time MoL"
+
+Time::dtfac = 0.25
+
+MoL::ODE_Method = ICN
+
+#MoL::ODE_Method = RK2
+#MoL::MoL_Intermediate_Steps = 2
+
+
+
+ActiveThorns = "IDScalarWaveMoL WaveToyMoL"
+
+IDScalarWaveMoL::wave_number[0] = 1.0
+
+
+
+#ActiveThorns = "Dissipation"
+
+#Dissipation::vars = "WaveToyMoL::scalarevolve"
+
+
+
+ActiveThorns = "IOBasic CarpetIOASCII"
+
+IO::out_dir = "sw1d_2l_0040"
+
+IOBasic::outInfo_every = 32
+IOBasic::outInfo_vars = "WaveToyMoL::phi IDScalarWaveMoL::phierror"
+ 
+IOBasic::outScalar_every = 32
+IOBasic::outScalar_reductions = "minimum maximum norm1 norm2 norm_inf"
+IOBasic::outScalar_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
+
+#IOASCII::output_all_timelevels = yes
+
+IOASCII::out1D_every = 32
+IO::out_yline_x = 0.5
+IO::out_zline_x = 0.5
+IOASCII::out1D_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
+
+#IOASCII::out3D_every = 32
+#IOASCII::out3D_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
diff --git a/CarpetExtra/WaveToyMoL/par/sw1d_2l_0080.par b/CarpetExtra/WaveToyMoL/par/sw1d_2l_0080.par
new file mode 100644
index 000000000..5bfcdf1cc
--- /dev/null
+++ b/CarpetExtra/WaveToyMoL/par/sw1d_2l_0080.par
@@ -0,0 +1,99 @@
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyMoL/par/sw1d_2l_0080.par,v 1.5 2004/07/07 16:15:31 tradke Exp $
+
+Cactus::cctk_initial_time = 0
+Cactus::cctk_itlast = 12800
+
+
+
+ActiveThorns = "IOUtil Carpet CarpetLib CarpetRegrid CarpetInterp CarpetReduce CarpetSlab"
+
+driver::global_nx = 84
+driver::global_ny = 5
+driver::global_nz = 5
+driver::ghost_size = 2
+
+#Carpet::verbose = yes
+#Carpet::veryverbose = yes
+Carpet::max_refinement_levels = 3
+Carpet::prolongation_order_space = 3
+Carpet::prolongation_order_time = 2
+#Carpet::init_3_timelevels = yes
+Carpet::init_each_timelevel = yes
+
+#CarpetLib::output_bboxes = yes
+
+CarpetRegrid::refinement_levels = 2
+CarpetRegrid::refined_regions = manual-gridpoint-list
+CarpetRegrid::gridpoints = "[[ ([90,4,4]:[246,12,12]:[2,2,2]) ]]"
+CarpetRegrid::outerbounds = "[[ [[0,0],[1,1],[1,1]] ]]"
+
+
+
+ActiveThorns = "LocalInterp NaNChecker Slab"
+
+NaNChecker::check_every = 64
+NaNChecker::check_vars = "all"
+NaNChecker::action_if_found = terminate
+
+
+
+ActiveThorns = "Boundary CoordBase SymBase CartGrid3D Periodic"
+
+grid::domain = full
+grid::type = byrange
+grid::xmin = -0.025
+grid::ymin = -0.025
+grid::zmin = -0.025
+grid::xmax = +1.0125
+grid::ymax = +0.025
+grid::zmax = +0.025
+
+Periodic::periodic = yes
+
+
+
+ActiveThorns = "Time MoL"
+
+Time::dtfac = 0.25
+
+MoL::ODE_Method = ICN
+
+#MoL::ODE_Method = RK2
+#MoL::MoL_Intermediate_Steps = 2
+
+
+
+ActiveThorns = "IDScalarWaveMoL WaveToyMoL"
+
+IDScalarWaveMoL::wave_number[0] = 1.0
+
+
+
+#ActiveThorns = "Dissipation"
+#
+#Dissipation::epsdis = 0.1
+#Dissipation::timelevel = "n+1"
+#Dissipation::vars = "WaveToyMoL::scalarevolve"
+
+
+
+ActiveThorns = "IOBasic CarpetIOASCII"
+
+IO::out_dir = "sw1d_2l_0080"
+
+IOBasic::outInfo_every = 64
+IOBasic::outInfo_vars = "WaveToyMoL::phi IDScalarWaveMoL::phierror"
+ 
+IOBasic::outScalar_every = 64
+IOBasic::outScalar_reductions = "minimum maximum norm1 norm2 norm_inf"
+IOBasic::outScalar_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
+
+#IOASCII::output_all_timelevels = yes
+
+IOASCII::out1D_every = 64
+IO::out_yline_x = 0.5
+IO::out_zline_x = 0.5
+IOASCII::out1D_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
+
+#IOASCII::out3D_every = 64
+#IOASCII::out3D_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
diff --git a/CarpetExtra/WaveToyMoL/par/sw1d_2l_0160.par b/CarpetExtra/WaveToyMoL/par/sw1d_2l_0160.par
new file mode 100644
index 000000000..c47ead004
--- /dev/null
+++ b/CarpetExtra/WaveToyMoL/par/sw1d_2l_0160.par
@@ -0,0 +1,99 @@
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyMoL/par/sw1d_2l_0160.par,v 1.5 2004/07/07 16:15:31 tradke Exp $
+
+Cactus::cctk_initial_time = 0
+Cactus::cctk_itlast = 25600
+
+
+
+ActiveThorns = "IOUtil Carpet CarpetLib CarpetRegrid CarpetInterp CarpetReduce CarpetSlab"
+
+driver::global_nx = 164
+driver::global_ny = 5
+driver::global_nz = 5
+driver::ghost_size = 2
+
+#Carpet::verbose = yes
+#Carpet::veryverbose = yes
+Carpet::max_refinement_levels = 3
+Carpet::prolongation_order_space = 3
+Carpet::prolongation_order_time = 2
+#Carpet::init_3_timelevels = yes
+Carpet::init_each_timelevel = yes
+
+#CarpetLib::output_bboxes = yes
+
+CarpetRegrid::refinement_levels = 2
+CarpetRegrid::refined_regions = manual-gridpoint-list
+CarpetRegrid::gridpoints = "[[ ([170,4,4]:[486,12,12]:[2,2,2]) ]]"
+CarpetRegrid::outerbounds = "[[ [[0,0],[1,1],[1,1]] ]]"
+
+
+
+ActiveThorns = "LocalInterp NaNChecker Slab"
+
+NaNChecker::check_every = 128
+NaNChecker::check_vars = "all"
+NaNChecker::action_if_found = terminate
+
+
+
+ActiveThorns = "Boundary CoordBase SymBase CartGrid3D Periodic"
+
+grid::domain = full
+grid::type = byrange
+grid::xmin = -0.0125
+grid::ymin = -0.0125
+grid::zmin = -0.0125
+grid::xmax = +1.00625
+grid::ymax = +0.0125
+grid::zmax = +0.0125
+
+Periodic::periodic = yes
+
+
+
+ActiveThorns = "Time MoL"
+
+Time::dtfac = 0.25
+
+MoL::ODE_Method = ICN
+
+#MoL::ODE_Method = RK2
+#MoL::MoL_Intermediate_Steps = 2
+
+
+
+ActiveThorns = "IDScalarWaveMoL WaveToyMoL"
+
+IDScalarWaveMoL::wave_number[0] = 1.0
+
+
+
+#ActiveThorns = "Dissipation"
+#
+#Dissipation::epsdis = 0.1
+#Dissipation::timelevel = "n+1"
+#Dissipation::vars = "WaveToyMoL::scalarevolve"
+
+
+
+ActiveThorns = "IOBasic CarpetIOASCII"
+
+IO::out_dir = "sw1d_2l_0160"
+
+IOBasic::outInfo_every = 128
+IOBasic::outInfo_vars = "WaveToyMoL::phi IDScalarWaveMoL::phierror"
+ 
+IOBasic::outScalar_every = 128
+IOBasic::outScalar_reductions = "minimum maximum norm1 norm2 norm_inf"
+IOBasic::outScalar_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
+
+#IOASCII::output_all_timelevels = yes
+
+IOASCII::out1D_every = 128
+IO::out_yline_x = 0.5
+IO::out_zline_x = 0.5
+IOASCII::out1D_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
+
+#IOASCII::out3D_every = 128
+#IOASCII::out3D_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
diff --git a/CarpetExtra/WaveToyMoL/par/sw1d_2l_0320.par b/CarpetExtra/WaveToyMoL/par/sw1d_2l_0320.par
new file mode 100644
index 000000000..4836c42b9
--- /dev/null
+++ b/CarpetExtra/WaveToyMoL/par/sw1d_2l_0320.par
@@ -0,0 +1,99 @@
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyMoL/par/sw1d_2l_0320.par,v 1.5 2004/07/07 16:15:31 tradke Exp $
+
+Cactus::cctk_initial_time = 0
+Cactus::cctk_itlast = 51200
+
+
+
+ActiveThorns = "IOUtil Carpet CarpetLib CarpetRegrid CarpetInterp CarpetReduce CarpetSlab"
+
+driver::global_nx = 324
+driver::global_ny = 5
+driver::global_nz = 5
+driver::ghost_size = 2
+
+#Carpet::verbose = yes
+#Carpet::veryverbose = yes
+Carpet::max_refinement_levels = 3
+Carpet::prolongation_order_space = 3
+Carpet::prolongation_order_time = 2
+#Carpet::init_3_timelevels = yes
+Carpet::init_each_timelevel = yes
+
+#CarpetLib::output_bboxes = yes
+
+CarpetRegrid::refinement_levels = 2
+CarpetRegrid::refined_regions = manual-gridpoint-list
+CarpetRegrid::gridpoints = "[[ ([330,4,4]:[966,12,12]:[2,2,2]) ]]"
+CarpetRegrid::outerbounds = "[[ [[0,0],[1,1],[1,1]] ]]"
+
+
+
+ActiveThorns = "LocalInterp NaNChecker Slab"
+
+NaNChecker::check_every = 256
+NaNChecker::check_vars = "all"
+NaNChecker::action_if_found = terminate
+
+
+
+ActiveThorns = "Boundary CoordBase SymBase CartGrid3D Periodic"
+
+grid::domain = full
+grid::type = byrange
+grid::xmin = -0.00625
+grid::ymin = -0.00625
+grid::zmin = -0.00625
+grid::xmax = +1.003125
+grid::ymax = +0.00625
+grid::zmax = +0.00625
+
+Periodic::periodic = yes
+
+
+
+ActiveThorns = "Time MoL"
+
+Time::dtfac = 0.25
+
+MoL::ODE_Method = ICN
+
+#MoL::ODE_Method = RK2
+#MoL::MoL_Intermediate_Steps = 2
+
+
+
+ActiveThorns = "IDScalarWaveMoL WaveToyMoL"
+
+IDScalarWaveMoL::wave_number[0] = 1.0
+
+
+
+#ActiveThorns = "Dissipation"
+#
+#Dissipation::epsdis = 0.1
+#Dissipation::timelevel = "n+1"
+#Dissipation::vars = "WaveToyMoL::scalarevolve"
+
+
+
+ActiveThorns = "IOBasic CarpetIOASCII"
+
+IO::out_dir = "sw1d_2l_0320"
+
+IOBasic::outInfo_every = 256
+IOBasic::outInfo_vars = "WaveToyMoL::phi IDScalarWaveMoL::phierror"
+ 
+IOBasic::outScalar_every = 256
+IOBasic::outScalar_reductions = "minimum maximum norm1 norm2 norm_inf"
+IOBasic::outScalar_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
+
+#IOASCII::output_all_timelevels = yes
+
+IOASCII::out1D_every = 256
+IO::out_yline_x = 0.5
+IO::out_zline_x = 0.5
+IOASCII::out1D_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
+
+#IOASCII::out3D_every = 256
+#IOASCII::out3D_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
diff --git a/CarpetExtra/WaveToyMoL/param.ccl b/CarpetExtra/WaveToyMoL/param.ccl
new file mode 100644
index 000000000..37879c4d8
--- /dev/null
+++ b/CarpetExtra/WaveToyMoL/param.ccl
@@ -0,0 +1,24 @@
+# Parameter definitions for thorn WaveToyMoL
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyMoL/param.ccl,v 1.2 2004/04/07 17:08:21 schnetter Exp $
+
+SHARES: MethodOfLines
+
+USES CCTK_INT MoL_Num_Evolved_Vars
+
+RESTRICTED:
+
+CCTK_INT WaveToyMoL_MaxNumEvolvedVars "The maximum number of evolved variables used by WaveMoL" ACCUMULATOR-BASE=MethodofLines::MoL_Num_Evolved_Vars
+{
+  2:2 :: ""
+} 2
+
+PRIVATE:
+
+STRING bound "Type of boundary condition to use"
+{
+  .* :: "must be a registered boundary condition"
+} "None"
+
+BOOLEAN estimate_error "Estimate the truncation error"
+{
+} "no"
diff --git a/CarpetExtra/WaveToyMoL/schedule.ccl b/CarpetExtra/WaveToyMoL/schedule.ccl
new file mode 100644
index 000000000..ffce79084
--- /dev/null
+++ b/CarpetExtra/WaveToyMoL/schedule.ccl
@@ -0,0 +1,69 @@
+# Schedule definitions for thorn WaveToyMoL
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyMoL/schedule.ccl,v 1.5 2004/04/07 17:08:21 schnetter Exp $
+
+STORAGE: scalarevolve[3]
+STORAGE: scalarevolvedot
+
+if (estimate_error) {
+  STORAGE: scalarevolveerrorestimate
+}
+
+
+
+SCHEDULE WaveToyMol_Startup AT startup
+{
+  LANG: Fortran
+  OPTIONS: meta
+} "Register banner"
+
+SCHEDULE WaveToyMoL_InitSymBound AT basegrid
+{
+  LANG: Fortran
+  OPTIONS: meta
+} "Schedule symmetries"
+
+SCHEDULE WaveToyMoL_RegisterVars IN MoL_Register
+{
+  LANG: Fortran
+  OPTIONS: meta
+} "Register variables for MoL"
+
+
+
+SCHEDULE WaveToyMoL_CalcRHS IN MoL_CalcRHS
+{
+  LANG: Fortran
+} "Calculate RHS for MoL"
+
+SCHEDULE WaveToyMoL_Boundaries IN MoL_PostStep
+{
+  LANG: Fortran
+  OPTIONS: level
+  SYNC: scalarevolve
+} "Select boundary conditions in MoL"
+
+SCHEDULE GROUP ApplyBCs IN MoL_PostStep AFTER WaveToyMoL_Boundaries
+{
+} "Apply boundary conditions in MoL"
+
+
+
+if (estimate_error) {
+  
+  SCHEDULE WaveToyMoL_EstimateError AT poststep
+  {
+    LANG: Fortran
+  } "Estimate the truncation error"
+  
+  SCHEDULE WaveToyMoL_EstimateErrorBoundaries AT poststep AFTER WaveToyMoL_EstimateError
+  {
+    LANG: Fortran
+    OPTIONS: level
+    SYNC: scalarevolveerrorestimate
+  } "Select boundary conditions for the truncation error"
+  
+  SCHEDULE GROUP ApplyBCs AT poststep AFTER WaveToyMoL_EstimateError
+  {
+  } "Apply boundary conditions"
+  
+}
diff --git a/CarpetExtra/WaveToyMoL/src/boundaries.F77 b/CarpetExtra/WaveToyMoL/src/boundaries.F77
new file mode 100644
index 000000000..6b86f6383
--- /dev/null
+++ b/CarpetExtra/WaveToyMoL/src/boundaries.F77
@@ -0,0 +1,23 @@
+c     $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyMoL/src/boundaries.F77,v 1.3 2004/04/07 17:08:21 schnetter Exp $
+
+#include "cctk.h"
+#include "cctk_Arguments.h"
+#include "cctk_Functions.h"
+#include "cctk_Parameters.h"
+
+      subroutine WaveToyMoL_Boundaries (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,
+     $     "WaveToyMoL::scalarevolve", fbound);
+      if (ierr.ne.0) call CCTK_WARN (0, "internal error")
+      end
diff --git a/CarpetExtra/WaveToyMoL/src/calcrhs.F77 b/CarpetExtra/WaveToyMoL/src/calcrhs.F77
new file mode 100644
index 000000000..b3df37eb2
--- /dev/null
+++ b/CarpetExtra/WaveToyMoL/src/calcrhs.F77
@@ -0,0 +1,37 @@
+c     $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyMoL/src/calcrhs.F77,v 1.3 2003/10/27 15:31:41 schnetter Exp $
+
+#include "cctk.h"
+#include "cctk_Arguments.h"
+#include "cctk_Functions.h"
+#include "cctk_Parameters.h"
+
+      subroutine WaveToyMoL_CalcRHS (CCTK_ARGUMENTS)
+      implicit none
+      DECLARE_CCTK_ARGUMENTS
+      DECLARE_CCTK_FUNCTIONS
+      DECLARE_CCTK_PARAMETERS
+      CCTK_INT bndwidth
+      parameter (bndwidth = 1)
+      integer imin(3), imax(3)
+      CCTK_REAL dx(3)
+      integer i, j, k
+      integer d
+      do d=1,3
+         imin(d) = 1+bndwidth
+         imax(d) = cctk_lsh(d)-bndwidth
+      end do
+      do d=1,3
+         dx(d) = CCTK_DELTA_SPACE(3)
+      end do
+      do k=imin(3),imax(3)
+         do j=imin(2),imax(2)
+            do i=imin(1),imax(1)
+               phidot(i,j,k) = psi(i,j,k)
+               psidot(i,j,k) =
+     $              ( (phi(i-1,j,k) - 2*phi(i,j,k) + phi(i+1,j,k)) / dx(1)**2
+     $              + (phi(i,j-1,k) - 2*phi(i,j,k) + phi(i,j+1,k)) / dx(2)**2
+     $              + (phi(i,j,k-1) - 2*phi(i,j,k) + phi(i,j,k+1)) / dx(3)**2)
+            end do
+         end do
+      end do
+      end
diff --git a/CarpetExtra/WaveToyMoL/src/estimateerror.F77 b/CarpetExtra/WaveToyMoL/src/estimateerror.F77
new file mode 100644
index 000000000..33e5b602f
--- /dev/null
+++ b/CarpetExtra/WaveToyMoL/src/estimateerror.F77
@@ -0,0 +1,49 @@
+c     $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyMoL/src/estimateerror.F77,v 1.1 2004/04/07 17:08:21 schnetter Exp $
+
+#include "cctk.h"
+#include "cctk_Arguments.h"
+#include "cctk_Functions.h"
+#include "cctk_Parameters.h"
+
+      subroutine WaveToyMoL_EstimateError (CCTK_ARGUMENTS)
+      implicit none
+      DECLARE_CCTK_ARGUMENTS
+      DECLARE_CCTK_FUNCTIONS
+      DECLARE_CCTK_PARAMETERS
+      CCTK_INT bndwidth
+      parameter (bndwidth = 1)
+      integer imin(3), imax(3)
+      CCTK_REAL dx(3)
+      integer i, j, k
+      integer d
+      do d=1,3
+         imin(d) = 1+bndwidth
+         imax(d) = cctk_lsh(d)-bndwidth
+      end do
+      do d=1,3
+         dx(d) = CCTK_DELTA_SPACE(3)
+      end do
+      do k=imin(3),imax(3)
+         do j=imin(2),imax(2)
+            do i=imin(1),imax(1)
+               error(i,j,k) = max (
+     $              abs(phi(i+1,j,k) - phi(i-1,j,k)),
+     $              abs(phi(i,j+1,k) - phi(i,j-1,k)),
+     $              abs(phi(i,j,k+1) - phi(i,j,k-1)))
+            end do
+         end do
+      end do
+      end
+
+      subroutine WaveToyMoL_EstimateErrorBoundaries (CCTK_ARGUMENTS)
+      implicit none
+      DECLARE_CCTK_ARGUMENTS
+      DECLARE_CCTK_FUNCTIONS
+      DECLARE_CCTK_PARAMETERS
+      CCTK_INT bndwidth
+      parameter (bndwidth = 1)
+      integer ierr
+      ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, bndwidth, -1,
+     $     "WaveToyMoL::scalarevolveerrorestimate", "scalar");
+      if (ierr.ne.0) call CCTK_WARN (0, "internal error")
+      end
diff --git a/CarpetExtra/WaveToyMoL/src/initsymbound.F77 b/CarpetExtra/WaveToyMoL/src/initsymbound.F77
new file mode 100644
index 000000000..554eb16f0
--- /dev/null
+++ b/CarpetExtra/WaveToyMoL/src/initsymbound.F77
@@ -0,0 +1,25 @@
+c     $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyMoL/src/initsymbound.F77,v 1.3 2004/04/07 17:08:21 schnetter Exp $
+
+#include "cctk.h"
+#include "cctk_Arguments.h"
+#include "cctk_Functions.h"
+#include "cctk_Parameters.h"
+
+      subroutine WaveToyMoL_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 SetCartSymVN (ierr, cctkGH, sym, "WaveToyMoL::phi")
+      if (ierr.ne.0) call CCTK_WARN (0, "internal error")
+      call SetCartSymVN (ierr, cctkGH, sym, "WaveToyMoL::phidot")
+      if (ierr.ne.0) call CCTK_WARN (0, "internal error")
+      call SetCartSymVN (ierr, cctkGH, sym, "WaveToyMoL::error")
+      if (ierr.ne.0) call CCTK_WARN (0, "internal error")
+      end
diff --git a/CarpetExtra/WaveToyMoL/src/make.code.defn b/CarpetExtra/WaveToyMoL/src/make.code.defn
new file mode 100644
index 000000000..59782e596
--- /dev/null
+++ b/CarpetExtra/WaveToyMoL/src/make.code.defn
@@ -0,0 +1,9 @@
+# Main make.code.defn file for thorn WaveToyMoL
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyMoL/src/make.code.defn,v 1.2 2004/04/07 17:08:21 schnetter Exp $
+
+# Source files in this directory
+SRCS = boundaries.F77 calcrhs.F77 estimateerror.F77 initsymbound.F77 registervars.F77 startup.F77
+
+# Subdirectories containing source files
+SUBDIRS = 
+
diff --git a/CarpetExtra/WaveToyMoL/src/registervars.F77 b/CarpetExtra/WaveToyMoL/src/registervars.F77
new file mode 100644
index 000000000..51710f60a
--- /dev/null
+++ b/CarpetExtra/WaveToyMoL/src/registervars.F77
@@ -0,0 +1,21 @@
+c     $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyMoL/src/registervars.F77,v 1.2 2003/10/27 15:31:41 schnetter Exp $
+
+#include "cctk.h"
+#include "cctk_Arguments.h"
+#include "cctk_Functions.h"
+#include "cctk_Parameters.h"
+
+      subroutine WaveToyMoL_RegisterVars (CCTK_ARGUMENTS)
+      implicit none
+      DECLARE_CCTK_ARGUMENTS
+      DECLARE_CCTK_FUNCTIONS
+      DECLARE_CCTK_PARAMETERS
+      integer group, groupdot
+      integer ierr
+      call CCTK_GroupIndex(group, "wavetoymol::scalarevolve")
+      if (group.lt.0) call CCTK_WARN (0, "internal error")
+      call CCTK_GroupIndex(groupdot, "wavetoymol::scalarevolvedot")
+      if (groupdot.lt.0) call CCTK_WARN (0, "internal error")
+      ierr = MoLRegisterEvolvedGroup (group, groupdot)
+      if (ierr.ne.0) call CCTK_WARN (0, "internal error")
+      end
diff --git a/CarpetExtra/WaveToyMoL/src/startup.F77 b/CarpetExtra/WaveToyMoL/src/startup.F77
new file mode 100644
index 000000000..703d0d922
--- /dev/null
+++ b/CarpetExtra/WaveToyMoL/src/startup.F77
@@ -0,0 +1,14 @@
+c     $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyMoL/src/startup.F77,v 1.2 2003/10/27 15:31:41 schnetter Exp $
+
+#include "cctk.h"
+#include "cctk_Functions.h"
+#include "cctk_Parameters.h"
+
+      subroutine WaveToyMoL_Startup
+      implicit none
+      DECLARE_CCTK_FUNCTIONS
+      DECLARE_CCTK_PARAMETERS
+      integer ierr
+      call CCTK_RegisterBanner (ierr, "WaveToyMoL: Evolving a scalar field with MoL")
+      if (ierr.ne.0) call CCTK_WARN (0, "internal error")
+      end