Z4co: copy ccl files from Z4c

Z4co/configuration.ccl

@@ -0,0 +1,3 @@
+# Configuration definitions for thorn Z4co
+
+REQUIRES Arith Loop NewRadX

Z4co/interface.ccl

@@ -0,0 +1,52 @@
+# Interface definition for thorn Z4co
+
+IMPLEMENTS: Z4co
+
+INHERITS: ADMBaseX TmunuBaseX
+
+USES INCLUDE HEADER: defs.hxx
+USES INCLUDE HEADER: div.hxx
+USES INCLUDE HEADER: dual.hxx
+USES INCLUDE HEADER: loop_device.hxx
+USES INCLUDE HEADER: mat.hxx
+USES INCLUDE HEADER: simd.hxx
+USES INCLUDE HEADER: sum.hxx
+USES INCLUDE HEADER: vec.hxx
+USES INCLUDE HEADER: vect.hxx
+
+USES INCLUDE HEADER: newradx.hxx
+
+
+CCTK_INT FUNCTION GetCallFunctionCount()
+REQUIRES FUNCTION GetCallFunctionCount
+
+
+
+# All variables have been shifted so that they tend to zero in flat space
+
+CCTK_REAL chi TYPE=gf TAGS='rhs="chi_rhs" dependents="ADMBaseX::metric"' "chi"
+CCTK_REAL gamma_tilde TYPE=gf TAGS='parities={+1 +1 +1   -1 -1 +1   -1 +1 -1   +1 +1 +1   +1 -1 -1   +1 +1 +1} rhs="gamma_tilde_rhs" dependents="ADMBaseX::metric"' { gammatxx gammatxy gammatxz gammatyy gammatyz gammatzz } "gamma-tilde"
+CCTK_REAL K_hat TYPE=gf TAGS='rhs="K_hat_rhs" dependents="ADMBaseX::curv"' { Kh } "K-hat"
+CCTK_REAL A_tilde TYPE=gf TAGS='parities={+1 +1 +1   -1 -1 +1   -1 +1 -1   +1 +1 +1   +1 -1 -1   +1 +1 +1} rhs="A_tilde_rhs" dependents="ADMBaseX::curv"' { Atxx Atxy Atxz Atyy Atyz Atzz } "A-tilde"
+CCTK_REAL Gam_tilde TYPE=gf TAGS='parities={-1 +1 +1   +1 -1 +1   +1 +1 -1} rhs="Gam_tilde_rhs"' { Gamtx Gamty Gamtz } "Gamma-tilde"
+CCTK_REAL Theta TYPE=gf TAGS='rhs="Theta_rhs" dependents="ADMBaseX::curv"' "Theta"
+CCTK_REAL alphaG TYPE=gf TAGS='rhs="alphaG_rhs" dependents="ADMBaseX::lapse ADMBaseX::dtlapse"' "alpha"
+CCTK_REAL betaG TYPE=gf TAGS='parities={-1 +1 +1   +1 -1 +1   +1 +1 -1} rhs="betaG_rhs" dependents="ADMBaseX::shift ADMBaseX::dtshift"' { betaGx betaGy betaGz } "beta"
+
+
+
+CCTK_REAL ZtC TYPE=gf TAGS='parities={-1 +1 +1   +1 -1 +1   +1 +1 -1} checkpoint="no"' { ZtCx ZtCy ZtCz } "Z-tilde"
+CCTK_REAL HC TYPE=gf TAGS='checkpoint="no"' "H"
+CCTK_REAL MtC TYPE=gf TAGS='parities={-1 +1 +1   +1 -1 +1   +1 +1 -1} checkpoint="no"' { MtCx MtCy MtCz } "M-tilde"
+CCTK_REAL allC TYPE=gf TAGS='checkpoint="no"' "constraint monitor"
+
+
+
+CCTK_REAL chi_rhs TYPE=gf TAGS='checkpoint="no"' "chi"
+CCTK_REAL gamma_tilde_rhs TYPE=gf TAGS='parities={+1 +1 +1   -1 -1 +1   -1 +1 -1   +1 +1 +1   +1 -1 -1   +1 +1 +1} checkpoint="no"' { gammatxx_rhs gammatxy_rhs gammatxz_rhs gammatyy_rhs gammatyz_rhs gammatzz_rhs } "gamma-tilde"
+CCTK_REAL K_hat_rhs TYPE=gf TAGS='checkpoint="no"' { Kh_rhs } "K-hat"
+CCTK_REAL A_tilde_rhs TYPE=gf TAGS='parities={+1 +1 +1   -1 -1 +1   -1 +1 -1   +1 +1 +1   +1 -1 -1   +1 +1 +1} checkpoint="no"' { Atxx_rhs Atxy_rhs Atxz_rhs Atyy_rhs Atyz_rhs Atzz_rhs } "A-tilde"
+CCTK_REAL Gam_tilde_rhs TYPE=gf TAGS='parities={-1 +1 +1   +1 -1 +1   +1 +1 -1} checkpoint="no"' { Gamtx_rhs Gamty_rhs Gamtz_rhs } "Gamma-tilde"
+CCTK_REAL Theta_rhs TYPE=gf TAGS='checkpoint="no"' "Theta"
+CCTK_REAL alphaG_rhs TYPE=gf TAGS='checkpoint="no"' "alpha"
+CCTK_REAL betaG_rhs TYPE=gf TAGS='parities={-1 +1 +1   +1 -1 +1   +1 +1 -1} checkpoint="no"' { betaGx_rhs betaGy_rhs betaGz_rhs } "beta"

Z4co/param.ccl

@@ -0,0 +1,103 @@
+# Parameter definitions for thorn Z4co
+
+BOOLEAN calc_ADM_vars "Calculate ADM variables" STEERABLE=recover
+{
+} yes
+
+BOOLEAN calc_ADMRHS_vars "Calculate RHS of ADM variables" STEERABLE=recover
+{
+} yes
+
+BOOLEAN calc_constraints "Calculate constraints" STEERABLE=recover
+{
+} yes
+
+
+BOOLEAN set_Theta_zero "set Theta to zero, which converts Z4c to BSSN"
+{
+} no
+
+CCTK_REAL kappa1 "kappa1" STEERABLE=always
+{
+  *:* :: ""
+} 0.02
+
+CCTK_REAL kappa2 "kappa2" STEERABLE=always
+{
+  *:* :: ""
+} 0.0
+
+CCTK_REAL f_mu_L "mu_L = f_mu_L / alpha" STEERABLE=always
+{
+  *:* :: ""
+} 2.0
+
+CCTK_REAL f_mu_S "mu_S = f_mu_S / alpha^2" STEERABLE=always
+{
+  *:* :: ""
+} 1.0
+
+CCTK_REAL eta "eta" STEERABLE=always
+{
+  *:* :: ""
+} 2.0
+
+CCTK_REAL chi_floor "Floor for chi" STEERABLE=always
+{
+  (0:* :: ""
+} 1.0e-10
+
+CCTK_REAL alphaG_floor "Floor for alphaG" STEERABLE=always
+{
+  (0:* :: ""
+} 1.0e-10
+
+CCTK_REAL epsdiss "Dissipation coefficient <arXiv:gr-qc/0610128>" STEERABLE=always
+{
+  0.0:* :: ""
+} 0.32
+
+KEYWORD boundary_conditions "boundary conditions"{
+  "NewRadX" :: "radiative boundary conditions using NewRadX"
+  "CarpetX" :: "use CarpetX default boundary conditions"
+} "CarpetX"
+
+CCTK_INT n_chi "n power of outgoing boundary r^n fall off rate for chi"
+{
+  0:2           :: "1 is reasonable"
+} 1
+
+CCTK_INT n_gammat "n power of outgoing boundary r^n fall off rate for gammat_ij"
+{
+  0:2           :: "1 is reasonable"
+} 1
+
+CCTK_INT n_Kh "n power of outgoing boundary r^n fall off rate for Kh"
+{
+  0:2           :: "Maybe 1?"
+} 1
+
+CCTK_INT n_At "n power of outgoing boundary r^n fall off rate for At_ij"
+{
+  0:2           :: "Maybe 1?"
+} 1
+
+CCTK_INT n_Gamt "n power of outgoing boundary r^n fall off rate for Gamt^i"
+{
+  0:2           :: "Maybe 1?"
+} 1
+
+CCTK_INT n_Theta "n power of outgoing boundary r^n fall off rate for Theta"
+{
+  0:2           :: "Maybe 1?"
+} 1
+
+CCTK_INT n_alphaG "n power of outgoing boundary r^n fall off rate for alpha"
+{
+  0:2           :: "1 is my guess"
+} 1
+
+CCTK_INT n_betaG "n power of outgoing boundary r^n fall off rate for beta"
+{
+  0:2           :: "1 is my guess"
+} 1