Reduce code duplication in jacobian functions

[juliasloan25]

Is your feature request related to a problem? Please describe. There is a significant amount of duplicated code between the update_jacobian functions for RichardsModel and for EnergyHydrology. It would be good to pull out the shared parts of this code into a shared location and reuse it in both of these functions.