Bugs in mtlkt.py

[gdefilippis]

1. Lines 167 - 174 should be replaced by:

self.coldlak = [Util2d(model, (nlkinit,), np.float32, coldlak,
                             name='coldlak', locat=self.unit_number[0],
                             array_free_format=False)]

This is because you would get an error at line 190, saying that the method append cannot be used if self.coldlak is an Util2d object instead of a list.

2. Line 227 should be replaced by:

for coldlak in self.coldlak:
    f_lkt.write(coldlak.get_file_entry())

This is because in the .lkt file you would get coldlak values for the first species only

3. Line 446 should be replaced by:

for icomp in range(1, ncomp): 

[langevin-usgs]

Thanks @gdefilippis, we are taking a look.

[emorway-usgs]

@gdefilippis thanks for bringing this to our attention. I'm going to work on putting together a multi-species test problem with LKT. If you have a small example already put together that builds the model from scratch, and wouldn't mind sharing, I will incorporate that instead. Thanks, Eric

…

On Mon, Oct 1, 2018 at 11:15 AM gdefilippis ***@***.***> wrote: 1. Lines 167 - 174 should be replaced by: self.coldlak = [Util2d(model, (nlkinit,), np.float32, coldlak, name='coldlak', locat=self.unit_number[0], array_free_format=False)] This is because you would get an error at line 190, saying that the method append cannot be used if self.coldlak is an Util2d object instead of a list. 1. Line 227 should be replaced by: for coldlak in self.coldlak: f_lkt.write(coldlak.get_file_entry()) This is because in the .lkt file you would get coldlak values for the first species only 1. Line 446 should be replaced by: for icomp in range(1, ncomp): — You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub <#402>, or mute the thread <https://github.com/notifications/unsubscribe-auth/ADFi4AvVR_vlmgEuFfL9ZaXI_Gedt8azks5uglu6gaJpZM4XCnN9> .

[gdefilippis]

I have a model I used to integrate the LKT package in the QGIS-based FREEWAT interface (www.freewat.eu), which uses flopy for writing codes input files.
The flow model, including 3 lakes, is reasonable, meaning that it runs successfully and it is based on a real-world case study. The transport model has been set with totally unreasonable concentration values, because I used it just to test if the .lkt file was written correctly.
I hope you can benefit from this anyway.

Thank you,
Giovanna
lkt_test.zip

[emorway-usgs]

@gdefilippis If you pull the latest changes on the develop branch, your problem should be working with the changes made under pull request #412 Thanks for your patience. Let us know if something doesn't look quite right.

[jdhughes-usgs]

@emorway-usgs should this issue be closed?

[gdefilippis]

Dear all,
I am sorry for delay.
I tested this new version and it works fine.
Thank you,
Giovanna