anchor.diff:    fostar = (mass*h/td**2)
   
     if (ncap > 0) then
  +     ! All objects described as spheres are re-posed as ellipsoids.
        if (nsph > 0) then
           rad = $rad$
           ellipa(1:nsph) = rad
  
positive.code:   ! at the +/- extremes in the x, y, z directions of the sphere before the flow
    ! is applied.
    integer, allocatable :: nvec_i(:,:)
    integer ncap
    integer nsph, nellip
    double precision ::  LCUBE,NU,MU,MASS,LENGTH
    !     Velocities are always expressed in program units
  
negative.code:         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
  