;$Id: muller_grappin.pro,v 1.3 2006/03/20 15:46:05 nilshau Exp $ if !d.name eq 'PS' then begin device,xsize=14,ysize=10,yoffset=3 !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3 end ; ; mv idl.ps ../fig/MullerGrappin.ps ; !p.charsize=1.4 !x.title='!8k!6/!8k!6!dd,eff!n' !y.title='!8E!6(!8k!6) !8k!6!u5/3!n' ; k=10^grange(-0.5,4.5,101) kK=1.5 & kM=5. & kd=5e3 & aT=1. & aR=2. cutoff=exp(-(k/kd)) subinK=1./(1.+(k/kK)^(-11./3.)) subinM=1./(1.+(k/kM)^(-11./3.)) ET=aT*k^(-1.66667)*cutoff ;*subinertial ER=aR*k^(-2.33333)*cutoff ;*subinertial ; EM=.5*(ET+ER)*subinM EK=(ET-EM)*subinK ; comp=k^1.66667 plot_oo,k/kd,comp*EK,xr=[3e-4,3e0],yr=[.01,1.] oplot,k/kd,comp*EM,li=2 ; xx=.002 & dx=3.*xx & siz=1.4 legend,xx,dx,.08,0,' !6Kinetic',siz=siz legend,xx,dx,.04,2,' !6Magnetic',siz=siz END