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

Merge branch 'feat-fst280' of bitbucket.org:k1jt/wsjtx into feat-fst280

Browse Source
This commit is contained in:
Bill Somerville 2020-07-27 19:54:08 +01:00
commit 5bad5dbd73
No known key found for this signature in database
GPG Key ID: D864B06D1E81618F
1 changed files with 17 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

17
lib/fst4_decode.f90
View File

@ -821,6 +821,7 @@ contains
complex, allocatable :: cwave(:) !Reconstructed complex signal complex, allocatable :: cwave(:) !Reconstructed complex signal
complex, allocatable :: g(:) !Channel gain, g(t) in QEX paper complex, allocatable :: g(:) !Channel gain, g(t) in QEX paper
real,allocatable :: ss(:) !Computed power spectrum of g(t) real,allocatable :: ss(:) !Computed power spectrum of g(t)
real,allocatable,save :: ssavg(:) !Computed power spectrum of g(t)
integer itone(160) !Tones for this message integer itone(160) !Tones for this message
integer*2 iwave(nmax) !Raw Rx data integer*2 iwave(nmax) !Raw Rx data
integer hmod !Modulation index integer hmod !Modulation index
@ -907,6 +908,22 @@ contains
1010 format(f12.6,f12.6) 1010 format(f12.6,f12.6)
enddo 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 return
end subroutine write_ref end subroutine write_ref