get_phi Subroutine

private subroutine get_phi(phi, dist, skip_fsa)

The 'phi' variable passed in is: sum_s int dv (J0 * g) This routine divides by the appropriate gamtot factor depending on if we have kinetic or adiabatic electrons, and also on whether we are using 'g' or 'h' as our distribution function that we are evolving. Note that this routine is only called in the Electrostatic, Fluxtube case.

Arrays Parameters Grids Geometry

Divide by sum_s (1 - Gamma_0) Z^2*n/T to get phi The kx = ky = 0.0 mode is not evolved by stella so make sure this term is set to zero. Now handle adiabatic electrons only if needed.

Arguments

Type	Intent	Optional	Attributes		Name
complex,	intent(inout),		dimension(:, :, -nzgrid:, :)	::	phi
character(len=*),	intent(in)			::	dist
logical,	intent(in),	optional		::	skip_fsa