From 310f0ea48d18866b773136aed11200b6eda6378b Mon Sep 17 00:00:00 2001 From: eschnett <> Date: Thu, 1 Mar 2001 11:40:00 +0000 Subject: Initial revision darcs-hash:20010301114010-f6438-12fb8a9ffcc80e86c0a97e37b5b0dae0dbc59b79.gz --- CarpetExtra/WaveToyMoL/README | 10 ++ CarpetExtra/WaveToyMoL/configuration.ccl | 4 + CarpetExtra/WaveToyMoL/doc/documentation.tex | 144 +++++++++++++++++++++++++++ CarpetExtra/WaveToyMoL/interface.ccl | 42 ++++++++ CarpetExtra/WaveToyMoL/par/sw1d.gnu | 55 ++++++++++ CarpetExtra/WaveToyMoL/par/sw1d_1l_0020.par | 99 ++++++++++++++++++ CarpetExtra/WaveToyMoL/par/sw1d_1l_0040.par | 99 ++++++++++++++++++ CarpetExtra/WaveToyMoL/par/sw1d_1l_0080.par | 99 ++++++++++++++++++ CarpetExtra/WaveToyMoL/par/sw1d_2l_0020.par | 99 ++++++++++++++++++ CarpetExtra/WaveToyMoL/par/sw1d_2l_0040.par | 99 ++++++++++++++++++ CarpetExtra/WaveToyMoL/par/sw1d_2l_0080.par | 99 ++++++++++++++++++ CarpetExtra/WaveToyMoL/par/sw1d_2l_0160.par | 99 ++++++++++++++++++ CarpetExtra/WaveToyMoL/par/sw1d_2l_0320.par | 99 ++++++++++++++++++ CarpetExtra/WaveToyMoL/param.ccl | 24 +++++ CarpetExtra/WaveToyMoL/schedule.ccl | 69 +++++++++++++ CarpetExtra/WaveToyMoL/src/boundaries.F77 | 23 +++++ CarpetExtra/WaveToyMoL/src/calcrhs.F77 | 37 +++++++ CarpetExtra/WaveToyMoL/src/estimateerror.F77 | 49 +++++++++ CarpetExtra/WaveToyMoL/src/initsymbound.F77 | 25 +++++ CarpetExtra/WaveToyMoL/src/make.code.defn | 9 ++ CarpetExtra/WaveToyMoL/src/registervars.F77 | 21 ++++ CarpetExtra/WaveToyMoL/src/startup.F77 | 14 +++ 22 files changed, 1318 insertions(+) create mode 100644 CarpetExtra/WaveToyMoL/README create mode 100644 CarpetExtra/WaveToyMoL/configuration.ccl create mode 100644 CarpetExtra/WaveToyMoL/doc/documentation.tex create mode 100644 CarpetExtra/WaveToyMoL/interface.ccl create mode 100644 CarpetExtra/WaveToyMoL/par/sw1d.gnu create mode 100644 CarpetExtra/WaveToyMoL/par/sw1d_1l_0020.par create mode 100644 CarpetExtra/WaveToyMoL/par/sw1d_1l_0040.par create mode 100644 CarpetExtra/WaveToyMoL/par/sw1d_1l_0080.par create mode 100644 CarpetExtra/WaveToyMoL/par/sw1d_2l_0020.par create mode 100644 CarpetExtra/WaveToyMoL/par/sw1d_2l_0040.par create mode 100644 CarpetExtra/WaveToyMoL/par/sw1d_2l_0080.par create mode 100644 CarpetExtra/WaveToyMoL/par/sw1d_2l_0160.par create mode 100644 CarpetExtra/WaveToyMoL/par/sw1d_2l_0320.par create mode 100644 CarpetExtra/WaveToyMoL/param.ccl create mode 100644 CarpetExtra/WaveToyMoL/schedule.ccl create mode 100644 CarpetExtra/WaveToyMoL/src/boundaries.F77 create mode 100644 CarpetExtra/WaveToyMoL/src/calcrhs.F77 create mode 100644 CarpetExtra/WaveToyMoL/src/estimateerror.F77 create mode 100644 CarpetExtra/WaveToyMoL/src/initsymbound.F77 create mode 100644 CarpetExtra/WaveToyMoL/src/make.code.defn create mode 100644 CarpetExtra/WaveToyMoL/src/registervars.F77 create mode 100644 CarpetExtra/WaveToyMoL/src/startup.F77 (limited to 'CarpetExtra/WaveToyMoL') 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 +Thorn Maintainer(s) : Erik Schnetter +-------------------------------------------------------------------------- + +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 -- cgit v1.2.3