Skip to content

The threed1 cross-section table reflects the stored half for lasym runs #18

Description

@CharlesCNorton

The cross-section table of the threed1 output (yinden, yellip, ytrian, yshift, eqfor.f90:579-588) scans itheta = 1, ntheta2 at the plane v and then, in the second pass of the icount loop (eqfor.f90:522-534), the reflected plane 2 pi - v with Z -> -Z to obtain the other half of the contour. That reconstruction is the stellarator-symmetric one. For lasym = T the stored surface already covers [0, 2 pi) and its second half is not the reflection of the first, so the ellipticity, triangularity, indentation and shift of an asymmetric cross-section are computed from a symmetrised contour: the lower half of every cross-section in the table is the upper half with Z negated.

The 2017 VMEC2000 8.52 source and PARVMEC carry the same loop. For lasym = T a single pass over itheta = 1, ntheta3 at the plane v gives the actual contour.

Activity

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

      Milestone

      No milestone

      Relationships

      None yet

      Development

      No branches or pull requests

      Issue actions