mirror of https://github.com/saitohirga/WSJT-X.git
51 lines
1010 B
Fortran
51 lines
1010 B
Fortran
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(int(id(1,i)))**2 + float(int(id(2,i)))**2 + &
|
|
float(int(id(3,i)))**2 + float(int(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
|