Fixed index error in calc_tau and added mean intensity stuff.

For rays with mu > 0, tau was *decreasing* along the ray, which allowed
Delta_tau() to be *negative* (which it should never be), which made lots
of exponentials explode.