Mathematical modelling of supraglacial meltwater production and drainage