WIP: Save BAYESTAR integrand at inclination quadrature points
Save the BAYESTAR posterior values at the inclination values that are used for N=10th order Gauss-Legendre quadrature over cos(i).
Each HEALPix pixel gets 40 extra columns: PROBDENSITY_i, DISTMU_i, DISTSIGMA_i, and DISTNORM_i, for i=0 to N-1.