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