nd electron temperature in eV</span>
    <span class="c">!REAL(rp) 				:: ne</span>
    <span class="c">! Background electron density in 1/m^3</span>
    <span class="c">!INTEGER 				:: num_impurity_species</span>
    <span class="c">!REAL(rp), DIMENSION(10) 		:: Zo</span>
    <span class="c">! Full nuclear charge of each impurity: Z=1 for D, Z=10 for Ne</span>
    <span class="c">!REAL(rp), DIMENSION(10) 		:: Zj</span>
    <span class="c">! Atomic number of each impurity: Z=1 for D, Z=10 for Ne</span>
    