if !d.name eq 'PS' then begin
  device,xsize=18,ysize=21,yoffset=3
  !p.charthick=4 & !p.thick=4 & !x.thick=4 & !y.thick=4
end
power,'_mag','hel_mag',k=k,spec1=MM,spec2=HH,i=n,tt=t,/noplot & print,n
power,'_kin','hel_kin',k=k,spec1=EE,spec2=FF,i=n,tt=t,/noplot & print,n
col1=122
;
;  example; must be set in local param.pro file
;yout1_M=.009 & yout1_C=.0004 & yout2_M=.010 & yout2_C=.0010
;
;  mv idl.ps ~/tex/mhd/alpha_B0dedt/fig/ppowerhel.ps
;  mv idl.ps ~/tex/mhd/alpha_B0dedt/fig/ppowerhelM512b3.ps
;  mv idl.ps ~/tex/mhd/alpha_B0dedt/fig/ppowerhelH512b3.ps
;
;  $Id: ppowerhel.pro,v 1.5 2005/04/09 03:51:32 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 param.pro'
@param
xr=[1,max(k)]
default,t1,.0
!p.multi=[0,1,2]
!p.charsize=2.
!p.title='!6'
!y.margin=[3.3,.8]
;
;  calculate dissipation
;
default,nu,2e-4
default,eta,2e-4
epsK=nu*total(k^2*EE)
epsM=eta*total(k^2*MM)
epsT=epsK+epsM
print,'epsK,epsM,epsT=',epsK,epsM,epsT
;
i=0
s=1./epsT^.6667
!x.title='!6'
half='!s!u 1!n!r!s-!r!d 2!n'
;
plot_oo,k,s*MM(*,i)*k^1.667,xr=xr,yr=[1e-2,5e-0]
oplot,k,.5*s*HH(*,i)*k^2.6667,col=col1
oplot,k,s*EE(*,i)*k^1.6667,li=1
oplot,k,.5*s*FF(*,i)*k^0.6667,li=1,col=col1
xyouts,20,yout1_M,'!8k!6!u5/3!n!7e!6!u-2/3!n!8M!6(!8k!6)'
xyouts,96,yout1_M*.5,'!8E!6(!8k!6)'
xyouts,1.6,yout1_C,half+'!8k!6!u2/3!n!7e!6!u-2/3!n!8C!6(!8k!6)',col=col1
xyouts,10,yout1_C*.5,'!8F!6(!8k!6)',col=col1
;
!x.title='!8k!6'
plot_oo,k,.5*s*HH(*,i)*k^3.667,xr=xr,yr=[1e-1,1e1]
oplot,k,.5*s*FF(*,i)*k^1.6667,li=1
xyouts,17,yout2_C,half+'!8k!6!u5/3!n!7e!6!u-2/3!n!8C!6(!8k!6)'
xyouts,2.3,yout2_F,half+'!8k!6!u5/3!n!7e!6!u-2/3!n!8F!6(!8k!6)'
;
END