WSJT-X/lib/test_q65.f90

192 lines
5.1 KiB
Fortran
Raw Normal View History

program test_q65
character*82 cmd1,cmd2,line
character*22 msg
character*8 arg
2020-10-03 08:25:17 -04:00
character*1 csubmode
2021-01-01 15:28:41 -05:00
integer naptype(0:5)
logical decok
nargs=iargc()
if(nargs.ne.12) then
print*,'Usage: test_q65 "msg" A-D depth freq DT fDop f1 Stp TRp Q nfiles SNR'
print*,'Example: test_q65 "K1ABC W9XYZ EN37" A 3 1500 0.0 5.0 0.0 1 60 3 100 -20'
print*,'Use SNR = 0 to loop over all relevant SNRs'
print*,'Use MyCall=K1ABC, HisCall=W9XYZ, HisGrid="EN37" for AP decodes'
2020-11-17 15:31:33 -05:00
print*,'Option Q sets QSOprogress (0-5) for AP decoding.'
print*,'Add 16 to requested depth to enable message averaging.'
go to 999
endif
call getarg(1,msg)
2020-10-03 08:25:17 -04:00
call getarg(2,csubmode)
call getarg(3,arg)
2020-10-03 08:25:17 -04:00
read(arg,*) ndepth
call getarg(4,arg)
2020-10-03 08:25:17 -04:00
read(arg,*) nf0
call getarg(5,arg)
2020-10-03 08:25:17 -04:00
read(arg,*) dt
call getarg(6,arg)
2020-10-03 08:25:17 -04:00
read(arg,*) fDop
call getarg(7,arg)
read(arg,*) f1
call getarg(8,arg)
read(arg,*) nstp
2020-10-03 08:25:17 -04:00
call getarg(9,arg)
read(arg,*) ntrperiod
call getarg(10,arg)
read(arg,*) nQSOprogress
call getarg(11,arg)
read(arg,*) nfiles
call getarg(12,arg)
read(arg,*) snr
if(ntrperiod.eq.15) then
nsps=1800
i50=-23
else if(ntrperiod.eq.30) then
nsps=3600
i50=-26
else if(ntrperiod.eq.60) then
nsps=7200
i50=-29
else if(ntrperiod.eq.120) then
nsps=16000
i50=-31
else if(ntrperiod.eq.300) then
nsps=41472
i50=-35
else
stop 'Invalid TR period'
endif
2020-12-06 08:12:16 -05:00
i50=i50 + 8.0*log(1.0+fDop)/log(240.0)
ia=i50 + 7
2020-11-19 13:50:59 -05:00
ib=i50 - 10
if(snr.ne.0.0) then
ia=99
ib=99
endif
baud=12000.0/nsps
tsym=1.0/baud
! 1 2 3 4 5 6 7
! 123456789012345678901234567890123456789012345678901234567890123456789012345'
cmd1='q65sim "K1ABC W9XYZ EN37 " A 1500 5.0 0.0 0.0 1 60 100 -10.0 > junk0'
2021-01-31 09:06:20 -05:00
cmd2='jt9 -3 -p 15 -L 300 -H 3000 -d 3 -b A -Q 3 -f 1500 *.wav > junk'
write(cmd1(10:33),'(a)') '"'//msg//'"'
2020-10-03 08:25:17 -04:00
cmd1(35:35)=csubmode
write(cmd1(37:40),'(i4)') nf0
write(cmd1(41:45),'(f5.0)') fDop
write(cmd1(46:50),'(f5.2)') dt
write(cmd1(51:54),'(f4.0)') f1
write(cmd1(55:57),'(i3)') nstp
write(cmd1(58:61),'(i4)') ntrperiod
write(cmd1(62:66),'(i5)') nfiles
write(cmd2(11:13),'(i3)') ntrperiod
2021-01-31 09:06:20 -05:00
write(cmd2(33:35),'(i3)') ndepth
write(cmd2(45:45),'(i1)') nQSOprogress
write(cmd2(50:53),'(i4)') nf0
cmd2(40:40)=csubmode
call system('rm -f *.wav')
write(*,1008) ntrperiod,csubmode,ndepth,fDop,f1,nstp
1008 format('Mode:',i4,a1,' Depth:',i3,' fDop:',f6.1,' Drift:',f6.1, &
' Steps:',i3)
2021-01-01 15:28:41 -05:00
write(*,1010) (j,j=0,5)
write(12,1010) (j,j=0,5)
1010 format(' SNR Sync Avg Dec Bad',6i4,' tdec avg rms'/64('-'))
dterr=tsym/4.0
nferr=max(1,nint(0.5*baud),nint(fdop/3.0))
2020-12-05 11:09:15 -05:00
ndec1z=nfiles
do nsnr=ia,ib,-1
snr1=nsnr
if(ia.eq.99) snr1=snr
nsync=0
ndec1=0
nfalse=0
2020-11-11 15:20:35 -05:00
naptype=0
ndecn=0
write(cmd1(70:74),'(f5.1)') snr1
call system(cmd1)
call sec0(0,tdec)
call system(cmd2)
call sec0(1,tdec)
open(10,file='junk',status='unknown')
n=0
snrsum=0.
snrsq=0.
nsum=0
do iline=1,9999
read(10,'(a71)',end=10) line
if(index(line,'<Decode').eq.1) cycle
read(line,*) itime,xsnr,xdt,nf
decok=index(line,trim(msg)).gt.0
if(decok) then
snrsum=snrsum + xsnr
snrsq=snrsq + xsnr*xsnr
nsum=nsum+1
endif
if((abs(xdt-dt).le.dterr .and. abs(nf-nf0).le.nferr) .or. decok) then
nsync=nsync+1
endif
2020-12-01 11:24:21 -05:00
idec=-1
iavg=0
i0=23
if(ntrperiod.le.30) i0=25
if(line(i0:i0).ne.' ') then
i1=index(line,'q')
idec=-1
2021-01-31 09:06:20 -05:00
read(line(i1+1:i1+1),*) idec
if(line(i1+2:i1+2).eq.'*') then
iavg=10
else
read(line(i1+2:i1+2),*,end=100) iavg
endif
endif
2021-01-31 09:06:20 -05:00
100 if(idec.lt.0) cycle
if(decok) then
ndecn=ndecn + 1
if(iavg.le.1) ndec1=ndec1 + 1
2020-12-01 11:24:21 -05:00
naptype(idec)=naptype(idec) + 1
else
nfalse=nfalse + 1
print*,'False: ',line
endif
enddo
10 close(10)
snr_avg=0.
snr_rms=0.
if(nsum.ge.1) then
snr_avg=snrsum/nsum
snr_rms=sqrt(snrsq/nsum - snr_avg**2)
endif
write(*,1100) snr1,nsync,ndecn,ndec1,nfalse,naptype,tdec/nfiles, &
snr_avg,snr_rms
write(12,1100) snr1,nsync,ndecn,ndec1,nfalse,naptype,tdec/nfiles, &
snr_avg,snr_rms
1100 format(f5.1,4i4,i5,5i4,f6.2,f6.1,f5.1)
2020-12-05 11:09:15 -05:00
if(ndec1.lt.nfiles/2 .and. ndec1z.ge.nfiles/2) then
snr_thresh=snr1 + float(nfiles/2 - ndec1)/(ndec1z-ndec1)
open(13,file='snr_thresh.out',status='unknown',position='append')
write(13,1200) ntrperiod,csubmode,ndepth,nQSOprogress,nfiles, &
fdop,f1,nstp,nfalse,snr_thresh,trim(msg)
1200 format(i3,a1,2i3,i5,2f7.1,2i3,f7.1,2x,a)
close(13)
endif
flush(6)
flush(12)
if(ndec1.eq.0 .and. ndecn.eq.0) exit !Bail out if no decodes at this SNR
2020-12-05 11:09:15 -05:00
ndec1z=ndec1
enddo ! nsnr
999 end program test_q65
include 'sec0.f90'