Fracton topological order at finite temperature