# PBHVlasov (GRTeclyn) — LTB e = 0 late start on particles, mu = 0.3, t_0 = 0.5 t_C(0) = 4.4086 t_H
# t_C(0) = 2159.782021, t_H = 244.948974, M_H = 183.7117; GRTeclyn time t = physical time - t_0
amr.verbose = 1
amr.n_cell = 64 64 64
geometry.prob_extent = 6194.411298 6194.411298 6194.411298
geometry.is_periodic = 1 1 1
amr.max_level = 6
amr.regrid_int = -1
amr.max_grid_size = 16
amr.blocking_factor = 8
amr.n_error_buf = 0
geometry.center = 3097.205649 3097.205649 3097.205649
evolution.dt_multiplier = 0.25
evolution.max_steps = 10000000
evolution.sigma = 0.3
# gauge: 0 = geodesic slicing (LTB comparison before the caustic), 1 = weak profile-referenced 1+log + Gamma driver
gauge.lapse_advec_coeff = 1.0
gauge.lapse_coeff = 0.3
gauge.lapse_power = 1.0
gauge.shift_advec_coeff = 0.0
gauge.shift_Gamma_coeff = 0.75
gauge.eta = 5.443311e-03
ccz4.formulation = 1
pbh_vlasov.init = table
pbh_vlasov.table = ltb_vlasov_table.dat
pbh_vlasov.n_per_dir = 2
pbh_vlasov.lattice_buffer = 4
# apparent-horizon search: every ah_interval coarse steps from t = ah_start (run time), guess radius in coordinates
pbh_vlasov.ah_interval = 4
pbh_vlasov.ah_start = 1727.8256
pbh_vlasov.ah_num_particles = 2000
pbh_vlasov.ah_max_iter = 3000
integration.type = RungeKutta
integration.rk.type = 4
ah_finder.eta = 3.0
ah_finder.c = 1.0
ah_finder.tolerance = 1e-4
ah_finder.r = 1.15
rho_line_extraction.enabled = true
rho_line_extraction.end = 6097.623621 3097.205649 3097.205649
chi_line_extraction.enabled = true
chi_line_extraction.end = 6097.623621 3097.205649 3097.205649
K_line_extraction.enabled = true
K_line_extraction.end = 6097.623621 3097.205649 3097.205649
lapse_line_extraction.enabled = true
lapse_line_extraction.end = 6097.623621 3097.205649 3097.205649

# --- rerun from t_0 (2 Oct 08:50 UTC): particle proper time, ray expansion diagnostic, guarded push
evolution.stop_time = 3671.6294
amr.plot_int = 10
amr.check_int = 20
amr.plot_vars = chi lapse K rho_p h11 h22 A11 A22 shift1
pbh_vlasov.ah_guess = 8.0
ah_finder.cfl_factor = 2.0
pbh_vlasov.tau0 = 1079.891010
pbh_vlasov.theta_rays = 1
pbh_vlasov.theta_ray_n = 256
pbh_vlasov.theta_ray_dr = 0.5
pbh_vlasov.theta_profile_interval = 4
pbh_vlasov.max_disp_cells = 0.8
# --- rerun with the fixed geodesic rhs (-Gamma d_i alpha), 2 Oct 22:xx UTC; stop at 2.2 t_C(0)