mirror of
				https://github.com/saitohirga/WSJT-X.git
				synced 2025-11-04 05:50:31 -05:00 
			
		
		
		
	
		
			
				
	
	
		
			91 lines
		
	
	
		
			1.9 KiB
		
	
	
	
		
			Fortran
		
	
	
	
	
	
			
		
		
	
	
			91 lines
		
	
	
		
			1.9 KiB
		
	
	
	
		
			Fortran
		
	
	
	
	
	
subroutine fast_decode(id2,narg,trperiod,line,mycall_12,   &
 | 
						|
     hiscall_12)
 | 
						|
 | 
						|
  parameter (NMAX=30*12000)
 | 
						|
  integer*2 id2(NMAX)
 | 
						|
  integer*2 id2a(NMAX)
 | 
						|
  integer*2 id2b(NMAX)
 | 
						|
  integer narg(0:14)
 | 
						|
  double precision trperiod
 | 
						|
  real dat(30*12000)
 | 
						|
  complex cdat(262145),cdat2(262145)
 | 
						|
  logical pick,first
 | 
						|
  character*6 cfile6
 | 
						|
  character*80 line(100)
 | 
						|
  character*12 mycall_12,hiscall_12
 | 
						|
  character*6 mycall,hiscall
 | 
						|
  data first/.true./,nutca/0/,nutcb/0/
 | 
						|
  save npts,cdat,cdat2,id2a,id2b,nutca,nutcb
 | 
						|
 | 
						|
  if(first) then
 | 
						|
     id2a=0
 | 
						|
     id2b=0
 | 
						|
     first=.false.
 | 
						|
  endif
 | 
						|
 | 
						|
  mycall=mycall_12(1:6)
 | 
						|
  hiscall=hiscall_12(1:6)
 | 
						|
  nutc=narg(0)
 | 
						|
  ndat0=narg(1)
 | 
						|
  nsubmode=narg(2)
 | 
						|
  newdat=narg(3)
 | 
						|
  minsync=narg(4)
 | 
						|
  npick=narg(5)
 | 
						|
  t0=0.001*narg(6)
 | 
						|
  t1=0.001*narg(7)
 | 
						|
  tmid=0.5*(t0+t1)
 | 
						|
  maxlines=narg(8)
 | 
						|
  nmode=narg(9)
 | 
						|
  nrxfreq=narg(10)
 | 
						|
  ntol=narg(11)
 | 
						|
  nhashcalls=narg(12)
 | 
						|
 | 
						|
  line(1:100)(1:1)=char(0)
 | 
						|
  if(t0.gt.trperiod) go to 900
 | 
						|
  if(t0.gt.t1) go to 900
 | 
						|
 | 
						|
  if(nmode.eq.102) then
 | 
						|
     call fast9(id2,narg,line)
 | 
						|
     go to 900
 | 
						|
  endif
 | 
						|
 | 
						|
  if(newdat.eq.1) then
 | 
						|
     cdat2=cdat
 | 
						|
     ndat=ndat0
 | 
						|
     call wav11(id2,ndat,dat)
 | 
						|
     nzz=11025*int(trperiod)    !beware if fractional T/R period ever used here
 | 
						|
     if(ndat.lt.nzz) dat(ndat+1:nzz)=0.0
 | 
						|
     ndat=min(ndat,30*11025)
 | 
						|
     call ana932(dat,ndat,cdat,npts)          !Make downsampled analytic signal
 | 
						|
  endif
 | 
						|
 | 
						|
! Now cdat() is the downsampled analytic signal.  
 | 
						|
! New sample rate = fsample = BW = 11025 * (9/32) = 3100.78125 Hz
 | 
						|
! NB: npts, nsps, etc., are all reduced by 9/32
 | 
						|
 | 
						|
  write(cfile6,'(i6.6)') nutc
 | 
						|
  nfreeze=1
 | 
						|
  mousedf=0
 | 
						|
  mousebutton=0
 | 
						|
  mode4=1
 | 
						|
  if(nsubmode.eq.1) mode4=2
 | 
						|
  nafc=0
 | 
						|
  ndebug=0
 | 
						|
  t2=0.
 | 
						|
  ia=1
 | 
						|
  ib=npts
 | 
						|
  pick=.false.
 | 
						|
 | 
						|
  if(npick.gt.0) then
 | 
						|
     pick=.true.
 | 
						|
     dt=1.0/11025.0 * (32.0/9.0)
 | 
						|
     ia=t0/dt + 1.
 | 
						|
     ib=t1/dt + 1.
 | 
						|
     t2=t0
 | 
						|
  endif
 | 
						|
  jz=ib-ia+1
 | 
						|
  line(1:100)(1:1)=char(0)
 | 
						|
 | 
						|
900 return
 | 
						|
end subroutine fast_decode
 |