我正在使用 Flopy (3.3.5) 和 MF6 (6.2.2) 运输来模拟从工业现场向小溪的羽流运动。我正在使用具有 6 层流量和运输同时运行的结构化规则网格，并且使用 RIV 包在流模型中表示小溪，沿岸有排水管。使用新的 ATS 包，模拟效果很好！并产生合理的浓度结果。我遇到的困难是试图从各种预算和通量输出选项中将质量通量提取到代表小溪的边界条件中。我已经建立了一个观察包来记录边界条件单元与模型其余部分共享的每个“JA-FACE”的质量通量，这似乎有效，
• 我注意到该模型模拟了非常小的质量通量，而不是在几次迭代后远离羽流的零。当具有三位指数的数字被写入列表文件和 obs 包的 csv 输出时，数字中的“E”被删除（例如“-2.90933775-241”，而两位指数数字被写成“- 2.90933775E-41”），我尝试更改输出的位数，但行为似乎始终如一。我已经为观察包的输出文件编写了一个修复脚本，但我想这可能是为什么当我尝试使用 flopy listBudget 实用程序读取列表文件时得到一个包含所有“无”条目的对象的原因。
• 我已将 SinkSourceMixing (SSM) 包设置为打印和保存预算信息。当我查看写入列表文件的输出时（在文本编辑器中而不是在 python 中），使用 SSM 信息写入列表文件的单元格（l，r，c 元组）与压力中的原始单元格不相关原始流量模型包的周期数据（井、河流和排水包之间的不一致）
• 我已经加载了二进制预算文件，并且可以使用“cbb.get_data(totem=t,text=”SOURCE-SINK MIX”) 来选择 SSM 数据以获得适当的recarray。我一直在根据“discretization/structuredgrid”包“node=L nrow ”的“get_node”方法中的公式计算节点数ncol+R*ncol+C” 与 flopy/mf6 复杂模型示例笔记本中使用的不同。使用这个节点号，我可以选择在“节点”列中包含它的recarray 的行。当我这样做时，我注意到“node2”列中的所有相关条目似乎都不是我能找到的任何尝试从节点号中退出 LRC 的方法的相邻单元格。当我查看预算中的“FLOW-JA-FACE”数据时，我使用从观察包中获得的计算节点号（用 L，R 定义,C 元组）所以看起来不错，但是从“FLOW-JA-FACE”数据进入单元的所有通量的总和似乎与预算中 SSM 数据中的值相加。
我认为出于我的目的，我从 obs 包中进入边界节点的流量就足够了，我只想发布一个问题来记录这些观察结果，希望我可以：
- 了解它是否是由于我的模型设置中的某些原因，以及
- 如果它们是由于需要解决的代码中的问题引起的，请根据我的经验为 flopy/mf6 项目做出一些贡献。
感谢此软件包的所有贡献者 Kevin