anchor.diff:       do m1 = 1, 3
           ev(m1)=elmv(m1,k)
        end do
  -     do i = 1, 2*n2
  +     do i = 1, nfsize2
           do m2 = 1, 3
              el(m2) = elmnew(m2, i)
           end do
  
positive.code: subroutine inmv(lcube,radx,h,xfn,elmnew,n1,n2,elmv,imic,irec, rlbs,ltb)
  !     Initialize microvilli
    implicit none
  
    integer, parameter :: nfsize=10242, nfsize2=20480
  ! mic - number of microvilli; mrec number of microvilli per receptor
    integer, parameter :: mic=252,mrec=50
    integer n1,n2,i,j,k
    real lcube,h,radx
    real xc1, yc1, zc1, xc2, yc2, zc2, dist, distm
    real x1,y1,z1,x2,y2,z2,x3,y3,z3
    integer j1,j2,j3,m1,m2,m3,m4
    integer imic(1:mic),irec(1:mrec,1:mic)
  !     elmv - element nearest a given microvillus
    integer elmnew(1:3,1:nfsize2),elmv(1:12,1:mic)
    integer el(1:3),ev(1:12),ltb(1:2,1:mrec,1:mic)
  
  !     cmv - coordinates of microvilli
    real xfn(1:3,1:nfsize),cmv(1:3,1:mic)
  !     edr does nothing. rlbs - 
    real edr(1:nfsize2),rlbs(1:3,1:mrec,1:mic)
  
  
negative.code:            endif
          endif
       end do
       do i=1,nfsize2
          do m3=1,3
             el(m3)=elmnew(m3,i)
          end do
  