A Physics-Constrained Neural Differential Equation Framework for Data-Driven Snowpack Simulation