  We report our preliminary progress in the calculation of the interquark
potential of bottomonium at non-zero temperature using the HAL QCD method. We
use NRQCD correlation functions of non-local mesonic S-wave states to obtain
the central potential as a function of temperature. These results have been
obtained using our anisotropic 2+1 flavour "Generation 2" FASTSUM ensembles.
