1 subroutine mkherm3(fun,x,nx,y,ny,z,nz,fherm)
23 real fherm(0:7,nx,ny,nz)
40 > fun(x(ix),y(iy),z(iz),dfdx,dfdy,dfdz,
41 > d2fdxdy,d2fdxdz,d2fdydz,d3fdxyz)
42 fherm(1,ix,iy,iz)=dfdx
43 fherm(2,ix,iy,iz)=dfdy
44 fherm(3,ix,iy,iz)=dfdz
45 fherm(4,ix,iy,iz)=d2fdxdy
46 fherm(5,ix,iy,iz)=d2fdxdz
47 fherm(6,ix,iy,iz)=d2fdydz
48 fherm(7,ix,iy,iz)=d3fdxyz