;
; look at ri data cband
file='/share/olda/datafile.24sep04.x101.1'
openr,lun,file,/get_lun

nrec=952L
istat=riget(lun,d,numrec=952)

nsmp=10000L
da=reform(d.d1[1,*],nsmp*nrec)
db=reform(d.d1[0,*],nsmp*nrec)
ntot=nsmp*nrec
nn=nsmp*nrec/2
spca=abs(fft(da))
spca=spca[0:nn-1L]
spcb=abs(fft(db))
spcb=spcb[0:nn-1l]
;
;
tm=findgen(ntot)*10e-6
tosmo=1000L
tms=findgen(ntot/tosmo)*10e-6*tosmo
das=total(reform(da,tosmo,ntot/tosmo),1)/tosmo
dbs=total(reform(db,tosmo,ntot/tosmo),1)/tosmo
hard=1
if hard then pscol,'24sep04_cboscillations.ps',/full
!p.multi=[0,1,2]
hor,0,10
ver,-.02,.2
plot,tms,das/median(das) - 1,$
	xtitle='time secs',ytitle='pwr [linear]',$
	title='24sep04 cband pwr vs time. oscillations in polA'
oplot,tms,dbs/median(dbs) - 1.,color=2
ln=2.
scl=.7
xp=.02
note,ln,'Smooth to 10 millseconds',xp=xp 
note,ln+1*scl,'polA',xp=xp 
note,ln+2*scl,'polB',xp=xp,color=2 
;
ver,0,10
hor,0,20
plot,frq[ind],spca[ind],$
	xtitle='hz',ytitle='spectra of time series',$
	title='24sep04 cb spectra of time series.'
ln=18
note,ln,'Horizontal lines ever 1.2 hz',xp=xp,color=4
oplot,frq[ind],spcb[ind],color=2
flag,findgen(20)*1.2001139,linestyle=1,color=4
if hard then hardcopy
x
;