  We study S-wave and P-wave bottomonium spectral functions at non-zero
temperature in 2+1 flavor QCD using the NRQCD formulation for bottom quarks. We
use a novel Bayesian approach to reconstruct the spectral functions and find
that $\chi_{b1}$ survives up to $T=249 {\rm MeV}$. We also study the effect of
different temporal discretizations of the NRQCD formalism on the bottomonium
correlation functions
