Reduced Histogram calculations by saving prob to sfh_hist and ir_hist and expanding later
 Replaced x**2 in main loop with x*x.
 Reversed 2D arrays for more sequential memory access.
 Added -malign-double to x64 compiler options.
 Tried real*8 savlog10 = log(10.0) had slight errors.
 Replaced with real** savelog10 = log(10.0D0), no errors.
 Converted to constant...
 Replaced 10**a with exp(a*2.30258509299404568402D0)
 60% gain in speed with gfortran.