source: PSPA/parmelaPSPA/trunk/spl2d1.f @ 465

      subroutine spl2d1 (nbrx,x,nbry,y,f,fx,fy,fxy,idm,ibd,t1,t2,t3)
c--------------------------------------------------------------------------
      dimension ibd(6),f(idm,nbrx),fx(idm,nbrx),fy(idm,nbrx),
     * fxy(idm,nbrx),iloc(2)
c--------------------------------------------------------------------------
c*
      data iloc/1,1/
      do 1 i=1,nbry
    1 call spl1d1 (nbrx,x,f(i,1),fx(i,1),ibd(1),idm,t1,t2,t3)
      do 2 j=1,nbrx
    2 call spl1d1 (nbry,y,f(1,j),fy(1,j),ibd(3),1,t1,t2,t3)
      call spl1d1 (nbry,y,fx(1,1),fxy(1,1),ibd(5),1,t1,t2,t3)
      call spl1d1 (nbry,y,fx(1,nbrx),fxy(1,nbrx),ibd(5),1,t1,t2,t3)
      do 3 i=1,nbry
    3 call spl1d1 (nbrx,x,fy(i,1),fxy(i,1),iloc,idm,t1,t2,t3)
      return
      end
c+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++