A new excision thorn.

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/NoExcision/trunk@2 4ec1db94-0e4f-0410-ada3-8bed251432c9

7 files changed, 313 insertions, 0 deletions

diff --git a/README b/README
new file mode 100644
index 0000000..551d895
--- /dev/null
+++ b/README
@@ -0,0 +1,17 @@
+CVS info   : $Header$
+
+Cactus Code Thorn NoExcision
+Thorn Author(s)     : Erik Schnetter <schnetter@aei.mpg.de>
+Thorn Maintainer(s) : Erik Schnetter <schnetter@aei.mpg.de>
+--------------------------------------------------------------------------
+
+Purpose of the thorn:
+
+Overwrite the regions near a singularity with "something else",
+possibly flat space, ignoring the constraints.  This is instead of
+excision.
+
+This method is modelled after an idea by Dennis Pollney, following an
+serendipitous bug by Ian Hawke.  The word "serendipity" was coined by
+Horace Walpole in the 18th century, from the Persian fairy tale "The
+Three Princes of Serendip" (www.wikipedia.org).
diff --git a/doc/documentation.tex b/doc/documentation.tex
new file mode 100644
index 0000000..cdcab77
--- /dev/null
+++ b/doc/documentation.tex
@@ -0,0 +1,144 @@
+% *======================================================================*
+%  Cactus Thorn template for ThornGuide documentation
+%  Author: Ian Kelley
+%  Date: Sun Jun 02, 2002
+%  $Header$
+%
+%  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
+%  relevant 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 separated 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 graphicx package.
+%     More specifically, with the "\includegraphics" command.  Do
+%     not specify any graphic file extensions in your .tex file. This
+%     will allow us 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$
+
+\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@aei.mpg.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$ $}
+\date{June 17 2004}
+
+\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/interface.ccl b/interface.ccl
new file mode 100644
index 0000000..0aa4916
--- /dev/null
+++ b/interface.ccl
@@ -0,0 +1,6 @@
+# Interface definition for thorn NoExcision
+# $Header$
+
+IMPLEMENTS: NoExcision
+
+INHERITS: ADMBase StaticConformal grid
diff --git a/param.ccl b/param.ccl
new file mode 100644
index 0000000..ce5c85d
--- /dev/null
+++ b/param.ccl
@@ -0,0 +1,43 @@
+# Parameter definitions for thorn NoExcision
+# $Header$
+
+BOOLEAN verbose "Produce some screen output"
+{
+} "no"
+
+BOOLEAN overwrite_geometry "Set the geometry to Minkowski"
+{
+} "yes"
+
+BOOLEAN overwrite_lapse "Set the lapse to one"
+{
+} "yes"
+
+BOOLEAN overwrite_shift "Set the shift to zero"
+{
+} "yes"
+
+INT num_regions "Number of no-excision regions"
+{
+  0:10 :: ""
+} 0
+
+REAL centre_x[10] "X-coordinate of the centre of the region"
+{
+  (*:*) :: ""
+} 0.0
+
+REAL centre_y[10] "Y-coordinate of the centre of the region"
+{
+  (*:*) :: ""
+} 0.0
+
+REAL centre_z[10] "Z-coordinate of the centre of the region"
+{
+  (*:*) :: ""
+} 0.0
+
+REAL radius[10] "Radius of the region"
+{
+  0.0:*) :: ""
+} 1.0
diff --git a/schedule.ccl b/schedule.ccl
new file mode 100644
index 0000000..ac7b7b1
--- /dev/null
+++ b/schedule.ccl
@@ -0,0 +1,7 @@
+# Schedule definitions for thorn NoExcision
+# $Header$
+
+SCHEDULE NoExcision_Overwrite AT initial AFTER (ADMBase_InitialData ADMBase_InitialGauge)
+{
+  LANG: Fortran
+} "Overwrite regions with Minkowski"
diff --git a/src/make.code.defn b/src/make.code.defn
new file mode 100644
index 0000000..7dd625c
--- /dev/null
+++ b/src/make.code.defn
@@ -0,0 +1,8 @@
+# Main make.code.defn file for thorn NoExcision
+# $Header$
+
+# Source files in this directory
+SRCS = overwrite.F90
+
+# Subdirectories containing source files
+SUBDIRS = 
diff --git a/src/overwrite.F90 b/src/overwrite.F90
new file mode 100644
index 0000000..c7949ae
--- /dev/null
+++ b/src/overwrite.F90
@@ -0,0 +1,88 @@
+! $Header$
+
+#include "cctk.h"
+#include "cctk_Arguments.h"
+#include "cctk_Functions.h"
+#include "cctk_Parameters.h"
+
+subroutine NoExcision_Overwrite (CCTK_ARGUMENTS)
+  implicit none
+  DECLARE_CCTK_ARGUMENTS
+  DECLARE_CCTK_FUNCTIONS
+  DECLARE_CCTK_PARAMETERS
+  
+  CCTK_REAL :: cx, cy, cz, rad
+  integer   :: n
+  
+  do n = 1, num_regions
+     
+     cx = centre_x(n)
+     cy = centre_y(n)
+     cz = centre_z(n)
+     rad = radius(n)
+     
+     if (overwrite_geometry /= 0) then
+        
+        if (conformal_state >= 1) then
+           where ((x - cx)**2 + (y - cy)**2 + (z - cz)**2 <= rad**2)
+              psi = 1
+           end where
+        end if
+        if (conformal_state >= 2) then
+           where ((x - cx)**2 + (y - cy)**2 + (z - cz)**2 <= rad**2)
+              psix = 0
+              psiy = 0
+              psiz = 0
+           end where
+        end if
+        if (conformal_state >= 3) then
+           where ((x - cx)**2 + (y - cy)**2 + (z - cz)**2 <= rad**2)
+              psixx = 0
+              psixy = 0
+              psixz = 0
+              psiyy = 0
+              psiyz = 0
+              psizz = 0
+           end where
+        end if
+        
+        where ((x - cx)**2 + (y - cy)**2 + (z - cz)**2 <= rad**2)
+           gxx = 1
+           gxy = 0
+           gxz = 0
+           gyy = 1
+           gyz = 0
+           gzz = 1
+           kxx = 0
+           kxy = 0
+           kxz = 0
+           kyy = 0
+           kyz = 0
+           kzz = 0
+        end where
+        
+     end if
+     
+     if (overwrite_lapse /= 0) then
+        
+        where ((x - cx)**2 + (y - cy)**2 + (z - cz)**2 <= rad**2)
+           alp = 1
+        end where
+        
+     end if
+        
+     if (overwrite_shift /= 0) then
+        
+        if (shift_state /= 0) then
+           where ((x - cx)**2 + (y - cy)**2 + (z - cz)**2 <= rad**2)
+              betax = 0
+              betay = 0
+              betaz = 0
+           end where
+        end if
+        
+     end if
+     
+  end do
+  
+end subroutine NoExcision_Overwrite