power,'u','_ro',k=k,spec1=spec1,spec2=spec2,i=n,tt=t,/noplot,datatopdir='../hd_256_4x8f/data'
power,'u','_ro',k=kb,spec1=spec1b,spec2=spec2b,i=nb,tt=tb,/noplot,datatopdir='../hd_512_8x4a/data'
power,'r2u','r3u',k=k,spec1=spec3,spec2=spec4,i=n,tt=t,/noplot,datatopdir='../hd_256_4x8f/data'
power,'r2u','r3u',k=kb,spec1=spec3b,spec2=spec4b,i=nb,tt=tb,/noplot,datatopdir='../hd_512_8x4a/data'
;
;  $Id: pspec.pro,v 1.1 2009/10/10 23:07:24 brandenb Exp $
;
;  This routine is used to view all power spectra
;  that were written during the run (controlled by dspec)
;  The plot range yr is set to a default value.
;  The time interval of plotting is "w" (default value)
;
;   5-oct-02/axel: written
;
;  mv idl.ps pspec1.ps
;  mv idl.ps pspec2.ps
;
spawn,'if (! -f param.pro) touch param.pro'
@param
xr=[1,max(k)]
yr=[4e-1,4e+2]
default,yr,[1e-7,3e+2]
default,w,.01
default,t1,.0
siz=1.8
;
nexp=2.
nexp=1.66667
comp=k^nexp
compb=kb^nexp
good=where(t ge t1)
xmax=max(k)
!x.range=[1,256]
!p.multi=[0,1,3]
!p.charsize=3.
;
;-----------------------------------------------------------------------------
!y.title='!8P!6(!8u!6)'
plot_oo,k,spec1(*,0)*comp,yr=yr,li=1,/nodata
oplot,kb,spec1b(*,0)*compb
for i=9,n-2,10 do begin
  if(t(i) ge t1) then begin
    print,t(i)
    ;oplot,k,spec1(*,i)*comp,li=1
    wait,w
  end
end
;
for i=9,nb-2,10 do begin
  if(t(i) ge t1) then begin
    print,tb(i)
    oplot,kb,spec1b(*,i)*compb
    wait,w
  end
end
;
;-----------------------------------------------------------------------------
!y.title='!8P!6(!7q!6!u1/2!n !8u!6)'
plot_oo,k,spec3(*,0)*comp,yr=yr,li=1,/nodata
oplot,kb,spec3b(*,0)*compb
for i=9,n-2,10 do begin
  if(t(i) ge t1) then begin
    print,t(i)
    ;oplot,k,spec3(*,i)*comp,li=1
    wait,w
  end
end
;
for i=9,nb-2,10 do begin
  if(t(i) ge t1) then begin
    print,tb(i)
    oplot,kb,spec3b(*,i)*compb
    wait,w
  end
end
;
;-----------------------------------------------------------------------------
!y.title='!8P!6(!7q!6!u1/3!n !8u!6)'
plot_oo,k,spec4(*,0)*comp,yr=yr,li=1,/nodata
oplot,kb,spec4b(*,0)*compb
for i=9,n-2,10 do begin
  if(t(i) ge t1) then begin
    print,t(i)
    ;oplot,k,spec4(*,i)*comp,li=1
    wait,w
  end
end
;
for i=9,nb-2,10 do begin
  if(t(i) ge t1) then begin
    print,tb(i)
    oplot,kb,spec4b(*,i)*compb
    wait,w
  end
end
;
;save,file='spec_aver.sav',k,Ek
END