All z/r values interpolated to lambdas, which means tiros and rees1 values can now be calculated. But first, must interpolate the s/r values for rees2, then they can all be swept through together...

all s/r values obtained

All three models now calculated for all ten energy regimes in all four pitch angle distributions at all 52 altitudes. I have the ion-pairs per electron that is (what a weird way of putting it) need to convert these into an erg per second per cm's (microwatt per metre squared) worth of energy flux.

I did have a goal when I started this, didn't I?

It is done. That's a hell of a lot of graphs to plot!