mirror of
https://github.com/saitohirga/WSJT-X.git
synced 2024-11-18 01:52:05 -05:00
6d6596500a
git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/map65@2434 ab8295b8-cf94-4d9e-aec4-7959e3be5d79
83 lines
2.2 KiB
Fortran
83 lines
2.2 KiB
Fortran
program pulsar
|
|
|
|
! Receives timf2 data from Linrad and saves it for pulsar processing.
|
|
|
|
parameter (NSPP=174) !Samples per UDP packet
|
|
integer*1 userx_no,iusb
|
|
integer*2 id
|
|
integer*2 nblock
|
|
integer*2 id2(1000)
|
|
integer nt(8)
|
|
character pname*12,cdate*8,ctime*12,czone*8
|
|
character*40 fname20,fname21
|
|
real*8 center_freq,tsec
|
|
common/plrscom/center_freq,msec,fselect,iptr,nblock,userx_no,iusb,id(4,NSPP)
|
|
! 8 4 4 4 2 1 1 1392
|
|
data nb0/0/,sqave/0.0/,u/0.001/,nw/0/
|
|
save
|
|
|
|
nargs=iargc()
|
|
if(nargs.ne.1) then
|
|
print*,'Usage: pulsar <pname>'
|
|
go to 999
|
|
endif
|
|
call getarg(1,pname)
|
|
|
|
! nt: year, month, day, time difference in minutes, hours, minutes,
|
|
! seconds, and milliseconds.
|
|
call date_and_time(cdate,ctime,czone,nt)
|
|
|
|
write(fname20,1001) nt(1)-2000,nt(2),nt(3),nt(5),nt(6),nt(7)
|
|
1001 format(3i2.2,'_',3i2.2,'.raw')
|
|
open(20,file=fname20,access='stream',status='unknown')
|
|
write(20) cdate,ctime,czone,nt,pname !Write header
|
|
fname21=fname20(1:13)//'.sq'
|
|
open(21,file=fname21,access='stream',status='unknown')
|
|
write(21) cdate,ctime,czone,nt,pname !Write header
|
|
|
|
call setup_rsocket !Prepare to receive UDP packets
|
|
k=0
|
|
|
|
10 call recv_pkt(center_freq)
|
|
nb=nblock
|
|
if(nb.lt.0) nb=nb+65536
|
|
lost=nb-nb0-1
|
|
if(lost.eq.-65536) lost=0
|
|
if(lost.ne.0 .and. nb0.ne.0) print*,'Lost packets:',lost,nb,nb0
|
|
|
|
write(20) id !Write raw data
|
|
sq=0.
|
|
do i=1,NSPP
|
|
sq=sq + float(int(id(1,i)))**2 + float(int(id(2,i)))**2 + &
|
|
float(int(id(3,i)))**2 + float(int(id(4,i)))**2
|
|
enddo
|
|
sqave=sqave + u*(sq-sqave)
|
|
rxnoise=10.0*log10(sqave) - 48.0
|
|
|
|
k=k+1
|
|
id2(k)=0.001*sq
|
|
if(k.eq.1000) then
|
|
tsec=0.001d0*msec
|
|
nsec=tsec
|
|
ih=nsec/3600
|
|
im=mod(nsec/60,60)
|
|
is=mod(nsec,60)
|
|
nw=nw+1
|
|
write(*,1000) nw,center_freq,ih,im,is,nblock,sqave,rxnoise,id2(1)
|
|
write(13,1000) nw,center_freq,ih,im,is,nblock,sqave,rxnoise,id2(1)
|
|
1000 format(i10,f7.3,i4,2i3.2,i7,f10.0,f8.2,i8)
|
|
write(21) id2
|
|
|
|
call flush(13)
|
|
call flush(20)
|
|
call flush(21)
|
|
k=0
|
|
endif
|
|
|
|
nb0=nb
|
|
go to 10
|
|
|
|
999 end program pulsar
|
|
|
|
! To compile: % gfortran -o pulsar pulsar.f90 plrr_subs.c
|