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. &
  lpmax=1,phi=0.0,.75,1.1,2.1 &
 &debugs           lvacdsk=.f
  ldelta=.true. , wall=.f., infwal=.true. , checki=.false. fast= .t. &
  aw=100.,bw=0.  &
  a=-10.0, b=0.5   &
  rho1=0.26 rho2=0.74  rpof1=1.2 rpof2 = 1.6 rpof3 = 1.0
  alphap=2.5,delp=1.0  &
  scale=1.0 ,varmin=-0.677,varmax=-0.677,nvar=1 &
 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
Defines up-down symmetry.

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

Defines minimum value for truncation of Fourier series

Typically: -5

Defines maximum value for truncation of Fourier series

Typically: n * q-edge + 15

The toroidal mode number n = 1,2,3 ...
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
Sets the subdivision of the finite elements for integration. Allowed values are 2 and 4. Typically mdiv = 2
The number of divisions in the poloidal direction Preferably with a value = 2**k, eg. 64, 128, 256
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.
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
The ratio of specific heats, can be reset to look at the effect of compressibility. GAMMA=1.666667
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.
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.