Text this: Numerical strategies for representing Richards' equation and its couplings in snowpack models