MF6 SFR perioddata Input Format?

[jklakovich-INTERA]

Hello,

I tried to input what I though to be "stress period data" to the perioddata parameter of flopy.mf6.ModflowGwfsfr - that is, a dict of lists of formatted strings. Here is the python input and output for perioddata (just one stress period here):

Input:

Output:

Here is the error that I get:

My format is wrong but I am having trouble deciphering from the documentation
(https://flopy.readthedocs.io/en/3.3.5/source/flopy.mf6.modflow.mfgwfsfr.html?highlight=sfr)
what is the correct input format. Could you please advise?

Thank you,
Jeff

[jklakovich-INTERA]

I also tried a dict of lists of tuples comprising an integer and a string - e.g., {0: [(102, 'status active'), ... ]} - and got the same error.

Thanks,
Jeff

[langevin-usgs]

How about something like this:

    perioddata = [
        [0, "status", "active"],
        [1, "status", "active"],
        [2, "status", "active"],
        [3, "status", "active"],
        [4, "status", "active"],
        [0, "diversion", 0, divflow],
        [1, "diversion", 0, divflow],
        [2, "diversion", 0, divflow2],
        [3, "diversion", 0, divflow3],
    ]

[jklakovich-INTERA]

Hi @langevin-usgs,

I input something similar -

And got the same error -

Also, perioddata must index the stress period, I think, so I'm not sure how a list of lists could have worked.

Jeff

[jklakovich-INTERA]

@langevin-usgs, I tried a variant of that - a dict of lists of lists:

Input:

Output:

Still got the same error....

[langevin-usgs]

There are a bunch of examples in this folder here: https://github.com/MODFLOW-USGS/modflow6/tree/develop/autotest

I was looking at the test_gwf_mvr01.py test, which has the following functioning code. Flopy should take a list and assume it is for stress period 0, or you can provide a dictionary.

[jklakovich-INTERA]

Hi @langevin-usgs, I tried to input the variable 'p_data' from autotest/test_gwf_buy_sfr01.py, shown below, into my perioddata variable:

... still get same error.

[jlarsen-usgs]

@jklakovich-INTERA

Here's what I've done in before to create perioddata for SFR.

    period_data = {
        i: [(0, "INFLOW", 50),] for i in range(12)
    }

Maybe the record needs to be in a tuple, even though FloPy should be able to handle any sort of iterable as input in these cases

[jklakovich-INTERA]

Hi @jlarsen-usgs,

I copied your code exactly -

And got the same error -

Maybe it's a bug?

I am using
FloPy 3.3.5
Microsoft Windows 10 Pro 10.0.19044 Build 19044

Strangely, I am not able to update FloPy using "conda update flopy" because I may not have the 'conda infrastructure' in my virtual environment (StackOverflow). I just discovered this while looking up my conda version.

Thanks for your help with this,

Jeff

[jlarsen-usgs]

@jklakovich-INTERA
Interesting, I haven't seen that before. Maybe try upgrading your flopy to the latest development version (3.3.7) using pip.

pip install -U https://github.com/modflowpy/flopy/zipball/develop

[jklakovich-INTERA]

Flopy 3.3.7 is handling the SFR perioddata with no issues.

Thanks for your help @jlarsen-usgs, @langevin-usgs .

Jeff