w4kek/WSJT-X
1
0
Fork 0
mirror of https://github.com/saitohirga/WSJT-X.git synced 2025-03-28 15:09:02 -04:00
Code Issues Packages Projects Releases Wiki Activity

File 'plotspec' in execution directory ==> save channel-gain spectrum to fort.52.

Browse Source
This commit is contained in:
Joe Taylor 2020-07-15 13:15:38 -04:00
parent fe86a6562c
commit ff0d31986f
1 changed files with 18 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

23
lib/fst240_decode.f90
View File

@ -809,11 +809,15 @@ contains
subroutine write_ref(itone,iwave,nsps,nmax,ndown,hmod,i0,fc) subroutine write_ref(itone,iwave,nsps,nmax,ndown,hmod,i0,fc)
complex cwave(nmax) complex cwave(nmax)
complex c(0:1440000-1) complex, allocatable :: c(:)
integer itone(160) integer itone(160)
integer*2 iwave(nmax) integer*2 iwave(nmax)
integer hmod integer hmod
data ncall/0/
save ncall
ncall=ncall+1
allocate( c(0:nmax-1) )
wave=0 wave=0
fsample=12000.0 fsample=12000.0
nsym=160 nsym=160
@ -829,17 +833,26 @@ contains
c=fac*float(iwave)*conjg(cwave) c=fac*float(iwave)*conjg(cwave)
call four2a(c,nmax,1,-1,1) !Forward c2c FFT call four2a(c,nmax,1,-1,1) !Forward c2c FFT
df=12000.0/nmax df=12000.0/nmax
ia=-10.0/df ia=-10.1/df
ib=10.0/df ib=10.1/df
smax=0.
do i=ia,ib do i=ia,ib
j=i j=i
if(j.lt.0) j=i+nmax if(j.lt.0) j=i+nmax
s=real(c(j))**2 + aimag(c(j))**2 s=real(c(j))**2 + aimag(c(j))**2
smax=max(s,smax)
enddo
do i=ia,ib
j=i
if(j.lt.0) j=i+nmax
s=(real(c(j))**2 + aimag(c(j))**2)/smax
s=s + ncall-1
f=i*df f=i*df
write(52,1010) f,s,db(s) write(52,1010) f,s,db(s)
1010 format(f10.3,e12.3,f10.3) 1010 format(f12.6,f12.6,f10.3)
enddo enddo
! close(52)
return return
end subroutine write_ref end subroutine write_ref