184       real*8 table(4,0:maxdata,0:maxdata,maxtypes)  local
185       save table
210          call readtable_2(itype,table,numdata,dx,
218          if( x .gt. table(1,0,0,itype) .or. 
219      &        x .lt. table(1,numdata(1,itype)-1,0,itype) ) then
222             write(*,*)  '#1  [ ',table(1,0,0,itype), ' , ',
223      &           table(1,numdata(1,itype)-1,0,itype), ' ] ',
229          if( y .gt. table(2,0,0,itype) .or. 
230      &        y .lt. table(2,0,numdata(2,itype)-1,itype) ) then
233             write(*,*)  '#2   [ ',table(2,0,0,itype), ' , ',
234      &           table(2,0,numdata(2,itype)-1,itype), ' ] ',
241       i = int ( (x - table(1,0,0,itype) ) / dx(1,itype))
242       j = int ( (y - table(2,0,0,itype) ) / dx(2,itype))
251          y = table(2,0,0,itype)
261          y = table(1,0,j+1,itype)
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)
286       subroutine readtable_2(islot,table,numdata,dx,maxdata,maxslots)  argument
300       real*8 table(4,0:maxdata,0:maxdata,maxslots), dx(2,maxslots)  local
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
357       dx(1,islot) = table(1,2,1,islot) - table(1,1,1,islot)
358       dx(2,islot) = table(2,1,2,islot) - table(2,1,1,islot)
373                write(*,*) i,j,(table(k,i,j,islot), k=1,4)
401       real*8 table(0:maxdata,2,maxtypes)  local
402       save table
428          call readtable_1(itype,table,numdata(itype),dx(itype),
450       i = int ( (conc_BC - table(0,1,itype) ) / dx(itype))
452       if( conc_BC .lt. table(0,1,itype))then
454          conc_BC =  table(i,1,itype)
455       elseif( conc_BC .gt. table(numdata(itype),1,itype)) then
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)
472       subroutine readtable_1(islot,table,numdata,dx,maxdata,maxslots)  argument
486       real*8 table(0:maxdata,2,maxslots),dx  local
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
526       dx = table(1,1,islot)-table(0,1,islot)
535             write(*,*) j,table(j,1,islot), table(j,2,islot)