Executable: EvolveGhSingleBlackHole

---

Parallelization:
  ElementDistribution: NumGridPointsAndGridSpacing


InitialData:
  NumericInitialData:
    FileGlob: "/central/groups/sxs/geoffrey/BBH/v14/Lev1/002_Inspiral/Segment_0000/BbhVolume*.h5"
    Subgroup: "ForContinuation"
    ObservationValue: &InitialTime 4981.840000000611 
    Interpolate: True
    SetPiPhiFromConstraints: False
    Variables:
      SpacetimeMetric: SpacetimeMetric
      Pi: Pi
      Phi: Phi

DomainCreator:
  Sphere:
    InnerRadius: &InnerRadius 1.45
    OuterRadius: "599.9999"
    Interior:
      ExciseWithBoundaryCondition:
        DemandOutgoingCharSpeeds:
    InitialRefinement:
      Shell0: [2, 2, 4]
      Shell1: [2, 2, 6]
    InitialGridPoints: 12
    UseEquiangularMap: False
    EquatorialCompression: None
    RadialPartitioning: [50.0]
    RadialDistribution: [Logarithmic, Linear]
    WhichWedges: All
    TimeDependentMaps:
      InitialTime: *InitialTime
      ShapeMap:
        LMax: &MaximumL 33
        InitialValues:
          H5Filename: "/central/groups/sxs/geoffrey/BBH/v14/Lev1/003_Ringdown/RingdownDistortedCoefs.h5"
          SubfileNames: ["DistortedAhC_Ylm",
                         "DistorteddtAhC_Ylm",
                         "Distorteddt2AhC_Ylm"]
          MatchTime: *InitialTime
          MatchTimeEpsilon: 1.e-10
          SetL1CoefsToZero: True
          CheckFrame: True
        SizeInitialValues: [0.0, -1.0, 0.0]
      RotationMap:
        InitialValues:
          - [-0.27190379957418187, -1.6255871298884641e-07, -1.7375965802234716e-07, 0.9623244378987091]
          - [-0.13368166660674433, -3.184216135788834e-07, -6.056766391604965e-08,
             -0.03777162010269037]
          - [-0.016430452313229812, -9.095918473367667e-07, -2.4976971660460986e-07,
             -0.024695398712742688]
        DecayTimescaleRotation: 1.0
      ExpansionMap:
        InitialValues: [1.0, 0.0, 0.0]
        InitialValuesOuterBoundary: [0.9950186617720754,
                                     -1.0001006999419623e-06,
                                     4.041459073536267e-14]
        DecayTimescaleExpansion: 1.0
        DecayTimescaleExpansionOuterBoundary: 1.0
      TranslationMap:
        InitialValues:
          - [0.0, 0.0, 0.0]
          - [0.0, 0.0, 0.0]
          - [0.0, 0.0, 0.0]
    OuterBoundaryCondition:
      ConstraintPreservingBjorhus:
        Type: ConstraintPreservingPhysical

Evolution:
  InitialTime: *InitialTime
  InitialTimeStep: 0.0001
  InitialSlabSize: 0.1
  StepChoosers:
    - LimitIncrease:
        Factor: 2
    - ErrorControl:
        AbsoluteTolerance: 1e-10
        RelativeTolerance: 1e-8
        MaxFactor: 2
        MinFactor: 0.25
        SafetyFactor: 0.95
  TimeStepper:
    AdamsMoultonPcMonotonic:
      Order: 4

EvolutionSystem:
  GeneralizedHarmonic:
    GaugeCondition:
      DampedHarmonic:
        SpatialDecayWidth: 17.0152695482514
        Amplitudes: [1.0, 0.0, 1.0]
        Exponents: [2, 2, 2]
    DampingFunctionGamma0:
      TimeDependentTripleGaussian:
        Constant: 0.01
        Gaussian1:
          Amplitude: 7.0
          Width: 2.5
          Center: [0.0, 0.0, 0.0]
        Gaussian2:
          Amplitude: 1.0
          Width: 100.0
          Center: [0.0, 0.0, 0.0]
        Gaussian3:
          Amplitude: 0.0
          Width: 100.0
          Center: [0.0, 0.0, 0.0]
    DampingFunctionGamma1:
      Constant:
        Value: -0.999999
    DampingFunctionGamma2:
      TimeDependentTripleGaussian:
        Constant: 0.001
        Gaussian1:
          Amplitude: 7.0
          Width: 2.5
          Center: [0.0, 0.0, 0.0]
        Gaussian2:
          Amplitude: 0.1
          Width: 100.0
          Center: [0.0, 0.0, 0.0]
        Gaussian3:
          Amplitude: 0.0
          Width: 100.0
          Center: [0.0, 0.0, 0.0]

Filtering:
  ExpFilter0:
    Alpha: 64.0
    HalfPower: 210
    Enable: true
    BlocksToFilter: All

SpatialDiscretization:
  BoundaryCorrection:
    UpwindPenalty:
  DiscontinuousGalerkin:
    Formulation: StrongInertial
    Quadrature: GaussLobatto

Observers:
  VolumeFileName: "BbhVolume"
  ReductionFileName: "BbhReductions"
  SurfaceFileName: "BbhSurfaces"

Amr:
  Criteria:
  Policies:
    EnforceTwoToOneBalanceInNormalDirection: true
    Isotropy: Anisotropic
    Limits:
      RefinementLevel: Auto
      NumGridPoints: Auto
  Verbosity: Quiet

PhaseChangeAndTriggers:
  - Trigger:
      Slabs:
       EvenlySpaced:
         Interval: 100
         Offset: 0
    PhaseChanges:
      - CheckpointAndExitAfterWallclock:
          WallclockHours: 165.0

EventsAndTriggers:
  - Trigger:
      Slabs:
        EvenlySpaced:
          Interval: 1
          Offset: 0
    Events:
      - ObserveTimeStep:
          SubfileName: TimeSteps
          PrintTimeToTerminal: True
          ObservePerCore: False
  - Trigger:
      Slabs:
        EvenlySpaced:
          Interval: 10
          Offset: 0
    Events:
      - ErrorIfDataTooBig:
          Threshold: 100
          VariablesToCheck: [SpacetimeMetric]
  - Trigger:
      TimeCompares:
        Comparison: GreaterThanOrEqualTo
        Value: 5581.8399
    Events:
      - Completion

EventsAndDenseTriggers:
  - Trigger:
      Times:
        EvenlySpaced:
          Interval: 0.1
          Offset: 0.0
    Events:
      - BondiSachsInterpolation
  - Trigger:
      Times:
        EvenlySpaced:
          Interval: 0.1
          Offset: 0.0
    Events:
      - ObserveNorms:
          SubfileName: Norms
          TensorsToObserve:
          - Name: Lapse
            NormType: L2Norm
            Components: Individual
          - Name: PointwiseL2Norm(GaugeConstraint)
            NormType: L2Norm
            Components: Sum
          - Name: PointwiseL2Norm(TwoIndexConstraint)
            NormType: L2Norm
            Components: Sum
          - Name: PointwiseL2Norm(ThreeIndexConstraint)
            NormType: L2Norm
            Components: Sum
          - Name: PointwiseL2Norm(FourIndexConstraint)
            NormType: L2Norm
            Components: Sum
          - Name: PointwiseL2Norm(FConstraint)
            NormType: L2Norm
            Components: Sum
          - Name: ConstraintEnergy
            NormType: L2Norm
            Components: Sum
  - Trigger:
      Times:
        EvenlySpaced:
          Interval: 0.1
          Offset: 0.0
    Events:
      - ApparentHorizon
      - ExcisionBoundary
  - Trigger:
      Times:
        EvenlySpaced:
          Interval: 1.0
          Offset: 0.0
    Events:
      - ObserveFields:
          SubfileName: VolumeData
          VariablesToObserve:
            - Lapse
            - ConstraintEnergy
            - Psi4Real
            - SpacetimeMetric
          InterpolateToMesh: None
          CoordinatesFloatingPointType: Double
          FloatingPointTypes: [Double]

Interpolator:
  DumpVolumeDataOnFailure: false

ApparentHorizons:
  ApparentHorizon: &Ah
    InitialGuess:
      LMax: *MaximumL
      Radius: 2.2
      Center: [0., 0., 0.]
    FastFlow:
      Flow: Fast
      Alpha: 1.0
      Beta: 0.5
      AbsTol: 1e-12
      TruncationTol: 1e-5
      DivergenceTol: 5
      DivergenceIter: 5
      MaxIts: 100
    Verbosity: Quiet
  ControlSystemSingleAh: *Ah
  ControlSystemCharSpeedAh: *Ah

InterpolationTargets:
  BondiSachsInterpolation:
    LMax: 16
    Radius: [100, 150, 200]
    Center: [0, 0, 0]
    AngularOrdering: Cce
  ExcisionBoundary: &ExBdry
    LMax: *MaximumL
    Center: [0., 0., 0.]
    Radius: *InnerRadius
    AngularOrdering: "Strahlkorper"
  ControlSystemCharSpeedExcision: *ExBdry

Cce:
  BondiSachsOutputFilePrefix: "BondiSachs"

ControlSystems:
  WriteDataToDisk: true
  MeasurementsPerUpdate: 4
  Verbosity: Silent
  Shape:
    IsActive: true
    Averager:
      AverageTimescaleFraction: 0.25
      Average0thDeriv: true
    Controller:
      UpdateFraction: 0.3
    TimescaleTuner:
      InitialTimescales: 0.04
      MinTimescale: 1.0e-3
      MaxTimescale: 10.0
      IncreaseThreshold: 2.5e-4
      DecreaseThreshold: 1.0e-3
      IncreaseFactor: 1.01
      DecreaseFactor: 0.98
    ControlError:
  Size:
    IsActive: true
    Averager:
      AverageTimescaleFraction: 0.25
      Average0thDeriv: true
    Controller:
      UpdateFraction: 0.1 
    TimescaleTuner:
      InitialTimescales: 0.008
      MinTimescale: 1.0e-5
      MaxTimescale: 20.0
      IncreaseThreshold: 2.5e-4
      IncreaseFactor: 1.01
    ControlError:
      MaxNumTimesForZeroCrossingPredictor: 4
      SmoothAvgTimescaleFraction: 0.25
      DeltaRDriftOutwardOptions: None
      InitialState: DeltaR
      SmootherTuner:
        InitialTimescales: [0.02]
        MinTimescale: 1.0e-4
        MaxTimescale: 20.0
        IncreaseThreshold: 2.5e-4
        DecreaseThreshold: 1.0e-3
        IncreaseFactor: 1.01
        DecreaseFactor: 0.9

ResourceInfo:
  AvoidGlobalProc0: false
  Singletons: Auto