Add Kranc generated WaveToy

1 files changed, 182 insertions, 0 deletions

diff --git a/ML_WaveToy/param.ccl b/ML_WaveToy/param.ccl
new file mode 100644
index 0000000..606ee34
--- /dev/null
+++ b/ML_WaveToy/param.ccl
@@ -0,0 +1,182 @@
+# File produced by user eschnett
+# Produced with Mathematica Version 6.0 for Mac OS X x86 (32-bit) (April 20, 2007)
+
+# Mathematica script written by Ian Hinder and Sascha Husa
+
+
+shares: GenericFD
+
+USES CCTK_INT stencil_width
+USES CCTK_INT stencil_width_x
+USES CCTK_INT stencil_width_y
+USES CCTK_INT stencil_width_z
+USES CCTK_INT boundary_width
+
+
+shares: MethodOfLines
+
+USES CCTK_INT MoL_Num_Evolved_Vars
+USES CCTK_INT MoL_Num_Constrained_Vars
+
+restricted:
+CCTK_REAL verbose "verbose"
+{
+  "*:*" :: "no restrictions"
+} 0
+
+restricted:
+CCTK_INT ML_WaveToy_MaxNumEvolvedVars "Number of evolved variables used by this thorn" ACCUMULATOR-BASE=MethodofLines::MoL_Num_Evolved_Vars
+{
+  2:2 :: "Number of evolved variables used by this thorn"
+} 2
+
+restricted:
+CCTK_INT ML_WaveToy_MaxNumConstrainedVars "Number of constrained variables used by this thorn" ACCUMULATOR-BASE=MethodofLines::MoL_Num_Constrained_Vars
+{
+  4:4 :: "Number of constrained variables used by this thorn"
+} 4
+
+restricted:
+CCTK_INT WT_Gaussian_calc_every "WT_Gaussian_calc_every"
+{
+  *:* :: "no restrictions"
+} 1
+
+restricted:
+CCTK_INT WT_RHS_calc_every "WT_RHS_calc_every"
+{
+  *:* :: "no restrictions"
+} 1
+
+restricted:
+CCTK_INT WT_Gaussian_calc_offset "WT_Gaussian_calc_offset"
+{
+  *:* :: "no restrictions"
+} 0
+
+restricted:
+CCTK_INT WT_RHS_calc_offset "WT_RHS_calc_offset"
+{
+  *:* :: "no restrictions"
+} 0
+
+private:
+KEYWORD rho_bound "Boundary condition to implement"
+{
+  "flat" :: "Flat boundary condition"
+  "none" :: "No boundary condition"
+  "static" :: "Boundaries held fixed"
+  "radiative" :: "Radiation boundary condition"
+  "scalar" :: "Dirichlet boundary condition"
+  "newrad" :: "Improved radiative boundary condition"
+  "skip" :: "skip boundary condition code"
+} "skip"
+
+private:
+KEYWORD u_bound "Boundary condition to implement"
+{
+  "flat" :: "Flat boundary condition"
+  "none" :: "No boundary condition"
+  "static" :: "Boundaries held fixed"
+  "radiative" :: "Radiation boundary condition"
+  "scalar" :: "Dirichlet boundary condition"
+  "newrad" :: "Improved radiative boundary condition"
+  "skip" :: "skip boundary condition code"
+} "skip"
+
+private:
+KEYWORD WT_rho_bound "Boundary condition to implement"
+{
+  "flat" :: "Flat boundary condition"
+  "none" :: "No boundary condition"
+  "static" :: "Boundaries held fixed"
+  "radiative" :: "Radiation boundary condition"
+  "scalar" :: "Dirichlet boundary condition"
+  "newrad" :: "Improved radiative boundary condition"
+  "skip" :: "skip boundary condition code"
+} "skip"
+
+private:
+KEYWORD WT_u_bound "Boundary condition to implement"
+{
+  "flat" :: "Flat boundary condition"
+  "none" :: "No boundary condition"
+  "static" :: "Boundaries held fixed"
+  "radiative" :: "Radiation boundary condition"
+  "scalar" :: "Dirichlet boundary condition"
+  "newrad" :: "Improved radiative boundary condition"
+  "skip" :: "skip boundary condition code"
+} "skip"
+
+private:
+CCTK_REAL rho_bound_speed "characteristic speed at boundary"
+{
+  "0:*" :: "outgoing characteristic speed > 0"
+} 1.
+
+private:
+CCTK_REAL u_bound_speed "characteristic speed at boundary"
+{
+  "0:*" :: "outgoing characteristic speed > 0"
+} 1.
+
+private:
+CCTK_REAL WT_rho_bound_speed "characteristic speed at boundary"
+{
+  "0:*" :: "outgoing characteristic speed > 0"
+} 1.
+
+private:
+CCTK_REAL WT_u_bound_speed "characteristic speed at boundary"
+{
+  "0:*" :: "outgoing characteristic speed > 0"
+} 1.
+
+private:
+CCTK_REAL rho_bound_limit "limit value for r -> infinity"
+{
+  "*:*" :: "value of limit value is unrestricted"
+} 0.
+
+private:
+CCTK_REAL u_bound_limit "limit value for r -> infinity"
+{
+  "*:*" :: "value of limit value is unrestricted"
+} 0.
+
+private:
+CCTK_REAL WT_rho_bound_limit "limit value for r -> infinity"
+{
+  "*:*" :: "value of limit value is unrestricted"
+} 0.
+
+private:
+CCTK_REAL WT_u_bound_limit "limit value for r -> infinity"
+{
+  "*:*" :: "value of limit value is unrestricted"
+} 0.
+
+private:
+CCTK_REAL rho_bound_scalar "Dirichlet boundary value"
+{
+  "*:*" :: "unrestricted"
+} 0.
+
+private:
+CCTK_REAL u_bound_scalar "Dirichlet boundary value"
+{
+  "*:*" :: "unrestricted"
+} 0.
+
+private:
+CCTK_REAL WT_rho_bound_scalar "Dirichlet boundary value"
+{
+  "*:*" :: "unrestricted"
+} 0.
+
+private:
+CCTK_REAL WT_u_bound_scalar "Dirichlet boundary value"
+{
+  "*:*" :: "unrestricted"
+} 0.
+