SFR package load issue

[John1205]

I am loading a groundwater model built with MODFLOW-NWT in Santa Rose Plain. The files can be downloaded here. The packages (SFR, HFB ) were not loaded successfully. When I load the package, the model shows the following error. Can anyone help me with this issue? Thanks.

sfr = flopy.modflow.ModflowSfr2.load(f, model = ml)
loading sfr2 package file...
Traceback (most recent call last):

File "", line 1, in
sfr = flopy.modflow.ModflowSfr2.load(f, model = ml)

File "C:\Users\liqia\AppData\Local\Continuum\anaconda3\lib\site-packages\flopy\modflow\mfsfr2.py", line 885, in load
_get_dataset(f.readline(), [0.0] * 8))

File "C:\Users\liqia\AppData\Local\Continuum\anaconda3\lib\site-packages\flopy\modflow\mfsfr2.py", line 2781, in _get_dataset
dataset[i] = n

IndexError: list assignment index out of range

[John1205]

I am here to upload the input files for this model. This model is GSFlOW, which includes MODFLOW for groundwater simulation. Thanks for your help!
elev33_icalc2_K2.zip

[andyrich]

This seems to be a bug. When reachinput is true, and isfropt ==3, and icalc =2, there is an issue with reading items 6d. The mrsfr2.py at line 881 is expecting [ XCPT1 XCPT2 ... XCPT8
ZCPT1 ZCPT2 ... ZCPT8 ]. But it should be expecting the following:
1 2 3 0 0 0 0 0 0.041 0.111

here is a summary of the sfr package referenced by @John1205:
OPTIONS
REACHINPUT
END
2902 324 0 0 128390 0.0001 119 0 3 10 1 30 0 4 0.75 91.54
1 39 138 1 1 525 1678.155 0.238999 2 0.007086 0.3 0.0202 3.5 0.007086
...
...
324 2 0
1 2 3 0 0 0 0 0 0.041 0.111 #this line causes the problem
0 3.64 7.28 10.92 14.55 18.19 21.83 25.47 #code is currently written to expect this
2.55 1.02 0.76 0 0.25 0.76 1.02 2.55
2 2 3 0 0 0 0 0 0.041 0.111
0 3.96 7.92 11.88 15.83 19.79 23.75 27.71
2.77 1.11 0.83 0 0.28 0.83 1.11 2.77
3 2 5 0 0 0 0 0 0.041 0.111

@aleaf-usgs

[John1205]

Hello Andrew, Thanks for identifying the issue for me. I am the freshman to Flopy. Do you know how to fix this issue? Thanks, John

…

On Tue, Oct 15, 2019 at 11:03 AM andyrich ***@***.***> wrote: This seems to be a bug. When reachinput is true, and isfropt ==3, and icalc =2, there is an issue with reading items 6d. The mrsfr2.py at line 881 is expecting [ XCPT1 XCPT2 ... XCPT8 ZCPT1 ZCPT2 ... ZCPT8 ]. But it should be expecting the following: 1 2 3 0 0 0 0 0 0.041 0.111 here is a summary of the sfr package referenced by @John1205 <https://github.com/John1205>: OPTIONS REACHINPUT END 2902 324 0 0 128390 0.0001 119 0 3 10 1 30 0 4 0.75 91.54 1 39 138 1 1 525 1678.155 0.238999 2 0.007086 0.3 0.0202 3.5 0.007086 ... ... 324 2 0 1 2 3 0 0 0 0 0 0.041 0.111 #this line causes the problem 0 3.64 7.28 10.92 14.55 18.19 21.83 25.47 #code is currently written to expect this 2.55 1.02 0.76 0 0.25 0.76 1.02 2.55 2 2 3 0 0 0 0 0 0.041 0.111 0 3.96 7.92 11.88 15.83 19.79 23.75 27.71 2.77 1.11 0.83 0 0.28 0.83 1.11 2.77 3 2 5 0 0 0 0 0 0.041 0.111 @aleaf-usgs <https://github.com/aleaf-usgs> — You are receiving this because you were mentioned. Reply to this email directly, view it on GitHub <#671?email_source=notifications&email_token=ABKEBATX52D3J2YDM3IVRULQOYAYTA5CNFSM4I7RZY22YY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOEBJV77Y#issuecomment-542334975>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/ABKEBARBV6TJ2EABM4G225LQOYAYTANCNFSM4I7RZY2Q> .

-- Qiang Li, Ph.D. Postdoctoral Fellow Department of Civil Engineering, University of Victoria, BC, Canada

[aleaf]

Hi @John1205, I'll take a look at this as soon as I can get to it.

[John1205]

Thanks @aleaf . This issue has been solved.

[fahadkhan469]

Hi @John1205

I am trying the same thing. Could you please share with me the earlier lines of your script defining the variable 'f' and 'ml'? My idea is to load the .sfr file and then modify some parameters.

Thanks!

Fahad

[John1205]

@fahadkhan469 You can download the files at https://ca.water.usgs.gov/sustainable-groundwater-management/california-groundwater-modeling.html . In this link, only download INPUT files for your purpose. You could found MODFLOW files in the folder "GSFLOW COUPLED MODEL". I only changed the code that developers corrected and the issue was solved. I hope this helps. When you have input files, you can extract features of streamS and write your own input files.

[fahadkhan469]

Hi @John1205

Thanks for replying quick. I may have mistranslated my intentions, actually I am running the modflow-nwt model in a field site at Ethiopia, which is part of my PhD research. I have all the files I need. I just needed to know in your script ("sfr = flopy.modflow.ModflowSfr2.load(f, model = ml)") what is "f" and "ml".

Is "f" the .sfr filename? And what about "ml"? How did you load the "ml" variable in python? This is it, thanks for your response.

Fahad

[John1205]

hi @fahadkhan469 . I am using the model built by USGS, which I mentioned earlier. "f " is the .sfr name. ml is the loaded MODFLOW by using . NAM file. Please see below.

ml = flopy.modflow.Modflow.load('SRPgsf_newHfb_0530132.nam',
verbose=True, check= False,
exe_name=path2mf, version = 'mfnwt')