1 subroutine r8pc2ev(xget,yget,x,nx,y,ny,ilinx,iliny,
17 INTEGER,
PARAMETER :: R8=selected_real_kind(12,100)
87 call r8herm2xy(xget,yget,x,nx,y,ny,ilinx,iliny,
88 > i,j,xparam,yparam,hx,hxi,hy,hyi,ier)
91 call r8pc2fcn(ict,1,1,
92 > fval,i,j,xparam,yparam,hx,hxi,hy,hyi,f,inf2,ny)
100 subroutine r8pc2fcn(ict,ivec,ivecd,
101 > fval,ii,jj,xparam,yparam,hx,hxi,hy,hyi,
107 INTEGER,
PARAMETER :: R8=selected_real_kind(12,100)
108 INTEGER ny,inf2,i,j,iadr
117 integer ii(ivec),jj(ivec)
118 real*8 xparam(ivec),yparam(ivec)
121 real*8 hx(ivec),hy(ivec)
122 real*8 hxi(ivec),hyi(ivec)
168 sum=ypi*(xpi*fin(i,j)+xp*fin(i+1,j))
169 > + yp*(xpi*fin(i,j+1)+xp*fin(i+1,j+1))
180 sum=ypi*(fin(i+1,j)-fin(i,j))
181 > + yp*(fin(i+1,j+1)-fin(i,j+1))
182 fval(v,iadr)=sum*hxi(v)
192 sum=xpi*(fin(i,j+1)-fin(i,j))
193 > + xp*(fin(i+1,j+1)-fin(i+1,j))
194 fval(v,iadr)=sum*hyi(v)
203 sum=fin(i+1,j+1)-fin(i,j+1)-fin(i+1,j)+fin(i,j)
204 fval(v,iadr)=sum*hxi(v)*hyi(v)
216 subroutine r8pc2fcnx(ict,ivec,ivecd,
217 > fval,ii,jj,xparam,yparam,hx,hxi,hy,hyi,
223 INTEGER,
PARAMETER :: R8=selected_real_kind(12,100)
224 INTEGER ny,inf2,j,i,iadr
234 real*8 xparam(ivec),yparam
285 sum=ypi*(xpi*fin(i,j)+xp*fin(i+1,j))
286 > + yp*(xpi*fin(i,j+1)+xp*fin(i+1,j+1))
297 sum=ypi*(fin(i+1,j)-fin(i,j))
298 > + yp*(fin(i+1,j+1)-fin(i,j+1))
299 fval(v,iadr)=sum*hxi(v)
309 sum=xpi*(fin(i,j+1)-fin(i,j))
310 > + xp*(fin(i+1,j+1)-fin(i+1,j))
320 sum=fin(i+1,j+1)-fin(i,j+1)-fin(i+1,j)+fin(i,j)
321 fval(v,iadr)=sum*hxi(v)*hyi