power,'_u.b','_bEP',k=k,spec1=spec1,spec2=spec2,i=n,tt=t,/noplot & print,n
power,'_kin','hel_kin',k=k,spec1=specK,spec2=specF,i=n,tt=t,/noplot & print,n
;
;  $Id: ppower_EP.pro,v 1.3 2009/08/30 13:35:22 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
;
spawn,'if (! -f param.pro) touch parameters.pro'
@parameters
xr=[1,max(k)]
yr=[1e-6,1e+1]
istride=1
default,w,.2
default,t1,.0
default,iover,0
default,istride,1
ampl15=3e-10 & ampl4=1e-13
ampl15=3e-6 & ampl4=1e-8
;
for i=0,n-2,istride do begin
  if(t(i) ge t1) then begin
    if iover eq 0 or i eq 0 then begin
      plot_oo,k,spec1(*,i),xr=xr,yr=yr,tit='t='+str(t(i))
    endif else begin
      oplot,k,spec1(*,i)
    endelse
    oplot,k,spec2(*,i),col=122
    oplot,k,specK(*,i),col=55
    oplot,k,specF(*,i),col=188,li=2
    kk=[1,20] & oplot,kk,ampl15*kk^1.5,col=188
    kk=[1,20] & oplot,kk,.01*ampl15*kk^1.5,col=188
    kk=[1,20] & oplot,kk,.001*ampl15*kk^2;,col=188
    ;kk=[1,5] & oplot,kk,ampl4*kk^4,col=55
    print,i
    wait,w
  end
end
;
print,'E_K (white), E_M (red), E_T (yellow)'
;
; averages
;
i=n-2
ts=t(i)
spec1s=spec1(*,i)
spec2s=spec2(*,i)
specKs=specK(*,i)
;
save,file='specs.sav',k,ts,spec1s,spec2s,specKs
;
END