anchor.diff:    do i = 1, mic
        read(25,21) cmv(1,i),cmv(2,i),cmv(3,i)
     end do
  +  close(25)
     
  +  ! Position roots of microvilli
     do i = 1, mic
        x = cmv(1,i)
        y = cmv(2,i)
  
positive.code: subroutine fmv(klok,mass,length,time,xfn,frc,elmv, &
       imic,irec,rlbs,rlbf,ihist,ltb)
    implicit none
  
    integer, parameter :: nfsize=10242
    integer, parameter :: mic=252, mrec=50
    double precision :: mass,length,time
  
negative.code: subroutine inmv(my_ellipa, my_ellipb, my_ellipc, my_cap_center, h, &
       xfn,elmnew, elmv,imic,irec, rlbs,ltb)
    implicit none
  !     Initialize microvilli
    interface
       subroutine saveallsolid(xfn, strfname)
         character(len=*) strfname
  