flux term for Dirichlet BC does not contain conductivity factor
The formulation of the Dirichlet BC residual mirrors the one of the skeleton terms, but the flux term lacks the conductivity factor
relCond. I would guess this cannot be on purpose.
Associated with this might be the time step instabilities encountered when applying Dirichlet BCs at the upper domain boundary (where conductivity is generally far lower than for full saturation).