How to find and plot the impurity cooling rate (carbon, oxygen, and iron) in a tokamak plasma. In fact, I want to compute the plasma emissivity profile, which is a function of plasma density, temperature, and impurity cooling rate.
Emissivity= Ne*Te*L(Te)
If anybody knows the formula or any program, please share it with me.