m65 noe mostly OK for single-pol mode. Needs more work for xpol.

This commit is contained in:
Joe Taylor 2021-04-28 10:13:42 -04:00
parent 57f48187e7
commit c7af9945d3
3 changed files with 34 additions and 45 deletions

View File

@ -80,16 +80,17 @@ contains
integer, intent(in) :: idec integer, intent(in) :: idec
integer, intent(in) :: nused integer, intent(in) :: nused
integer, intent(in) :: ntrperiod integer, intent(in) :: ntrperiod
character*3 cq
if(nutc+snr1+nsnr+dt+freq+idec+nused+ntrperiod.eq.-999) stop if(nutc+snr1+nsnr+dt+freq+idec+nused+ntrperiod.eq.-999) stop
if(decoded.eq.'-999') stop if(decoded.eq.'-999') stop
cq='q ' cq0='q '
write(cq(2:2),'(i1)') idec write(cq0(2:2),'(i1)') idec
if(nused.ge.2) write(cq(3:3),'(i1)') nused if(nused.ge.2) write(cq0(3:3),'(i1)') nused
write(24,1010) nutc,nsnr,dt,nint(freq),decoded,cq nsnr0=nsnr
1010 format(i4.4,i4,f5.1,i5,4x,a37,1x,a3) xdt0=dt
nfreq0=nint(freq)
msg0=decoded
select type(this) select type(this)
type is (counting_q65_decoder) type is (counting_q65_decoder)

View File

@ -1,5 +1,9 @@
module q65_decode module q65_decode
integer nsnr0,nfreq0
real xdt0
character msg0*37,cq0*3
type :: q65_decoder type :: q65_decoder
procedure(q65_decode_callback), pointer :: callback procedure(q65_decode_callback), pointer :: callback
contains contains

View File

@ -1,7 +1,7 @@
subroutine q65b(nutc,fcenter,nfcal,nfsample,ikhz,mousedf,ntol,xpol, & subroutine q65b(nutc,fcenter,nfcal,nfsample,ikhz,mousedf,ntol,xpol, &
mycall0,hiscall0,hisgrid,mode_q65) mycall0,hiscall0,hisgrid,mode_q65)
use wavhdr use q65_decode
parameter (MAXFFT1=5376000) !56*96000 parameter (MAXFFT1=5376000) !56*96000
parameter (MAXFFT2=336000) !56*6000 (downsampled by 1/16) parameter (MAXFFT2=336000) !56*6000 (downsampled by 1/16)
parameter (NMAX=60*12000) parameter (NMAX=60*12000)
@ -14,9 +14,9 @@ subroutine q65b(nutc,fcenter,nfcal,nfsample,ikhz,mousedf,ntol,xpol, &
character*12 mycall,hiscall character*12 mycall,hiscall
character*6 hisgrid character*6 hisgrid
character*4 grid4 character*4 grid4
character*62 line character*80 line
character*80 line2 character*37 msg1
character*40 msg40 character*3 cq1
character*80 wsjtx_dir character*80 wsjtx_dir
common/cacb/ca,cb common/cacb/ca,cb
save save
@ -24,7 +24,6 @@ subroutine q65b(nutc,fcenter,nfcal,nfsample,ikhz,mousedf,ntol,xpol, &
open(9,file='wsjtx_dir.txt',status='old') open(9,file='wsjtx_dir.txt',status='old')
read(9,'(a)') wsjtx_dir read(9,'(a)') wsjtx_dir
close(9) close(9)
open(24,file='q65_decodes.txt',status='unknown')
mycall='K1JT' mycall='K1JT'
hiscall='IV3NWV' hiscall='IV3NWV'
@ -67,6 +66,7 @@ subroutine q65b(nutc,fcenter,nfcal,nfsample,ikhz,mousedf,ntol,xpol, &
! 96000 5376000 0.017857143 336000 6000.000 ! 96000 5376000 0.017857143 336000 6000.000
! 95238 5120000 0.018601172 322560 5999.994 ! 95238 5120000 0.018601172 322560 5999.994
nsnr1=-99
npol=1 npol=1
if(xpol) npol=4 if(xpol) npol=4
do ipol=1,npol do ipol=1,npol
@ -90,38 +90,23 @@ subroutine q65b(nutc,fcenter,nfcal,nfsample,ikhz,mousedf,ntol,xpol, &
newdat=1 newdat=1
nagain=0 nagain=0
call map65_mmdec(nutc,iwave,nsubmode,nfa,nfb,nfqso,ntol,newdat,nagain, & call map65_mmdec(nutc,iwave,nsubmode,nfa,nfb,nfqso,ntol,newdat,nagain, &
mycall,hiscall,hisgrid) mycall,hiscall,hisgrid)
enddo ! print*,'bbb',ipol,nsnr0,xdt0,nfreq0,msg0,cq0
if(nsnr0.gt.nsnr1) then
! 1 2 3 4 5 6 nsnr1=nsnr0
! 1234567890123456789012345678901234567890123456789012345678901234567 xdt1=xdt0
! 0001 -22 2.9 1081 : EA2AGZ IK4WLV -16 q0 nfreq1=nfreq0
! 110 101 2 1814 2.9 -11 # QRZ HB9Q JN47 1 0 30 H msg1=msg0
nsnr0=-99 cq1=cq0
line2=' ' ipol1=45*ipol
rewind 24
do i=1,8
read(24,1002,end=100) line
1002 format(a62)
if(line(1:4).eq.'<Dec') cycle
read(line,1010) nsnr,xdt,nfreq,msg40
1010 format(4x,i4,f5.1,i5,4x,a40)
nfreq=nfreq+mousedf
if(nsnr.gt.nsnr0) then
ipol0=(i/2)*45
nsnr0=nsnr
write(line2,1020) ikhz,nfreq-1000,ipol0,nutc,xdt,nsnr0,msg40(1:27),msg40(39:40)
1020 format('!',i3.3,i5,i4,i6.4,f5.1,i5,' : ',a27,a2)
nsnr0=nsnr
freq0=144.0 + 0.001*ikhz
nfreq0=nfreq-1000
xdt0=xdt
endif endif
enddo enddo
100 if(nsnr0.gt.-40) then nfreq=nfreq1+mousedf-1000
write(*,1100) trim(line2) write(line,1020) ikhz,nfreq,ipol1,nutc,xdt1,nsnr1,msg1(1:27),cq1
1020 format('!',i3.3,i5,i4,i6.4,f5.1,i5,' : ',a27,a3)
freq0=144.0 + 0.001*ikhz
write(*,1100) trim(line)
1100 format(a) 1100 format(a)
! Should write to lu 26 here, for Messages and Band Map windows ? ! Should write to lu 26 here, for Messages and Band Map windows ?
@ -129,10 +114,9 @@ subroutine q65b(nutc,fcenter,nfcal,nfsample,ikhz,mousedf,ntol,xpol, &
! nsnr0,nutc,msg40(1:22),' ',char(ichar('A') + mode_q65-1) ! nsnr0,nutc,msg40(1:22),' ',char(ichar('A') + mode_q65-1)
!1014 format(f8.3,i5,3i3,f5.1,i4,i3,i4,i5.4,4x,a22,2x,a1,3x,':',a1) !1014 format(f8.3,i5,3i3,f5.1,i4,i3,i4,i5.4,4x,a22,2x,a1,3x,':',a1)
write(21,1110) freq0,nfreq0,xdt0,ipol0,nsnr0,nutc,msg40(1:28),msg40(39:40) ! Write to file map65_rx.log:
1110 format(f8.3,i5,f5.1,2i4,i5.4,2x,a28,': A',2x,a2) write(21,1110) freq0,nfreq,xdt1,ipol1,nsnr1,nutc,msg1(1:28),cq1
endif 1110 format(f8.3,i5,f5.1,2i4,i5.4,2x,a28,': A',2x,a3)
close(24,status='delete')
900 return 900 return
end subroutine q65b end subroutine q65b