mirror of
https://github.com/saitohirga/WSJT-X.git
synced 2024-11-23 04:38:37 -05:00
560c7550db
git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/map65@418 ab8295b8-cf94-4d9e-aec4-7959e3be5d79
51 lines
993 B
Fortran
51 lines
993 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(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
|