anchor.diff:       end do
     end do
   
  +  ! Record the actual microvillus root positions
     call saveallsolid(cmv, 'mv-xyz-start.txt')
   
   !       do i = 1,12
  
positive.code: subroutine inmv(my_ellipa, my_ellipb, my_ellipc, my_cap_center, h, &
       xfn,elmnew, elmv,imic,irec, rlbs,ltb)
  !     Initialize microvilli
    implicit none
    interface
       subroutine saveallsolid(xfn, strfname)
         character(len=*) strfname
  
negative.code: subroutine fmv(klok,mass,length,time,xfn,frc,elmv, &
       imic,irec,rlbs,rlbf,ihist,ltb)
    implicit none
    ! Determine force due to microvillus/bond physics
    integer, parameter :: nfsize=10242
    integer, parameter :: mic=252, mrec=50
    double precision :: mass,length,time
  