;$Id: ppower_ro.pro,v 1.1 2021/04/09 14:32:48 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=12,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
cwd,run
print,' $mv idl.ps ~/tex/tina/rospec/fig/ppower_ro_'+run+'.ps'
;
!p.charsize=1.7
!x.margin=[7.6,.5]
!y.margin=[3.2,.5]
;
;  sample to compare time series form different directories
;
@parameters
;
default,kin,0
default,iread,0
if iread eq 0 then begin
power,'_kin','_ro',k=k,spec1=spec1 ,spec2=spec2 ,i=n,tt=t,/noplot,/lks
;
if kin then begin
  print,'kin'
endif else begin
  ;power,'_mag','hel_mag',k=k,spec1=magn1,spec2=magn2,i=n,tt=t,/noplot,/lks
  print,'mag'
endelse
endif
;
yr=[1e-8,2e-3]
nt=n_elements(t)
!x.title='!8k!6'
!y.title='!8E!6(!8k!6)'
plot_oo,k[1:*],spec2[1:*,nt-1],/nodata,yr=yr
;
s=.14
s=2.2
for it=80,100,5 do begin
  oplot,k[1:*],spec2[1:*,it]*s,col=122
  oplot,k[1:*],spec1[1:*,it],col=55
  wait,.1
endfor
;
END