diff --git a/map65/libm65/wideband_sync.f90 b/map65/libm65/wideband_sync.f90 index f358d01d7..c71c5ebf5 100644 --- a/map65/libm65/wideband_sync.f90 +++ b/map65/libm65/wideband_sync.f90 @@ -7,7 +7,7 @@ module wideband_sync contains -subroutine wb_sync(ss,savg) +subroutine wb_sync(ss,savg,ntone_spacing) ! Compute "orange sync curve" using the Q65 sync pattern @@ -15,6 +15,8 @@ subroutine wb_sync(ss,savg) real ss(4,322,NFFT) real savg(4,NFFT) logical first + character*1 c1 +! integer hist(0:20) integer isync0(22) ! Q65 sync symbols data isync0/1,9,12,13,15,22,23,26,27,33,35,38,46,50,55,60,62,66,69,74,76,85/ @@ -60,12 +62,36 @@ subroutine wb_sync(ss,savg) sync_dat(i,2)=ccfmax sync_dat(i,3)=xdt sync_dat(i,4)=ipolbest -! f0=0.001*i*df -! write(70,3010) f0,sync_dat(i,2:3),nint(sync_dat(i,4)),0.001*i*df,nkhz_center -!3010 format(3f10.3,i5,f10.3,i5) - enddo + call pctile(sync_dat(ia:ib,2),ib-ia+1,50,base) + sync_dat(ia:ib,2)=sync_dat(ia:ib,2)/base +! hist=0 +! s2_avg=63.5*ntone_spacing +! do i=ia,ib +! f0=0.001*i*df +! write(70,3010) f0,sync_dat(i,2:3),nint(sync_dat(i,4)), & +! 0.001*i*df,nkhz_center +!3010 format(3f10.3,i5,f10.3,i5) +! x=min(sync_dat(i,2),20.0) +! nx=x +! hist(nx)=hist(nx)+1 +! if(x.gt.2.5) then +! c1=' ' +! s0=sync_dat(i,2) - 1.0 +! s1=(sync_dat(i-1,2) + sync_dat(i+1,2) - 2.0)/s0 +! s2=(sum(sync_dat(i+ntone_spacing/2+1:i+64*ntone_spacing,2)) - s2_avg)/s0 +! if(s2.ge.0.5) c1='*' +! write(72,3072) f0,s0,s1,s2,c1 +!3072 format(f12.6,3f10.3,2x,a1) +! endif +! enddo + +! do i=0,20 +! write(71,3071) i,hist(i) +!3071 format(i2,i8) +! enddo + return end subroutine wb_sync