WSJT-X/savetf2.F90

subroutine savetf2(id,nsave,nutc)
  parameter (NZ=60*96000)
  parameter (NSPP=174)
  parameter (NPKTS=NZ/NSPP)
  integer*2 id(4,NZ)
  real*4 ss(NPKTS),ss2(60)
  real*8 dt,t,t2

  dt=NSPP/96000.d0
  t=0.d0
  nh=nutc/100
  nm=mod(nutc,100)
  t2=3600*nh + 60*nm
  fac=1.0/(4.0*NSPP)

  do i=1,NPKTS
     s=0.
     do n=1,NSPP
        s=s + float(id(1,i))**2 + float(id(2,i))**2 + float(id(3,i))**2 + &
             float(id(4,i))**2 
     enddo
     ss(i)=fac*s
     t=t+dt
     t2=t2+dt
     if(nsave.eq.3) write(24,1010) t,t2,ss(i)
1010 format(f9.6,f15.6,f10.3)
  enddo

  if(nsave.eq.2) then
     dt2=551*dt
     t=0.d0
     t2=3600*nh + 60*nm
     k=0
     do i=1,60
        s=0.
        ns=0
        do n=1,551
           k=k+1
           s=s + ss(k)
           if(ss(k).gt.0.0) ns=ns+1
        enddo
        ss2(i)=s/ns
        t=t+dt
        t2=t2+dt2
        write(25,1010) t,t2,ss2(i)
     enddo
  endif

  return
end subroutine savetf2