您所在位置：主题：Error with curfit
Belinda124
[专家分：0] 发布于 2011-03-16 02:19:00
OPEN(11,file='test.dat',status="old",action='read',iostat=error)
if(error/=0)then
write(*,*)"open test.dat fail"
stop
end if
do l=-ni,ni
read(11,'(2(f9.6,3x))',END=20)xi(l),yi(l),s
end do
20 continue
call LSQUnivariateSpline(xi,yi,xn,xo,yo,wi,k,ni,nn,no,ier)
OPEN(10,file='testo.dat')
do l=-no,no
WRITE(10,'(2(f9.6,3x))')xo(l),yo(l)
end do
CLOSE(10)
if(error/=0)then
write(*,*)"open test.dat fail"
stop
end if
do l=-ni,ni
read(11,'(2(f9.6,3x))',END=20)xi(l),yi(l),s
end do
20 continue
call LSQUnivariateSpline(xi,yi,xn,xo,yo,wi,k,ni,nn,no,ier)
OPEN(10,file='testo.dat')
do l=-no,no
WRITE(10,'(2(f9.6,3x))')xo(l),yo(l)
end do
CLOSE(10)