1 subroutine r8dnherm2(x,nx,y,ny,fherm,nf2,ilinx,iliny,ier)
16 INTEGER,
PARAMETER :: R8=selected_real_kind(12,100)
17 INTEGER iliny,ier,ilinx,ierx,iery,iy,iyp,iym,ix,ixp,ixm
25 real*8 fherm(0:3,nf2,ny)
44 call r8splinck(x,nx,ilinx,1.0e-3_r8,ierx)
48 write(6,
'('' ?dnherm2: x axis not strict ascending'')')
51 call r8splinck(y,ny,iliny,1.0e-3_r8,iery)
55 write(6,
'('' ?dnherm2: y axis not strict ascending'')')
60 write(6,*)
'?dnherm2: fherm array dimension too small.'
77 zd=(fherm(0,ixp,iy)-fherm(0,ixm,iy))/(x(ixp)-x(ixm))
83 zd=(fherm(0,ix,iyp)-fherm(0,ix,iym))/(y(iyp)-y(iym))
89 fherm(3,ix,iy)=(fherm(0,ixp,iyp)-fherm(0,ixm,iyp)
90 > -fherm(0,ixp,iym)+fherm(0,ixm,iym))/
91 > ((x(ixp)-x(ixm))*(y(iyp)-y(iym)))