Add savetf2.F90

git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/map65@417 ab8295b8-cf94-4d9e-aec4-7959e3be5d79
This commit is contained in:
Joe Taylor 2007-06-25 23:51:09 +00:00
parent 55dfeabf15
commit 74dd138074
2 changed files with 54 additions and 1 deletions

View File

@ -1096,7 +1096,10 @@ def update():
bmtext.insert(END,'Freq DF Pol UTC\n')
bmtext.insert(END,'----------------------------------------\n')
for i in range(len(lines)):
try:
nage=int(lines[i][41:])
except:
nage=0
lines[i]=lines[i][:41]
if nage==0: attr='age0'
if nage==1: attr='age1'

50
savetf2.F90 Normal file
View File

@ -0,0 +1,50 @@
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(f6.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