There are several ways to categorize this work. One way is to see it inserted in the GEOtop 3.0 project that aims to rebuild GEOtop physics and informatics. Another way is to see it inserted in the work to get better cryospheric parameterisations. A third way is to see it as a part of the process delineate in my professorship talk. More pragmatically, it represents a new implementation of the coupled Richards equation and energy budget by using sound numerics and OMS3.
No easy to read. It requires application.