babywade/pystark

Lineshape models for hydrogen emission from fusion plasmas (Stark, Zeeman and Doppler broadened)

pystark

Lineshape models for spectral emission from fusion plasma. Four line models are implemented. Not optimised!

Authorship:

• Joseph Allcock ([email protected] for contact)
• James Harrison

Prerequisites

Python 3 only.

• numpy
• scipy
• matplotlib
• pandas
• f2py (used to wrap fortran90 code in python. This comes with numpy, but if there are any problems, check you are able to run 'f2py' from the terminal. When i run locally on macOS, i use the GNU Compiler Collection for the Fortran backend, which can be installed easily using homebrew: 'brew install gcc')

Installation

In the top-level cloned pystark directory run terminal command:

pip install -e . --user

Then from python try:

import pystark
pystark.demo()

to see if it is all working.

Objects

• pystark.demo(): Run this to see what is and isn't working. Demonstrates how to use each of the lineshape models, produces an example plot.
• pystark.BalmerLineshape(): Main class, see docstring and demo() for usage.

Line models

• Voigt: Using Griem's empirical scaling for Stark profile FWHM, assuming a Lorentzian lineshape. Quick, but inaccurate in the wings.
• Rosato: interpolates Rosato tabulated profiles (Stark-Zeeman), convolved with Doppler profile. Narrow parameter range but fine tabulated grid and self-consistent treatment of Stark-Zeeman effects. paper.
• Stehle: interpolates Stehle tabulated profiles (Stark-Doppler), convolved with Zeeman split lines (simple, strong-Zeeman approximation). This is SLOW (~0.5 sec to evaluate) and the tabulated grid is coarse -- included for its wide input parameter range. paper
• Stehle param: B. Lomanowski's parameterised model of the tabulated Stehle data. Quick AND accurate, thanks Bart!