1 subroutine dnherm3(x,nx,y,ny,z,nz,fherm,nf2,nf3,
2 > ilinx,iliny,ilinz,ier)
14 integer nx,ny,nz,nf2,nf3
18 real fherm(0:7,nf2,nf3,nz)
48 call splinck(x,nx,ilinx,1.0e-3,ierx)
52 write(6,
'('' ?dnherm3: x axis not strict ascending'')')
55 call splinck(y,ny,iliny,1.0e-3,iery)
59 write(6,
'('' ?dnherm3: y axis not strict ascending'')')
62 call splinck(z,nz,ilinz,1.0e-3,ierz)
66 write(6,
'('' ?dnherm3: z axis not strict ascending'')')
71 write(6,*)
'?dnherm3: fherm (x) array dimension too small.'
76 write(6,*)
'?dnherm3: fherm (y) array dimension too small.'
100 zd=(fherm(0,ixp,iy,iz)-fherm(0,ixm,iy,iz))/
107 zd=(fherm(0,ix,iyp,iz)-fherm(0,ix,iym,iz))/
114 zd=(fherm(0,ix,iy,izp)-fherm(0,ix,iy,izm))/
122 > (fherm(0,ixp,iyp,iz)-fherm(0,ixm,iyp,iz)
123 > -fherm(0,ixp,iym,iz)+fherm(0,ixm,iym,iz))/
124 > ((x(ixp)-x(ixm))*(y(iyp)-y(iym)))
129 > (fherm(0,ixp,iy,izp)-fherm(0,ixm,iy,izp)
130 > -fherm(0,ixp,iy,izm)+fherm(0,ixm,iy,izm))/
131 > ((x(ixp)-x(ixm))*(z(izp)-z(izm)))
136 > (fherm(0,ix,iyp,izp)-fherm(0,ix,iym,izp)
137 > -fherm(0,ix,iyp,izm)+fherm(0,ix,iym,izm))/
138 > ((y(iyp)-y(iym))*(z(izp)-z(izm)))
143 > ((fherm(0,ixp,iyp,izp)-fherm(0,ixp,iym,izp)
144 > -fherm(0,ixp,iyp,izm)+fherm(0,ixp,iym,izm))-
145 > (fherm(0,ixm,iyp,izp)-fherm(0,ixm,iym,izp)
146 > -fherm(0,ixm,iyp,izm)+fherm(0,ixm,iym,izm)))/
147 > ((x(ixp)-x(ixm))*(y(iyp)-y(iym))*(z(izp)-z(izm)))