Modelling the evolution of a comet subsurface: implications for 67P/Churyumov–Gerasimenko