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 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:      end do
    end do
  
    ! Record the actual microvillus root positions
    call saveallsolid(cmv, 'mv-xyz-start.txt')
  
  !       do i = 1,12
  