Function adas.run_adas310

def run_adas310(adf01=None, adf18=None, izimp=[0], amimp=[0.0], frimp=[0.0], te=[0.0], tp=[0.0], dens=[0.0], denp=[0.0], bmeng=[0.0], files=None, tref=0.0, dref=0.0, bref=0.0, iz0=1, iz1=1, bsim=1.0, ts=8617.0, wdil=0.0, cion=1.0, cpy=1.0, wdpi=1e-08, nip=2, intd=3, iprs=1, ilow=True, ionip=True, nionip=2, ilprs=1, ivdisp=True, nmin=1, nmax=110, nrep=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 15, 20, 30, 40, 50, 60, 70, 80, 90, 100, 110], bdens=10000000000.0, no_ion=False, no_cx=False)

Runs the adas310 beam stopping generation code.

Parameters

adf01 : str
full name of ADAS adf01 charge exchange cross section file
adf18 : str
full name of ADAS adf18 expansion file
izimp : int, array
impurity atomic number
amimp : float, array
impurity atomic mass number
frimp : float, array
impurity fraction - values between 0-1.0
te : float, array
requested electron temperatures (eV)
tp : float, array
requested ion temperatures (eV)
dens : float, array
requested electron densities (cm-3)
denp : float, array
requested ion densities (cm-3)
bmeng : float, array
requested beam energies (eV/amu)
files : dictionary
names of the output files
  'adf21'   - stopping coefficient
  'adf26'   - population dataset
  'crmat'   - collisional-radiative matrix
  'repmap'  - mapping between n-representations
  'logfile' - summary of the runtref : float
dref : float, optional
reference density (electron) (cm-3), defaults to middle value of te
bref : float, optional
reference beam energy (eV/amu), defaults to middle value of bmeng
iz0 : int, optional
nuclear charge of beam species (default =1)
iz1 : int, optional
beam species recombining ion charge, ie beam charge + 1 (default = 1)
bsim : float, optional
beam species isotope mass (default = 1.0)
ts : float, optional
radiation field temperature (default = 8617eV)
wdil : float, optional
external radiation field dilution factor (default = 0.0)
cion : float, optional
multiplier of ground level e-impact ionisation rate (default = 1.0)
cpy : float, optional
multiplier of ground level e-impact excitation rate (default = 1.0)
wdpi : float, optional
external radiation field dilution factor for photo-ionisation (default = 1.0e-8)
nip : int, optional
range of delta n for impact parameter cross sections (default = 2)
intd : int, optional
order of Maxwell quadrature for cross sections (default = 3)
iprs : int, optional
above nip, iprs=0 use van Regemorter or Percival-Richards for iprs=1
ilow : boolean, optional
if True access special low level data (default = True)
ionip : boolean, optional
if True include ion impact collisions (default = True)
nionip : int, optional
Delta n for ion impact excitation cross sections (default = 2)
ilprs : int, optional
for ilprs=0 Vainshtein for ion impact (default) but use Lodge-Percival-Richards when ilprs=1
ivdisp : boolean, optional
if True use beam energy in calculation of cross section (default = True)
nmin : int, optional
minimum representative shell (default = 1)
nmax : int, optional
maximunm representative shell (default = 110)
nrep : int, array, optional
representative n-shells (default as above)
bdens : float, optional
beam density (cm-3)
no_ion : boolean, optional
if True switch OFF the ion impact contribution (default is False)
no_cx : boolean, optional
if True switch OFF CX contribution (default is False)

Returns

code : int
success code

Notes

There are a lot of inputs, with most of the control parameters having sensible defaults. The choice of target elements and plasma conditions are set to 0.0 and must be supplied for the routine to work properly.

The output of the routine is the set of files.

Calls a fortran executable code and communicates via a bi-directional pipe connected to stdout.

References

ADAS manual description of adas310: http://www.adas.ac.uk/man/chap3-10.pdf

Version History

  • Martin O'Mullane, 16-04-2019
    • First version