search

MODFLOW-USGS / modflow6

USGS Modular Hydrologic Model
https://modflow6.readthedocs.io/
Other
253 stars 116 forks source link

Specifying auxiliary variables in the Lake Package connectiondata table #432

Open aleaf opened 4 years ago

aleaf commented 4 years ago

From the IO doc and the ModflowGwflak constructor in flopy, it looks like auxilary variables aren't supported for the Lake Package connectiondata block. Would it be possible to implement this? @mnfienen-usgs and I were thinking this would be useful for parameter estimation. Specifically, we want to parameterize bedleak within two zones representing the littoral and profundal zones of the lake. It would be nice to flag these separately in another column of connectiondata, in case we want to start the two zones at the same values. Thanks for your consideration.

langevin-usgs commented 4 years ago

Hey @aleaf, auxiliary variables are already implemented by lake, and not by connection. Thinking about the structure of the code, it might be tricky to support two different types of auxiliary variables; some of the code we're working on now requires auxiliary variables by lake (such as concentration). Is there an alternative concept that would work?

ghost commented 4 years ago

Hey @aleaf and @langevin-usgs --- the issue is, currently bedleak is only specified in the connectiondata block. So, for parameterization, it would be ideal to be able to identify different values (like Andy said in this example profundal and littoral areas in the lake). So, Aux column identifying them would be the most straightforward way to facilitate that. Not sure how to reconcile that with the need for aux variables lake by lake as you say.....