contour

pfsspy.plot.contour(phi, costheta, field, levels, ax=None)