Skip to content

NPF and TVK -- allocation of K22 and K33 arrays #1000

@aprovost-usgs

Description

@aprovost-usgs

Describe the bug
If K22 or K33 (or both) is not specified explicitly in the NPF input, it is not allocated its own array but instead points to K11. This is done to save on storage and is fine unless TVK changes the value of that conductivity (K22 and/or K33). In that case, the value TVK specifies for that conductivity overwrites the value of K11 because that conductivity points to K11. In general terms, the fix would be to take the needs of TVK into account when (re)allocating storage for K22 and K33.

To Reproduce
Steps to reproduce the behavior:

  1. Start with the tvk01 test problem, in which K11 and K33 are explicitly specified in the NPF input. K22 is not and so points to K11.
  2. Modify the TVK input so that K11 (i.e., K) is set to a value other than 1.0000000 at one of the cells in period 1. For example, set K to 5.0000000 instead of 1.00000000 at cell (3, 3, 3). If rerun at this point, the heads and flows will obviously be somewhat different than in the original problem.
  3. Now add a line to the bottom of the period 1 input that sets K22 to 1.00000000. This will cause the value of K to be set back to the original value of 1.00000000. Rerunning the problem now gives the original heads and flows.
  4. Starting fresh with the tvk01 test problem, modify the NPF input to specify K22 explicitly to some value. This causes K22 to be allocated its own array, i.e., it does not point to K11. Also, add a line to the bottom of the period 1 input that sets K22 at cell (3, 3, 3) to a value other than 1.0000000, say, 5.0000000. Rerunning the problem still gives the original heads and flows because K22 does not overwrite K11 (i.e., K), and the value of K22 itself has no effect on the solution in this problem.

Expected behavior
The expected behavior when TVK changes K22 or K33 is that only that conductivity changes, without overwriting K11.

Screenshots
N/A

Desktop (please complete the following information):

  • OS: Windows 10 Enterprise
  • Version 10.0.19044 Build 19044

Metadata

Metadata

Assignees

No one assigned

    Labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions