When running ve_run I still get divide by zero warnings emitted from the soil.uptake submodule, even though these cases have been explicitly handled.
I should go through this submodule systematic and switch uses of np.where to np.divide with a where argument provided (which stops the warning from being emitted)