w4kek/WSJT-X
1
0
Fork 0
mirror of https://github.com/saitohirga/WSJT-X.git synced 2024-11-22 04:11:16 -05:00
Code Issues Packages Projects Releases Wiki Activity

Protect against bounds error in dopspread.f90. Remove the temporary ssavg code.

Browse Source
This commit is contained in:
Joe Taylor 2020-07-29 09:48:19 -04:00
parent 2ed0dafd0d
commit 902c0b7a59
1 changed files with 7 additions and 22 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

29
lib/fst4_decode.f90
View File

@ -545,7 +545,7 @@ contains
inquire(file='plotspec',exist=ex)
fmid=-999.0
if(ex) then
call write_ref(itone,iwave,nsps,nmax,ndown,hmod, &
call dopspread(itone,iwave,nsps,nmax,ndown,hmod, &
isbest,fc_synced,fmid,w50)
endif
xsig=0
@ -813,7 +813,7 @@ contains
return
end subroutine get_candidates_fst4
subroutine write_ref(itone,iwave,nsps,nmax,ndown,hmod,i0,fc,fmid,w50)
subroutine dopspread(itone,iwave,nsps,nmax,ndown,hmod,i0,fc,fmid,w50)
! On "plotspec" special request, compute Doppler spread for a decoded signal
@ -902,29 +902,14 @@ contains
fmid=xi2*df !Frequency midpoint of signal powere
do i=-ia,ia !Save the spectrum for plotting
f=i*df
y=0.99*ss(i+nint(xi2)) + ncall-1
write(52,1010) f,y
y=ncall-1
j=i+nint(xi2)
if(abs(j*df).lt.10.0) y=0.99*ss(i+nint(xi2)) + ncall-1
write(52,1010) i*df,y
1010 format(f12.6,f12.6)
enddo
if(nsps.eq.720) then
ia=101.0/df
if(ncall.eq.1) then
allocate(ssavg(-ia:ia))
ssavg=0.
endif
rewind 53
do i=-ia,ia !Find smax in +/- 1 Hz around 0.
j=i
if(j.lt.0) j=i+nfft
ssavg(i)=ssavg(i) + real(g(j))**2 + aimag(g(j))**2
write(53,1020) i*df,ssavg(i)
1020 format(f12.6,e12.3)
enddo
endif
return
end subroutine write_ref
end subroutine dopspread
end module fst4_decode