1 subroutine genxpkg(nx,x,xpkg,iper,imsg,itol,ztol,ialg,ier)
83 write(6,*) .ge.
' %genxpkg: nx2 required!'
91 if(ialgu.eq.0) ialgu=3
92 if(iabs(ialgu).gt.3) ialgu=3
102 ztola=max(abs(x(1)),abs(x(nx)))*ztolr
134 xpkg(nx,2)=(x(nx)-x(1))/(nx-1)
138 if((ier.eq.0).and.(ix.lt.nx))
then
139 if(x(ix+1).le.x(ix))
then
141 write(6,*)
' %genxpkg: x axis not strict ascending!'
148 if(abs(zh-xpkg(nx,2)).gt.ztola) xpkg(3,4)=1.0
162 xpkg(nx,3)=1.0/xpkg(nx,2)
169 if(xpkg(3,4).eq.0.0)
then
172 xpkg(ixp,1)=xpkg(1,1)+ix*xpkg(nx,2)
185 if(xpkg(3,4).ne.0.0)
then
195 if(abs(ialgu).eq.3)
then
203 xtest=xtest+xpkg(nx,2)
205 if((xpkg(itest,1).le.xtest).and.
206 > (xtest.le.xpkg(itest+1,1)))
then
207 xpkg(i,2)=itest+(xtest-xpkg(itest,1))/
208 > (xpkg(itest+1,1)-xpkg(itest,1))