WSJT-X/lib/tweak1.f90
Steven Franke 5b5915a024 New short ping decoding algorithm for msk144 realtime decoder. Not yet enabled.
git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/wsjtx@7108 ab8295b8-cf94-4d9e-aec4-7959e3be5d79
2016-09-23 22:28:11 +00:00

24 lines
406 B
Fortran

subroutine tweak1(ca,jz,f0,cb)
! Shift frequency of analytic signal ca, with output to cb
complex ca(jz),cb(jz)
real*8 twopi
complex*16 w,wstep
complex w4
data twopi/0.d0/
save twopi
if(twopi.eq.0.d0) twopi=8.d0*atan(1.d0)
w=1.d0
dphi=twopi*f0/12000.d0
wstep=cmplx(cos(dphi),sin(dphi))
do i=1,jz
w=w*wstep
w4=w
cb(i)=w4*ca(i)
enddo
return
end subroutine tweak1