2015-11-17 20:28:12 -05:00
|
|
|
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
|
2016-09-23 18:28:11 -04:00
|
|
|
complex w4
|
2015-11-17 20:28:12 -05:00
|
|
|
data twopi/0.d0/
|
|
|
|
|
save twopi
|
|
|
|
|
|
|
|
|
|
if(twopi.eq.0.d0) twopi=8.d0*atan(1.d0)
|
|
|
|
|
w=1.d0
|
2016-01-27 13:51:55 -05:00
|
|
|
dphi=twopi*f0/12000.d0
|
2015-11-17 20:28:12 -05:00
|
|
|
wstep=cmplx(cos(dphi),sin(dphi))
|
|
|
|
|
do i=1,jz
|
|
|
|
|
w=w*wstep
|
2016-09-20 12:25:25 -04:00
|
|
|
w4=w
|
|
|
|
|
cb(i)=w4*ca(i)
|
2015-11-17 20:28:12 -05:00
|
|
|
enddo
|
|
|
|
|
|
|
|
|
|
return
|
|
|
|
|
end subroutine tweak1
|
