The also work in variably saturated flow models, but: Discrete features connecting along the well pipe are switched off above the water table (at p<0) and pumping well nodes (placed at the botttom of the multi-layer well) are also switched off if they are on dry nodes (p<0).