file='kpm.txt'
a=rtable(file,6,head=1)
t=reform(a[1,*])
EEM=reform(a[2,*])
xiM=reform(a[3,*])
;
;plot_io,t,EEM*t/xiM^2;,yr=[1e-2,1e-0]
plot_oo,t,EEM*t/xiM^2
;
END