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