Try to get c_bigfft sorted out.

This commit is contained in:
Steven Franke 2020-07-11 16:39:31 -05:00
parent 5b804d90a3
commit 4cb902a4cc

View File

@ -211,7 +211,7 @@ contains
nfft1=nfft2*ndown nfft1=nfft2*ndown
nh1=nfft1/2 nh1=nfft1/2
allocate( c_bigfft(0:nfft1/2+1) ) allocate( c_bigfft(0:nfft1/2) )
allocate( c2(0:nfft2-1) ) allocate( c2(0:nfft2-1) )
allocate( cframe(0:160*nss-1) ) allocate( cframe(0:160*nss-1) )
@ -232,10 +232,9 @@ contains
! The big fft is done once and is used for calculating the smoothed spectrum ! The big fft is done once and is used for calculating the smoothed spectrum
! and also for downconverting/downsampling each candidate. ! and also for downconverting/downsampling each candidate.
do i=1,nfft1/2 do i=0,nfft1/2
c_bigfft(i)=cmplx(float(iwave(2*i-1)),float(iwave(2*i))) c_bigfft(i)=cmplx(float(iwave(2*i+1)),float(iwave(2*i+2)))
enddo enddo
c_bigfft(nfft1/2+1)=0.
call four2a(c_bigfft,nfft1,1,-1,0) call four2a(c_bigfft,nfft1,1,-1,0)
if(iwspr.eq.0) then if(iwspr.eq.0) then
itype1=1 itype1=1
@ -727,7 +726,7 @@ contains
subroutine get_candidates_fst240(c_bigfft,nfft1,nsps,hmod,fs,fa,fb, & subroutine get_candidates_fst240(c_bigfft,nfft1,nsps,hmod,fs,fa,fb, &
minsync,ncand,candidates,base) minsync,ncand,candidates,base)
complex c_bigfft(0:nfft1/2+1) !Full length FFT of raw data complex c_bigfft(0:nfft1/2) !Full length FFT of raw data
integer hmod !Modulation index (submode) integer hmod !Modulation index (submode)
integer im(1) !For maxloc integer im(1) !For maxloc
real candidates(100,4) !Candidate list real candidates(100,4) !Candidate list