Branch divergenceSnowDrift: Mismatch between WINDEROSIONDEPOSITION and delta SWE

[EricKeenan]

I'm noticing some strange differences between WINDEROSIONDEPOSITION and delta SWE.

For example, after a few time steps, the time sum of WINDEROSIONDEPOSITION looks like this. Sensible right?

But then delta SWE looks like this.

Somehow I think the different variables in SnowpackInterface.cc are not being updated properly.

[EricKeenan]

Variables to look at include:

mns, winderosiondeposition, tmp_ErodedMass, and ErodedMass.

[EricKeenan]

In the driftingsnow branch winderosiondeposition is defined twice in calcExplicitSnowDrift. Is this relevant for my problem? E.g.

L1541

winderosiondeposition(ix, iy) = tmp_ErodedMass(ix, iy);

and L1665

winderosiondeposition(ix, iy) = dM(ix, iy);

[EricKeenan]

After discussion with @nwever, we determined that the mismatch between delta SWE and WINDEROSIONDEPOSITION can be explained by the MS_WIND variable. MS_WIND represents the mass of eroded snow (kg/m^2). This mass, while suspended above the snow surface is not included in the model SWE. This is why delta SWE has a much larger magnitude than WINDEROSIONDEPOSITION.