search

MODFLOW-USGS / mt3d-usgs

MT3D-USGS Repository
23 stars 12 forks source link

Multi-species not working in SFT #27

Closed emorway-usgs closed 6 years ago

emorway-usgs commented 6 years ago

SFT not properly accounting for multi-species. For example, compare the list files (*.m3d) from the two attached models, 1 and 2 species, respectively.

The following call in gnt1.f (line 39):

CALL xmdsolv(AMATSF, RHSSF, CNEWSF, hclosexmd, rrctol, IASF, 
&            JASF,NJASF,numactive, north, iter, iacl, ierr)

sends in numactive equal to the number of SFT cells. However, CNEWSF gets picked-up as x and numactive gets picked up as n. In xmdsolv.f, on line 1961, x is dimensioned as x(n), but this is a problem because CNEWSF is dimensioned as CNEWSF(1:650,1:2) when it gets sent in, and so the second species is not properly accounted for.

CrnkNic.zip

emorway-usgs commented 6 years ago

@jtwhite79 This issue should be resolved with commit d916dd4

Tested with multi-species version of the 'CrnkNic' test problem. Please confirm all looks OK with your real-world model and reopen the issue if something is still amiss.

Screen shot showing budget for second species in SFT:

sft_multi_spec