2024-01-03 17:36:46 -05:00
|
|
|
subroutine save_qm(fname,nutc,dd,ntx30a,ntx30b)
|
2024-01-03 16:17:57 -05:00
|
|
|
|
|
|
|
|
parameter(NMAX=60*96000)
|
2024-01-03 17:36:46 -05:00
|
|
|
character*(*) fname
|
2024-01-03 16:17:57 -05:00
|
|
|
real*4 dd(2,NMAX)
|
|
|
|
|
integer*1 id1(2,NMAX)
|
2024-01-03 17:36:46 -05:00
|
|
|
|
2024-01-03 16:17:57 -05:00
|
|
|
ia=1
|
|
|
|
|
ib=NMAX
|
|
|
|
|
if(ntx30a.gt.5) ia=NMAX/2+1
|
|
|
|
|
if(ntx30b.gt.5) ib=NMAX/2
|
|
|
|
|
|
|
|
|
|
sq=0.
|
|
|
|
|
do i=ia,ib
|
|
|
|
|
x=dd(1,i)
|
2024-01-03 17:36:46 -05:00
|
|
|
y=dd(2,i)
|
2024-01-03 16:17:57 -05:00
|
|
|
sq=sq + x*x + y*y
|
|
|
|
|
enddo
|
2024-01-03 17:36:46 -05:00
|
|
|
nsum=2*(ib-ia+1)
|
2024-01-03 16:17:57 -05:00
|
|
|
rms=sqrt(sq/nsum)
|
|
|
|
|
|
2024-01-03 17:36:46 -05:00
|
|
|
jz=len(fname)
|
|
|
|
|
fname(jz-1:jz)="qm"
|
|
|
|
|
write(*,3001) fname(jz-13:jz),nutc,rms,db(rms*rms),ia,ib, &
|
|
|
|
|
nsum/(2*96000),ntx30a,ntx30b
|
|
|
|
|
3001 format(a14,2x,i4.4,2f7.1,2i9,3i5)
|
2024-01-03 16:17:57 -05:00
|
|
|
|
|
|
|
|
return
|
|
|
|
|
end subroutine save_qm
|
