search

MODFLOW-USGS / mt3d-usgs

MT3D-USGS Repository
23 stars 12 forks source link

Check mass balance in SFR network when negative values used for FLOW in SFR input #20

Closed emorway-usgs closed 6 years ago

emorway-usgs commented 6 years ago

Noticed that the From-To node listing wasn't as expected in the linker file. Note the 3rd line (in context) in some example lines from linker file below:

36 37 1 55692.22 1.208313
37 38 1 59427.54 1.283066
38 39 1 65349.79 1.370677
-999 40 0 -4.108000 6.5532312E-02 40 41 1 502.3974 9.3378410E-02 41 42 1 2249.491 0.2373254

As currently constructed, the line looks like it might work, but need to verify. Based on the node order, it is suggesting flow is from the sink (-999) to the first reach of iseg 4 (the 40th item in the REACHINPUT list), but because the flow is negative the order is actually reversed. Just want to make sure this is OK.

emorway-usgs commented 6 years ago

Followed up on this awhile back. Turned out to be related to issue #21 . MT3D-USGS was already written to account for negative flows, they weren't being passed correctly by MF-NWT. Fix was provided for NZ model awhile back and all seems to be working OK