desc.vmec.VMECIO.compute_coord_surfaces

classmethod VMECIO.compute_coord_surfaces(equil, vmec_data, Nr=10, Nt=8, Nz=None, **kwargs)[source]

Compute points on surfaces of constant rho, vartheta for both DESC and VMEC.

Parameters:

  • equil (Equilibrium) – desc equilibrium to compare

  • vmec_data (str or path-like or dict) – path to VMEC output file, or dictionary of vmec outputs

  • Nr (int, optional) – number of rho contours

  • Nt (int, optional) – number of vartheta contours

  • Nz (int, optional) – Number of zeta planes to compare. If None, use 1 plane for axisymmetric cases or 6 for non-axisymmetric.

Returns:

coords (dict of ndarray) – dictionary of coordinate arrays with keys Xy_code where X is R or Z, y is r for rho contours, or v for vartheta contours, and code is vmec or desc