Version v2.2.1bn18 (Jan. 2014)

New atomic data for K shell absorption by neutral and once- and twice- ionized stages of Ne and Mg from Gorczyca.

Fix errors in the routine binemis which puts out binned emission lines. These errors led to spurious features in models with very high spectral resolution.

Change to the value of the constant hc used in conversion from ev to A and back, to reflect more accurate values for constants. Old value was 12398.54, new value is 12398.41. Also change to Rydberg constant; old value was 13.598, new value is 13.605. Adoption of consistent value of proton thermal speed as 1.29e+6 cm/s at 10$^4$K.

Adoption of routine which calculates photoionization integrals (phint53) which uses interpolation and smoothing.

Inclusion of code for calculating aped rates for collisional excitation (not yet fully implemented).

Inclusion of Bryans rates for collisional ionization.

Added feature which allows an array of densities to be read in. This is described in the chapter on inputs. It requires that the 'radexp' input variable be set to a number more negative than -100. Then ordered pairs of (radius, density) are read in from a file called 'density.dat'. Reading continues until the end of the file is reached. The density and radius values override the values derived from the ordinary input parameters. But execution will stop if other ending criteria are satisfied, i.e. if the model column density exceeds the input value, or the electron fraction falls below the specified minimum. The code will stop with an error if the density.dat file does not exist, or if the radius values are not monotonically increasing.

Another new feature allows reading in of table spectra in units of log10(F$_\varepsilon$). This requires that the spectrum_units input parameter be set to 2.