From: VAXK::LAO
Date: 18-NOV-1993  11:17:24
Description: Kinetic EFIT                                            

The following are outlines of steps for preparing the input
file to run kinetic EFIT using  magnetic data, kinetic
profile data, SXR inversion information, and MSE.  A number
of these steps were documented by politzer@gav.gat.com (Pete Politzer).
Some of these steps are implemented in the command file
EFIT:DO12.COM on Vaxs.
 
1.  Running ENERGY code to prepare profile data for ONETWO in order
to get the beam profile.  As a first guess of magnetic topology, use
a G EQDSK from EFIT based on an equilibrium fit using magnetics and MSE.
See command file EFIT:DOENERGY.COM for ENERGY code run.
 
2. After running ONETWO, three tables are pulled out of QIKONE file into a
temporary file: the table following "values on rho grid", and the two
tables following "densities, temperatures, currents, and related
quantities".  This file will then be used in the next step to form
the total pressure for use in the fitting.
 
3. The first line of data in EFIT:COPY.DAT is put at the head of this file.
This contains npress (26), ndata (51), ndel (2), scale (1.0), sigmin (1000.0),
sigmap (0.1), sigmab (0.0), psimin (0.0), psimax (0.0), aion (2.), abim
(2.), aimp (16.).  These may be changed to suit for the following
procedure.
 
4. Read this file (with EFIT:MAKE12U.EXE).  This looks for input tables with
ndata items, and generates tables with npress items by picking every
ndel'th item.  From the input tables it picks up (ONETWO names) psi,
enbeam, p ions, i ions, (beam, primary and impurity densities), ti, we, wi,
and wbeam.  Note: this only works if all 51 (and only 51) values are in the
tables.  Also, stick to one primary and one impurity species.
The output is a namelist-like file headed by $kinetic containing:
 -- tables of normalized flux (0 -> -1) (names are rpress,
    sibeam (0 -> +1), and rion),
 -- the number of items in the tables (respectively npress,
    nbeam, and nion; all equal)
 -- total pressure table: pressr = (2/3)*10^6*(we+wi+wbeam)*scale
    sigma: sigpre = Max[sigmin, (sigmap*pressr)]  (with sigmab
    instead of sigmap at edge)
 -- beam pressure: pbeam = (2/3)*10^6*wbeam*scale
 -- beam ion density: dnbeam = 10^6*enbeam
 -- mass density: dmass = 10^6 * m_proton * (A_p*n_p + A_i*n_i + A_b*n_b)
 
5. Make a new K0 file by deleting the $kinetic line from the MAKE12 output
and appending the rest at the end of the $IN1 namelist in a copy of the old
K0 file.
 
6. Finally, the items in the file NAME.KIN2 are appended at the end of the
$IN1 namelist.  Note that some of these may supersede values that have been
set earlier in the K0 file.
 
7. If using SXR data, set ITEK=-4 or -5, IQPLOT=12.  Append namelist
   EFIT:INSXR.DAT to K0 file.  Update inversion channels KSXR0 and KSXR2
   for the side and top SXR channels.
 
8. Edit K) file mannually to check all input variables.  Update edge
   values of SIGPRE array if needed.