25       do i=1,nshg
26          itype=ifeature(i)
27          if(btest(iBC(i),13)) then
29                tmp(j)=y(i,5+j)
39                BC(i,6+j)=-tmp(j)
66       integer i  local
78       do i=1,nscalar 
79          tmp(i) = 0.0d0
90       do i=1,nscalar 
91          tmp(i) = D(i) * ( tmp(i) - 0.0 ) / distance  * units
92          if(itype.eq.2) tmp(i) = tmp(i) / 2.0d+00
125       do i=1,npro
127             if(iBCB(i,2).lt.0) then  !this is a face with dtn
129                   BCB(i,k,6+j)=tmpBCB(i,k,j)
187       integer i,j,n  local
190       data (readfile(i),i=1,maxtypes) / maxtypes*.false./
241       i = int ( (x - table(1,0,0,itype) ) / dx(1,itype))
244       if(i .lt. 0) then
245          i = 0
254       if(i .ge. numdata(1,itype)) then
255          i = numdata(1,itype)-2
267          z(1,1) = table(n,i,j,itype)
268          z(3,1) = table(n,i+1,j,itype)
269          z(1,2) = table(n,i,j+1,itype)
270          z(3,2) = table(n,i+1,j+1,itype)
273      &        *(x-table(1,i,j,itype)) + z(1,1)
275      &        *(x-table(1,i,j,itype)) + z(1,2)
277      &        *(y-table(2,i,j,itype)) + z(2,1)
305       integer i,j,k  local
310       i=0
327       if (i.gt.maxdata**2+maxdata+1) then
335       if(i .gt. 0 .and. x2 .ne. x2old) then
339          i=0
342       table(1,i,j,islot) = x1*1.d-0
343       table(2,i,j,islot) = x2*1.d-0
344       table(3,i,j,islot) = y1*1.d-0
345       table(4,i,j,islot) = y2*1.d-0
347       i=i+1
361          write(*,*) 'Table is ',i,' by ',j+1
363          write(*,*) 'there are ',i*(j+1),' flux data points'
371          do i=0,numdata(1,islot)-1
373                write(*,*) i,j,(table(k,i,j,islot), k=1,4)
404       integer i  local
407       data (readfile(i),i=1,maxtypes) / maxtypes*.false./
450       i = int ( (conc_BC - table(0,1,itype) ) / dx(itype))
453          i = 0
454          conc_BC =  table(i,1,itype)
456          i = numdata(itype)
457          conc_BC =  table(i,1,itype)
461       dt = conc_BC - table(i,1,itype)
462       flux_BC = dt * (table(i+1,2,itype) - table(i,2,itype)) +
463      &     table(i,2,itype)
490       integer i,j  local
494       i=-1
509       i=i+1
510       if (i.ge.maxdata) then
517       read (linein,*) table(i,1,islot), table(i,2,islot)
518       table(i,1,islot)= table(i,1,islot)*1.0d-0
519       table(i,2,islot)= table(i,2,islot)*1.0d-0
525       numdata = i
534          do 101 j=0,i