mirror of
https://github.com/saitohirga/WSJT-X.git
synced 2025-11-03 13:30:52 -05:00
Re-organize some of the FT8 a7 code.
This commit is contained in:
Joe Taylor
parent
4a68d3c4cf
commit
9aa4786f11
@ -336,7 +336,7 @@ set (wsjt_FSRCS
|
|||||||
lib/timer_module.f90
|
lib/timer_module.f90
|
||||||
lib/wavhdr.f90
|
lib/wavhdr.f90
|
||||||
lib/qra/q65/q65_encoding_modules.f90
|
lib/qra/q65/q65_encoding_modules.f90
|
||||||
lib/ft8_a7.f90
|
lib/ft8/ft8_a7.f90
|
||||||
|
|
||||||
# remaining non-module sources
|
# remaining non-module sources
|
||||||
lib/addit.f90
|
lib/addit.f90
|
||||||
@ -422,7 +422,6 @@ set (wsjt_FSRCS
|
|||||||
lib/freqcal.f90
|
lib/freqcal.f90
|
||||||
lib/ft8/ft8apset.f90
|
lib/ft8/ft8apset.f90
|
||||||
lib/ft8/ft8b.f90
|
lib/ft8/ft8b.f90
|
||||||
lib/ft8/ft8c.f90
|
|
||||||
lib/ft8/ft8code.f90
|
lib/ft8/ft8code.f90
|
||||||
lib/ft8/ft8_downsample.f90
|
lib/ft8/ft8_downsample.f90
|
||||||
lib/ft8/ft8sim.f90
|
lib/ft8/ft8sim.f90
|
||||||
|
|||||||
339
lib/ft8/ft8_a7.f90
Normal file
339
lib/ft8/ft8_a7.f90
Normal file
@ -0,0 +1,339 @@
|
|||||||
|
module ft8_a7
|
||||||
|
|
||||||
|
parameter(MAXDEC=100)
|
||||||
|
|
||||||
|
! For the following three arrays
|
||||||
|
! First index i=decode number in this sequence
|
||||||
|
! Second index j=0 or 1 for even or odd sequence
|
||||||
|
! Third index k=0 or 1 for previous or current tally for this j
|
||||||
|
real dt0(MAXDEC,0:1,0:1) !dt0(i,j,k)
|
||||||
|
real f0(MAXDEC,0:1,0:1) !f0(i,j,k)
|
||||||
|
character*37 msg0(MAXDEC,0:1,0:1) !msg0(i,j,k)
|
||||||
|
|
||||||
|
integer itone_a7(79)
|
||||||
|
integer jseq !even=0, odd=1
|
||||||
|
integer ndec(0:1,0:1) !ndec(j,k)
|
||||||
|
data ndec/4*0/,jseq/0/
|
||||||
|
|
||||||
|
contains
|
||||||
|
|
||||||
|
subroutine ft8_a7_save(nutc,dt,f,msg)
|
||||||
|
|
||||||
|
use packjt77
|
||||||
|
character*37 msg,msg1
|
||||||
|
character*13 w(19)
|
||||||
|
character*4 g4
|
||||||
|
integer nw(19)
|
||||||
|
logical isgrid4
|
||||||
|
|
||||||
|
! Statement function:
|
||||||
|
isgrid4(g4)=(len_trim(g4).eq.4 .and. &
|
||||||
|
ichar(g4(1:1)).ge.ichar('A') .and. ichar(g4(1:1)).le.ichar('R') .and. &
|
||||||
|
ichar(g4(2:2)).ge.ichar('A') .and. ichar(g4(2:2)).le.ichar('R') .and. &
|
||||||
|
ichar(g4(3:3)).ge.ichar('0') .and. ichar(g4(3:3)).le.ichar('9') .and. &
|
||||||
|
ichar(g4(4:4)).ge.ichar('0') .and. ichar(g4(4:4)).le.ichar('9'))
|
||||||
|
|
||||||
|
j=mod(nutc/5,2) !j is 0 or 1 for odd/even sequence
|
||||||
|
jseq=j
|
||||||
|
|
||||||
|
! Add this decode to current table for this sequence
|
||||||
|
ndec(j,1)=ndec(j,1)+1 !Number of decodes in this sequence
|
||||||
|
i=ndec(j,1) !i is pointer to new table entry
|
||||||
|
if(i.ge.MAXDEC-1) return !Prevent table overflow
|
||||||
|
if(index(msg,'<...>').ge.1) return !Don't save an unknown hashcall
|
||||||
|
|
||||||
|
dt0(i,j,1)=dt !Save dt in table
|
||||||
|
f0(i,j,1)=f !Save f in table
|
||||||
|
f0(i+1,j,1)=-99.0 !Flag after last entry in current table
|
||||||
|
call split77(msg,nwords,nw,w) !Parse msg into words
|
||||||
|
msg0(i,j,1)=trim(w(1))//' '//trim(w(2))
|
||||||
|
if(w(1)(1:3).eq.'CQ ' .and. nw(2).le.2) then
|
||||||
|
msg0(i,j,1)='CQ '//trim(w(2))//' '//trim(w(3))
|
||||||
|
endif
|
||||||
|
msg1=msg0(i,j,1) !Message without grid
|
||||||
|
nn=len(trim(msg1)) !Message length without grid
|
||||||
|
if(isgrid4(w(nwords))) msg0(i,j,1)=trim(msg0(i,j,1))//' '//trim(w(nwords))
|
||||||
|
|
||||||
|
! If a transmission at this frequency with this message fragment
|
||||||
|
! was decoded in the previous sequence, flag it as "DO NOT USE" because
|
||||||
|
! we have already decoded that station's next transmission.
|
||||||
|
|
||||||
|
call split77(msg1,nwords,nw,w) !Parse msg into words
|
||||||
|
do i=1,ndec(j,0)
|
||||||
|
if(f0(i,j,0).le.-98.0) cycle
|
||||||
|
i2=index(msg0(i,j,0),' '//trim(w(2)))
|
||||||
|
if(abs(f-f0(i,j,0)).lt.2.0 .and. i2.ge.3) then
|
||||||
|
f0(i,j,0)=-98.0 !Remove from list of to-be-tried a7 decodes
|
||||||
|
endif
|
||||||
|
enddo
|
||||||
|
|
||||||
|
return
|
||||||
|
end subroutine ft8_a7_save
|
||||||
|
|
||||||
|
subroutine ft8_a7d(dd0,newdat,call_1,call_2,grid4,xdt,f1,nharderrors,dmin, &
|
||||||
|
msg37,xsnr)
|
||||||
|
|
||||||
|
use crc
|
||||||
|
use timer_module, only: timer
|
||||||
|
use packjt77
|
||||||
|
include 'ft8_params.f90'
|
||||||
|
parameter(NP2=2812)
|
||||||
|
character*37 msg37,msg,msgsent,msgbest
|
||||||
|
character*12 call_1,call_2
|
||||||
|
character*4 grid4
|
||||||
|
real a(5)
|
||||||
|
real s8(0:7,NN)
|
||||||
|
real s2(0:511)
|
||||||
|
real bmeta(174),bmetb(174),bmetc(174),bmetd(174)
|
||||||
|
real llra(174),llrb(174),llrc(174),llrd(174),llrbest(174) !Soft symbols
|
||||||
|
real dd0(15*12000)
|
||||||
|
real ss(9)
|
||||||
|
real rcw(174)
|
||||||
|
integer*1 cw(174)
|
||||||
|
integer*1 msgbits(77)
|
||||||
|
integer*1 nxor(174),hdec(174)
|
||||||
|
integer itone(NN)
|
||||||
|
integer icos7(0:6),ip(1)
|
||||||
|
logical one(0:511,0:8)
|
||||||
|
integer graymap(0:7)
|
||||||
|
integer iloc(1)
|
||||||
|
complex cd0(0:3199)
|
||||||
|
complex ctwk(32)
|
||||||
|
complex csymb(32)
|
||||||
|
complex cs(0:7,NN)
|
||||||
|
logical std_1,std_2
|
||||||
|
logical first,newdat
|
||||||
|
data icos7/3,1,4,0,6,5,2/ ! Flipped w.r.t. original FT8 sync array
|
||||||
|
data first/.true./
|
||||||
|
data graymap/0,1,3,2,5,6,4,7/
|
||||||
|
save one
|
||||||
|
|
||||||
|
if(first) then
|
||||||
|
one=.false.
|
||||||
|
do i=0,511
|
||||||
|
do j=0,8
|
||||||
|
if(iand(i,2**j).ne.0) one(i,j)=.true.
|
||||||
|
enddo
|
||||||
|
enddo
|
||||||
|
first=.false.
|
||||||
|
endif
|
||||||
|
|
||||||
|
call stdcall(call_1,std_1)
|
||||||
|
if(call_1(1:3).eq.'CQ ') std_1=.true.
|
||||||
|
call stdcall(call_2,std_2)
|
||||||
|
|
||||||
|
nharderrors=-1
|
||||||
|
fs2=12000.0/NDOWN
|
||||||
|
dt2=1.0/fs2
|
||||||
|
twopi=8.0*atan(1.0)
|
||||||
|
delfbest=0.
|
||||||
|
ibest=0
|
||||||
|
|
||||||
|
call timer('ft8_down',0)
|
||||||
|
call ft8_downsample(dd0,newdat,f1,cd0) !Mix f1 to baseband and downsample
|
||||||
|
call timer('ft8_down',1)
|
||||||
|
|
||||||
|
i0=nint((xdt+0.5)*fs2) !Initial guess for start of signal
|
||||||
|
smax=0.0
|
||||||
|
do idt=i0-10,i0+10 !Search over +/- one quarter symbol
|
||||||
|
call sync8d(cd0,idt,ctwk,0,sync)
|
||||||
|
if(sync.gt.smax) then
|
||||||
|
smax=sync
|
||||||
|
ibest=idt
|
||||||
|
endif
|
||||||
|
enddo
|
||||||
|
|
||||||
|
! Now peak up in frequency
|
||||||
|
smax=0.0
|
||||||
|
do ifr=-5,5 !Search over +/- 2.5 Hz
|
||||||
|
delf=ifr*0.5
|
||||||
|
dphi=twopi*delf*dt2
|
||||||
|
phi=0.0
|
||||||
|
do i=1,32
|
||||||
|
ctwk(i)=cmplx(cos(phi),sin(phi))
|
||||||
|
phi=mod(phi+dphi,twopi)
|
||||||
|
enddo
|
||||||
|
call sync8d(cd0,ibest,ctwk,1,sync)
|
||||||
|
if( sync .gt. smax ) then
|
||||||
|
smax=sync
|
||||||
|
delfbest=delf
|
||||||
|
endif
|
||||||
|
enddo
|
||||||
|
a=0.0
|
||||||
|
a(1)=-delfbest
|
||||||
|
call twkfreq1(cd0,NP2,fs2,a,cd0)
|
||||||
|
f1=f1+delfbest !Improved estimate of DF
|
||||||
|
|
||||||
|
call timer('ft8_down',0)
|
||||||
|
call ft8_downsample(dd0,.false.,f1,cd0) !Mix f1 to baseband and downsample
|
||||||
|
call timer('ft8_down',1)
|
||||||
|
|
||||||
|
smax=0.0
|
||||||
|
do idt=-4,4 !Search over +/- one quarter symbol
|
||||||
|
call sync8d(cd0,ibest+idt,ctwk,0,sync)
|
||||||
|
ss(idt+5)=sync
|
||||||
|
enddo
|
||||||
|
smax=maxval(ss)
|
||||||
|
iloc=maxloc(ss)
|
||||||
|
ibest=iloc(1)-5+ibest
|
||||||
|
xdt=(ibest-1)*dt2 - 0.5
|
||||||
|
sync=smax
|
||||||
|
|
||||||
|
do k=1,NN
|
||||||
|
i1=ibest+(k-1)*32
|
||||||
|
csymb=cmplx(0.0,0.0)
|
||||||
|
if( i1.ge.0 .and. i1+31 .le. NP2-1 ) csymb=cd0(i1:i1+31)
|
||||||
|
call four2a(csymb,32,1,-1,1)
|
||||||
|
cs(0:7,k)=csymb(1:8)/1e3
|
||||||
|
s8(0:7,k)=abs(csymb(1:8))
|
||||||
|
enddo
|
||||||
|
|
||||||
|
! sync quality check
|
||||||
|
is1=0
|
||||||
|
is2=0
|
||||||
|
is3=0
|
||||||
|
do k=1,7
|
||||||
|
ip=maxloc(s8(:,k))
|
||||||
|
if(icos7(k-1).eq.(ip(1)-1)) is1=is1+1
|
||||||
|
ip=maxloc(s8(:,k+36))
|
||||||
|
if(icos7(k-1).eq.(ip(1)-1)) is2=is2+1
|
||||||
|
ip=maxloc(s8(:,k+72))
|
||||||
|
if(icos7(k-1).eq.(ip(1)-1)) is3=is3+1
|
||||||
|
enddo
|
||||||
|
! hard sync sum - max is 21
|
||||||
|
nsync=is1+is2+is3
|
||||||
|
! if(nsync .le. 6) return ! bail out
|
||||||
|
|
||||||
|
do nsym=1,3
|
||||||
|
nt=2**(3*nsym)
|
||||||
|
do ihalf=1,2
|
||||||
|
do k=1,29,nsym
|
||||||
|
if(ihalf.eq.1) ks=k+7
|
||||||
|
if(ihalf.eq.2) ks=k+43
|
||||||
|
amax=-1.0
|
||||||
|
do i=0,nt-1
|
||||||
|
i1=i/64
|
||||||
|
i2=iand(i,63)/8
|
||||||
|
i3=iand(i,7)
|
||||||
|
if(nsym.eq.1) then
|
||||||
|
s2(i)=abs(cs(graymap(i3),ks))
|
||||||
|
elseif(nsym.eq.2) then
|
||||||
|
s2(i)=abs(cs(graymap(i2),ks)+cs(graymap(i3),ks+1))
|
||||||
|
elseif(nsym.eq.3) then
|
||||||
|
s2(i)=abs(cs(graymap(i1),ks)+cs(graymap(i2),ks+1)+cs(graymap(i3),ks+2))
|
||||||
|
else
|
||||||
|
print*,"Error - nsym must be 1, 2, or 3."
|
||||||
|
endif
|
||||||
|
enddo
|
||||||
|
i32=1+(k-1)*3+(ihalf-1)*87
|
||||||
|
if(nsym.eq.1) ibmax=2
|
||||||
|
if(nsym.eq.2) ibmax=5
|
||||||
|
if(nsym.eq.3) ibmax=8
|
||||||
|
do ib=0,ibmax
|
||||||
|
bm=maxval(s2(0:nt-1),one(0:nt-1,ibmax-ib)) - &
|
||||||
|
maxval(s2(0:nt-1),.not.one(0:nt-1,ibmax-ib))
|
||||||
|
if(i32+ib .gt.174) cycle
|
||||||
|
if(nsym.eq.1) then
|
||||||
|
bmeta(i32+ib)=bm
|
||||||
|
den=max(maxval(s2(0:nt-1),one(0:nt-1,ibmax-ib)), &
|
||||||
|
maxval(s2(0:nt-1),.not.one(0:nt-1,ibmax-ib)))
|
||||||
|
if(den.gt.0.0) then
|
||||||
|
cm=bm/den
|
||||||
|
else ! erase it
|
||||||
|
cm=0.0
|
||||||
|
endif
|
||||||
|
bmetd(i32+ib)=cm
|
||||||
|
elseif(nsym.eq.2) then
|
||||||
|
bmetb(i32+ib)=bm
|
||||||
|
elseif(nsym.eq.3) then
|
||||||
|
bmetc(i32+ib)=bm
|
||||||
|
endif
|
||||||
|
enddo
|
||||||
|
enddo
|
||||||
|
enddo
|
||||||
|
enddo
|
||||||
|
call normalizebmet(bmeta,174)
|
||||||
|
call normalizebmet(bmetb,174)
|
||||||
|
call normalizebmet(bmetc,174)
|
||||||
|
call normalizebmet(bmetd,174)
|
||||||
|
|
||||||
|
scalefac=2.83
|
||||||
|
llra=scalefac*bmeta
|
||||||
|
llrb=scalefac*bmetb
|
||||||
|
llrc=scalefac*bmetc
|
||||||
|
llrd=scalefac*bmetd
|
||||||
|
|
||||||
|
! apmag=maxval(abs(llra))*1.01
|
||||||
|
|
||||||
|
MAXMSG=206
|
||||||
|
pbest=0.
|
||||||
|
dmin=1.e30
|
||||||
|
do imsg=1,MAXMSG
|
||||||
|
msg=trim(call_1)//' '//trim(call_2)
|
||||||
|
i=imsg
|
||||||
|
if(call_1(1:3).eq.'CQ ' .and. i.ne.5) msg='QQ0XYZ '//trim(call_2)
|
||||||
|
if(.not.std_1) then
|
||||||
|
if(i.eq.1 .or. i.ge.6) msg='<'//trim(call_1)//'> '//trim(call_2)
|
||||||
|
if(i.ge.2 .and. i.le.4) msg=trim(call_1)//' <'//trim(call_2)//'>'
|
||||||
|
else if(.not.std_2) then
|
||||||
|
if(i.le.4 .or. i.eq.6) msg='<'//trim(call_1)//'> '//trim(call_2)
|
||||||
|
if(i.ge.7) msg=trim(call_1)//' <'//trim(call_2)//'>'
|
||||||
|
endif
|
||||||
|
j0=len(trim(msg))+2
|
||||||
|
if(i.eq.2) msg(j0:j0+2)='RRR'
|
||||||
|
if(i.eq.3) msg(j0:j0+3)='RR73'
|
||||||
|
if(i.eq.4) msg(j0:j0+1)='73'
|
||||||
|
if(i.eq.5) then
|
||||||
|
if(std_2) then
|
||||||
|
msg='CQ '//trim(call_2)
|
||||||
|
if(call_1(3:3).eq.'_') msg=trim(call_1)//' '//trim(call_2)
|
||||||
|
if(grid4.ne.'RR73') msg=trim(msg)//' '//grid4
|
||||||
|
endif
|
||||||
|
if(.not.std_2) msg='CQ '//trim(call_2)
|
||||||
|
endif
|
||||||
|
if(i.eq.6 .and. std_2) msg(j0:j0+3)=grid4
|
||||||
|
if(i.ge.7) then
|
||||||
|
isnr = -50 + (i-7)/2
|
||||||
|
if(iand(i,1).eq.1) then
|
||||||
|
write(msg(j0:j0+2),'(i3.2)') isnr
|
||||||
|
if(msg(j0:j0).eq.' ') msg(j0:j0)='+'
|
||||||
|
else
|
||||||
|
write(msg(j0:j0+3),'("R",i3.2)') isnr
|
||||||
|
if(msg(j0+1:j0+1).eq.' ') msg(j0+1:j0+1)='+'
|
||||||
|
endif
|
||||||
|
endif
|
||||||
|
|
||||||
|
! Source-encode, then get codeword
|
||||||
|
i3=-1
|
||||||
|
n3=-1
|
||||||
|
call genft8(msg,i3,n3,msgsent,msgbits,itone)
|
||||||
|
call encode174_91(msgbits,cw)
|
||||||
|
rcw=2*cw-1
|
||||||
|
pa=sum(llra*rcw)
|
||||||
|
pb=sum(llrb*rcw)
|
||||||
|
pc=sum(llrc*rcw)
|
||||||
|
pd=sum(llrd*rcw)
|
||||||
|
|
||||||
|
hdec=0
|
||||||
|
where(llra.ge.0.0) hdec=1
|
||||||
|
nxor=ieor(hdec,cw)
|
||||||
|
da=sum(nxor*abs(llra))
|
||||||
|
|
||||||
|
if(da.lt.dmin) then
|
||||||
|
dmin=da
|
||||||
|
pbest=pa
|
||||||
|
msgbest=msgsent
|
||||||
|
llrbest=llra
|
||||||
|
nharderrors=count((2*cw-1)*llra.lt.0.0)
|
||||||
|
endif
|
||||||
|
|
||||||
|
enddo ! imsg
|
||||||
|
|
||||||
|
write(41,3041) nharderrors,pbest,dmin,trim(msgbest)
|
||||||
|
3041 format(i5,2f10.3,2x,a)
|
||||||
|
msg37=msgbest
|
||||||
|
|
||||||
|
return
|
||||||
|
end subroutine ft8_a7d
|
||||||
|
|
||||||
|
end module ft8_a7
|
||||||
211
lib/ft8_a7.f90
211
lib/ft8_a7.f90
@ -1,211 +0,0 @@
|
|||||||
module ft8_a7
|
|
||||||
|
|
||||||
parameter(MAXDEC=100)
|
|
||||||
|
|
||||||
! For the following three arrays
|
|
||||||
! First index i=decode number in this sequence
|
|
||||||
! Second index j=0 or 1 for even or odd sequence
|
|
||||||
! Third index k=0 or 1 for previous or current tally for this j
|
|
||||||
real dt0(MAXDEC,0:1,0:1) !dt0(i,j,k)
|
|
||||||
real f0(MAXDEC,0:1,0:1) !f0(i,j,k)
|
|
||||||
character*37 msg0(MAXDEC,0:1,0:1) !msg0(i,j,k)
|
|
||||||
|
|
||||||
integer itone_a7(79)
|
|
||||||
integer jseq !even=0, odd=1
|
|
||||||
integer ndec(0:1,0:1) !ndec(j,k)
|
|
||||||
data ndec/4*0/,jseq/0/
|
|
||||||
|
|
||||||
contains
|
|
||||||
|
|
||||||
subroutine ft8_a7_save(nutc,dt,f,msg)
|
|
||||||
|
|
||||||
use packjt77
|
|
||||||
character*37 msg,msg1
|
|
||||||
character*13 w(19)
|
|
||||||
character*4 g4
|
|
||||||
integer nw(19)
|
|
||||||
logical isgrid4
|
|
||||||
|
|
||||||
! Statement function:
|
|
||||||
isgrid4(g4)=(len_trim(g4).eq.4 .and. &
|
|
||||||
ichar(g4(1:1)).ge.ichar('A') .and. ichar(g4(1:1)).le.ichar('R') .and. &
|
|
||||||
ichar(g4(2:2)).ge.ichar('A') .and. ichar(g4(2:2)).le.ichar('R') .and. &
|
|
||||||
ichar(g4(3:3)).ge.ichar('0') .and. ichar(g4(3:3)).le.ichar('9') .and. &
|
|
||||||
ichar(g4(4:4)).ge.ichar('0') .and. ichar(g4(4:4)).le.ichar('9'))
|
|
||||||
|
|
||||||
j=mod(nutc/5,2) !j is 0 or 1 for odd/even sequence
|
|
||||||
jseq=j
|
|
||||||
|
|
||||||
! Add this decode to current table for this sequence
|
|
||||||
ndec(j,1)=ndec(j,1)+1 !Number of decodes in this sequence
|
|
||||||
i=ndec(j,1) !i is pointer to new table entry
|
|
||||||
if(i.ge.MAXDEC-1) return !Prevent table overflow
|
|
||||||
if(index(msg,'<...>').ge.1) return !Don't save an unknown hashcall
|
|
||||||
|
|
||||||
dt0(i,j,1)=dt !Save dt in table
|
|
||||||
f0(i,j,1)=f !Save f in table
|
|
||||||
f0(i+1,j,1)=-99.0 !Flag after last entry in current table
|
|
||||||
call split77(msg,nwords,nw,w) !Parse msg into words
|
|
||||||
msg0(i,j,1)=trim(w(1))//' '//trim(w(2))
|
|
||||||
if(w(1)(1:3).eq.'CQ ' .and. nw(2).le.2) then
|
|
||||||
msg0(i,j,1)='CQ '//trim(w(2))//' '//trim(w(3))
|
|
||||||
endif
|
|
||||||
msg1=msg0(i,j,1) !Message without grid
|
|
||||||
nn=len(trim(msg1)) !Message length without grid
|
|
||||||
if(isgrid4(w(nwords))) msg0(i,j,1)=trim(msg0(i,j,1))//' '//trim(w(nwords))
|
|
||||||
|
|
||||||
! If a transmission at this frequency with this message fragment
|
|
||||||
! was decoded in the previous sequence, flag it as "DO NOT USE" because
|
|
||||||
! we have already decoded that station's next transmission.
|
|
||||||
|
|
||||||
call split77(msg1,nwords,nw,w) !Parse msg into words
|
|
||||||
do i=1,ndec(j,0)
|
|
||||||
if(f0(i,j,0).le.-98.0) cycle
|
|
||||||
i2=index(msg0(i,j,0),' '//trim(w(2)))
|
|
||||||
if(abs(f-f0(i,j,0)).lt.2.0 .and. i2.ge.3) then
|
|
||||||
f0(i,j,0)=-98.0 !Remove from list of to-be-tried a7 decodes
|
|
||||||
endif
|
|
||||||
enddo
|
|
||||||
|
|
||||||
return
|
|
||||||
end subroutine ft8_a7_save
|
|
||||||
|
|
||||||
subroutine ft8_dec7(cd,xdt0,f0,msg0,xdt,xsnr,msgbest,snr7,snr7b)
|
|
||||||
|
|
||||||
! Get a7 (q3-style) decodes for FT8.
|
|
||||||
|
|
||||||
use packjt77
|
|
||||||
parameter(NN=79,NSPS=32)
|
|
||||||
parameter(NWAVE=NN*NSPS) !2528
|
|
||||||
parameter(NZ=3200,NLAGS=NZ-NWAVE)
|
|
||||||
parameter(MAXMSG=206)
|
|
||||||
character*12 call_1,call_2
|
|
||||||
character*13 w(19)
|
|
||||||
character*4 grid4
|
|
||||||
character*37 msg0,msg,msgbest,msgsent
|
|
||||||
character c77*77
|
|
||||||
complex cwave(0:NWAVE-1)
|
|
||||||
complex cd(0:NZ-1)
|
|
||||||
complex z
|
|
||||||
real xjunk(NWAVE)
|
|
||||||
real ccf(0:NLAGS-1)
|
|
||||||
real ccfmsg(MAXMSG)
|
|
||||||
integer itone(NN)
|
|
||||||
integer nw(19)
|
|
||||||
integer*1 msgbits(77)
|
|
||||||
logical std_1,std_2
|
|
||||||
|
|
||||||
if(xdt0.eq.-999.0) return !Silence compiler warning
|
|
||||||
|
|
||||||
snr7=0.
|
|
||||||
ccfmsg=0.
|
|
||||||
call split77(msg0,nwords,nw,w) !Parse msg0 into words
|
|
||||||
call_1=w(1)(1:12)
|
|
||||||
call_2=w(2)(1:12)
|
|
||||||
grid4=w(3)(1:4)
|
|
||||||
if(call_1(1:3).eq.'CQ_') call_1(3:3)=' '
|
|
||||||
|
|
||||||
call stdcall(call_1,std_1)
|
|
||||||
if(call_1(1:3).eq.'CQ ') std_1=.true.
|
|
||||||
call stdcall(call_2,std_2)
|
|
||||||
|
|
||||||
fs=200.0 !Sample rate (Hz)
|
|
||||||
bt=2.0
|
|
||||||
ccfbest=0.
|
|
||||||
lagbest=-1
|
|
||||||
imsgbest=1
|
|
||||||
|
|
||||||
do imsg=1,MAXMSG
|
|
||||||
msg=trim(call_1)//' '//trim(call_2)
|
|
||||||
i=imsg
|
|
||||||
if(call_1(1:3).eq.'CQ ' .and. i.ne.5) msg='CQ0XYZ '//trim(call_2)
|
|
||||||
if(.not.std_1) then
|
|
||||||
if(i.eq.1 .or. i.ge.6) msg='<'//trim(call_1)//'> '//trim(call_2)
|
|
||||||
if(i.ge.2 .and. i.le.4) msg=trim(call_1)//' <'//trim(call_2)//'>'
|
|
||||||
else if(.not.std_2) then
|
|
||||||
if(i.le.4 .or. i.eq.6) msg='<'//trim(call_1)//'> '//trim(call_2)
|
|
||||||
if(i.ge.7) msg=trim(call_1)//' <'//trim(call_2)//'>'
|
|
||||||
endif
|
|
||||||
j0=len(trim(msg))+2
|
|
||||||
if(i.eq.2) msg(j0:j0+2)='RRR'
|
|
||||||
if(i.eq.3) msg(j0:j0+3)='RR73'
|
|
||||||
if(i.eq.4) msg(j0:j0+1)='73'
|
|
||||||
if(i.eq.5) then
|
|
||||||
if(std_2) then
|
|
||||||
msg='CQ '//trim(call_2)
|
|
||||||
if(call_1(3:3).eq.'_') msg=trim(call_1)//' '//trim(call_2)
|
|
||||||
if(grid4.ne.'RR73') msg=trim(msg)//' '//grid4
|
|
||||||
endif
|
|
||||||
if(.not.std_2) msg='CQ '//trim(call_2)
|
|
||||||
endif
|
|
||||||
if(i.eq.6 .and. std_2) msg(j0:j0+3)=grid4
|
|
||||||
if(i.ge.7) then
|
|
||||||
isnr = -50 + (i-7)/2
|
|
||||||
if(iand(i,1).eq.1) then
|
|
||||||
write(msg(j0:j0+2),'(i3.2)') isnr
|
|
||||||
if(msg(j0:j0).eq.' ') msg(j0:j0)='+'
|
|
||||||
else
|
|
||||||
write(msg(j0:j0+3),'("R",i3.2)') isnr
|
|
||||||
if(msg(j0+1:j0+1).eq.' ') msg(j0+1:j0+1)='+'
|
|
||||||
endif
|
|
||||||
endif
|
|
||||||
|
|
||||||
! Source-encode, then get itone()
|
|
||||||
i3=-1
|
|
||||||
n3=-1
|
|
||||||
call pack77(msg,i3,n3,c77)
|
|
||||||
call genft8(msg,i3,n3,msgsent,msgbits,itone)
|
|
||||||
! Generate complex cwave
|
|
||||||
f00=0.0
|
|
||||||
call gen_ft8wave(itone,NN,NSPS,bt,fs,f00,cwave,xjunk,1,NWAVE)
|
|
||||||
|
|
||||||
lagmax=-1
|
|
||||||
ccfmax=0.
|
|
||||||
nsum=32*2
|
|
||||||
lag0=200.0*(xdt0+0.5)
|
|
||||||
lag1=max(0,lag0-20)
|
|
||||||
lag2=min(nlags-1,lag0+20)
|
|
||||||
do lag=lag1,lag2
|
|
||||||
z=0.
|
|
||||||
s=0.
|
|
||||||
do i=0,NWAVE-1
|
|
||||||
z=z + cd(i+lag)*conjg(cwave(i))
|
|
||||||
if(mod(i,nsum).eq.nsum-1 .or. i.eq.NWAVE-1) then
|
|
||||||
s=s + abs(z)
|
|
||||||
z=0.
|
|
||||||
endif
|
|
||||||
enddo
|
|
||||||
ccf(lag)=s
|
|
||||||
if(ccf(lag).gt.ccfmax) then
|
|
||||||
ccfmax=ccf(lag)
|
|
||||||
lagmax=lag
|
|
||||||
endif
|
|
||||||
enddo ! lag
|
|
||||||
ccfmsg(imsg)=ccfmax
|
|
||||||
if(ccfmax.gt.ccfbest) then
|
|
||||||
ccfbest=ccfmax
|
|
||||||
lagbest=lagmax
|
|
||||||
msgbest=msg
|
|
||||||
imsgbest=imsg
|
|
||||||
itone_a7=itone
|
|
||||||
endif
|
|
||||||
enddo ! imsg
|
|
||||||
|
|
||||||
call pctile(ccfmsg,MAXMSG,50,base)
|
|
||||||
call pctile(ccfmsg,MAXMSG,84,sigma)
|
|
||||||
sigma=sigma-base
|
|
||||||
if(sigma.eq.0.0) sigma=1.0
|
|
||||||
ccfmsg=(ccfmsg-base)/sigma
|
|
||||||
xdt=lagbest/200.0 - 0.5
|
|
||||||
snr7=maxval(ccfmsg)
|
|
||||||
ccfmsg(imsgbest)=0.
|
|
||||||
snr7b=snr7/maxval(ccfmsg)
|
|
||||||
if(index(msgbest,'CQ0XYZ').ge.1) snr7=0.
|
|
||||||
xsnr=-99.0
|
|
||||||
if(snr7.gt.4.0) xsnr=db(snr7)-24.0
|
|
||||||
|
|
||||||
return
|
|
||||||
end subroutine ft8_dec7
|
|
||||||
|
|
||||||
|
|
||||||
end module ft8_a7
|
|
||||||
@ -254,10 +254,10 @@ contains
|
|||||||
xdt=dt0(i,jseq,0)
|
xdt=dt0(i,jseq,0)
|
||||||
f1=f0(i,jseq,0)
|
f1=f0(i,jseq,0)
|
||||||
msg37=' '
|
msg37=' '
|
||||||
call timer('ft8c ',0)
|
call timer('ft8_a7d ',0)
|
||||||
call ft8c(dd,newdat,call_1,call_2,grid4,xdt,f1,nharderrors, &
|
call ft8_a7d(dd,newdat,call_1,call_2,grid4,xdt,f1,nharderrors, &
|
||||||
dmin,msg37,xsnr)
|
dmin,msg37,xsnr)
|
||||||
call timer('ft8c ',1)
|
call timer('ft8_a7d ',1)
|
||||||
! write(51,3051) i,xdt,nint(f1),nharderrors,dmin,call_1,call_2,grid4
|
! write(51,3051) i,xdt,nint(f1),nharderrors,dmin,call_1,call_2,grid4
|
||||||
!3051 format(i3,f7.2,2i5,f7.1,1x,a12,a12,1x,a4)
|
!3051 format(i3,f7.2,2i5,f7.1,1x,a12,a12,1x,a4)
|
||||||
if(nharderrors.ge.0 .and. dmin.le.80.0) then
|
if(nharderrors.ge.0 .and. dmin.le.80.0) then
|
||||||
|
|||||||
Loading…
x
Reference in New Issue
Block a user