Fixed an error in the calculation of flux uncertainties for the output arrays. Covariances are now only calculated for TC regions (23,23) and not for the full state vector.