Function adas.run_adas316

def run_adas316(adf25=None, adf26=None, adf12=None, adf40=None, npix=[300], wvmin=[100.0], wvmax=[9999.0], thresh=1e-12, tr_upper=[0], tr_lower=[0])

Prepare effective emission coefficient data for bundle-n charge exchange and output via adf12 and adf40 datasets.

Parameters

adf25 : str
driver filename
adf26 : str, optional
filename for detailed information on the reference case
adf12 : str, optional
filename for effective emission for each transition within all (input) wavelength ranges
adf12 : str, optional
filename for feature emission for input wavelength ranges
npix : int, array, optional
number of pixels in range, up to 5 ranges are allowed but defaults to one with 300 pixels
wvmin : float, array, optional
minimum wavelength (A) with default of 100.0A
wvmax : float, array, optional
maximum wavelength (A) with default of 9999.0A
thresh : float
emissivity threshold (default is 1e-12)
tr_upper : int, array
upper n for transitions in adf12 output file, this over-rides wmin/wmax choice
tr_lower : int, array
lower n for transitions in adf12 output file, this over-rides wmin/wmax choice

Returns

nothing - the output is aset of files written to disk
 

Notes

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

References

ADAS manual description of adas316: http://www.ac.uk/man/chap3-16.pdf ADAS manual description of adf25: http://www.adas.ac.uk/man/appxa-25.pdf

Version History

  • Martin O'Mullane, 08-06-2021
    • First version

Examples

Generate an adf12 file for C5+ and print the first line.

>>> import adas as adas
>>> adf25='/home/adas/adas/adf25/p316_a25/bn_cxs_c5_example.dat'
>>> adas.run_adas316(adf25=adf25, adf12='/tmp/c5_adf12.pass')
>>> with open('/tmp/c5_adf12.pass') as f: first_line = f.readline()
>>> first_line.strip()
'Ne+ 9     H + 1 (1)  / receiver, donor (donor state index)    / LPARMS /''