anchor.diff:          ellipn3(nsph*3+1:ncap*3) = $ellipn3$
        end if
   
  -
        cap_center(1,:)=$xc_cap$
        cap_center(2,:)=$yc_cap$
        cap_center(3,:)=$zc_cap$
        a_prestress = $a_prestress$
     end if
   
  +!todo take these lines out once ellipsoid code is validated
     print *, "fineness", fineness
     print *, "ellipa", ellipa
     print *, "ellipb", ellipb
  
positive.code:   fostar = (mass*h/td**2)
  
    if (ncap > 0) then
       if (nsph > 0) then
          rad = $rad$
          ellipa(1:nsph) = rad
  
negative.code:   ! at the +/- extremes in the x, y, z directions of the sphere before the flow
    ! is applied.
    integer, allocatable :: nvec_i(:,:)
    ! ncap is the total number of elastic capsules: all elastic capsules are
    ! modeled as ellipsoids
    integer ncap
    ! nsph is the number of spheres input from the bp files, nellip is the number
    ! of ellipsoids input from the bp files. ncap=nsph+nellip, and all of
    ! the arrays that talk about the parameters of the capsules have entries for
    ! the spheres, then for the ellipsoids.
    integer nsph, nellip
    double precision ::  LCUBE,NU,MU,MASS,LENGTH
    !     Velocities are always expressed in program units
  