anchor.diff:       enddo
   
   
  -     ! read in model for the secondary WD
  +     ! Read in model for the secondary WD
        call read_model_file(model_S_name)
        
  -     ! copy the data over to B's arrays
  +     ! Copy the data over to B's arrays
        allocate(model_S_r(npts_model))
        allocate(model_S_state(npts_model, nvars_model))
   
  
positive.code:         call bl_error("ERROR: Stars are touching!")
       endif
  
       ! make sure the domain is big enough
  
       ! for simplicity, make sure the x and y sizes are 2x the greatest a + R
       length = max(a_P_initial + radius_P_initial, a_S_initial + radius_S_initial)
       if (length > HALF*(probhi(1) - problo(1)) .or. &
           length > HALF*(probhi(2) - problo(2))) then
          call bl_error("ERROR: The domain width is too small to include the stars (along their axis).")
       endif
  
       ! for the height, let's take 2x the radius
       if (TWO*max(radius_P_initial, radius_S_initial) > HALF*(probhi(3) - problo(3))) then
          call bl_error("ERROR: The domain height is too small to include the stars (perpendicular to their axis).")
       endif
  
       
       ! star center positions -- we'll put them in the midplane on the
       ! x-axis, with the CM at the center of the domain
       center_P_initial(1) = center(1) - a_P_initial
       center_P_initial(2) = center(2)
  
negative.code: 
       dens_ambient = model_S_state(npts_model_S,idens_model)
  
       ! Check that the cutoff temperature match
       if (model_S_state(npts_model_S,itemp_model) /= &
           model_P_state(npts_model_P,itemp_model)) then
          call bl_error("ERROR: Cutoff temperatures for primary and secondary models do not agree.")
  