anchor.diff:  
       ! Establish binary parameters and create initial models.
   
  -    if (init == 1) then 
  -       call binary_setup
  -    endif
  +    call binary_setup
   
     end subroutine
   
  
positive.code:     call eos(eos_input_rt, eos_state, .false.)
  
      small_pres = eos_state % p
  !    small_ener = eos_state % e
  
    end subroutine set_small
  
  
negative.code: 
    subroutine set_small
  
      use meth_params_module, only: small_temp, small_pres, small_dens, small_ener
  
      implicit none
  
  