diff --git a/lib/afc9.f90 b/lib/afc9.f90 index bcf291ae8..dd106aa48 100644 --- a/lib/afc9.f90 +++ b/lib/afc9.f90 @@ -54,13 +54,14 @@ subroutine afc9(c3a,npts,fsample,a,syncpk) enddo if(a(3).ne.a3) call shft(c3a,a(3),a3,c3) chisqr=fchisq(c3,npts,fsample,a) +! write(*,4000) 0,0,a,-chisqr if(chisqr/chisqr0.gt.0.99) exit chisqr0=chisqr enddo syncpk=-chisqr c3a=c3 -! write(*,4001) a,-chisq2 +! write(*,4001) a,syncpk !4001 format(3x,3f10.4,f11.3) return diff --git a/lib/fchisq.f90 b/lib/fchisq.f90 index 6eb8b6b68..1aa2b0bb7 100644 --- a/lib/fchisq.f90 +++ b/lib/fchisq.f90 @@ -34,8 +34,11 @@ real function fchisq(c3,npts,fsample,a) sum0=sum0+pp endif enddo - sync=sum1/10000.0 - fchisq=-sync + sync_4992=(sum1/16.0)/(sum0/69.0) - 1.0 !r4992 + sync_4993=sum1/10000.0 !r4993+ +! write(80,3001) 1.e-5*sum1,1.e-5*sum0,sync_4992,sync_4993,sync +!3001 format(5f11.4) + fchisq=-sync_4993 return end function fchisq diff --git a/lib/softsym.f90 b/lib/softsym.f90 index 120fb6df9..6475b0e46 100644 --- a/lib/softsym.f90 +++ b/lib/softsym.f90 @@ -28,10 +28,12 @@ subroutine softsym(id2,npts8,nsps8,newdat,fpk,syncpk,snrdb,xdt, & fsample=1500.0/ndown a=0. call timer('afc9 ',0) - call afc9(c3,nz3,fsample,a,syncpk) !Find deltaF, fDot, fDDot + call afc9(c3,nz3,fsample,a,syncpk) !Find deltaF, fDot, extra DT call timer('afc9 ',1) freq=fpk - a(1) drift=-2.0*a(2) +! write(*,3301) fpk,freq,a +!3301 format(2f9.3,3f10.4) a3=a(3) a(3)=0.