anchor.diff:                   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
  +        print *, ""
  +        write (*,1001), mass_P_initial, central_density_P, radius_P_initial
  +        1001 format ("Generated initial model for primary WD of mass ", f4.2, &
  +                     " solar masses, central density ", ES8.2, " g cm**-3, and radius ", ES8.2, " cm.")
  +        print *, ""
       endif
   
       roche_rad_P = radius_P_initial
  
positive.code:            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
             endif      
  
             ! Direction of initial motion -- it is the position axis
  
negative.code:                     stellar_temp, stellar_comp, ambient_state)
  
         if (ioproc == 1 .and. init == 1) then
            write (*,1002) mass_S_initial, central_density_S, radius_S_initial
            1002 format ("Generated initial model for secondary WD of mass ", f4.2, &
                         " solar masses, central density ", ES8.2, " g cm**-3, and radius ", ES8.2, " cm.")
            print *, ""
         endif
  
         roche_rad_S = radius_S_initial
  