&binary_job

   evolve_both_stars = .true.
   inlist_names(1) = 'inlist1'  ! read the star properties from 'inlist1'
   inlist_names(2) = 'inlist2'  ! idem for the secondary

   pgbinary_flag = .true.

/ ! end of binary_job namelist

&binary_controls

   ! initial conditions


   ! tidal sync setting


   ! physics switches
   do_jdot_gr = .true.  ! gravitational radiation
   do_jdot_mb = .false. ! magnetic braking, off for massive stars
   do_jdot_ls = .true.  ! spin-orbit coupling (through tides)
   do_jdot_ml = .true.  ! angular momentum loss due to mass loss
   do_jdot_missing_wind = .true.  ! account for angular momentum loss from an accretor's wind
   do_j_accretion = .true.  ! compute infall stream angular momentum

   ! timestep controls
   fr = 5d-2
   fr_limit = 5d-2
   fr_dt_limit = 5d3
   fm = -1d0
   fm_hard = -1d0
   fdm = 5d-3
   fdm_hard = 5d-2
   fj = 1d-2
   fj_hard = -1d0

   ! implicit mdot solver settings
   mdot_scheme = 'contact'  ! extension of the roche_lobe scheme to allow for contact configurations
   report_rlo_solver_progress = .true.
   max_implicit_abs_mdot = 1d-1  ! max mdot; if reached, we terminate simulation
   implicit_scheme_tolerance = 1d-2
   implicit_scheme_tiny_factor = 1d-4
   starting_mdot = 1d-5  ! lower this for lower mass stars!
   roche_min_mdot = 1d-7  ! any MT computed to be lower than this is put to 0
   max_change_factor = 5d0
   num_tries_for_increase_change_factor = 5
   change_factor_increase = 1.5d0
   min_mdot_for_implicit = 1d-6

   ! don't terminate if the accretor overflows
   accretor_overflow_terminate = 1d99

   ! output settings
   photo_interval = 100
   history_interval = 1
   terminal_interval = 1
   write_header_frequency = 10
   history_dbl_format = '(1pes32.5e3, 1x)'  ! only 5 digits

/ ! end of binary_controls namelist
