From e87b82bcc302e4740688d72577725f86ff1fe10a Mon Sep 17 00:00:00 2001 From: Joe Taylor Date: Sat, 29 Jul 2006 16:10:39 +0000 Subject: [PATCH] Added missing routine symsync65.f. git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/trunk@213 ab8295b8-cf94-4d9e-aec4-7959e3be5d79 --- symsync65.f | 35 +++++++++++++++++++++++++++++++++++ wsjt.py | 2 +- 2 files changed, 36 insertions(+), 1 deletion(-) create mode 100644 symsync65.f diff --git a/symsync65.f b/symsync65.f new file mode 100644 index 000000000..dbf04fe1c --- /dev/null +++ b/symsync65.f @@ -0,0 +1,35 @@ + subroutine symsync65(c5,n5,k0,s,flip,pr,kmax,kpk,ccf,smax) + + complex c5(n5) + real s(n5),pr(126),ccf(-128:128) + complex z + + z=0. + do i=1,32 + z=z + c5(i) + enddo + s(1)=real(z)*real(z) + aimag(z)*aimag(z) + smax=s(1) + do i=33,n5 + z=z + c5(i) - c5(i-32) + s(i-31)=real(z)*real(z) + aimag(z)*aimag(z) + smax=max(s(i-31),smax) + enddo + iz=n5-31 + + smax=0. + do k=-kmax,kmax + sum=0. + do i=1,126 + j=32*(i-1)+k+k0 + if(j.ge.1 .and. j.le.iz) sum=sum + flip*pr(i)*s(j) + enddo + ccf(k)=sum + if(sum.gt.smax) then + smax=sum + kpk=k + endif + enddo + + return + end diff --git a/wsjt.py b/wsjt.py index f9bacd87a..1e25b0839 100644 --- a/wsjt.py +++ b/wsjt.py @@ -1,4 +1,4 @@ -#------------------------------------------------------------------- WSJT +#------------------------------------------------------------------ WSJT from Tkinter import * from tkFileDialog import * import Pmw