Encorporated non matom estimators for ionization scheme into estimators....

In order to fix bug #43 we need ionization mode 7 to work properly in macro atom mode.

This change means that we now store values of average frequency and J for macro atom mode by calculating them in bf_estimators_increment, as macro mode never enters the routine radiation().

This means that Python now at least manages to get a model for the spectrum in a cell, but there are still some unsolved problems. This is the first step to getting QSO macro models working.

I'm also testing out the pull & fork model to see how it works!