LHM conversion: Water balance of wells

[JoerivanEngelen]

In my conversion of the LHM, I notice that the total water balance of wells differs a lot between the iMOD5 conversion and iMOD Python conversion.

In the iMOD5 conversion 1.4e6 m3/d is injected and 4.3e6 m3/d is extracted, whereas in iMOD Python's 8e6 m3/d is injected and 1.5e6 m3/d is extracted. I expect more to be extracted than injected, so iMOD5's seems to behave better.

I first thought this was because iMOD5 and MODFLOW6 use different sign conventions, but this doesn't appear to be the case:

From iMOD 5 manual:

Extraction rate, negative rates are extractions, positive rates injections

From MODFLOW6 manual:

Negative values of QWEL are used to indicate well discharge (pumping), whereas positive values of QWEL indicate
an injection well

Things to investigate:

• Are rates similarly distributed?
• Are the same wells filtered?
• Differences between text files?

[JoerivanEngelen]

I've verified with a visual check that for a few individual LayeredWells at least: the same indices + flux combination occur. So I can confirm there is no sign difference between iMOD Python and iMOD5. For these wells the multiplication factors were also correctly applied.