anchor.diff:  
       ! Establish binary parameters and create initial models.
   
  -    if (init == 1) then 
  -       call binary_setup
  -    endif
  +    call binary_setup
   
     end subroutine
   
  
positive.code:   double precision :: ambient_density
  
    ! Bulk system motion
  
    double precision :: bulk_velx, bulk_vely, bulk_velz
  
  contains
  
    ! This routine calls all of the other subroutines at the beginning
    ! of a simulation to fill in the basic problem data.
  
    subroutine initialize(name, namlen, init)
  
      use bl_constants_module, only: ZERO
      use bl_error_module, only: bl_error
  
      implicit none
  
      integer :: namlen, i, init
      integer :: name(namlen)
  
      ! Build "probin" filename -- the name of the file containing the fortin namelist.
      allocate(character(len=namlen) :: probin)
  
negative.code:     call eos(eos_input_rt, eos_state, .false.)
  
      small_pres = eos_state % p
      small_ener = eos_state % e
  
    end subroutine set_small
  
  