From 3ce9963f38406ee6352f456cb456b08ed200d8c6 Mon Sep 17 00:00:00 2001 From: Joe Taylor Date: Sat, 14 Jan 2006 14:59:30 +0000 Subject: [PATCH] git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/trunk@98 ab8295b8-cf94-4d9e-aec4-7959e3be5d79 --- four2.f | 6 +++--- pix2d.f90 | 2 ++ wsjt.py | 2 +- 3 files changed, 6 insertions(+), 4 deletions(-) diff --git a/four2.f b/four2.f index 54a7e1422..0a2e78c25 100755 --- a/four2.f +++ b/four2.f @@ -36,7 +36,7 @@ C the naive NTOT**2. Furthermore, less error is built up. C Written by Norman Brenner of MIT Lincoln Laboratory, January 1969. C See IEEE Audio Transactions (June 1967), Special issue on FFT. - DIMENSION DATA(1), N(1) + DIMENSION DATA(1048576), N(1) NTOT=1 DO 10 IDIM=1,NDIM 10 NTOT=NTOT*N(IDIM) @@ -78,7 +78,7 @@ C EXCHANGE DATA(J1,J4REV,J5) WITH DATA(J1,J4,J5) FOR ALL J1 FROM 1 C TO NPREV, ALL J4 FROM 1 TO N (WHICH MUST BE A POWER OF TWO), AND C ALL J5 FROM 1 TO NREM. J4REV-1 IS THE BIT REVERSAL OF J4-1. E.G. C SUPPOSE N = 32. THEN FOR J4-1 = 10011, J4REV-1 = 11001, ETC. - DIMENSION DATA(1) + DIMENSION DATA(1048576) IP0=2 IP1=IP0*NPREV IP4=IP1*N @@ -135,7 +135,7 @@ C IT IS NOT NECESSARY TO REWRITE THIS SUBROUTINE INTO COMPLEX C NOTATION SO LONG AS THE FORTRAN COMPILER USED STORES REAL AND C IMAGINARY PARTS IN ADJACENT STORAGE LOCATIONS. IT MUST ALSO C STORE ARRAYS WITH THE FIRST SUBSCRIPT INCREASING FASTEST. - DIMENSION DATA(1) + DIMENSION DATA(1048576) real*8 twopi,wstpr,wstpi,wr,wi,w2r,w2i,w3r,w3i,wtempr diff --git a/pix2d.f90 b/pix2d.f90 index f9fd05447..e723416b9 100644 --- a/pix2d.f90 +++ b/pix2d.f90 @@ -6,6 +6,8 @@ subroutine pix2d(d2,jz,mousebutton,s2,nchan,nz,b) integer*2 d2(jz) !Raw input data real s2(nchan,nz) !2-d spectrum integer*2 b(60000) !Pixels corresponding to 2-d spectrum + data nx0/0/ + save tbest=s2(2,1) s2(1,1)=s2(3,1) diff --git a/wsjt.py b/wsjt.py index 4fd987e84..89c940248 100644 --- a/wsjt.py +++ b/wsjt.py @@ -1,4 +1,4 @@ -#------------------------------------------------------------- WSJT +#-------------------------------------------------------------- WSJT from Tkinter import * from tkFileDialog import * import Pmw