mf.change_model_ws(etx_model_ws, reset_external=True) not writting external for well package

[rosskush]

Howdy again, I have run into an issue when using the new MfUsg object.

In the 3.3.4 release this writes out a .ref file for each stress period for the well file however with 3.3.5 the well file looks like a normal well file with the data stored internally.

ext_model_ws = os.path.join('model_ws','freyberg.usg.ext')
mf.external_path = "."

if os.path.exists(ext_model_ws):
    shutil.rmtree(ext_model_ws)

# change dir and write
mf.change_model_ws(ext_model_ws, reset_external=True)
mf.write_input()

results in:

It seems to work fine for static packages like the bas, disu and lpf but it also seems to work for drn, and ghb

Here is a link to a working example: https://github.com/rosskush/usg_freyberg

I know there has been some changes to how flopy is handling USG so please let me know if I just missed something.

Thanks!

[cnicol-gwlogic]

@swfwmd - any chance you could have a look at this?

[jdhughes-dev]

@cnicol-gwlogic Any thoughts on how to proceed?

[cnicol-gwlogic]

@jdhughes-usgs I'll try take a quick look today. If I can get up a quick solution I will. I'm a bit pushed for time is all.

[jdhughes-dev]

@cnicol-gwlogic no rush. I was just working on clearing issues today. Thanks.

[cnicol-gwlogic]

@jdhughes-usgs @rosskush - fyi yep this will take a bit of work. The MfusgWel load and write stuff is a bit hokey - doesn't take advantage of MFList capability really, and I can see a few other issues with AUX variables. I'll put on my list but it will be a while I'd say. Cheers.

[cnicol-gwlogic]

@jdhughes-usgs @rosskush this should be sorted via PR #1551

[rosskush]

Awesome looking forward to trying it out!