I am using the open-source Python package, pygfunction, to model a BTES system to meet heating demands in a district heating network. Apart from obtaining fluid temperature profiles inside the borehole and inlet/outlet temperatures, I am interested in investigating the development of temperature outside the bore field in the surrounding soil.