From 76d4393537e1e66fe34351af4c629078963e9b4e Mon Sep 17 00:00:00 2001 From: Joe Taylor Date: Fri, 8 Jan 2016 20:55:17 +0000 Subject: [PATCH] Commit several tentative attempts (now commented out) addressing a possible JT9 decoding issue. git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/wsjtx@6366 ab8295b8-cf94-4d9e-aec4-7959e3be5d79 --- lib/afc9.f90 | 3 ++- lib/fchisq.f90 | 7 +++++-- lib/softsym.f90 | 4 +++- 3 files changed, 10 insertions(+), 4 deletions(-) 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.