============================================================================== Castro Job Information ============================================================================== job name: Castro inputs file: inputs_2d_collision.test number of MPI processes: 4 number of threads: 4 hydro tile size: (1024,16) CPU time used since start of simulation (CPU-hours): 0.806802 ============================================================================== Plotfile Information ============================================================================== output date / time: Wed Mar 25 02:45:52 2026 output dir: /raid/testing/castro-gfortran/Castro-SBU-tests/2026-03-25/wdmerger_collision-nosubcycle I/O time (s): 0.00800733 ============================================================================== Build Information ============================================================================== build date: 2026-03-25 02:42:46.055589 build machine: Linux groot.astro.sunysb.edu 6.14.6-200.fc41.x86_64 #1 SMP PREEMPT_DYNAMIC Fri May 9 19:55:50 UTC 2025 x86_64 GNU/Linux build dir: /raid/testing/castro-gfortran/Castro/Exec/science/wdmerger/tests/wdmerger_collision AMReX dir: /raid/testing/castro-gfortran/amrex/ make flags: -j30 -- AMREX_CCACHE_ENV=CCACHE_MAXSIZE=20G USE_CCACHE=TRUE TINY_PROFILE=TRUE TEST=TRUE COMP=gnu CASTRO_HOME=/raid/testing/castro-gfortran/Castro/ DIM=2 USE_OMP=TRUE USE_MPI=TRUE USE_ACC=FALSE DEBUG=TRUE MICROPHYSICS_HOME=/raid/testing/castro-gfortran/Microphysics/ AMREX_HOME=/raid/testing/castro-gfortran/amrex/ COMP: gnu COMP version: 15.2.1 C++ compiler: mpicxx C++ flags: -Werror=return-type -g -O0 -ggdb -ftrapv -Wall -Wextra -Wlogical-op -Wfloat-conversion -Wnull-dereference -Wmisleading-indentation -Wduplicated-cond -Wduplicated-branches -Wmissing-include-dirs -Wpedantic -Wshadow -Woverloaded-virtual -Wnon-virtual-dtor -std=c++20 -fopenmp -pthread -DAMREX_DEBUG -DAMREX_TINY_PROFILING -DAMREX_TESTING -DBL_LAZY -DAMREX_LAZY -DBL_USE_MPI -DAMREX_USE_MPI -DAMREX_MPI_THREAD_MULTIPLE -DBL_USE_OMP -DAMREX_USE_OMP -DBL_NO_FORT -DAMREX_GPU_MAX_THREADS=0 -DBL_SPACEDIM=2 -DAMREX_SPACEDIM=2 -DBL_FORT_USE_UNDERSCORE -DAMREX_FORT_USE_UNDERSCORE -DBL_Linux -DAMREX_Linux -DAMREX_DIMENSION_AGNOSTIC -DBL_USE_ASSERTION -DAMREX_USE_ASSERTION -DAMREX_NO_PROBINIT -DOMPI_SKIP_MPICXX -DSPONGE -DMODEL_PARSER -DNPTS_MODEL=10000 -DNUM_MODELS=2 -DGRAVITY -DREACTIONS -DSHOCK_VAR -DROTATION -DSTRANG -DNET_LOOP_UNROLL_LEN=4 -DSCREEN_METHOD=SCREEN_METHOD_screen5 -DNEUTRINO_METHOD=NEUTRINO_METHOD_sneut5 -DAUTODIFF_DEVICE_FUNC=AMREX_GPU_HOST_DEVICE -DALLOW_JACOBIAN_CACHING -DRATES -DSCREENING -DNEUTRINOS -DNAUX_NET=0 -Itmp_build_dir/s/2d.gnu.DEBUG.TPROF.MPI.OMP.EXE -I. -I/raid/testing/castro-gfortran/amrex//Src/Base -I/raid/testing/castro-gfortran/amrex//Src/Base/Parser -I/raid/testing/castro-gfortran/amrex//Src/AmrCore -I/raid/testing/castro-gfortran/amrex//Src/Amr -I/raid/testing/castro-gfortran/amrex//Src/Boundary -I/raid/testing/castro-gfortran/amrex//Src/LinearSolvers/MLMG -I/raid/testing/castro-gfortran/Microphysics//util -I/raid/testing/castro-gfortran/Microphysics//util/autodiff -I/raid/testing/castro-gfortran/Microphysics//util/gcem/include -I/raid/testing/castro-gfortran/Microphysics//integration/VODE -I/raid/testing/castro-gfortran/Microphysics//integration/utils -I/raid/testing/castro-gfortran/Microphysics//integration -I/raid/testing/castro-gfortran/Microphysics//rates -I/raid/testing/castro-gfortran/Microphysics//screening -I/raid/testing/castro-gfortran/Microphysics//neutrinos -I../../ -I/raid/testing/castro-gfortran/Castro//Source/driver -I/raid/testing/castro-gfortran/Castro//Source/hydro -I/raid/testing/castro-gfortran/Castro//Source/problems -I/raid/testing/castro-gfortran/Castro//Source/sources -I/raid/testing/castro-gfortran/Castro//Util/model_parser -I/raid/testing/castro-gfortran/Castro//Source/gravity -I/raid/testing/castro-gfortran/Castro//Source/scf -I/raid/testing/castro-gfortran/Castro//Source/reactions -I/raid/testing/castro-gfortran/Castro//Source/rotation -I/raid/testing/castro-gfortran/Microphysics//EOS -I/raid/testing/castro-gfortran/Microphysics//EOS/helmholtz -I/raid/testing/castro-gfortran/Microphysics//networks/aprox13 -I/raid/testing/castro-gfortran/Microphysics//EOS -I/raid/testing/castro-gfortran/Microphysics//networks -I/raid/testing/castro-gfortran/Microphysics//interfaces -I/raid/testing/castro-gfortran/Microphysics//constants -I/raid/testing/castro-gfortran/Microphysics//util/approx_math -I/raid/testing/castro-gfortran/Microphysics//integration/VODE -I/raid/testing/castro-gfortran/Microphysics//integration/utils -I/raid/testing/castro-gfortran/Microphysics//integration -I/raid/testing/castro-gfortran/Microphysics//rates -I/raid/testing/castro-gfortran/Microphysics//screening -I/raid/testing/castro-gfortran/Microphysics//neutrinos -Itmp_build_dir/castro_sources/2d.gnu.DEBUG.TPROF.MPI.OMP.EXE -I/raid/testing/castro-gfortran/amrex//Tools/C_scripts Fortran comp: mpif90 Fortran flags: -g -O0 -ggdb -fcheck=bounds -fbacktrace -Wuninitialized -Wunused -ffpe-trap=invalid,zero -finit-real=snan -finit-integer=2147483647 -ftrapv -ffree-line-length-none -fno-range-check -fno-second-underscore -fimplicit-none Link flags: -L. Libraries: -I/usr/include/mpich-x86_64 -L/usr/lib64/mpich/lib -lmpicxx -Wl,-rpath -Wl,/usr/lib64/mpich/lib -Wl,--enable-new-dtags -lmpi EOS: /raid/testing/castro-gfortran/Microphysics//EOS/helmholtz NETWORK: /raid/testing/castro-gfortran/Microphysics//networks/aprox13 INTEGRATOR: VODE SCREENING: screen5 Castro git describe: 26.03-27-gb3f9c35c8 AMReX git describe: 26.03-99-g9d8951b79 Microphysics git describe: 26.03-7-g736a456b ============================================================================== Grid Information ============================================================================== level: 0 number of boxes = 4 maximum zones = 64 128 level: 1 number of boxes = 6 maximum zones = 128 256 level: 2 number of boxes = 9 maximum zones = 256 512 Boundary conditions -x: symmetry +x: outflow -y: outflow +y: outflow Domain geometry info center: 0 , 0 , 0 geometry.is_periodic: 0 0 geometry.coord_sys: 1 geometry.prob_lo: 0 -5.12e+09 geometry.prob_hi: 5.12e+09 5.12e+09 amr.n_cell: 64 128 amr.max_level: 2 amr.ref_ratio: 2 2 amr.n_error_buf: 2 2 amr.regrid_int: 2 2 amr.blocking_factor: 16 16 amr.max_grid_size: 256 256 amr.subcycling_mode: None ============================================================================== Species Information ============================================================================== index name A Z ------------------------------------------------------------------------------ 0 He4 4 2 1 C12 12 6 2 O16 16 8 3 Ne20 20 10 4 Mg24 24 12 5 Si28 28 14 6 S32 32 16 7 Ar36 36 18 8 Ca40 40 20 9 Ti44 44 22 10 Cr48 48 24 11 Fe52 52 26 12 Ni56 56 28 ============================================================================== Inputs File Parameters ============================================================================== castro.use_point_mass = 0 castro.point_mass = 0 castro.point_mass_fix_solution = 0 castro.gw_dist = 0 castro.point_mass_offset_is_true = 0 castro.point_mass_location_offset = 0 castro.rotational_period = -1e+200 castro.rotation_include_centrifugal = 1 castro.rotation_include_coriolis = 1 castro.rot_source_type = 4 castro.implicit_rotation_update = 1 [*] castro.rot_axis = 1 castro.sponge_lower_radius = -1 castro.sponge_upper_radius = -1 [*] castro.sponge_lower_density = 10 [*] castro.sponge_upper_density = 10 castro.sponge_lower_pressure = -1 castro.sponge_upper_pressure = -1 castro.sponge_lower_factor = 0 castro.sponge_upper_factor = 1 castro.sponge_target_x_velocity = 0 castro.sponge_target_y_velocity = 0 castro.sponge_target_z_velocity = 0 [*] castro.sponge_timescale = 0.01 [*] castro.state_interp_order = 0 [*] castro.lin_limit_state_interp = 1 castro.do_reflux = 1 castro.update_sources_after_reflux = 1 castro.allow_non_unit_aspect_zones = 0 castro.difmag = 0.1 [*] castro.small_dens = 1e-05 [*] castro.small_temp = 100000 [*] castro.small_pres = 4.65799e+07 [*] castro.small_ener = 6.8548e+12 castro.do_hydro = 1 castro.time_integration_method = 0 castro.limit_fourth_order = 1 castro.initialization_is_cell_average = 0 castro.use_reconstructed_gamma1 = 0 castro.add_ext_src = 0 castro.hybrid_hydro = 0 castro.ppm_type = 1 castro.ppm_do_limiting = 1 castro.mhd_limit_characteristic = 1 castro.ppm_temp_fix = 0 castro.plm_iorder = 2 castro.plm_limiter = 2 castro.hybrid_riemann = 0 castro.riemann_solver = 0 castro.riemann_shock_maxiter = 12 castro.riemann_pstar_tol = 1e-05 castro.riemann_cg_blend = 2 castro.use_flattening = 1 castro.transverse_use_eos = 0 castro.transverse_reset_density = 1 [*] castro.transverse_reset_rhoe = 1 [*] castro.dual_energy_eta1 = 0.001 castro.dual_energy_eta2 = 0.0001 castro.use_pslope = 0 castro.ppm_well_balanced = 0 castro.pslope_cutoff_density = -1e+20 castro.limit_fluxes_on_small_dens = 0 castro.speed_limit = 0 [*] castro.do_sponge = 1 castro.sponge_implicit = 1 [*] castro.ext_src_implicit = 1 [*] castro.source_term_predictor = 1 castro.first_order_hydro = 0 castro.xl_ext_bc_type = -1 castro.xr_ext_bc_type = -1 castro.yl_ext_bc_type = -1 castro.yr_ext_bc_type = -1 castro.zl_ext_bc_type = -1 castro.zr_ext_bc_type = -1 castro.hse_zero_vels = 0 castro.hse_interp_temp = 0 castro.hse_fixed_temp = -1e+200 castro.hse_reflect_vels = 0 castro.fill_ambient_bc = 0 castro.ambient_fill_dir = -1 castro.ambient_outflow_vel = 0 castro.clamp_ambient_temp = 0 castro.ambient_safety_factor = 1.1 [*] castro.ambient_density = 0.0001 [*] castro.ambient_temp = 1e+07 castro.ambient_energy = -1e+200 castro.sdc_order = 2 castro.sdc_quadrature = 0 castro.sdc_extra = 0 castro.sdc_solver = 1 castro.use_geom_source = 1 castro.add_sdc_react_source_to_advection = 1 castro.hydro_memory_footprint_ratio = -1 castro.fixed_dt = -1 castro.initial_dt = -1 castro.dt_cutoff = 1e-12 castro.max_dt = 1e+200 [*] castro.cfl = 0.5 castro.init_shrink = 1 [*] castro.change_max = 1.05 castro.check_dt_before_advance = 1 castro.check_dt_after_advance = 1 castro.plot_per_is_exact = 0 castro.small_plot_per_is_exact = 0 castro.use_retry = 1 castro.retry_subcycle_factor = 0.5 castro.retry_small_density_cutoff = -1e+200 castro.abundance_failure_tolerance = 0.01 castro.abundance_failure_rho_cutoff = -1e+200 castro.use_post_step_regrid = 0 castro.max_subcycles = 10 castro.sdc_iters = 2 castro.stopping_criterion_field = castro.stopping_criterion_value = 1e+200 castro.dtnuc_e = 1e+200 castro.dtnuc_X = 1e+200 castro.dtnuc_X_threshold = 0.001 castro.do_react = 1 [*] castro.react_T_min = 1e+08 [*] castro.react_T_max = 1e+12 [*] castro.react_rho_min = 1e+06 [*] castro.react_rho_max = 1e+12 [*] castro.disable_shock_burning = 1 castro.shock_detection_threshold = 0.666667 castro.shock_detection_include_sources = 1 castro.T_guess = 1e+08 castro.drive_initial_convection = 0 castro.drive_initial_convection_tmax = 1e+200 castro.drive_initial_convection_reinit_period = 1e+200 castro.do_grav = 1 castro.moving_center = 0 castro.grav_source_type = 4 [*] castro.do_rotation = 0 castro.bndry_func_thread_safe = 1 castro.grown_factor = 1 castro.star_at_center = 1 castro.do_scf_initial_model = 0 castro.scf_maximum_density = -1e+06 castro.scf_equatorial_radius = -1e+09 castro.scf_polar_radius = -1e+09 castro.scf_relax_tol = 0.001 castro.scf_max_iterations = 30 castro.do_special_tagging = 0 [*] castro.max_tagging_radius = 0.75 [*] castro.verbose = 1 [*] castro.dump_old = 1 castro.domain_is_plane_parallel = 0 [*] castro.print_update_diagnostics = 1 [*] castro.sum_interval = 1 castro.sum_per = -1 castro.job_name = Castro castro.output_at_completion = 1 castro.reset_checkpoint_time = -1e+200 castro.reset_checkpoint_step = -1 [*] castro.store_omegadot = 1 [*] castro.store_burn_weights = 1 castro.abort_on_invalid_params = 0 castro.do_radiation = 1 [*] gravity.gravity_type = PoissonGrav gravity.const_grav = 0 gravity.direct_sum_bcs = 0 gravity.drdxfac = 1 [*] gravity.lnum = 6 [*] gravity.verbose = 1 gravity.no_sync = 0 [*] gravity.do_composite_phi_correction = 0 gravity.max_solve_level = 14 gravity.get_g_from_phi = 0 gravity.mlmg_max_fmg_iter = 0 gravity.mlmg_agglomeration = 1 gravity.mlmg_consolidation = 1 gravity.mlmg_nsolve = 0 [*] problem.mass_P = 1.25795e+33 [*] problem.mass_S = 1.25795e+33 [*] problem.central_density_P = 4.51328e+06 [*] problem.central_density_S = 4.51328e+06 problem.stellar_temp = 1e+07 problem.orbital_eccentricity = 0 problem.orbital_angle = 0 [*] problem.axis_1 = 2 [*] problem.axis_2 = 3 problem.axis_3 = 3 problem.center_fracx = 0.5 problem.center_fracy = 0.5 problem.center_fracz = 0.5 problem.initial_model_dx = 625000 [*] problem.max_he_wd_mass = 8.9478e+32 [*] problem.max_hybrid_wd_mass = 1.19304e+33 [*] problem.hybrid_wd_he_shell_mass = 1.9884e+32 [*] problem.max_co_wd_mass = 2.08782e+33 problem.co_wd_he_shell_mass = 0 problem.hybrid_wd_c_frac = 0.5 problem.hybrid_wd_o_frac = 0.5 problem.co_wd_c_frac = 0.4 problem.co_wd_o_frac = 0.6 problem.onemg_wd_o_frac = 0.6 problem.onemg_wd_ne_frac = 0.35 problem.onemg_wd_mg_frac = 0.05 problem.interp_temp = 1 problem.nsub = 1 [*] problem.problem = 0 problem.roche_radius_factor = 1 [*] problem.collision_separation = 1.09023e+09 problem.collision_impact_parameter = 0 problem.collision_velocity = -1 problem.tde_separation = 8 problem.tde_beta = 6 problem.tde_initial_velocity = 1 [*] problem.stellar_density_threshold = 100 problem.relaxation_damping_factor = -0.1 problem.relaxation_density_cutoff = 1000 problem.relaxation_cutoff_time = -1 problem.radial_damping_velocity_factor = 0 problem.radial_damping_roche_factor = 1 problem.B_x = 0 problem.B_y = 0 problem.B_z = 0 problem.m_0 = 0 integrator.X_reject_buffer = 1 integrator.call_eos_in_rhs = 1 integrator.integrate_energy = 1 integrator.jacobian = 1 integrator.burner_verbose = 0 [*] integrator.rtol_spec = 1e-06 integrator.rtol_enuc = 1e-06 [*] integrator.atol_spec = 1e-06 integrator.atol_enuc = 1e-06 [*] integrator.renormalize_abundances = 1 integrator.SMALL_X_SAFE = 1e-30 [*] integrator.MAX_TEMP = 1e+10 integrator.react_boost = -1 integrator.ode_max_steps = 150000 integrator.ode_max_dt = 1e+30 integrator.use_jacobian_caching = 1 integrator.nonaka_i = 0 integrator.nonaka_j = 0 integrator.nonaka_k = 0 integrator.nonaka_level = 0 integrator.nonaka_file = nonaka_plot.dat integrator.use_burn_retry = 0 integrator.retry_swap_jacobian = 1 integrator.retry_rtol_spec = -1 integrator.retry_rtol_enuc = -1 integrator.retry_atol_spec = -1 integrator.retry_atol_enuc = -1 integrator.do_species_clip = 0 integrator.do_corrector_validation = 1 integrator.use_number_densities = 0 integrator.scale_system = 0 integrator.nse_deriv_dt_factor = 0.05 integrator.nse_include_enu_weak = 1 integrator.linalg_do_pivoting = 1 screening.enable_chabrier1998_quantum_corr = 0 screening.enable_debye_huckel_skip = 0 screening.debye_huckel_skip_threshold = 0.01 neutrino_cooling.include_recomb = 0 integrator.subtract_internal_energy = 0 eos.use_eos_coulomb = 1 eos.eos_input_is_constant = 1 eos.eos_ttol = 1e-08 eos.eos_dtol = 1e-08 eos.prad_limiter_rho_c = -1 eos.prad_limiter_delta_rho = -1 [*] network.small_x = 1e-12 [*] network.use_tables = 1 network.use_c12ag_deboer17 = 0