 Magnetohydrodynamic mass loading

When $\rho {\mathbf{u}}_p=k{\mathbf{B}}_p$, the flux-surface function $k$ is the mass flux per unit poloidal magnetic flux.