2012-11-26 16:06:41 -05:00
|
|
|
real function fchisq(c3,npts,fsample,a)
|
2012-11-21 12:42:53 -05:00
|
|
|
|
2012-11-26 16:06:41 -05:00
|
|
|
parameter (NMAX=85*16)
|
|
|
|
|
complex c3(npts)
|
|
|
|
|
complex c4(NMAX)
|
2012-11-21 12:42:53 -05:00
|
|
|
real a(3)
|
2012-11-26 16:06:41 -05:00
|
|
|
complex z
|
|
|
|
|
complex w,wstep
|
|
|
|
|
data a1,a2,a3/99.,99.,99./
|
|
|
|
|
include 'jt9sync.f90'
|
2012-11-21 12:42:53 -05:00
|
|
|
save
|
|
|
|
|
|
|
|
|
|
if(a(1).ne.a1 .or. a(2).ne.a2 .or. a(3).ne.a3) then
|
|
|
|
|
a1=a(1)
|
|
|
|
|
a2=a(2)
|
|
|
|
|
a3=a(3)
|
2012-11-26 16:06:41 -05:00
|
|
|
call twkfreq(c3,c4,npts,fsample,a)
|
2012-11-21 12:42:53 -05:00
|
|
|
endif
|
|
|
|
|
|
2012-11-26 16:06:41 -05:00
|
|
|
! Get sync power.
|
|
|
|
|
nspsd=16
|
|
|
|
|
sum1=0.
|
|
|
|
|
sum0=0.
|
|
|
|
|
k=-1
|
|
|
|
|
do i=1,85
|
|
|
|
|
z=0.
|
|
|
|
|
do j=1,nspsd
|
|
|
|
|
k=k+1
|
|
|
|
|
z=z+c4(k+1)
|
|
|
|
|
enddo
|
|
|
|
|
pp=real(z)**2 + aimag(z)**2
|
|
|
|
|
if(isync(i).eq.1) then
|
|
|
|
|
sum1=sum1+pp
|
|
|
|
|
else
|
|
|
|
|
sum0=sum0+pp
|
2012-11-21 12:42:53 -05:00
|
|
|
endif
|
|
|
|
|
enddo
|
2012-11-26 16:06:41 -05:00
|
|
|
sync=(sum1/16.0)/(sum0/69.0) - 1.0
|
|
|
|
|
fchisq=-sync
|
2012-11-21 12:42:53 -05:00
|
|
|
|
|
|
|
|
return
|
|
|
|
|
end function fchisq
|
