1 subroutine pc3ev(xget,yget,zget,x,nx,y,ny,z,nz,
3 > f,inf2,inf3,ict,fval,ier)
78 real xparam,yparam,zparam
92 call herm3xyz(xget,yget,zget,x,nx,y,ny,z,nz,ilinx,iliny,ilinz,
93 > i,j,k,xparam,yparam,zparam,hx,hxi,hy,hyi,hz,hzi,ier)
97 > fval,i,j,k,xparam,yparam,zparam,
98 > hx,hxi,hy,hyi,hz,hzi,
107 subroutine pc3fcn(ict,ivec,ivecd,
108 > fval,ii,jj,kk,xparam,yparam,zparam,
109 > hx,hxi,hy,hyi,hz,hzi,
116 integer ii(ivec),jj(ivec),kk(ivec)
117 real xparam(ivec),yparam(ivec),zparam(ivec)
120 real hx(ivec),hy(ivec),hz(ivec)
121 real hxi(ivec),hyi(ivec),hzi(ivec)
124 real fin(inf2,inf3,nz)
174 > xpi*(ypi*fin(i,j,k) +yp*fin(i,j+1,k))+
175 > xp*(ypi*fin(i+1,j,k)+yp*fin(i+1,j+1,k)))
177 > xpi*(ypi*fin(i,j,k+1) +yp*fin(i,j+1,k+1))+
178 > xp*(ypi*fin(i+1,j,k+1)+yp*fin(i+1,j+1,k+1)))
189 > -(ypi*fin(i,j,k) +yp*fin(i,j+1,k))
190 > +(ypi*fin(i+1,j,k)+yp*fin(i+1,j+1,k)))
192 > -(ypi*fin(i,j,k+1) +yp*fin(i,j+1,k+1))
193 > +(ypi*fin(i+1,j,k+1)+yp*fin(i+1,j+1,k+1)))
195 fval(v,iadr)=sum*hxi(v)
204 > xpi*(-fin(i,j,k) +fin(i,j+1,k))+
205 > xp*(-fin(i+1,j,k)+fin(i+1,j+1,k)))
207 > xpi*(-fin(i,j,k+1) +fin(i,j+1,k+1))+
208 > xp*(-fin(i+1,j,k+1)+fin(i+1,j+1,k+1)))
210 fval(v,iadr)=sum*hyi(v)
219 > xpi*(ypi*fin(i,j,k) +yp*fin(i,j+1,k))+
220 > xp*(ypi*fin(i+1,j,k)+yp*fin(i+1,j+1,k)))
222 > xpi*(ypi*fin(i,j,k+1) +yp*fin(i,j+1,k+1))+
223 > xp*(ypi*fin(i+1,j,k+1)+yp*fin(i+1,j+1,k+1)))
225 fval(v,iadr)=sum*hzi(v)
234 > -(-fin(i,j,k) +fin(i,j+1,k))
235 > +(-fin(i+1,j,k)+fin(i+1,j+1,k)))
237 > -(-fin(i,j,k+1) +fin(i,j+1,k+1))
238 > +(-fin(i+1,j,k+1)+fin(i+1,j+1,k+1)))
240 fval(v,iadr)=sum*hxi(v)*hyi(v)
249 > -(ypi*fin(i,j,k) +yp*fin(i,j+1,k))
250 > +(ypi*fin(i+1,j,k)+yp*fin(i+1,j+1,k)))
252 > -(ypi*fin(i,j,k+1) +yp*fin(i,j+1,k+1))
253 > +(ypi*fin(i+1,j,k+1)+yp*fin(i+1,j+1,k+1)))
255 fval(v,iadr)=sum*hxi(v)*hzi(v)
264 > xpi*(-fin(i,j,k) +fin(i,j+1,k))+
265 > xp*(-fin(i+1,j,k)+fin(i+1,j+1,k)))
267 > xpi*(-fin(i,j,k+1) +fin(i,j+1,k+1))+
268 > xp*(-fin(i+1,j,k+1)+fin(i+1,j+1,k+1)))
270 fval(v,iadr)=sum*hyi(v)*hzi(v)
279 > -(-fin(i,j,k) +fin(i,j+1,k))
280 > +(-fin(i+1,j,k)+fin(i+1,j+1,k)))
282 > -(-fin(i,j,k+1) +fin(i,j+1,k+1))
283 > +(-fin(i+1,j,k+1)+fin(i+1,j+1,k+1)))
285 fval(v,iadr)=sum*hxi(v)*hyi(v)*hzi(v)
297 subroutine pc3fcnx(ict,ivec,ivecd,
298 > fval,ii,jj,kk,xparam,yparam,zparam,
299 > hx,hxi,hy,hyi,hz,hzi,
306 integer ii(ivec),jj,kk
307 real xparam(ivec),yparam,zparam
311 real hxi(ivec),hyi,hzi
314 real fin(inf2,inf3,nz)
365 > xpi*(ypi*fin(i,j,k) +yp*fin(i,j+1,k))+
366 > xp*(ypi*fin(i+1,j,k)+yp*fin(i+1,j+1,k)))
368 > xpi*(ypi*fin(i,j,k+1) +yp*fin(i,j+1,k+1))+
369 > xp*(ypi*fin(i+1,j,k+1)+yp*fin(i+1,j+1,k+1)))
380 > -(ypi*fin(i,j,k) +yp*fin(i,j+1,k))
381 > +(ypi*fin(i+1,j,k)+yp*fin(i+1,j+1,k)))
383 > -(ypi*fin(i,j,k+1) +yp*fin(i,j+1,k+1))
384 > +(ypi*fin(i+1,j,k+1)+yp*fin(i+1,j+1,k+1)))
386 fval(v,iadr)=sum*hxi(v)
395 > xpi*(-fin(i,j,k) +fin(i,j+1,k))+
396 > xp*(-fin(i+1,j,k)+fin(i+1,j+1,k)))
398 > xpi*(-fin(i,j,k+1) +fin(i,j+1,k+1))+
399 > xp*(-fin(i+1,j,k+1)+fin(i+1,j+1,k+1)))
410 > xpi*(ypi*fin(i,j,k) +yp*fin(i,j+1,k))+
411 > xp*(ypi*fin(i+1,j,k)+yp*fin(i+1,j+1,k)))
413 > xpi*(ypi*fin(i,j,k+1) +yp*fin(i,j+1,k+1))+
414 > xp*(ypi*fin(i+1,j,k+1)+yp*fin(i+1,j+1,k+1)))
425 > -(-fin(i,j,k) +fin(i,j+1,k))
426 > +(-fin(i+1,j,k)+fin(i+1,j+1,k)))
428 > -(-fin(i,j,k+1) +fin(i,j+1,k+1))
429 > +(-fin(i+1,j,k+1)+fin(i+1,j+1,k+1)))
431 fval(v,iadr)=sum*hxi(v)*hyi
440 > -(ypi*fin(i,j,k) +yp*fin(i,j+1,k))
441 > +(ypi*fin(i+1,j,k)+yp*fin(i+1,j+1,k)))
443 > -(ypi*fin(i,j,k+1) +yp*fin(i,j+1,k+1))
444 > +(ypi*fin(i+1,j,k+1)+yp*fin(i+1,j+1,k+1)))
446 fval(v,iadr)=sum*hxi(v)*hzi
455 > xpi*(-fin(i,j,k) +fin(i,j+1,k))+
456 > xp*(-fin(i+1,j,k)+fin(i+1,j+1,k)))
458 > xpi*(-fin(i,j,k+1) +fin(i,j+1,k+1))+
459 > xp*(-fin(i+1,j,k+1)+fin(i+1,j+1,k+1)))
461 fval(v,iadr)=sum*hyi*hzi
470 > -(-fin(i,j,k) +fin(i,j+1,k))
471 > +(-fin(i+1,j,k)+fin(i+1,j+1,k)))
473 > -(-fin(i,j,k+1) +fin(i,j+1,k+1))
474 > +(-fin(i+1,j,k+1)+fin(i+1,j+1,k+1)))
476 fval(v,iadr)=sum*hxi(v)*hyi*hzi