anchor.diff:  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
  
positive.code:   do i = 1, mic
       read(25,21) cmv(1,i),cmv(2,i),cmv(3,i)
    end do
    
    do i = 1, mic
       x = cmv(1,i)
       y = cmv(2,i)
  
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
  