freegs.critical.find_separatrix

freegs.critical.find_separatrix(eq, opoint=None, xpoint=None, ntheta=20, psi=None, axis=None, psival=1.0)

Find the R, Z coordinates of the separatrix for equilbrium eq. Returns a tuple of (R, Z, R_X, Z_X), where R_X, Z_X are the coordinates of the X-point on the separatrix. Points are equally spaced in geometric poloidal angle.

If opoint, xpoint or psi are not given, they are calculated from eq

eq - Equilibrium object opoint - List of O-point tuples of (R, Z, psi) xpoint - List of X-point tuples of (R, Z, psi) ntheta - Number of points to find psi - Grid of psi on (R, Z) axis - A matplotlib axis object to plot points on