pleiades.mirror_equilibrium

pleiades.mirror_equilibrium(brb, gplas, p0, alpha1, alpha2, tol=1e-10, maxiter=100, relax=0.0, plotit=False)

Compute jxB=gradP equilibrium for given P as function of psi/psi_lim given by P(psi) = p0*(1-(psi/psi_lim)**alpha1)**alpha2