1 subroutine r8mkherm3(fun,x,nx,y,ny,z,nz,fherm)
11 INTEGER,
PARAMETER :: R8=selected_real_kind(12,100)
12 INTEGER ny,nz,nx,iz,iy,ix
15 real*8 fun,dfdx,dfdy,dfdz,d2fdxdy,d2fdxdz,d2fdydz,d3fdxyz
32 real*8 fherm(0:7,nx,ny,nz)
49 > fun(x(ix),y(iy),z(iz),dfdx,dfdy,dfdz,
50 > d2fdxdy,d2fdxdz,d2fdydz,d3fdxyz)
51 fherm(1,ix,iy,iz)=dfdx
52 fherm(2,ix,iy,iz)=dfdy
53 fherm(3,ix,iy,iz)=dfdz
54 fherm(4,ix,iy,iz)=d2fdxdy
55 fherm(5,ix,iy,iz)=d2fdxdz
56 fherm(6,ix,iy,iz)=d2fdydz
57 fherm(7,ix,iy,iz)=d3fdxyz