将现有的mf2005模型加载到软盘中并更改参数

Question

我有一个在Processing Modflow gui中创建的现有MODFLOW2005模型。我想将此模型导入flopy，以便能够对模型参数进行敏感分析，我相信使用flopy可以更快地进行。

我可以使用以下方法加载现有的modflow模型：

ml = flopy.modflow.Modflow.load（“ modelnamw.nam”，model_ws = model_ws，verbose = True，check = False）

并且可以使用以下方式重命名模型以创建新的输出：

ml.name ='新模型'

ml.write_input（）

有没有办法让我保持整个模型不变，而只是更改水力传导率（hy）参数（将其余的bcf输入保持不变）？

谢谢

Answer 1

最简单的方法可能是创建模型的副本（通过更改model_ws或为其命名），然后使用修改后的参数创建新的BCF包。确保将所有未更改的参数传递给新的BCF软件包。

# get the BCF package
bcf = ml.get_package("BCF6")  

# new hy
new_hy = 2.

# don't forget to pass all the unchanged parameters from the old BCF
new_bcf = flopy.modflow.ModflowBcf(ml, laycon=bcf.laycon, hy=new_hy, vcont=bcf.vcont) 

new_bcf.write_file()  # write file
ml.run_model()  # run model with new BCF

也可以仅更改现有对象上的参数。为此，用新的bcf.hy对象替换现有的Util3d对象。注意：在这种情况下，它是一个Util3d，但对于其他参数，它可能是1D或2D。

# get the BCF package
bcf = ml.get_package("BCF6")

# create new util3d object
new_hy_util3d = flopy.utils.Util3d(ml, bcf.hy.array.shape, np.float32, new_hy, "hy")

# replace the old hy with the new object
bcf.hy= new_hy_util3d  

bcf.write_file()  # write file
ml.run_model()  # run model with new hy