WSJT-X/qmap/libqmap/chkstat.f90

29 lines
476 B
Fortran

subroutine chkstat(dd,nhsym,pdb)
real dd(2,5760000)
real pdb(4)
integer ia(4),ib(4)
logical*1 btx0,btx1
btx0=.false.
btx1=.false.
ia(1)=23*96000+1
ib(1)=24*96000
ia(2)=27*96000+1
ib(2)=28*96000
ia(3)=53*96000+1
ib(3)=54*96000
ia(4)=57*96000+1
ib(4)=58*96000
do j=1,4
sq=0.
do i=ia(j),ib(j)
sq=sq + dd(1,i)*dd(1,i) + dd(2,i)*dd(2,i)
enddo
pdb(j)=db(1.0 + sq/(2.0*96000.0))
enddo
return
end subroutine chkstat