&fortin 

  mass_P = 0.64
  mass_S = 0.64

  central_density_P = -1.0d0
  central_density_S = -1.0d0

  nsub = 4

  problem = 0

  roche_radius_factor = 2.0d0

  collision_separation = 4.0
  collision_impact_parameter = 0.0

  interp_temp = T

  ambient_density = 1.0d-4

  stellar_temp = 1.0d7
  ambient_temp = 1.0d7

  orbital_eccentricity = 0.0d0
  orbital_angle = 0.0d0

  axis_1 = 1
  axis_2 = 2 

  max_stellar_tagging_level = 0
  max_temperature_tagging_level = 0
  max_center_tagging_level = 20
  stellar_density_threshold = 1.0d0
  temperature_tagging_threshold = 5.0d8
  center_tagging_radius = 2.0d8
  max_tagging_radius = 0.75d0
  roche_tagging_factor = 2.0d0

  bulk_velx = 0.0d0
  bulk_vely = 0.0d0
  bulk_velz = 0.0d0

  smallu = 0.0d0

  center_fracx = 0.5d0
  center_fracy = 0.5d0
  center_fracz = 0.5d0

  initial_model_dx = 6.25d5	
  initial_model_npts = 4096
  initial_model_mass_tol = 1.d-6
  initial_model_hse_tol = 1.d-10

  gw_dist = 10.0

  max_he_wd_mass = 0.45d0
  max_hybrid_wd_mass = 0.6d0
  hybrid_wd_he_shell_mass = 0.1d0
  max_co_wd_mass = 1.05d0
  co_wd_he_shell_mass = 0.0d0

  hybrid_wd_c_frac = 0.50d0
  hybrid_wd_o_frac = 0.50d0

  co_wd_c_frac = 0.5d0
  co_wd_o_frac = 0.5d0

  onemg_wd_o_frac  = 0.60d0
  onemg_wd_ne_frac = 0.35d0
  onemg_wd_mg_frac = 0.05d0

/

&tagging 

/

&sponge

  sponge_lower_radius = 3.840d9
  sponge_upper_radius = 4.352d9
  sponge_timescale    = 0.01d0

/

&extern

  ! Note that some of the parameters in this
  ! namelist are specific to the default EOS,
  ! network, and/or integrator used in the 
  ! makefile. If you try a different set of 
  ! microphysics routines be sure to check that
  ! the parameters in here are consistent, as
  ! Fortran does not like seeing unknown variables
  ! in the namelist.

  small_x = 1.d-12

  use_eos_coulomb = T
  eos_input_is_constant = T

  burning_mode = 1

  do_constant_volume_burn = T
  call_eos_in_rhs = T

  rtol_spec = 1.d-10
  atol_spec = 1.d-10

  rtol_temp = 1.d-6
  atol_temp = 1.d-6

  rtol_enuc = 1.d-6
  atol_enuc = 1.d-6

  retry_burn = T

  renormalize_abundances = T

/
