How to plot the electron velocity distribution? about irfu-matlab HOT 2 CLOSED

Hi,

There are a couple of reasons why the distribution is different. mms.plot_projection plots a slice of the 3D distribution, whereas in the Torbert paper they are plotting a 2D reduced distribution, meaning the third velocity axis is integrated over (note that the units are different between the two methods). It can look different because you are plotting something different to the Torbert paper. Also in your script you are using MMS1 data, whereas Torbert et al. is using MMS3 data, so this may also account for a different looking distribution.

Here is an example of the 2D reduced distribution seen in Torbert et al. (see also Example_MMS_reduced_ele_dist_2D.m):
vg = linspace(-40e3,40e3,100); % km/s
nMC = 2e3;
f2D = dist.reduce('2D',hatB0,hatExB0,'base','cart','vg',vg,'nMC',nMC,'lowerelim',100);

fn=figure;
set(fn,'Position',[10 10 400 350])
h(1)=axes('position',[0.12 0.13 0.8 0.8]); % [x y dx dy]
ud=get(fn,'userdata');
ud.subplot_handles=h;
set(fn,'userdata',ud);
set(fn,'defaultLineLineWidth',1);

f2D.plot_plane(h(1),'docolorbar',0);
caxis(h(1),[-14 -10]);
hcb = colorbar(h(1));
ylabel(hcb,'log_{10} f_e (s^2 m^{-5})','interpreter','tex','fontsize',12)
hold(h(1),'on')
plot(h(1),[0 0 ],max(abs(vg))*[-1 1],'color','k','linewidth',1)
plot(h(1),max(abs(vg))*[-1 1],[0 0 ],'color','k','linewidth',1)
hold(h(1),'off')
grid(h(1),'off')
axis(h(1),'equal')
colormap(h(1),'jet')
irf_legend(h(1),'(c)',[0.98 0.85],'color','w','fontsize',12)
axis(h(1),1e-3*[min(vg) max(vg) min(vg) max(vg)])
xlabel(h(1),'V_{||} (10^3 km s^{-1})','interpreter','tex','fontsize',12)
ylabel(h(1),'V_{\perp} (10^3 km s^{-1})','interpreter','tex','fontsize',12)
set(gcf,'color','w')

from irfu-matlab.

thanks a lot!

from irfu-matlab.