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