问题描述
我正在使用 FloPy 加载现有的 MODFLOW-USG 模型。
load_model = flopy.modflow.Modflow.load('HTHModel',model_ws='model_ws',version='mfusg',exe_name='exe_name',verbose = True,check = False)
在加载LPF包的过程中,python显示hk和hani已经加载成功,然后报如下错误:
loading bas6 package file...
adding Package: BAS6
BAS6 package load...success
loading lpf package file...
loading IBCFCB,HDRY,NPLPF...
loading LAYTYP...
loading LAYAVG...
loading CHANI...
loading LAYVKA...
loading LAYWET...
loading hk layer 1...
loading hani layer 1...
D:\Anaconda\program\lib\site-packages\flopy\utils\util_array.py in parse_control_record(line,current_unit,dtype,ext_unit_dict,array_format)
3215 locat = int(line[0:10].strip())
ValueError: invalid literal for int() with base 10: '-877.0
我该如何解决这种问题。
顺便说一下,我通过使用 GMS 中的“保存原生文本副本”功能创建了这个模型。 Flopy可以正常读取LPF包中的其他内容,报错位置出现在读取[ANGLEX(NJAG)]数据的部分。
我将LFP文件与MODFLOW-USG的输入输出描述进行了对比,符合输入文件的格式要求。
我是 pyhton 和 flopy 的新手,这个问题让我很困惑。非常感谢你给我提供了一些参考资料,无论是关于Python、FloPy、MODFLOW-USG还是GMS。
解决方法
你能上传你的 lpf 文件吗?那我可以看看这个。但乍一看,-877.0 前面的“'”看起来很可疑——是在 lpf 文件中吗?