2007-07-09 16:36:44 -04:00
|
|
|
subroutine savetf2(id,nsave,ntime,nutc,savedir)
|
|
|
|
|
|
|
|
#ifdef Win32
|
|
|
|
use dfport
|
|
|
|
#else
|
|
|
|
external gmtime_r
|
|
|
|
#endif
|
|
|
|
|
2007-06-25 19:51:09 -04:00
|
|
|
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
|
2007-07-09 16:36:44 -04:00
|
|
|
character*80 savedir,fname
|
|
|
|
integer it(9)
|
|
|
|
|
|
|
|
if(nsave.eq.1) then
|
|
|
|
n2=ntime/60
|
|
|
|
n3=(n2-1)*60
|
|
|
|
#ifdef Win32
|
|
|
|
call gmtime(n3,it)
|
|
|
|
#else
|
|
|
|
call gmtime_r(n3,it)
|
|
|
|
#endif
|
|
|
|
it(5)=it(5)+1
|
|
|
|
it(6)=mod(it(6),100)
|
|
|
|
nhr=nutc/100
|
|
|
|
nmin=mod(nutc,100)
|
|
|
|
write(fname,1000) it(6),it(5),it(4),nhr,nmin
|
|
|
|
1000 format('/',3i2.2,'_',2i2.2,'.tf2')
|
|
|
|
do i=80,1,-1
|
|
|
|
if(savedir(i:i).ne.' ') go to 1
|
2007-06-25 19:51:09 -04:00
|
|
|
enddo
|
2007-07-09 16:36:44 -04:00
|
|
|
|
|
|
|
1 iz=i
|
|
|
|
fname=savedir(1:iz)//fname
|
|
|
|
#ifdef Win32
|
|
|
|
open(17,file=fname,status='unknown',form='binary',err=998)
|
|
|
|
write(17) id
|
|
|
|
close(17)
|
|
|
|
#else
|
|
|
|
open(17,file=fname,status='unknown',form='unformatted', &
|
|
|
|
access='direct',recl=nbytes,err=998)
|
|
|
|
write(17,rec=1) id
|
|
|
|
close(17)
|
|
|
|
#endif
|
|
|
|
|
|
|
|
|
|
|
|
else
|
|
|
|
dt=NSPP/96000.d0
|
2007-06-25 19:51:09 -04:00
|
|
|
t=0.d0
|
2007-07-09 16:36:44 -04:00
|
|
|
nh=nutc/100
|
|
|
|
nm=mod(nutc,100)
|
2007-06-25 19:51:09 -04:00
|
|
|
t2=3600*nh + 60*nm
|
2007-07-09 16:36:44 -04:00
|
|
|
fac=1.0/(4.0*NSPP)
|
|
|
|
do i=1,NPKTS
|
2007-06-25 19:51:09 -04:00
|
|
|
s=0.
|
2007-07-09 16:36:44 -04:00
|
|
|
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
|
2007-06-25 19:51:09 -04:00
|
|
|
enddo
|
2007-07-09 16:36:44 -04:00
|
|
|
ss(i)=fac*s
|
2007-06-25 19:51:09 -04:00
|
|
|
t=t+dt
|
2007-07-09 16:36:44 -04:00
|
|
|
t2=t2+dt
|
|
|
|
if(nsave.eq.3) write(24,1010) t,t2,ss(i)
|
|
|
|
1010 format(f9.6,f15.6,f10.3)
|
2007-06-25 19:51:09 -04:00
|
|
|
enddo
|
2007-07-09 16:36:44 -04:00
|
|
|
|
|
|
|
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
|
2007-06-25 19:51:09 -04:00
|
|
|
endif
|
2007-07-09 16:36:44 -04:00
|
|
|
go to 999
|
|
|
|
|
|
|
|
998 print*,'Cannot open file:'
|
|
|
|
print*,fname
|
2007-06-25 19:51:09 -04:00
|
|
|
|
2007-07-09 16:36:44 -04:00
|
|
|
999 return
|
2007-06-25 19:51:09 -04:00
|
|
|
end subroutine savetf2
|