Analysis of interfacial error in saturated–unsaturated flow models