force_balance
moment_kinetics.force_balance — Modulemoment_kinetics.force_balance.force_balance! — Methoduse the force balance equation d(nu)/dt + d(ppar + nuparupar)/dz = -(dens/2)dphi/dz + Rdensi*densn(uparn-upari) to update the parallel particle flux densupar for each species