35 subroutine bicubic(z,glamn,glomn,dla,dlo,
36 *nla,nlo,maxla,maxlo,xla,xlo,val)
39 implicit real*8 (a-h,o-z)
43 real*4 x,y,fx0,fx1,fx2,fx3,cubterp
52 if(jla.gt.nla-3)jla = nla-3
61 if(jlo.gt.nlo-3)jlo = nlo-3
70 x=(xlo-dlo*(jlo-1)-glomn)/dlo
77 *
'FATAL in bicubic: x<0 : ',f20.10)
84 *
'FATAL in bicubic: x>3 : ',f20.10)
93 y=(xla-dla*(jla-1)-glamn)/dla
100 *
'FATAL in bicubic: y<0 : ',f20.10)
107 *
'FATAL in bicubic: y>3 : ',f20.10)
120 fx0=cubterp(x,z(jla ,jlo ),z(jla ,jlo+1),
121 * z(jla ,jlo+2),z(jla ,jlo+3))
123 fx1=cubterp(x,z(jla+1,jlo ),z(jla+1,jlo+1),
124 * z(jla+1,jlo+2),z(jla+1,jlo+3))
126 fx2=cubterp(x,z(jla+2,jlo ),z(jla+2,jlo+1),
127 * z(jla+2,jlo+2),z(jla+2,jlo+3))
129 fx3=cubterp(x,z(jla+3,jlo ),z(jla+3,jlo+1),
130 * z(jla+3,jlo+2),z(jla+3,jlo+3))
132 val=dble(cubterp(y,fx0,fx1,fx2,fx3))
148 1001
format(f14.10,1x,f14.10,1x,i8,1x,f8.5)
149 1002
format(f14.10,1x,f14.10,1x,i8,1x,f8.5)
150 1003
format(4(f15.8,1x))
subroutine bicubic(z, glamn, glomn, dla, dlo, nla, nlo, maxla, maxlo, xla, xlo, val)
Subroutine to perform a 2-D cubic ("bicubic") interpolation.