PRO spec_all_hydro,del_itt=del_itt,comp=comp
;
; This routine plots all the kinetic energy spectra with a time 
; spacing of del_itt
;
; Define some default values
;
default,del_itt,10
;
; Read in data from file 
;
power,v1='u',spec1=spec,k=k,i=i,/noplot
;
; Find dimensions of the data
;
kmax=(size(spec))[1]
ns=(size(spec))[2]
print,'ns=',ns
;
;  factor (possibility of compensated spectra)
;
fac=1.
if keyword_set(comp) then fac=k^1.667
;
; Plot data
;
plot_oo,k,fac*spec(*,0),xr=[1,kmax]
;oplot,k,k^(+1.67)*spec(1,0),li=2
for i=0,ns-1,del_itt do begin
    oplot,k,fac*spec(*,i)
end
oplot,k,k^(-1.67)*spec(1,ns-1),li=2
;
END