WSJT-X/lib/ft8/genft8refsig.f90

subroutine genft8refsig(itone,cref,f0)
  complex cref(79*1920)
  integer itone(79)
!  real*8 twopi,phi,dphi,dt,xnsps
  real twopi,phi,dphi,dt,xnsps
  data twopi/0.d0/
  save twopi
  if( twopi .lt. 0.1 ) twopi=8.d0*atan(1.d0)

  xnsps=1920.d0
  dt=1.d0/12000.d0
  phi=0.d0
  k=1
  do i=1,79
    dphi=twopi*(f0*dt+itone(i)/xnsps)
    do is=1,1920
      cref(k)=cmplx(cos(phi),sin(phi))
      phi=mod(phi+dphi,twopi)
      k=k+1
    enddo
  enddo
  return
end subroutine genft8refsig
Add subtraction routines - not fully tested or tweaked yet. git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/wsjtx@7787 ab8295b8-cf94-4d9e-aec4-7959e3be5d79 2017-07-06 10:44:17 -04:00			`subroutine genft8refsig(itone,cref,f0)`
			`complex cref(79*1920)`
			`integer itone(79)`
Fix an off-by-one problem. 2018-06-30 15:56:19 -04:00			`! real*8 twopi,phi,dphi,dt,xnsps`
			`real twopi,phi,dphi,dt,xnsps`
Add subtraction routines - not fully tested or tweaked yet. git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/wsjtx@7787 ab8295b8-cf94-4d9e-aec4-7959e3be5d79 2017-07-06 10:44:17 -04:00			`data twopi/0.d0/`
			`save twopi`
			`if( twopi .lt. 0.1 ) twopi=8.d0*atan(1.d0)`

			`xnsps=1920.d0`
			`dt=1.d0/12000.d0`
			`phi=0.d0`
			`k=1`
			`do i=1,79`
			`dphi=twopi(f0dt+itone(i)/xnsps)`
			`do is=1,1920`
			`cref(k)=cmplx(cos(phi),sin(phi))`
			`phi=mod(phi+dphi,twopi)`
			`k=k+1`
			`enddo`
			`enddo`
			`return`
			`end subroutine genft8refsig`