Fix an out-of-bounds error.

This commit is contained in:
Steven Franke 2020-06-26 07:24:11 -05:00
parent 3d8ea2df5f
commit 37d5af9311

View File

@ -520,6 +520,9 @@ contains
s(i)=s(i) + real(c_bigfft(j))**2 + aimag(c_bigfft(j))**2
enddo
enddo
if(ina-hmod*3.lt.1) ina=1+hmod*3
if(inb+hmod*3.gt.18000) inb=18000-hmod*3
s2=0.
do i=ina,inb
s2(i)=s(i-hmod*3) + s(i-hmod) +s(i+hmod) +s(i+hmod*3)
enddo