correct mass for PDist.macroparticles

irfu · Dec 9, 2024 · c43a0c6 · c43a0c6
1 parent afc3d8e
commit c43a0c6
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/PDist.m b/PDist.m
@@ -1672,7 +1672,7 @@
                   if doScpot, tmp_energy = tmp_energy - scpot.data(it); end
                   tmp_azim   = azim_edges(iAzim)     + dazim*rand(N,1);   % deg
                   tmp_polar  = polar_edges(iPolar)   + dpolar*rand(N,1);  % deg
-                  tmp_v = sqrt(tmp_energy*units.eV*2/units.me)/1000; % km/s
+                  tmp_v = sqrt(tmp_energy*units.eV*2/mass)/1000; % km/s
 
                   if doScpot % check if energy is negative, then skip
                     if iEnergy>6;
@@ -3982,6 +3982,7 @@
 
       % Data from PDist in spherical coordinate system
       sizedata = obj.datasize;
+      mass = obj.mass;
 
       switch method
         case 'center' % add particles to the center of the bin