Correct flawed initial values for deltaa(1:2) in polfit.f90.

This commit is contained in:
Joe Taylor 2021-07-05 15:11:35 -04:00
parent 808b5f556d
commit 84e07c346a

View File

@ -17,7 +17,7 @@ subroutine polfit(y,npts,a)
ipk=maxloc(y) ipk=maxloc(y)
a(3)=(ipk(1)-1)*45.0 a(3)=(ipk(1)-1)*45.0
deltaa(1:2)=0.1*(a(2)-a(1)) deltaa(1:2)=0.1*a(2)
deltaa(3)=10.0 deltaa(3)=10.0
nterms=3 nterms=3