Hi All,
I've been reading through the FEM file format documentation and found something particularly interesting.
In the decoding/enconding bc type section (pages 13 and 14) one of the parameters (the third last) mention the following:
unsigned_h_bc_loc : 3; /*Save entry location for the changed boundary value :
...
...
4 - B.C. value is replaced by bc.v[n][5], zero flux */
Does it mean that if we use this option 4 (100 in binary), a zero flux constraint is applied when the boundary exceeds the constraint?
If so, that would be able to sort out the problem with well constraints when the hydraulic heads goes below the specified level, assigning zero flux instead of a constant head equating to the minimum constraint and thus injecting water in to the aquifer.
Thanks in advance,
Eduardo
