Fri Nov. 18, 1994

MODIN input deck for the PEST-1 stability code.

 TITLE card
 &gloctl lanal=.false. &
 &modes   lsymz=.t.  lreal=.true. limag=.true.
   llmin=-5,  llmax=20, n=1.0, m=200, mdiv=2, mth=64,lfunin=.false. &
 &cplots
  lpmax=1,phi=0.0,.75,1.1,2.1 &
 &debugs           lvacdsk=.f
  ldelta=.true. , wall=.f., infwal=.true. , checki=.false. fast= .t. &
 &vacdat
  aw=100.,bw=0.  &
 &shape
  a=-10.0, b=0.5   &
 &cprofl
 gamma=1.66667
  rho1=0.26 rho2=0.74  rpof1=1.2 rpof2 = 1.6 rpof3 = 1.0
  alphap=2.5,delp=1.0  &
 &vardat
  scale=1.0 ,varmin=-0.677,varmax=-0.677,nvar=1 &
 &eigdt
 alam=0.0, dtry=.1, nsteps=1, nitmax=10, epscon=1.0e-3 &

Index of INPUT variables

TITLE card
80 character ASCII string used to identify the run e.g. Analysis of TFTR Shot 76778 at 3.92 secs. Nov. 11, 1994 JM
LSYMZ
Defines up-down symmetry.

.true. for up-down symmetry .false. for non up-down symmetry

LLMIN
Defines minimum value for truncation of Fourier series

Typically: -5

LLMAX
Defines maximum value for truncation of Fourier series

Typically: n * q-edge + 15

n
The toroidal mode number n = 1,2,3 ...
M
The number of radial finite elements. This is related to the number of surfaces, NOSURF, as well as MDIV, and requires that nosurf = m * mdiv + 1
MDIV
Sets the subdivision of the finite elements for integration. Allowed values are 2 and 4. Typically mdiv = 2
MTH
The number of divisions in the poloidal direction Preferably with a value = 2**k, eg. 64, 128, 256
LDELTA
Switch to compute delta-W(PLASMA). Switch must be true atleast once to compute the plasma contribution. It may be turned off when studying different boundary/wall conditions, or if you are searching for an eigenvalue that has not converged.See ALAM below. .true. Computes delta-w plasma. Must be set true at least once .false. Expects to use previous calculation of delta-W plasma.
WALL
One of several variables which determines the boundary conditions. In hierarachy, these are:
  • WALL TRUE implies a perfectly conducting wall at the plasma edge
  • LVACDSK TRUE will read in the vacuum delta-w from a file FORT.36 pre-computed using Chances, VACUUM code. FALSE will compute the vacuum delta-w according to the following boundary conditions
  • INFWAL TRUE will use wall at infinity boundary conditions FALSE will use a closed wall whose shape is determined by the variables a,b,aw,bw. The key variable is a
    • a > 10.0 Set wall at infinity, equivalent to INFWAL=.TRUE>
    • -10.0 < a < 10 Set wall according to a shape formula using a,b,aw and bw
    • a <-10 or a= -10 Use a conforming wall at a distance b measured in units of minor radius. Note if a = -10 and the plasma is bean shaped, use a straight line segment on the indented section
GAMMA
The ratio of specific heats, can be reset to look at the effect of compressibility. GAMMA=1.666667
RHOs
Defines the density profile. rho = rho1 + rho2*(1-x**rpof2)**rpof1 + (1-rho1-rho2)*(1-x)**rpof3 Note that rho is normalized to be 1 on axis, and x refers to the Toroidal flux and not the Poloidal flux.
SCALE
The toroidal field scale factor, can be used to look at nearby equilibria with different q-axis or q-edge scale = 1.0 ==> No change scale = -qval ==> B-field is scaled to make q-axis = |qval| scale = +qval ==> B-field is scaled to make q-edge = qval
Eigenvalue Search
One of several variables used in determining the eigenvalue in the inverse iteration scheme. ALAM Initial guess for the eigenvalue. Note that if ALAM=0 then only the number of unstables eigenvalues is determined. If NSTEPS > 0 , then ALAM is reset to be ALAM - DTRY, NSTEPS times. NITMAX iterations are attempted for each value of ALAM, and EPSCON represents the convergence criteria for each vector element.

manickam@pppl.pppl.gov