A finite element adapted matricial programming language (FEMPL) for modelling variably saturated flow in porous media with Galerkin-type schemes.