A semi-classical field method for the equilibrium Bose gas and application to thermal vortices in two dimensions