Correct a flaw in qra64sim.

This commit is contained in:
Joe Taylor 2020-07-31 10:58:30 -04:00
parent ee24f03d09
commit 9146ce9667

View File

@ -107,7 +107,7 @@ program qra64sim
twopi=8*atan(1.0) twopi=8*atan(1.0)
cspread(0)=1.0 cspread(0)=1.0
cspread(NH)=0. cspread(NH)=0.
b=6.0 !Lorenzian 3/28 onward b=6.0 !Use truncated Lorenzian shape for fspread
do i=1,NH do i=1,NH
f=i*df f=i*df
x=b*f/fspread x=b*f/fspread
@ -129,13 +129,13 @@ program qra64sim
cspread(NFFT-i)=z cspread(NFFT-i)=z
enddo enddo
do i=0,NFFT-1 ! do i=0,NFFT-1
f=i*df ! f=i*df
if(i.gt.NH) f=(i-nfft)*df ! if(i.gt.NH) f=(i-nfft)*df
s=real(cspread(i))**2 + aimag(cspread(i))**2 ! s=real(cspread(i))**2 + aimag(cspread(i))**2
! write(13,3000) i,f,s,cspread(i) ! write(13,3000) i,f,s,cspread(i)
!3000 format(i5,f10.3,3f12.6) !3000 format(i5,f10.3,3f12.6)
enddo ! enddo
! s=real(cspread(0))**2 + aimag(cspread(0))**2 ! s=real(cspread(0))**2 + aimag(cspread(0))**2
! write(13,3000) 1024,0.0,s,cspread(0) ! write(13,3000) 1024,0.0,s,cspread(0)