w4kek/WSJT-X
1
0
Fork 0
mirror of https://github.com/saitohirga/WSJT-X.git synced 2025-10-27 02:50:39 -04:00
Code Issues Packages Projects Releases Wiki Activity

Deal with the FreqCal edge case of FTol>DF

Browse Source 
Reduce the tolerance to +-DF in this case to avoid array bounds error.

git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/wsjtx@7470 ab8295b8-cf94-4d9e-aec4-7959e3be5d79
This commit is contained in:
Bill Somerville 2017-01-09 22:54:47 +00:00
parent fb2592453d
commit b43567c38c
1 changed files with 8 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
lib/freqcal.f90
View File

@ -3,7 +3,7 @@ subroutine freqcal(id2,k,nkhz,noffset,ntol,line)
parameter (NZ=30*12000,NFFT=55296,NH=NFFT/2)
integer*2 id2(0:NZ-1)
real x(0:NFFT-1)
real s(NH)
real s(0:NH)
character line*80,cflag*1
complex cx(0:NH)
equivalence (x,cx)
@ -17,8 +17,13 @@ subroutine freqcal(id2,k,nkhz,noffset,ntol,line)
x=0.001*id2(k-NFFT:k-1)
call four2a(x,NFFT,1,-1,0) !Compute spectrum, r2c
df=12000.0/NFFT
ia=nint((noffset-ntol)/df)
ib=nint((noffset+ntol)/df)
if (ntol.gt.noffset) then
ia=0
ib=nint((noffset*2)/df)
else
ia=nint((noffset-ntol)/df)
ib=nint((noffset+ntol)/df)
endif
smax=0.
s=0.
do i=ia,ib