如何在FloPy中添加河流？

Question

我正在处理一个稳态问题，我有一个光栅，其中包含所有河流的头部数据值。我必须将其添加到我的FloPy代码（MODFLOW 6）中，因此我看到了RIV包和DRN，希望找到一个变量，我可以使用它添加光栅中的数据（使用GDAL将数据值读取为数组），但还没有答案。目前还没有找到任何示例作为指导

我一直在考虑将其添加为IC包中的STRT，在第一层中为河流添加STRT（我为另一层的头部添加了STRT）。但我不太确定，可以吗？请问这是正确的方法吗？有什么建议吗

另外，应该用什么值来表示层中没有头

很明显，我是一个新手：p事先谢谢你：）