looking at the code and comparing to the reference it looks like slightly different expressions for the fast ion density and pressure are used
around line 343
pressa = i4 .* taus .* 2.0 ./ 3.0 .* (sion_particles .* particle_energy .* mks.e)
cion.density_fast .+= sion_particles .* taut
In the reference the density is
and pressure
If I sub density into pressure this gives
which is similar to the expression given in the code, but the code doesn't look to include the (v_c/v_alpha) term.
I also looked at the density expression and it looks slightly different as well. The reference uses tau_slow and (v_c/v_alpha) but the code looked to use tau_therm defined as tau_s * log((vf3 + vc3) / vc3) / 3
I haven't double checked this but thought I'd post it before I forget!
looking at the code and comparing to the reference it looks like slightly different expressions for the fast ion density and pressure are used
around line 343
In the reference the density is
and pressure
If I sub density into pressure this gives
which is similar to the expression given in the code, but the code doesn't look to include the (v_c/v_alpha) term.
I also looked at the density expression and it looks slightly different as well. The reference uses tau_slow and (v_c/v_alpha) but the code looked to use tau_therm defined as
tau_s * log((vf3 + vc3) / vc3) / 3I haven't double checked this but thought I'd post it before I forget!