Initial commit.

2 files changed, 554 insertions, 0 deletions

diff --git a/contrib/ThornList b/contrib/ThornList
new file mode 100644
index 0000000..4501a1b
--- /dev/null
+++ b/contrib/ThornList
@@ -0,0 +1,204 @@
+# arrangement/thorn                # implements (inherits) [friend] {shares}
+#
+AEIThorns/ADMMass                  # ADMMass (ADMBase,ADMMacros,StaticConformal,SpaceMask) [ ] { }
+AEIThorns/AEILocalInterp           # AEILocalInterp ( ) [ ] { }
+AEIThorns/GSLFortran               # GSLFortran (GSL) [ ] { }
+AEIThorns/ManualTermination        # ManualTermination ( ) [ ] { }
+AEIThorns/PunctureTracker          # PunctureTracker (ADMBase,CarpetRegrid2,SphericalSurface) [ ] {SphericalSurface}
+CactusBase/Boundary                # boundary ( ) [ ] { }
+CactusBase/CartGrid3D              # grid (coordbase) [ ] {driver}
+CactusBase/CoordBase               # CoordBase ( ) [ ] { }
+CactusBase/Fortran                 # Fortran ( ) [ ] { }
+CactusBase/IOASCII                 # IOASCII ( ) [ ] {IO}
+CactusBase/IOBasic                 # IOBasic (IO) [ ] {IO}
+CactusBase/IOUtil                  # IO ( ) [ ] { }
+CactusBase/InitBase                # InitBase ( ) [ ] { }
+CactusBase/SymBase                 # SymBase ( ) [ ] { }
+CactusBase/Time                    # time ( ) [ ] {cactus}
+#CactusConnect/HTTPD                # HTTPD (Socket) [ ] {Cactus}
+#CactusConnect/HTTPDExtra           # http_utils ( ) [ ] { }
+#CactusConnect/Socket               # Socket ( ) [ ] { }
+#CactusElliptic/EllBase             # ellbase ( ) [ ] { }
+#CactusElliptic/EllPETSc            # ellpetsc (ellbase,driver) [ ] {grid,ellbase}
+#CactusElliptic/EllSOR              # ellsor (ellbase,boundary) [ ] {ellbase}
+#CactusIO/IOJpeg                    # IOJpeg (IO) [ ] {IO}
+#CactusIO/IOSDF                     # IOSDF (IO) [ ] {IO}
+CactusNumerical/Cartoon2D          # cartoon2d ( ) [ ] { }
+CactusNumerical/Dissipation        # Dissipation (grid,SphericalSurface,SpaceMask) [ ] {SpaceMask}
+CactusNumerical/InterpToArray      # InterpToArray ( ) [ ] { }
+CactusNumerical/LocalInterp        # LocalInterp ( ) [ ] { }
+CactusNumerical/LocalReduce        # LocalReduce ( ) [ ] { }
+CactusNumerical/MoL                # MethodOfLines ( ) [ ] { }
+CactusNumerical/Noise              # Noise (grid) [ ] { }
+CactusNumerical/Norms              # Norms ( ) [ ] { }
+CactusNumerical/Periodic           # Periodic ( ) [ ] { }
+CactusNumerical/ReflectionSymmetry # ReflectionSymmetry ( ) [ ] { }
+CactusNumerical/RotatingSymmetry180 # RotatingSymmetry180 ( ) [ ] { }
+CactusNumerical/RotatingSymmetry90 # RotatingSymmetry90 ( ) [ ] { }
+CactusNumerical/Slab               # Slab ( ) [ ] { }
+CactusNumerical/SlabTest           # SlabTest (Slab) [ ] { }
+CactusNumerical/SpaceMask          # SpaceMask (grid) [ ] { }
+CactusNumerical/SphericalSurface   # SphericalSurface (grid) [ ] { }
+CactusPUGH/PUGH                    # Driver ( ) [ ] {cactus}
+CactusPUGH/PUGHInterp              # Interp ( ) [ ] { }
+CactusPUGH/PUGHReduce              # Reduce ( ) [ ] { }
+CactusPUGH/PUGHSlab                # Hyperslab ( ) [ ] { }
+#CactusPUGHIO/IOFlexIO              # IOFlexIO (IO) [ ] {IO}
+#CactusPUGHIO/IOHDF5                # IOHDF5 ( ) [ ] {IO}
+#CactusPUGHIO/IOHDF5Util            # IOHDF5Util (IO) [ ] {IO}
+#CactusPUGHIO/IOPanda               # IOPanda (IO) [ ] {IO}
+#CactusPUGHIO/IOStreamedHDF5        # IOStreamedHDF5 (IO) [ ] {IO}
+#CactusPUGHIO/IsoSurfacer           # isosurfacer (Grid,IO) [ ] {IO}
+#CactusUtils/Formaline              # Formaline ( ) [ ] {IO}
+CactusUtils/NaNChecker             # NaNChecker (Reduce) [ ] { }
+CactusUtils/Nice                   # nice ( ) [ ] { }
+#CactusUtils/NoMPI                  # NoMPI ( ) [ ] { }
+CactusUtils/TerminationTrigger     # TerminationTrigger ( ) [ ] {IO}
+CactusUtils/TimerReport            # timerreport ( ) [ ] {IO}
+Carpet/Carpet                      # Driver ( ) [ ] {Cactus,InitBase,IO,CarpetLib}
+Carpet/CarpetEvolutionMask         # CarpetEvolutionMask ( ) [ ] { }
+Carpet/CarpetIOASCII               # IOASCII ( ) [ ] {IO}
+Carpet/CarpetIOBasic               # IOBasic (IO) [ ] {IO}
+Carpet/CarpetIOHDF5                # IOHDF5 ( ) [ ] {IO,Cactus}
+Carpet/CarpetIOScalar              # IOScalar ( ) [ ] {IO}
+#Carpet/CarpetIOStreamedHDF5        # IOStreamedHDF5 (IO) [ ] {IO}
+Carpet/CarpetInterp                # interp ( ) [ ] { }
+Carpet/CarpetInterp2               # interp2 ( ) [ ] { }
+Carpet/CarpetLib                   # CarpetLib ( ) [ ] {IO}
+Carpet/CarpetMask                  # CarpetMask (grid,SphericalSurface) [ ] {SphericalSurface}
+Carpet/CarpetReduce                # reduce ( ) [ ] { }
+Carpet/CarpetRegrid                # CarpetRegrid ( ) [ ] { }
+Carpet/CarpetRegrid2               # CarpetRegrid2 ( ) [ ] { }
+Carpet/CarpetSlab                  # Hyperslab ( ) [ ] { }
+#Carpet/CarpetTest                  # CarpetTest ( ) [ ] { }
+Carpet/CarpetTracker               # CarpetTracker (SphericalSurface,CarpetRegrid2) [ ] {SphericalSurface}
+Carpet/LoopControl                 # LoopControl ( ) [ ] { }
+Carpet/Timers
+Carpet/CycleClock
+#CarpetAttic/CarpetIOFlexIO         # IOFlexIO ( ) [ ] {IO}
+#CarpetAttic/CarpetIOFlexIOCheckpoint # IOFlexIO ( ) [ ] {IO,Carpet,IO}
+#CarpetAttic/Cart3dTest             # Cart3dTest (grid) [ ] { }
+#CarpetDev/CarpetAdaptiveRegrid     # CarpetAdaptiveRegrid ( ) [ ] { }
+#CarpetDev/CarpetCG                 # CarpetCG (TATelliptic) [ ] { }
+#CarpetDev/CarpetIOF5               # CarpetIOF5 (grid) [ ] {IO}
+#CarpetDev/CarpetIOF5_standalone    # IOF5 ( ) [ ] {IO}
+#CarpetDev/CarpetIONirvana          # IONirvana ( ) [ ] {IO,Cactus}
+#CarpetDev/CarpetJacobi             # CarpetJacobi (TATelliptic) [ ] { }
+#CarpetDev/CarpetMG                 # CarpetMG (TATelliptic) [ ] { }
+#CarpetDev/CarpetPETSc              # CarpetPETSc (grid) [ ] { }
+#CarpetExtra/AMRToy                 # AMRToy (grid,CarpetRegrid2) [ ] { }
+#CarpetExtra/CarpetCCTest           # CarpetCCTest ( ) [ ] { }
+#CarpetExtra/CarpetIOSer            # IOSer (Hyperslab) [ ] {IO}
+#CarpetExtra/CarpetIntegrateTest    # CarpetIntegrateTest (Reduce,Grid,Driver) [ ] { }
+#CarpetExtra/CarpetProlongateTest   # CarpetProlongateTest (grid) [ ] { }
+#CarpetExtra/CarpetRegridTest       # CarpetRegridTest (grid) [ ] { }
+#CarpetExtra/FOWaveToyF77           # FOwavetoy (grid) [ ] { }
+#CarpetExtra/HydroToy               # hydrotoy (boundary,grid) [ ] { }
+#CarpetExtra/IDFOScalarWave         # idfoscalarwave (fowavetoy,grid) [ ] { }
+#CarpetExtra/IDHydroToy             # idhydrotoy (hydrotoy,grid) [ ] { }
+#CarpetExtra/IDScalarWave           # idscalarwave (wavetoy,grid) [ ] { }
+#CarpetExtra/IDScalarWaveExpl       # IDScalarWaveExpl (Grid,WaveToyExpl) [ ] { }
+#CarpetExtra/IDScalarWaveFO         # IDScalarWaveFO (Grid,WaveToyFO) [ ] { }
+#CarpetExtra/IDScalarWaveMoL        # IDScalarWaveMoL (Grid,WaveToyMoL) [ ] { }
+#CarpetExtra/IDSpaceTimeToy         # idspacetimetoy (spacetimetoy,hydrotoy,grid) [ ] {hydrotoy}
+#CarpetExtra/Nirvana                # Nirvana ( ) [ ] { }
+#CarpetExtra/ReductionTest          # ReductionTest (GRHydro,grid) [ ] {grid}
+#CarpetExtra/ReductionTest2         # ReductionTest2 ( ) [ ] {grid}
+#CarpetExtra/ReductionTest3         # ReductionTest3 (grid) [ ] {grid}
+#CarpetExtra/RegridSyncTest         # RegridSyncTest ( ) [ ] { }
+#CarpetExtra/SpaceTimeToy           # spacetimetoy (hydrotoy,grid) [ ] {hydrotoy}
+#CarpetExtra/TestCarpetGridInfo     # TestCarpetGridInfo ( ) [ ] { }
+CarpetExtra/WaveToyExpl            # WaveToyExpl (grid) [ ] { }
+#CarpetExtra/WaveToyF77             # wavetoy (grid) [ ] { }
+CarpetExtra/WaveToyFO              # WaveToyFO (grid) [ ] {MethodOfLines}
+CarpetExtra/WaveToyMoL             # WaveToyMoL (grid) [ ] {MethodOfLines}
+EinsteinAnalysis/ADMAnalysis       # ADMAnalysis (ADMBase,StaticConformal,Grid,ADMMacros) [ ] {ADMBase}
+EinsteinAnalysis/ADMConstraints    # admconstraints (ADMBase,StaticConformal,SpaceMask,grid,ADMMacros) [ADMCoupling] {spacemask,ADMBase,StaticConformal,ADMMacros}
+EinsteinAnalysis/AHFinder          # AHFinder (ADMBase,StaticConformal,SpaceMask,Grid,IO,ADMMacros) [ ] {IO,grid,admbase,ADMMacros,spacemask}
+EinsteinAnalysis/AHFinderDirect    # AHFinderDirect (Grid,ADMBase,StaticConformal,SpaceMask,SphericalSurface,IO) [ ] {grid,ADMBase,SpaceMask,IO,SphericalSurface}
+EinsteinAnalysis/CalcK             # CalcK (ADMBase,StaticConformal) [ ] {ADMBase}
+EinsteinAnalysis/EHFinder          # ehfinder (grid,admbase,coordgauge,staticconformal,spacemask,boundary) [ ] {grid,admbase,spacemask,MethodOfLines}
+EinsteinAnalysis/Extract           # extract (grid,ADMBase,staticconformal,io) [ ] {grid,io,ADMBase}
+EinsteinAnalysis/Hydro_Analysis    # Hydro_Analysis (Grid,HydroBase) [ ] { }
+EinsteinAnalysis/Multipole         # multipole (Grid) [ ] {IO}
+EinsteinAnalysis/WeylScal4         # WeylScal4 (admbase,methodoflines,Grid,GenericFD,Boundary) [ ] {GenericFD,MethodOfLines}
+EinsteinBase/ADMBase               # ADMBase (grid) [ ] { }
+EinsteinBase/ADMCoupling           # ADMCoupling ( ) [ ] { }
+EinsteinBase/ADMMacros             # ADMMacros ( ) [ ] { }
+EinsteinBase/Constants             # Constants ( ) [ ] { }
+EinsteinBase/CoordGauge            # CoordGauge ( ) [ ] {ADMBase}
+EinsteinBase/EOSG_Base             # EOSBase ( ) [ ] { }
+EinsteinBase/EOS_Base              # EOS_Base ( ) [ ] { }
+EinsteinBase/HydroBase             # HydroBase (InitBase) [ ] {InitBase}
+EinsteinBase/StaticConformal       # StaticConformal (grid) [ ] {ADMBase}
+EinsteinBase/TmunuBase             # TmunuBase (ADMBase,StaticConformal) [ADMCoupling] { }
+EinsteinEOS/EOSG_Hybrid            # EOS_GeneralHybrid (HydroBase) [EOS_GeneralPolytrope] {EOS_GeneralPolytrope}
+EinsteinEOS/EOSG_IdealFluid        # EOS_IdealFluid ( ) [ ] { }
+EinsteinEOS/EOSG_Polytrope         # EOS_GeneralPolytrope ( ) [ ] { }
+EinsteinEOS/EOS_Hybrid             # EOS_Hybrid (HydroBase,EOS_2d_Polytrope) [ ] {EOS_2d_Polytrope}
+EinsteinEOS/EOS_IdealFluid         # EOS_Ideal_Fluid (EOS_Base) [ ] { }
+EinsteinEOS/EOS_Omni               # EOS_Omni ( ) [ ] {IO}
+EinsteinEOS/EOS_Polytrope          # EOS_2d_Polytrope (EOS_Base) [ ] { }
+#EinsteinEvolve/GRHydro             # GRHydro (ADMBase,Boundary,SpaceMask,ADMMacros,Tmunubase,HydroBase) [ ] {ADMBase,HydroBase,MethodOfLines,SpaceMask,ADMMacros}
+#EinsteinEvolve/LegoExcision        # legoexcision (spacemask,grid) [ ] { }
+EinsteinEvolve/NewRad              # NewRad ( ) [ ] { }
+#EinsteinExact/EinsteinExact_Test   # EinsteinExact_Test ( ) [ ] { }
+#EinsteinExact/GaugeWave            # GaugeWave (admbase,Grid,GenericFD,Boundary) [ ] {ADMBase,GenericFD,MethodOfLines}
+#EinsteinExact/KerrSchild           # KerrSchild (admbase,Grid,GenericFD,Boundary) [ ] {ADMBase,GenericFD,MethodOfLines}
+#EinsteinExact/Minkowski            # Minkowski (admbase,Grid,GenericFD,Boundary) [ ] {ADMBase,GenericFD,MethodOfLines}
+#EinsteinExact/ModifiedSchwarzschildBL # ModifiedSchwarzschildBL (admbase,Grid,GenericFD,Boundary) [ ] {ADMBase,GenericFD,MethodOfLines}
+#EinsteinExact/ShiftedGaugeWave     # ShiftedGaugeWave (admbase,Grid,GenericFD,Boundary) [ ] {ADMBase,GenericFD,MethodOfLines}
+#EinsteinExact/Vaidya2              # Vaidya2 (admbase,Grid,GenericFD,Boundary) [ ] {ADMBase,GenericFD,MethodOfLines}
+#EinsteinInitialData/DistortedBHIVP # distortedbhivp (ADMBase,grid,StaticConformal) [ ] {StaticConformal,admbase}
+EinsteinInitialData/Exact          # exact (admbase,grid,coordgauge,staticconformal) [ADMCoupling] {admbase,StaticConformal}
+#EinsteinInitialData/GRHydro_InitData # GRHydro_init_data (GRHydro,grid) [ ] {HydroBase,ADMBase,GRHydro}
+#EinsteinInitialData/Hydro_InitExcision # Hydro_InitExcision (HydroBase,Grid,SpaceMask) [ ] {grid}
+#EinsteinInitialData/IDAnalyticBH   # idanalyticbh (ADMBase,StaticConformal,Grid) [ ] {ADMBase,StaticConformal}
+#EinsteinInitialData/IDAxiBrillBH   # idaxibrillbh (ADMBase,grid,StaticConformal) [ ] {ADMBase}
+#EinsteinInitialData/IDAxiOddBrillBH # axiodd (ADMBase,grid,StaticConformal) [ ] {StaticConformal,admbase}
+#EinsteinInitialData/IDBrillData    # IDBrillData (grid,ADMBase,StaticConformal,ellbase) [ ] {ADMBase,StaticConformal}
+#EinsteinInitialData/IDConstraintViolate # idconstraintviolate (admbase,grid) [ ] {ADMBase}
+#EinsteinInitialData/IDFileADM      # IDFileADM (ADMBase,StaticConformal,IO) [ ] {ADMBase,StaticConformal,IO}
+#EinsteinInitialData/IDLinearWaves  # idlinearwaves (ADMBase,StaticConformal,Grid) [ ] {ADMBase,StaticConformal}
+#EinsteinInitialData/Meudon_Bin_BH  # ID_Bin_BH (grid,SummationByParts,ADMBase) [ ] {ADMBase}
+#EinsteinInitialData/Meudon_Bin_NS  # Meudon_Bin_NS (grid,SummationByParts,ADMBase,HydroBase) [ ] {ADMBase,HydroBase}
+#EinsteinInitialData/Meudon_Mag_NS  # ID_Mag_NS (grid,ADMBase,HydroBase) [ ] {ADMBase,HydroBase}
+EinsteinInitialData/NoExcision     # NoExcision (ADMBase,StaticConformal,grid,SphericalSurface,Boundary) [ ] { }
+#EinsteinInitialData/Puncture       #  ( ) [ ] { }
+#EinsteinInitialData/RotatingDBHIVP # rotatingdbhivp (ADMBase,grid,StaticConformal) [ ] {StaticConformal,ADMBase}
+#EinsteinInitialData/TOVSolver      # TOVSolver (GRHydro,Constants) [ ] {HydroBase,admbase,StaticConformal,GRHydro}
+EinsteinInitialData/Trumpet        # Trumpet (ADMBase,grid,CoordBase) [ ] {ADMBase}
+EinsteinInitialData/TwoPunctures   # TwoPunctures (ADMBase,StaticConformal,grid) [ ] {ADMBase,StaticConformal,IO}
+#EinsteinInitialData/TwoTrumpets    # TwoTrumpets (ADMBase,StaticConformal,grid) [ ] {ADMBase,StaticConformal}
+EinsteinUtils/SetMask_SphericalSurface # SetMask_SphericalSurface (HydroBase,SphericalSurface) [ ] { }
+EinsteinUtils/TGRtensor            # TGRtensor ( ) [ ] { }
+ExternalLibraries/BLAS             # BLAS ( ) [ ] { }
+ExternalLibraries/FFTW3            # FFTW3 ( ) [ ] { }
+ExternalLibraries/GSL              # GSL ( ) [ ] { }
+ExternalLibraries/HDF5             # HDF5 ( ) [ ] { }
+ExternalLibraries/LAPACK           # LAPACK ( ) [ ] { }
+#ExternalLibraries/LORENE           # LORENE ( ) [ ] { }
+ExternalLibraries/MPI              # MPI ( ) [ ] { }
+ExternalLibraries/OpenCL           # OpenCL ( ) [ ] { }
+ExternalLibraries/OpenSSL          # OpenSSL ( ) [ ] { }
+ExternalLibraries/libjpeg          # libjpeg ( ) [ ] { }
+ExternalLibraries/zlib             # zlib ( ) [ ] { }
+KrancNumericalTools/GenericFD      # GenericFD ( ) [ ] { }
+LSUThorns/QuasiLocalMeasures       # QuasiLocalMeasures (ADMBase,SphericalSurface,TmunuBase) [ ] {ADMBase,SphericalSurface}
+LSUThorns/SummationByParts         # SummationByParts ( ) [ ] { }
+LSUThorns/Vectors                  # Vectors ( ) [ ] { }
+McLachlan/ML_ADM                   # ML_ADM (ADMBase,TmunuBase,Grid,GenericFD,Boundary) [ ] {GenericFD,MethodOfLines}
+McLachlan/ML_ADMConstraints        # ML_ADMConstraints (ADMBase,TmunuBase,Grid,GenericFD,Boundary) [ ] {GenericFD,MethodOfLines}
+McLachlan/ML_ADMQuantities         # ML_ADMQuantities (ADMBase,ML_BSSN,TmunuBase,Grid,GenericFD,Boundary) [ ] {GenericFD,MethodOfLines}
+McLachlan/ML_BSSN                  # ML_BSSN (ADMBase,TmunuBase,Grid,GenericFD,Boundary) [ ] {ADMBase,GenericFD,MethodOfLines}
+McLachlan/ML_BSSN_Helper           # ML_BSSN_Helper (ADMBase,CoordGauge,ML_BSSN) [ ] {Boundary,MethodofLines,ADMBase,GenericFD,ML_BSSN}
+McLachlan/ML_Kretschmann
+McLachlan/ML_Kretschmann_Helper
+#McLachlan/ML_BSSN_Test             # ML_BSSN_Test ( ) [ ] { }
+McLachlan/ML_CCZ4                  # ML_CCZ4 (ADMBase,TmunuBase,Grid,GenericFD,Boundary) [ ] {ADMBase,GenericFD,MethodOfLines}
+McLachlan/ML_CCZ4_Helper           # ML_CCZ4_Helper (ADMBase,CoordGauge,ML_CCZ4) [ ] {Boundary,MethodofLines,ADMBase,GenericFD,ML_CCZ4}
+#McLachlan/ML_CCZ4_Test             # ML_BSSN_Test ( ) [ ] { }
+#McLachlan/ML_WaveToy               # ML_WaveToy (Grid,GenericFD,Boundary) [ ] {GenericFD,MethodOfLines}
+#McLachlan/ML_WaveToy_CL            # ML_WaveToy_CL (Grid,GenericFD,Boundary) [ ] {GenericFD,MethodOfLines}
+CactusWave/IDScalarWaveC
+CactusWave/WaveToyC
diff --git a/contrib/skel.par b/contrib/skel.par
new file mode 100644
index 0000000..59793b5
--- /dev/null
+++ b/contrib/skel.par
@@ -0,0 +1,350 @@
+##############################################################################
+#                              GRID SETUP                                    #
+##############################################################################
+ActiveThorns = "CoordBase"
+
+CoordBase::domainsize = "minmax"
+CoordBase::xmin = -256
+CoordBase::xmax = 256
+CoordBase::ymin = 0
+CoordBase::ymax = 256
+CoordBase::zmin = 0
+CoordBase::zmax = 256
+
+CoordBase::spacing = "numcells"
+CoordBase::ncells_x = 64
+CoordBase::ncells_y = 32
+CoordBase::ncells_z = 32
+
+CoordBase::boundary_size_x_lower     = 3
+CoordBase::boundary_size_y_lower     = 3
+CoordBase::boundary_size_z_lower     = 3
+CoordBase::boundary_size_x_upper     = 3
+CoordBase::boundary_size_y_upper     = 3
+CoordBase::boundary_size_z_upper     = 3
+
+CoordBase::boundary_shiftout_x_lower = 1
+CoordBase::boundary_shiftout_y_lower = 1
+CoordBase::boundary_shiftout_z_lower = 1
+
+ActiveThorns = "CartGrid3d"
+
+CartGrid3d::type                 = "coordbase"
+
+ActiveThorns = "Carpet CarpetLib CarpetRegrid2 CarpetTracker SphericalSurface IOUtil"
+
+Carpet::verbose                  = 0
+Carpet::domain_from_coordbase    = 1
+Carpet::init_fill_timelevels     = 1
+Carpet::poison_new_timelevels    = 1
+Carpet::max_refinement_levels    = 10
+Carpet::refinement_factor        = 2       # anything else not implemented
+Carpet::refinement_centering     = "vertex"
+Carpet::ghost_size               = 3
+Carpet::use_buffer_zones         = 1
+Carpet::prolongation_order_space = 5
+Carpet::prolongation_order_time  = 2
+Carpet::num_threads              = 4
+CarpetLib::poison_new_memory     = 1
+
+CarpetRegrid2::verbose                 = 1
+CarpetRegrid2::regrid_every            = 256
+CarpetRegrid2::num_centres             = 1
+CarpetRegrid2::movement_threshold_1    = 0.1
+CarpetRegrid2::position_x_1            = 0
+CarpetRegrid2::num_levels_1            = 8
+CarpetRegrid2::radius_1[1]             = 128
+CarpetRegrid2::radius_1[2]             = 64
+CarpetRegrid2::radius_1[3]             = 32
+CarpetRegrid2::radius_1[4]             = 16
+CarpetRegrid2::radius_1[5]             = 8
+CarpetRegrid2::radius_1[6]             = 4
+CarpetRegrid2::radius_1[7]             = 2
+CarpetRegrid2::radius_1[8]             = 1
+CarpetRegrid2::radius_1[9]             = 0.5
+
+CarpetTracker::verbose                 = 0
+CarpetTracker::surface[0]              = 1
+
+ActiveThorns = "ReflectionSymmetry"
+
+ReflectionSymmetry::reflection_z   = 1
+ReflectionSymmetry::avoid_origin_z = 0
+ReflectionSymmetry::reflection_y   = 1
+ReflectionSymmetry::avoid_origin_y = 0
+#ReflectionSymmetry::reflection_x   = 1
+#ReflectionSymmetry::avoid_origin_x = 0
+
+
+##############################################################################
+#                            INITIAL DATA                                    #
+##############################################################################
+
+ActiveThorns = "ADMBase ADMCoupling ADMMacros StaticConformal Trumpet"
+
+ADMBase::metric_timelevels        = 3
+ADMBase::initial_data             = "trumpet"
+ADMBase::initial_lapse            = "trumpet"
+ADMBase::initial_dtlapse          = "zero"
+ADMBase::initial_shift            = "trumpet"
+ADMBase::initial_dtshift          = "zero"
+
+ADMMacros::spatial_order = 4
+
+Trumpet::boost_velocity = 0.5
+
+##############################################################################
+#                             EVOLUTION                                      #
+##############################################################################
+
+Cactus::terminate       = "time"
+Cactus::cctk_final_time = 100
+
+ActiveThorns = "MoL Time"
+
+MoL::initial_data_is_crap   = 0
+MoL::ODE_Method             = "RK3"
+MoL::MoL_Intermediate_Steps = 3
+MoL::MoL_Num_Scratch_Levels = 1
+MoL::MoL_NaN_Check          = 1
+
+Time::dtfac = 0.125
+
+ActiveThorns = "ML_BSSN ML_BSSN_Helper CoordGauge GenericFD NewRad Boundary StaticConformal TmunuBase
+                SymBase LoopControl"
+
+ADMBase::evolution_method         = "ML_BSSN"
+ADMBase::lapse_evolution_method   = "ML_BSSN"
+ADMBase::dtlapse_evolution_method = "ML_BSSN"
+ADMBase::shift_evolution_method   = "ML_BSSN"
+ADMBase::dtshift_evolution_method = "ML_BSSN"
+
+ML_BSSN::verbose                       = 1
+ML_BSSN::timelevels                    = 3
+ML_BSSN::fdOrder                       = 4
+ML_BSSN::harmonicF                     = 2      # \dot{α} = - f α^n K
+ML_BSSN::harmonicN                     = 1      # 1+log
+ML_BSSN::conformalMethod               = 1      # χ method
+ML_BSSN::ShiftGammaCoeff               = 0.75
+ML_BSSN::HarmonicShift                 = 0
+ML_BSSN::LapseACoeff                   = 0.0
+ML_BSSN::ShiftBCoeff                   = 1.0
+ML_BSSN::AlphaDriver                   = 0.0
+ML_BSSN::BetaDriver                    = 2.0
+ML_BSSN::LapseAdvectionCoeff           = 1.0
+ML_BSSN::ShiftAdvectionCoeff           = 1.0
+ML_BSSN::MinimumLapse                  = 1e-08
+ML_BSSN::dt_lapse_shift_method         = "correct"
+ML_BSSN::EpsDiss                       = 0
+ML_BSSN::my_initial_boundary_condition = "extrapolate-gammas"
+ML_BSSN::my_rhs_boundary_condition     = "NewRad"
+Boundary::radpower                     = 2
+ML_BSSN::ML_log_confac_bound           = "none"
+ML_BSSN::ML_metric_bound               = "none"
+ML_BSSN::ML_Gamma_bound                = "none"
+ML_BSSN::ML_trace_curv_bound           = "none"
+ML_BSSN::ML_curv_bound                 = "none"
+ML_BSSN::ML_lapse_bound                = "none"
+ML_BSSN::ML_dtlapse_bound              = "none"
+ML_BSSN::ML_shift_bound                = "none"
+ML_BSSN::ML_dtshift_bound              = "none"
+
+ActiveThorns = "Dissipation SpaceMask"
+
+Dissipation::order = 5
+Dissipation::vars  = "
+        ML_BSSN::ML_log_confac
+        ML_BSSN::ML_metric
+        ML_BSSN::ML_trace_curv
+        ML_BSSN::ML_curv
+        ML_BSSN::ML_Gamma
+        ML_BSSN::ML_lapse
+        ML_BSSN::ML_shift
+        ML_BSSN::ML_dtlapse
+        ML_BSSN::ML_dtshift
+"
+ActiveThorns = "PunctureTracker"
+
+PunctureTracker::track[0] = 1
+PunctureTracker::initial_x[0] = 0
+PunctureTracker::which_surface_to_store_info[0] = 1
+
+###############################################################################
+#                               ANALYSIS                                      #
+###############################################################################
+
+ActiveThorns = "NaNChecker CarpetReduce"
+
+NaNChecker::check_every     = 1024
+NaNChecker::action_if_found = "just warn"
+NaNChecker::check_vars      = "
+        ML_BSSN::ML_log_confac
+        ML_BSSN::ML_metric
+        ML_BSSN::ML_trace_curv
+        ML_BSSN::ML_curv
+        ML_BSSN::ML_Gamma
+        ML_BSSN::ML_lapse
+        ML_BSSN::ML_shift
+        ML_BSSN::ML_dtlapse
+        ML_BSSN::ML_dtshift
+        ADMBase::metric
+        ADMBase::curv
+        ADMBase::lapse
+        ADMBase::shift
+        ADMBase::dtlapse
+        ADMBase::dtshift
+"
+
+SphericalSurface::nsurfaces = 6
+SphericalSurface::maxntheta = 39
+SphericalSurface::maxnphi   = 76
+
+SphericalSurface::ntheta      [0] = 39
+SphericalSurface::nphi        [0] = 76
+SphericalSurface::nghoststheta[0] = 2
+SphericalSurface::nghostsphi  [0] = 2
+SphericalSurface::ntheta      [1] = 39
+SphericalSurface::nphi        [1] = 76
+SphericalSurface::nghoststheta[1] = 2
+SphericalSurface::nghostsphi  [1] = 2
+SphericalSurface::ntheta      [2] = 39
+SphericalSurface::nphi        [2] = 76
+SphericalSurface::nghoststheta[2] = 2
+SphericalSurface::nghostsphi  [2] = 2
+SphericalSurface::ntheta      [3] = 39
+SphericalSurface::nphi        [3] = 76
+SphericalSurface::nghoststheta[3] = 2
+SphericalSurface::nghostsphi  [3] = 2
+SphericalSurface::ntheta      [4] = 39
+SphericalSurface::nphi        [4] = 76
+SphericalSurface::nghoststheta[4] = 2
+SphericalSurface::nghostsphi  [4] = 2
+SphericalSurface::ntheta      [5] = 39
+SphericalSurface::nphi        [5] = 76
+SphericalSurface::nghoststheta[5] = 2
+SphericalSurface::nghostsphi  [5] = 2
+
+ActiveThorns = "CarpetInterp AEILocalInterp"
+ActiveThorns = "AHFinderDirect CarpetInterp AEILocalInterp"
+
+AHFinderDirect::find_every   = 256
+AHFinderDirect::move_origins = 1
+AHFinderDirect::predict_origin_movement = 0
+AHFinderDirect::N_horizons   = 5
+
+AHFinderDirect::origin_x[1]                              = 0
+AHFinderDirect::initial_guess__coord_sphere__x_center[1] = 0
+AHFinderDirect::initial_guess__coord_sphere__radius[1]   = 1
+AHFinderDirect::which_surface_to_store_info[1]           = 0
+
+AHFinderDirect::surface_selection[2]                     = "areal radius"
+AHFinderDirect::origin_x[2]                              = 0
+AHFinderDirect::initial_guess__coord_sphere__x_center[2] = 0
+AHFinderDirect::initial_guess__coord_sphere__radius[2]   = 50
+AHFinderDirect::desired_value[2]                         = 50
+AHFinderDirect::which_surface_to_store_info[2]           = 2
+
+AHFinderDirect::surface_selection[3]                     = "areal radius"
+AHFinderDirect::origin_x[3]                              = 0
+AHFinderDirect::initial_guess__coord_sphere__x_center[3] = 0
+AHFinderDirect::initial_guess__coord_sphere__radius[3]   = 100
+AHFinderDirect::desired_value[3]                         = 100
+AHFinderDirect::which_surface_to_store_info[3]           = 3
+
+AHFinderDirect::surface_selection[4]                     = "areal radius"
+AHFinderDirect::origin_x[4]                              = 0
+AHFinderDirect::initial_guess__coord_sphere__x_center[4] = 0
+AHFinderDirect::initial_guess__coord_sphere__radius[4]   = 150
+AHFinderDirect::desired_value[4]                         = 150
+AHFinderDirect::which_surface_to_store_info[4]           = 4
+
+AHFinderDirect::surface_selection[5]                     = "areal radius"
+AHFinderDirect::origin_x[5]                              = 0
+AHFinderDirect::initial_guess__coord_sphere__x_center[5] = 0
+AHFinderDirect::initial_guess__coord_sphere__radius[5]   = 250
+AHFinderDirect::desired_value[5]                         = 250
+AHFinderDirect::which_surface_to_store_info[5]           = 5
+
+ActiveThorns = "Multipole WeylScal4"
+
+WeylScal4::verbose = 1
+Multipole::integration_method = "simpson"
+Multipole::nradii = 4
+Multipole::radius[0] = 50
+Multipole::radius[1] = 100
+Multipole::radius[2] = 150
+Multipole::radius[3] = 250
+Multipole::variables = "WeylScal4::Psi4r{sw=-2 cmplx='WeylScal4::Psi4i'}"
+Multipole::l_max = 4
+Multipole::interpolator_name = "Lagrange polynomial interpolation"
+
+ActiveThorns = "ML_Kretschmann ML_Kretschmann_Helper"
+ML_Kretschmann::conformalMethod = 1
+ML_Kretschmann::ML_Kretschmann_kretschmann_calc_every = 128
+
+ActiveThorns = "QuasiLocalMeasures"
+QuasiLocalMeasures::num_surfaces = 5
+QuasiLocalMeasures::surface_index[0] = 0
+QuasiLocalMeasures::surface_index[1] = 2
+QuasiLocalMeasures::surface_index[2] = 3
+QuasiLocalMeasures::surface_index[3] = 4
+QuasiLocalMeasures::surface_index[4] = 5
+
+###############################################################################
+#                             OUTPUT                                          #
+###############################################################################
+
+ActiveThorns = "CarpetIOBasic CarpetIOHDF5"
+
+IO::out_dir       = $parfile
+IO::out_criterion = "time"
+IO::out_dt        = 0.25
+IO::out_mode      = "onefile"
+IO::out_unchunked = 1
+IOHDF5::checkpoint = 1
+IO::checkpoint_every = 4096
+IO::checkpoint_on_terminate = 1
+IO::checkpoint_dir = $parfile
+IO::recover = "autoprobe"
+IO::recover_dir = $parfile
+
+IOBasic::outInfo_criterion  = "iteration"
+IOBasic::outInfo_every      = 128
+IOBasic::outInfo_reductions = "minimum maximum average norm1 norm2"
+IOBasic::outInfo_vars       = "
+        ML_BSSN::H
+        PunctureTracker::pt_loc_x[0]
+        PunctureTracker::pt_vel_x[0]
+        QuasiLocalMeasures::qlm_adm_energy[0]
+        QuasiLocalMeasures::qlm_adm_energy[1]
+        QuasiLocalMeasures::qlm_adm_energy[2]
+        QuasiLocalMeasures::qlm_adm_energy[3]
+        QuasiLocalMeasures::qlm_adm_energy[4]
+"
+
+IOHDF5::compression_level      = 5
+IOHDF5::output_symmetry_points = 0
+IOHDF5::out3D_ghosts        = 0
+IOHDF5::out3D_outer_ghosts  = 0
+IOHDF5::out_vars           = "
+    ML_BSSN::ML_metric
+    ML_BSSN::ML_curv
+    ML_BSSN::ML_trace_curv
+    ML_BSSN::ML_log_confac
+    ML_BSSN::ML_lapse
+    ML_BSSN::ML_shift
+    ML_BSSN::ML_dtshift
+    ML_BSSN::ML_Ham
+    ML_BSSN::ML_mom
+    ML_BSSN::ML_cons_detg
+    ML_BSSN::ML_cons_traceA
+    ML_BSSN::ML_Gamma
+    AHFinderDirect::ah_radius
+    AHFinderDirect::ah_origin
+    ML_Kretschmann::Kretsch
+    PunctureTracker::pt_loc_x[0]
+    PunctureTracker::pt_vel_x[0]
+    SphericalSurface::sf_info
+    SphericalSurface::sf_radius
+    QuasiLocalMeasures::qlm_scalars
+    "