26          itype=ifeature(i)
31             call Dirichlet2Neumann(nsclr, itype, tmp)
46       subroutine Dirichlet2Neumann_faux(nscalar, itype, tmp)  argument
63       integer itype, nscalar  local
92          if(itype.eq.2) tmp(i) = tmp(i) / 2.0d+00
141       subroutine Dirichlet2Neumann(nscalar, itype, tmp)  argument
142       integer nscalar, itype  local
151          call Dirichlet2Neumann_1(nscalar,itype,tmp)
155          call Dirichlet2Neumann_2(nscalar,itype,tmp)
164       subroutine Dirichlet2Neumann_2(nscalar, itype, tmp)  argument
180       integer itype, nscalar  local
208       if( .not. readfile(itype)) then
209          readfile(itype) = .true.
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), ' ] ',
224      &           numdata(1,itype)-1
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), ' ] ',
235      &           numdata(2,itype)-1
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)
254       if(i .ge. numdata(1,itype)) then
255          i = numdata(1,itype)-2
259       if(j .ge. numdata(2,itype)) then
260          j = numdata(2,itype)-2
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)
272          z(2,1) = (z(3,1) - z(1,1)) / dx(1,itype)
273      &        *(x-table(1,i,j,itype)) + z(1,1)
274          z(2,2) = (z(3,2) - z(1,2)) / dx(1,itype)
275      &        *(x-table(1,i,j,itype)) + z(1,2)
276          tmp(n-2) = (z(2,2) - z(2,1))/dx(2,itype)
277      &        *(y-table(2,i,j,itype)) + z(2,1)
383       subroutine Dirichlet2Neumann_1(nscalar, itype, tmp)  argument
399       integer itype, nscalar  local
426       if( .not. readfile(itype)) then
427          readfile(itype) = .true.
428          call readtable_1(itype,table,numdata(itype),dx(itype),
431          if(dx(itype) .eq. 0.0d0) then
432             write(*,*) 'DX for table ',itype,' is zero (I think!)'
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
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)