**ALPHA****exponent in pressure profile when****ALPHA2****exponent in pressure profile when****ARAD****minor radius of the plasma. Half width at****BEAN****indentation parameter. This value is approximately of the indentation defined as Chance et. al. PRL.****BETA****exponent in the pressure profile when****BETA2****exponent in the pressure profile when****DEITY****triangularity in units of the minor radius****DELP****increment in the pressure on axis -****DLAY(I)****layer width for flatteing the and profiles****ELLIPT****Ellipticity : max. height / minor radius****FACIMP****error criterion for convergence. When is greater than 1, then the criterion is . Thus is equivalent to . This criterion is ignored if the inner loop converges in one iteration.****FACM****factor ranging between 0 and 1, which allows the addition of the second term proportional to in the -profile so that the gradients can be set independently in different regions of the plasma.****GBETA****Primary exponent of the -function when****GBETA2****Exponent of -function for supplemental term****GBETA3****Exponent of -function for supplemental term****GBETA4****Exponent of -function for supplemental term****GP0****Factor defining the relative weight of the supplemental term in the -function. Is fixed externally. When then it is redefined internally to match the slope of the pressure term near the axis.****GP1****Factor defining the relative weight of a supplemental term in****IFUNC****Switch to choose the functional form of and/or****IGDISK****Switch to obtain an array for from the disk file will set the current from uses the value of from the disk file****IPDISK****Switch to obtain an array for the pressure profile from****IPRINT****switch to invoke output to the terminal , on each iteration. The output quantities are . Note that in the -solver mode is irrelevant****IRESET****Will cause the code to pause after iterations to allow for resetting the numerical relaxation parameters etc.****IQFUNC****Switch to choose between the -solver mode( = 2) and the -solver mode ( = 0)****IQSAVE****Switch to choose the method of defining the -profile. Only applies in the -solver mode.**- uses the functional form
- Expects an array of coefficients for Chebyshev polynomials in the file QPOLY
- Expects an array of and values in a file QDATA

**IVARY****Number of equilibria to generate. If , then it will use the -profile of the first equilibrium to create more equilibria. These are written into seperate disk files ( Max. )****LFOURY****Switch to use fourier coefficients from to determine the boundary shape****LRATS****Number of rational surfaces at which the and profiles will have the slopes reduced (Max. 3)****NIMAX****Number of iterations in the inner loop of the solver. Must be a multiple of 50****NPSI****Number of radial grid points. These will be spaced equally in . The number should be even and at least 4 less than the max. dimension****NPTS****Number of data points specified in the boundary shape****NTHE****Number of intervals. The number should be even and at least 4 less than the max. dimension****NUMIT****Number of outer iterations. Should be a multiple of 50****OMEGA****Under-relaxation parameter for the inner loop. It ranges between 0 and 1 in the -solver mode and between 1 and 2 in the -solver mode. In general the code is less sensitive to this and it may be set at its optimal value,****P0****Pressure at the magnetic axis in units of****P1****Coefficient of the auxiliary term in the pressure profile when . It is normalized to .****PPFAC****Parameter controlling the flattening of the -profile near the rational surface when .****Q0****Desired value of on axis in the -solver mode. Desired value of in the -solver mode.****QLIM****Coefficient in the -profile****QLIM2****Coefficient in the -profile****QLIM3****Coefficient in the -profile****QPFAC****Parameter to determine the degree of flatteing of the and profiles****QPOF****Exponent of the -profile****QPOF2****Exponent of the -profile****QPOF3****Exponent of the -profile****QPOF4****Exponent of the -profile****QRAT(i)****value at which the slope is reduced with****SHAPE****Shape coefficient which controls the of the bean tip****TSF****Under-relaxation parameter for the outer loop. The code is quite sensitive to this, and if convergence is not obtained , it should be reduced, sometimes to as low as****XZERO****Major radius of the geometric center. Midpoint at****XGUESS****Initial guess of the location of the magnetic axis. As the coordinate system is constructed about this point, it is crucial in obtaining convergence. In general at low choose a value slightly greater and at high a value slightly less than the expected magnetic axis location**

**
**

Tue Nov 22 17:17:59 EST 1994