Geometry

The code works in a flux coordinate system as a fixed boundary equilibrium solver. The outer boundary shape may be specified as:

• Conforming to the numerical expressions

  

  where

  • xzero - Major radius, Geometric center of the plasma

  • arad - minor radius, half width at

  • ellipt - Ellipticity

  • delta - triangularity- shift of from

  • - indentation parameter ( range from )

  • - controls of bean tip
    • crescent shape
    • rounded tip

• As an array of Fourier coefficients,

  

  

  The coefficients and are obtained from a DATA file:

  
        call link('unit23=(boundy.open)//') 
        read(23,2000)nbdyf
  2000  format(i3)
        read(23,2100)(rb(i),yb(i),i=1,nbdyf)
  2100  format(2e20.7)

  For this option, the only values in INEQU that are still relevant are

  • - switch to invoke this option, and

  • - the initial guess for the location of the magnetic axis. Other parameters, such as ...are deduced from the implied boundary shape and the input file values are ignored