Next: About this document Up: Introductory manual for the Previous: Input

INDEX

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



Next: About this document Up: Introductory manual for the Previous: Input


manickam@pppl.pppl.gov
Tue Nov 22 17:17:59 EST 1994