Minor cleanup: debugging a problem with FT4 decoding.

This commit is contained in:
Joe Taylor 2019-02-14 14:09:22 -05:00
parent c2cb55a68a
commit 073a22125a
4 changed files with 12 additions and 15 deletions

View File

@ -512,7 +512,6 @@ set (wsjt_FSRCS
lib/mskrtd.f90 lib/mskrtd.f90
lib/ft4/ft4sim.f90 lib/ft4/ft4sim.f90
lib/ft4/ft4sim_mult.f90 lib/ft4/ft4sim_mult.f90
lib/ft4/ft4d.f90
lib/ft4/ft4_decode.f90 lib/ft4/ft4_decode.f90
lib/ft4/ft4_downsample.f90 lib/ft4/ft4_downsample.f90
lib/77bit/my_hash.f90 lib/77bit/my_hash.f90

View File

@ -14,6 +14,7 @@ subroutine clockit(dname,k)
data first/.true./,eps/0.000001/,ntrace/0/ data first/.true./,eps/0.000001/,ntrace/0/
data level/0/,nmax/0/,space/' '/ data level/0/,nmax/0/,space/' '/
data limtrace/0/,lu/29/,ntimer/1/ data limtrace/0/,lu/29/,ntimer/1/
! data limtrace/1000000/,lu/29/,ntimer/1/
save save
if(ntimer.eq.0) return if(ntimer.eq.0) return
@ -56,8 +57,8 @@ subroutine clockit(dname,k)
endif endif
ntrace=ntrace+1 ntrace=ntrace+1
if(ntrace.lt.limtrace) write(lu,1020) ntrace,dname,k,level,nparent(n) if(ntrace.lt.limtrace) write(28,1020) ntrace,dname,k,level,nparent(n)
1020 format(i5,': ',a8,3i5) 1020 format(i8,': ',a8,3i5)
return return
! Write out the timer statistics ! Write out the timer statistics

View File

@ -16,7 +16,7 @@ subroutine ft4_decode(cdatetime0,tbuf,nfa,nfb,nQSOProgress,ncontest,nfqso, &
character*6 hhmmss character*6 hhmmss
complex cd2(0:NMAX/NDOWN-1) !Complex waveform complex cd2(0:NMAX/NDOWN-1) !Complex waveform
complex cds(0:NMAX/NDOWN-1) !Complex waveform ! complex cds(0:NMAX/NDOWN-1) !Complex waveform
complex cb(0:NMAX/NDOWN-1) complex cb(0:NMAX/NDOWN-1)
complex cd(0:NN*NSS-1) !Complex waveform complex cd(0:NN*NSS-1) !Complex waveform
complex ctwk(4*NSS),ctwk2(4*NSS) complex ctwk(4*NSS),ctwk2(4*NSS)
@ -66,8 +66,8 @@ subroutine ft4_decode(cdatetime0,tbuf,nfa,nfb,nQSOProgress,ncontest,nfqso, &
1,0,0,1,0,1,1,0,0,0,0,1,0,0,0,1,0,1,0,0,1,1,1,1,0,0,1,0,1, & 1,0,0,1,0,1,1,0,0,0,0,1,0,0,0,1,0,1,0,0,1,1,1,1,0,0,1,0,1, &
0,1,0,1,0,1,1,0,1,1,1,1,1,0,0,0,1,0,1/ 0,1,0,1,0,1,1,0,1,1,1,1,1,0,0,0,1,0,1/
save fs,dt,tt,txt,twopi,h,one,first,nrxx,linex,apbits,nappasses,naptypes, & save fs,dt,tt,txt,twopi,h,one,first,nrxx,linex,apbits,nappasses,naptypes, &
mycall0,hiscall0,ncontest0,msg0 mycall0,hiscall0,msg0
call clockit('ft4_deco',0) call clockit('ft4_deco',0)
hhmmss=cdatetime0(8:13) hhmmss=cdatetime0(8:13)
@ -333,7 +333,6 @@ subroutine ft4_decode(cdatetime0,tbuf,nfa,nfb,nQSOProgress,ncontest,nfqso, &
apmag=maxval(abs(llra))*1.1 apmag=maxval(abs(llra))*1.1
npasses=3+nappasses(nQSOProgress) npasses=3+nappasses(nQSOProgress)
if(ncontest.ge.5) npasses=3 ! Don't support Fox and Hound if(ncontest.ge.5) npasses=3 ! Don't support Fox and Hound
do ipass=1,npasses do ipass=1,npasses
if(ipass.eq.1) llr=llra if(ipass.eq.1) llr=llra
if(ipass.eq.2) llr=llrb if(ipass.eq.2) llr=llrb
@ -480,7 +479,7 @@ subroutine ft4_decode(cdatetime0,tbuf,nfa,nfb,nQSOProgress,ncontest,nfqso, &
nrx=-1 nrx=-1
if(index(message,'CQ ').eq.1) nrx=1 if(index(message,'CQ ').eq.1) nrx=1
if((index(message,trim(mycall)//' ').eq.1) .and. & if((index(message,trim(mycall)//' ').eq.1) .and. &
(index(message,' '//trim(hiscall)//' ').ge.4)) then (index(message,' '//trim(hiscall)//' ').ge.4)) then
if(index(message,' 559 ').gt.8) nrx=2 !### Not right ! if(index(message,' 559 ').gt.8) nrx=2 !### Not right !
if(index(message,' R 559 ').gt.8) nrx=3 !### Not right ! if(index(message,' R 559 ').gt.8) nrx=3 !### Not right !
if(index(message,' RR73 ').gt.8) nrx=4 if(index(message,' RR73 ').gt.8) nrx=4

View File

@ -5,7 +5,6 @@ program ft4d
character*8 arg character*8 arg
character*17 cdatetime character*17 cdatetime
character*512 data_dir character*512 data_dir
character*11 datetime
character*37 decodes(100) character*37 decodes(100)
character*16 fname character*16 fname
character*12 mycall character*12 mycall
@ -48,9 +47,10 @@ program ft4d
read(arg,*) nQSOProgress read(arg,*) nQSOProgress
iarg=iarg+2 iarg=iarg+2
endif endif
nfa=200 nfa=0
nfb=3000 nfb=4200
ncontest=0 ncontest=4
ndecodes=0
nfqso=1500 nfqso=1500
mycall="K9AN" mycall="K9AN"
hiscall="K1JT" hiscall="K1JT"
@ -60,10 +60,8 @@ program ft4d
j2=index(infile,'.wav') j2=index(infile,'.wav')
open(10,file=infile,status='old',access='stream') open(10,file=infile,status='old',access='stream')
read(10,end=999) ihdr,iwave read(10,end=999) ihdr,iwave
read(infile(j2-4:j2-1),*) nutc
datetime=infile(j2-11:j2-1)
cdatetime=' '//datetime
close(10) close(10)
cdatetime=infile(1:13)//'.000'
call ft4_decode(cdatetime,0.0,nfa,nfb,nQSOProgress,ncontest,nfqso,iwave, & call ft4_decode(cdatetime,0.0,nfa,nfb,nQSOProgress,ncontest,nfqso,iwave, &
ndecodes,mycall,hiscall,nrx,line,data_dir) ndecodes,mycall,hiscall,nrx,line,data_dir)