anchor.diff:             endif
   
              if (ioproc == 1 .and. init == 1) then
  -              print *, "Generated binary orbit of distance ", a, &
  -                       ", primary distance ", a_P_initial, ", and secondary distance", a_S_initial
  +              write (*,1003) a, a / AU
  +              write (*,1004) a_P_initial, a_P_initial / AU
  +              write (*,1005) a_S_initial, a_S_initial / AU
  +              1003 format ("Generated binary orbit of distance ", ES8.2, " cm = ", ES8.2, " AU.")
  +              1004 format ("The primary orbits the center of mass at distance ", ES8.2, " cm = ", ES8.2, " AU.")
  +              1005 format ("The secondary orbits the center of mass at distance ", ES8.2, " cm = ", ES8.2, " AU.")
  +              print *, ""
              endif      
   
              ! Direction of initial motion -- it is the position axis
  
positive.code:                  stellar_temp, stellar_comp, ambient_state)
  
      if (ioproc == 1 .and. init == 1) then
          print *, "Generated initial model for primary WD of mass", mass_P_initial, &
                   ", central density", central_density_P, ", and radius", radius_P_initial
      endif
  
      roche_rad_P = radius_P_initial
  
negative.code:   subroutine binary_setup
  
      use bl_constants_module, only: ZERO, ONE
      use fundamental_constants_module, only: M_solar, AU
      use meth_params_module, only: do_rotation, rot_period
      use initial_model_module, only: init_1d
      use prob_params_module, only: center, problo, probhi
  