2

我想为我的井文件分配一个注释掉的井 ID 号（也对 ghb 单元执行相同的操作），但我找不到任何关于如何执行此操作的信息。

我写了一些东西来创建我自己的 ghb 文件，但是如果我尝试将它加载回我的软盘 mf 类，然后用 mf.write_input() 将它写出我的其他包，注释不会留下，它会被覆盖。

我知道在 mf.wrtie_input() 中我可以指定要写出的包，如果我拿走我之前制作的 ghb 文件（或井文件），那么原始文件不会被写入哪个好。

但我想知道是否有办法直接为每个包的压力周期数据添加注释，以便我可以将它全部包含在软盘类中。

谢谢

4

2 回答 2

2

我知道从现有包中转移注释的唯一方法是以读取模式打开文件，使用列数据创建一个 pandas DataFrame，然后从中构建一个新包。这是一个例子：

import os
import pandas as pd
import flopy.modflow as fpm
from collections import OrderedDict


pak_nam = 'drn'
mf_version = 'mfnwt'

# the model from which the DRN package will be copied
inmod = fpm.Modflow.load('10kTDS.nam',
                     model_ws=r'..\10kTDS',
                     version=mf_version, 
                     load_only=['drn'], 
                     check=False)

# the model where the new DRN package will be attached
mf = fpm.Modflow.load('ss2010.nam',
                      model_ws=os.path.join('..', 'ss2010'),
                      version=mf_version,
                      load_only=['dis', 'bas6'],
                      check=False)

# read the contents of the DRN package
with open(inmod.drn.fn_path, 'r') as f:
    lines = f.readlines()

# create pandas DataFrame
data = []
for line in lines[3:]:
    pieces = line.strip().split('#')
    t = pieces[0].strip().split()
    remark = pieces[-1]
    if t[0] == '-1':
        break
    else:
        data.append([int(t[0]),
                     int(t[1]),
                     int(t[2]),
                     float(t[3]),
                     float(t[4]),
                     '# ' + remark.strip()])
pak_df = pd.DataFrame(data,
                      columns=['k', 'i', 'j', 'alt_va', 'cond', 'remark'])
pak_df.loc[:, ['k', 'i', 'j']] -= 1

# specify data format
formats = OrderedDict([('k', '{:>10d}'.format),
                       ('i', '{:>10d}'.format),
                       ('j', '{:>10d}'.format),
                       ('alt_va', '{:>.2F}'.format),
                       ('cond', '{:>15.6E}'.format),
                       ('remark', '{>:50}'.format)])

# create new stress period data: for numpy record array use DataFrame.to_records()
pak_spd = {0: pak_df[list(formats.keys())].to_records(index=False)}

# attach DRN package to new model
pak = fpm.ModflowDrn(mf,
                     stress_period_data=pak_spd,
                     ipakcb=53,
                     options=['NOPRINT'],
                     filenames=os.path.join('..', 'ss2010', 'ss2010.{}'.format(pak_nam)),
                     dtype=pak_spd[0].dtype)

pak.write_file(check=False)
于 2018-11-13T14:06:19.860 回答
2

就像在这个例子中一样，您可以扩展默认值dtype以包含MfList实例将携带到写入的额外属性：

well_dtype = [('k', '<i8'), ('i', '<i8'), ('j', '<i8'),('flux', '<f4'), ('wel_id', object)]
stress_period_data = np.zeros((3), dtype=well_dtype)
wel = flopy.modflow.ModflowWel(m, stress_period_data=stress_period_data, dtype=well_dtype)

我不确定加载wel具有额外属性的现有包的简单方法 - 仅供参考

于 2018-11-11T17:51:18.537 回答