157 lines
4.7 KiB
Fortran
157 lines
4.7 KiB
Fortran
|
module js8d_decode
|
||
|
|
|
||
|
|
type :: js8d_decoder
|
||
|
|
procedure(js8d_decode_callback), pointer :: callback
|
||
|
|
contains
|
||
|
|
procedure :: decode
|
||
|
|
end type js8d_decoder
|
||
|
|
|
||
|
|
abstract interface
|
||
|
|
subroutine js8d_decode_callback (this,sync,snr,dt,freq,decoded,nap,qual)
|
||
|
|
import js8d_decoder
|
||
|
|
implicit none
|
||
|
|
class(js8d_decoder), intent(inout) :: this
|
||
|
|
real, intent(in) :: sync
|
||
|
|
integer, intent(in) :: snr
|
||
|
|
real, intent(in) :: dt
|
||
|
|
real, intent(in) :: freq
|
||
|
|
character(len=37), intent(in) :: decoded
|
||
|
|
integer, intent(in) :: nap
|
||
|
|
real, intent(in) :: qual
|
||
|
|
end subroutine js8d_decode_callback
|
||
|
|
end interface
|
||
|
|
|
||
|
|
contains
|
||
|
|
|
||
|
|
subroutine decode(this,callback,iwave,nQSOProgress,nfqso,nftx,newdat, &
|
||
|
|
nutc,nfa,nfb,nexp_decode,ndepth,nagain,lft8apon,lapcqonly,napwid, &
|
||
|
|
mycall12,mygrid6,hiscall12,hisgrid6)
|
||
|
|
! use wavhdr
|
||
|
|
use timer_module, only: timer
|
||
|
|
! type(hdr) h
|
||
|
|
use js8d_module
|
||
|
|
|
||
|
|
class(js8d_decoder), intent(inout) :: this
|
||
|
|
procedure(js8d_decode_callback) :: callback
|
||
|
|
real s(NH1,NHSYM)
|
||
|
|
real sbase(NH1)
|
||
|
|
real candidate(3,200)
|
||
|
|
real dd(NMAX)
|
||
|
|
logical, intent(in) :: lft8apon,lapcqonly,nagain
|
||
|
|
logical newdat,lsubtract,ldupe,bcontest
|
||
|
|
character*12 mycall12, hiscall12
|
||
|
|
character*6 mygrid6,hisgrid6
|
||
|
|
integer*2 iwave(NMAX)
|
||
|
|
integer apsym(KK)
|
||
|
|
character datetime*13,message*22,msg37*37
|
||
|
|
character*22 allmessages(100)
|
||
|
|
integer allsnrs(100)
|
||
|
|
save s,dd
|
||
|
|
|
||
|
|
bcontest=iand(nexp_decode,128).ne.0
|
||
|
|
this%callback => callback
|
||
|
|
write(datetime,1001) nutc !### TEMPORARY ###
|
||
|
|
1001 format("000000_",i6.6)
|
||
|
|
|
||
|
|
call ft8apset(mycall12,mygrid6,hiscall12,hisgrid6,bcontest,apsym,iaptype)
|
||
|
|
|
||
|
|
dd=iwave
|
||
|
|
ndecodes=0
|
||
|
|
allmessages=' '
|
||
|
|
allsnrs=0
|
||
|
|
ifa=nfa
|
||
|
|
ifb=nfb
|
||
|
|
if(nagain) then
|
||
|
|
ifa=nfqso-10
|
||
|
|
ifb=nfqso+10
|
||
|
|
endif
|
||
|
|
|
||
|
|
! For now:
|
||
|
|
! ndepth=1: no subtraction, 1 pass, belief propagation only
|
||
|
|
! ndepth=2: subtraction, 3 passes, belief propagation only
|
||
|
|
! ndepth=3: subtraction, 3 passes, bp+osd
|
||
|
|
if(ndepth.eq.1) npass=1
|
||
|
|
if(ndepth.ge.2) npass=3
|
||
|
|
do ipass=1,npass
|
||
|
|
newdat=.true. ! Is this a problem? I hijacked newdat.
|
||
|
|
syncmin=ASYNCMIN
|
||
|
|
if(ipass.eq.1) then
|
||
|
|
lsubtract=.true.
|
||
|
|
if(ndepth.eq.1) lsubtract=.false.
|
||
|
|
elseif(ipass.eq.2) then
|
||
|
|
n2=ndecodes
|
||
|
|
if(ndecodes.eq.0) cycle
|
||
|
|
lsubtract=.true.
|
||
|
|
elseif(ipass.eq.3) then
|
||
|
|
if((ndecodes-n2).eq.0) cycle
|
||
|
|
lsubtract=.false.
|
||
|
|
endif
|
||
|
|
|
||
|
|
call timer('syncjs8 ',0)
|
||
|
|
call syncjs8(dd,ifa,ifb,syncmin,nfqso,s,candidate,ncand,sbase)
|
||
|
|
call timer('syncjs8 ',1)
|
||
|
|
|
||
|
|
if(NWRITELOG.eq.1) then
|
||
|
|
open(99, file="./js8.log", status="old", position="append", action="write")
|
||
|
|
write(99,*) ncand, "candidates"
|
||
|
|
close(99)
|
||
|
|
endif
|
||
|
|
|
||
|
|
do icand=1,ncand
|
||
|
|
sync=candidate(3,icand)
|
||
|
|
f1=candidate(1,icand)
|
||
|
|
xdt=candidate(2,icand)
|
||
|
|
xbase=10.0**(0.1*(sbase(nint(f1/(12000.0/NFFT1)))-40.0)) ! 3.125Hz
|
||
|
|
nsnr0=min(99,nint(10.0*log10(sync) - 25.5)) !### empirical ###
|
||
|
|
|
||
|
|
if(NWRITELOG.eq.1) then
|
||
|
|
open(99, file="./js8.log", status="old", position="append", action="write")
|
||
|
|
write(99,*) 'icand', icand, f1, sync, xdt, xbase
|
||
|
|
close(99)
|
||
|
|
endif
|
||
|
|
|
||
|
|
call timer('js8dec ',0)
|
||
|
|
call js8dec(dd,newdat,nQSOProgress,nfqso,nftx,ndepth,lft8apon, &
|
||
|
|
lapcqonly,napwid,lsubtract,nagain,iaptype,mycall12,mygrid6, &
|
||
|
|
hiscall12,bcontest,sync,f1,xdt,xbase,apsym,nharderrors,dmin, &
|
||
|
|
nbadcrc,iappass,iera,msg37,xsnr)
|
||
|
|
message=msg37(1:22) !###
|
||
|
|
nsnr=nint(xsnr)
|
||
|
|
xdt=xdt-ASTART
|
||
|
|
hd=nharderrors+dmin
|
||
|
|
|
||
|
|
if(NWRITELOG.eq.1) then
|
||
|
|
open(99, file="./js8.log", status="old", position="append", action="write")
|
||
|
|
write(99,*) 'icand', icand, hd, nbadcrc
|
||
|
|
close(99)
|
||
|
|
endif
|
||
|
|
|
||
|
|
call timer('js8dec ',1)
|
||
|
|
if(nbadcrc.eq.0) then
|
||
|
|
ldupe=.false.
|
||
|
|
do id=1,ndecodes
|
||
|
|
if(message.eq.allmessages(id).and.nsnr.le.allsnrs(id)) ldupe=.true.
|
||
|
|
enddo
|
||
|
|
if(.not.ldupe) then
|
||
|
|
ndecodes=ndecodes+1
|
||
|
|
allmessages(ndecodes)=message
|
||
|
|
allsnrs(ndecodes)=nsnr
|
||
|
|
endif
|
||
|
|
if(.not.ldupe .and. associated(this%callback)) then
|
||
|
|
qual=1.0-(nharderrors+dmin)/60.0 ! scale qual to [0.0,1.0]
|
||
|
|
call this%callback(sync,nsnr,xdt,f1,msg37,iaptype,qual)
|
||
|
|
endif
|
||
|
|
endif
|
||
|
|
|
||
|
|
if(NWRITELOG.eq.1) then
|
||
|
|
open(99, file="./js8.log", status="old", position="append", action="write")
|
||
|
|
write(99,*) '---'
|
||
|
|
close(99)
|
||
|
|
endif
|
||
|
|
enddo
|
||
|
|
enddo
|
||
|
|
return
|
||
|
|
end subroutine decode
|
||
|
|
|
||
|
|
end module js8d_decode
|