w4kek/WSJT-X
1
0
Fork 0
mirror of https://github.com/saitohirga/WSJT-X.git synced 2025-05-29 21:02:26 -04:00
Code Issues Packages Projects Releases Wiki Activity

Another minor tweak to q65sim.

Browse Source
This commit is contained in:
Joe Taylor 2023-02-15 15:58:20 -05:00
parent b5f9e0ff03
commit f9d2ca736d
1 changed files with 4 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
lib/qra/q65/q65sim.f90
View File

@ -18,7 +18,7 @@ program q65sim
complex cdat(NMAX) !Generated complex waveform complex cdat(NMAX) !Generated complex waveform
complex cspread(0:NMAX-1) !Complex amplitude for Rayleigh fading complex cspread(0:NMAX-1) !Complex amplitude for Rayleigh fading
complex z complex z
real*8 f0,dt,twopi,phi,dphi,baud,fsample,freq real*8 f00,f0,dt,twopi,phi,dphi,baud,fsample,freq
character fname*17,csubmode*1,arg*12,c2*2 character fname*17,csubmode*1,arg*12,c2*2
character*37 msg,msgsent,imsg(10) character*37 msg,msgsent,imsg(10)
@ -40,7 +40,7 @@ program q65sim
call getarg(2,csubmode) call getarg(2,csubmode)
mode65=2**(ichar(csubmode)-ichar('A')) mode65=2**(ichar(csubmode)-ichar('A'))
call getarg(3,arg) call getarg(3,arg)
read(arg,*) f0 read(arg,*) f00
call getarg(4,arg) call getarg(4,arg)
read(arg,*) fspread read(arg,*) fspread
call getarg(5,arg) call getarg(5,arg)
@ -147,7 +147,7 @@ program q65sim
bandwidth_ratio=2500.0/6000.0 bandwidth_ratio=2500.0/6000.0
sig=sqrt(2*bandwidth_ratio)*10.0**(0.05*snrdb) sig=sqrt(2*bandwidth_ratio)*10.0**(0.05*snrdb)
if(snrdb.gt.90.0) sig=1.0 if(snrdb.gt.90.0) sig=1.0
write(*,1020) ifile,ntrperiod,f0,csubmode,snrdb,fspread,xdt,f1,nstp,trim(msgsent) write(*,1020) ifile,ntrperiod,f00,csubmode,snrdb,fspread,xdt,f1,nstp,trim(msgsent)
1020 format(i4,i6,f7.1,2x,a1,2x,f5.1,1x,f6.2,2f6.1,i4,2x,a) 1020 format(i4,i6,f7.1,2x,a1,2x,f5.1,1x,f6.2,2f6.1,i4,2x,a)
n=65.0*baud*mode65/100.0 + 0.9999 n=65.0*baud*mode65/100.0 + 0.9999
@ -155,6 +155,7 @@ program q65sim
nf1=1500 - nfstep*(nsig-1)/2 nf1=1500 - nfstep*(nsig-1)/2
do n=1,nsig do n=1,nsig
if(nsig.ge.2) then if(nsig.ge.2) then
f0=f00
f0=nf1 + (n-1)*nfstep f0=nf1 + (n-1)*nfstep
itone=ntone(:,n) itone=ntone(:,n)
endif endif