Analytic functions with blackbody

See #3077 
You can use [`numpy.nan_to_num`](http://docs.scipy.org/doc/numpy/reference/generated/numpy.nan_to_num.html) for this--doesn't require the overhead of a mask or anything.