diff --git a/lib/decode65a.f90 b/lib/decode65a.f90 index ec9fbd311..3b82178b2 100644 --- a/lib/decode65a.f90 +++ b/lib/decode65a.f90 @@ -41,9 +41,7 @@ subroutine decode65a(dd,npts,newdat,nqd,f0,nflip,mode65,ntrials, & if(nflip.eq.0) go to 900 ! Find best DF, drift, curvature, and DT. Start by downsampling to 344.53125 Hz - call timer('fil6521 ',0) call fil6521(cx,n5,c5x,n6) - call timer('fil6521 ',1) fsample=1378.125/4. @@ -60,10 +58,8 @@ subroutine decode65a(dd,npts,newdat,nqd,f0,nflip,mode65,ntrials, & ! Apply AFC corrections to the time-domain signal ! Now we are back to using the 1378.125 Hz sample rate, enough to ! accommodate the full JT65C bandwidth. - a(3)=0 - call timer('twkfreq ',0) + a(3)=0 call twkfreq65(cx,n5,a) - call timer('twkfreq ',1) ! Compute spectrum for each symbol. nsym=126 @@ -71,7 +67,6 @@ subroutine decode65a(dd,npts,newdat,nqd,f0,nflip,mode65,ntrials, & df=1378.125/nfft j=int(dtbest*1378.125) - call timer('sh_ffts ',0) c5a=cmplx(0.0,0.0) do k=1,nsym do i=1,nfft @@ -89,7 +84,6 @@ subroutine decode65a(dd,npts,newdat,nqd,f0,nflip,mode65,ntrials, & s1(jj,k)=real(c5a(i))**2 + aimag(c5a(i))**2 enddo enddo - call timer('sh_ffts ',1) call timer('dec65b ',0) qualbest=0. diff --git a/lib/qra/q65/q65_loops.f90 b/lib/qra/q65/q65_loops.f90 index 8b1df6a9c..ebec128c8 100644 --- a/lib/qra/q65/q65_loops.f90 +++ b/lib/qra/q65/q65_loops.f90 @@ -58,9 +58,7 @@ subroutine q65_loops(c00,npts2,nsps2,nsubmode,ndepth,jpk0, & if(mod(idt,2).eq.0) ndt=-ndt jpk=jpk0 + nsps2*ndt/16 !tsym/16 if(jpk.lt.0) jpk=0 - call timer('spec64 ',0) call spec64(c0,nsps2,mode_q65,jpk,s3,LL,NN) - call timer('spec64 ',1) call pctile(s3,LL*NN,40,base) s3=s3/base where(s3(1:LL*NN)>s3lim) s3(1:LL*NN)=s3lim diff --git a/map65/libm65/decode1a.f90 b/map65/libm65/decode1a.f90 index bf106080a..648a5e803 100644 --- a/map65/libm65/decode1a.f90 +++ b/map65/libm65/decode1a.f90 @@ -43,7 +43,6 @@ subroutine decode1a(dd,newdat,f0,nflip,mode65,nfsample,xpol, & z=cmplx(cos(dphi),sin(dphi)) cy(:n5)=z*cy(:n5) !Adjust for cable length difference endif - call timer('fil6521 ',0) ! Add some zeros at start of c5 arrays -- empirical fix for negative DT's nadd=1089 c5x(:nadd)=0. @@ -53,7 +52,6 @@ subroutine decode1a(dd,newdat,f0,nflip,mode65,nfsample,xpol, & call fil6521(cy,n5,c5y(nadd+1),n6) endif n6=n6+nadd - call timer('fil6521 ',1) fsample=1378.125/4. a(5)=dt00 @@ -80,9 +78,7 @@ subroutine decode1a(dd,newdat,f0,nflip,mode65,nfsample,xpol, & ! Now we are back to using the 1378.125 Hz sample rate, enough to ! accommodate the full JT65C bandwidth. - call timer('twkfreq_',0) call twkfreq_xy(cx,cy,n5,a) - call timer('twkfreq_',1) ! Compute spectrum at best polarization for each half symbol. ! Adding or subtracting a small number (e.g., 5) to j may make it decode.\ @@ -94,7 +90,6 @@ subroutine decode1a(dd,newdat,f0,nflip,mode65,nfsample,xpol, & j=(dt00+dtbest+2.685)*1378.125 if(j.lt.0) j=0 - call timer('sh_ffts ',0) ! Perhaps should try full-symbol-length FFTs even in B, C sub-modes? ! (Tried this, found no significant difference in decodes.) @@ -123,8 +118,6 @@ subroutine decode1a(dd,newdat,f0,nflip,mode65,nfsample,xpol, & enddo enddo - call timer('sh_ffts ',1) - flip=nflip call timer('dec65b ',0) call decode65b(s2,flip,mycall,hiscall,hisgrid,mode65,neme,ndepth, &