mf2005 and flopy: Using the Flow and Head Boundary package - negative head

[Gal Bartfeld]

Hey,

I am trying to build a model where one side of a coastal unconfined aquifer has a constant flux boundary condition, and the other a constant head of 0. The top and bottom edge have no flow.

I tried to use the flow and boundary head package, and set -1 to the left edge to force the head to stay 0. Head is 0, but at the right, where the constant flux is, I get a small negative head, regardless of the values of k and q (the flux).

I am not sure if I am using the FHB package correctly, I found little examples. I was wondering if someone can look at my code and see if they can tell what went wrong? I expected the flux to create a head gradient.

The code:

https://pastebin.com/JnbvurMr

Thank you very much.