我正在研究一个稳态问题，我有一个包含所有河流头部数据值的栅格。我必须将它添加到我的 FloPy 代码（MODFLOW 6）中，所以我看到了 RIV 包和 DRN，希望能找到我可以使用一个变量从我的栅格中添加数据（使用 GDAL 将数据值作为数组读取）但还没有答案。目前还没有找到任何示例作为指南。
我一直在考虑将它作为 STRT 添加到 IC 封装中，在第一层有河流的 STRT（我为另一层的头部做了这个）。但我不太确定，这样好吗？这是正确的方法吗，请问有什么建议吗？
另外，应该用什么值来表示层中没有正面？
显然，我是新手 :p 提前谢谢你:)。