Unable to load .NAM file using Flopy

[fahadkhan469]

Hi everyone,

I am facing a problem while loading modflow input files with Flopy. I am trying to load the .NAM file and also the .SFR file, and in both cases, I am facing the following errors. All my model input files are uploaded in this Google Drive folder.

1. When I try loading the .NAM file, I get the following error

my code,

ml = flopy.modflow.Modflow.load("Tana_GW_climate_4_5.nam", model_ws=model_ws, verbose=False,
                               check=True, exe_name="mfnwt") 

error I receive,

IndexError: list index out of range

2. When I try loading the .SFR file, I get the following error

my code,

m = flopy.modflow.Modflow()
f = 'C:/Users/fkk17001/Desktop/FKK_UCONN/Tana_GW_climate/runs/4_5/test.sfr'
sfr = flopy.modflow.ModflowSfr2.load(f, m, nper=1465)

error I receive,

ValueError: invalid literal for int() with base 10: 'OPTIONS'

The model files are all uploaded in the following google drive folder (https://drive.google.com/drive/folders/1Q2NcZHeyFRsrrzFx5weTRw_hmwRaHUfh?usp=sharing).

I would highly appreciate if any of you could kindly explore this and let me know what is going on. I have been able to read other models with Flopy (which has different sets of packages), but not this one. The version I am using is, Flopy v 3.3.4!

Fahad

[fahadkhan469]

Hi @jlarsen-usgs, do you have some time to look into this?

Fahad

[jlarsen-usgs]

@fahadkhan469

It looks like you're using the modflow-nwt options block within your SFR file. Because Modflow-2005 does not support this block, try setting version="mfnwt" like this when you load your model:

ml = flopy.modflow.Modflow.load("Tana_GW_climate_4_5.nam", model_ws=model_ws, verbose=False,
                               check=True, exe_name="mfnwt", version="mfnwt")

Alternatively if you are just loading the SFR file, set version="mfnwt" in your modflow object, like this.

ml = flopy.modflow.Modflow("my_model", version="mfnwt")
sfr = flopy.modflow.ModflowSfr.load(.....)

Let me know if this works for you.

[fahadkhan469]

Hi @jlarsen-usgs

No it does not work. Below is the warning and error copied from Jupyter (also, the screenshot is available in:
https://drive.google.com/file/d/1qCfp01r8ZCYOj0_Lz5PdbvwRgjYTe5Uh/view?usp=sharing)

Warning

C:\Users\fkk17001\AppData\Local\Continuum\anaconda3\lib\site-packages\ipykernel\ipkernel.py:287: DeprecationWarning: should_run_async will not call transform_cell automatically in the future. Please pass the result to transformed_cell argument and any exception that happen during thetransform in preprocessing_exc_tuple in IPython 7.17 and above.
and should_run_async(code)

Additional text

dataset 2 will not be preserved in the created hfb object.
dataset 3 will not be preserved in the created hfb object.

Error

IndexError Traceback (most recent call last)
in
2
3 ml = flopy.modflow.Modflow.load("Tana_GW_climate_4_5.nam", model_ws=model_ws, verbose=False, check=True,
----> 4 exe_name="mfnwt", version="mfnwt")

~\AppData\Local\Continuum\anaconda3\lib\site-packages\flopy\modflow\mf.py in load(cls, f, version, exe_name, verbose, model_ws, load_only, forgive, check)
917 item.filehandle,
918 ml,
--> 919 ext_unit_dict=ext_unit_dict,
920 )
921 files_successfully_loaded.append(item.filename)

~\AppData\Local\Continuum\anaconda3\lib\site-packages\flopy\modflow\mffhb.py in load(cls, f, model, nper, ext_unit_dict)
551 bdtime = []
552 for n in range(nbdtim):
--> 553 bdtime.append(float(raw[n]))
554
555 # Dataset 5 and 6

IndexError: list index out of range

[fahadkhan469]

I want to also mention that I have another model (with the UZF package), which loads fine with FLOPY and the same code. But that model also does not have the LAK package. Could it be an issue with any of the packages? If so, could you please suggest any change in existing codes?

[jlarsen-usgs]

@fahadkhan469

I misunderstood the issue when I responded the first time. It looks like there are a few issues with the Fhb loading code in FloPy. The main issue is that it expects datasets 4b, 5b, 6b, 7b, and 8b to each be on a only a single line (even though spaning these datasets over multiple liness is valid in MODFLOW). I'm working on a fix to this issue and will put a PR in either this afternoon or sometime tomorrow morning.

[christianlangevin]

@jlarsen-usgs, if you haven't noticed it already, perhaps the read1d() function is a quick and easy way to handle the values spanning multiple lines.

[jlarsen-usgs]

@fahadkhan469, @langevin-usgs

I notice that there is no write support for datasets 2, 3, 6, and 8. Are these something that we should be supporting, at least with very basic functionality (if ModflowFhb reads it, it writes it back out)?

[fahadkhan469]

Hi @jlarsen-usgs, once the issue is resolved, could you please tell me what I need to do on my end? Do I need to update the flopy package?

[jlarsen-usgs]

@fahadkhan469

Yes. Once this issue is resolved update your flopy to the most recent development version using:

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

[fahadkhan469]

Hi @jlarsen-usgs I modified the mffhb.py file from my local drive, based on the changes you made in the mffhb.py file from pull #1264. The model now loads fine! I am not sure if I would run into any other error with the SFR file, but for now I am good. Thanks a lot!