173 implicit double precision(a-h,o-z)
177 common/gstuff/glamn,dgla,glamx,glomn,dglo,glomx,nla,nlo,nclip
179 if(gla.lt.glamn.or.gla.gt.glamx.or.
180 * glo.lt.glomn.or.glo.gt.glomx)
then 187 ix=(glo-glomn)/dglo+1.d0
200 x=(glo-dglo*(ix-1)-glomn)/dglo
204 if(x.lt.0.5.and.ix.gt.1)
then 210 if(x.lt.0..or.x.gt.2.) stop 55555
214 jy=(gla-glamn)/dgla+1.d0
223 y=(gla-dgla*(jy-1)-glamn)/dgla
224 if(y.lt.0.5.and.jy.gt.1)
then 230 if(y.lt.0..or.y.gt.2.) stop 33333
232 fx0=
qterp1(x, h(jy ,ix ), h(jy ,ix1), h(jy ,ix2))
233 fx1=
qterp1(x, h(jy1,ix ), h(jy1,ix1), h(jy1,ix2))
234 fx2=
qterp1(x, h(jy2,ix ), h(jy2,ix1), h(jy2,ix2))
235 val=
qterp2(y, fx0 , fx1 , fx2 )
real function qterp2(x, f0, f1, f2)
real function qterp1(x, if0, if1, if2)