diff --git a/qmap/libqmap/chkstat.f90 b/qmap/libqmap/chkstat.f90 index 760aa47df..7f3121892 100644 --- a/qmap/libqmap/chkstat.f90 +++ b/qmap/libqmap/chkstat.f90 @@ -1,22 +1,28 @@ -subroutine chkstat(dd,nhsym,dbdiff) +subroutine chkstat(dd,nhsym,pdb) real dd(2,5760000) + real pdb(4) + integer ia(4),ib(4) + logical*1 btx0,btx1 - sq0=0. - sq1=0. - k=0 - do i=1,60 + btx0=.false. + btx1=.false. + ia(1)=23*96000+1 + ib(1)=24*96000 + ia(2)=27*96000+1 + ib(2)=28*96000 + ia(3)=53*96000+1 + ib(3)=54*96000 + ia(4)=57*96000+1 + ib(4)=58*96000 + + do j=1,4 sq=0. - do j=1,96000 - k=k+1 - sq=sq + dd(1,k)*dd(1,k) + dd(2,k)*dd(2,k) - enddo - if(i.ge.12 .and. i.le.24) sq0=sq0+sq - if(i.ge.42 .and. i.le.54) sq1=sq1+sq + do i=ia(j),ib(j) + sq=sq + dd(1,i)*dd(1,i) + dd(2,i)*dd(2,i) + enddo + pdb(j)=db(1.0 + sq/(2.0*96000.0)) enddo - db0=db(1.0+sq0) - db1=db(1.0+sq1) - dbdiff=db0-db1 return end subroutine chkstat diff --git a/qmap/libqmap/q65c.f90 b/qmap/libqmap/q65c.f90 index e2f700e2d..bcc1f2867 100644 --- a/qmap/libqmap/q65c.f90 +++ b/qmap/libqmap/q65c.f90 @@ -10,6 +10,7 @@ subroutine q65c(itimer) parameter (NFFT=32768) include 'njunk.f90' real*8 fcenter + real*4 pdb(4) integer nparams0(NJUNK+3),nparams(NJUNK+3) logical first logical*1 bAlso30 @@ -41,12 +42,12 @@ subroutine q65c(itimer) npatience=1 newdat=1 !Always on ?? - call chkstat(dd,nhsym,dbdiff) -! if(nhsym.eq.390) write(*,3001) nutc,nhsym,dbdiff -!3001 format(i4.4,i6,f7.1) - - if(dbdiff.gt.5.0) ntx30b=30 - if(dbdiff.lt.-5.0) ntx30a=30 + call chkstat(dd,nhsym,pdb) + if(abs(pdb(1)-pdb(2)).gt.3.0) ntx30a=20 !We transmitted in first half + if(abs(pdb(3)-pdb(4)).gt.3.0) ntx30b=20 !We transmitted in second half + +! write(*,3001) nutc,nhsym,pdb,ntx30a,ntx30b +!3001 format(i4.4,i6,4f7.1,2i6) if(ntx30a.gt.5) then dd(1:2,1:30*96000)=0.