10.10 翻译设置

转换步骤允许您将模型转换为适用引擎将运行的输入文件。 每个引擎还有一些设置，允许您调整模型的运行方式（例如解算器设置、替代初始头、输出控制和单独的包设置）。

平移步骤由两个垂直面板组成。 左侧面板包括一个树结构，其中列出了在“选择引擎”步骤中选择的引擎以及每个适用引擎的关联设置。 右侧面板显示翻译树上所选节点的设置（和控件）。 文档本节中的每个子主题描述了翻译树上每个节点的可用设置。

Tools

10.10.1.1 General

10.10.1.2 MODFLOW Settings

10.10.1.3 LGR settings

10.10.1.4 Solvers

10.10.1.4.1  PCG

PCG 求解器（MODFLOW-SURFACT 中的 PCG2）使用预条件共轭梯度法来求解模型产生的联立方程。 可以模拟线性和非线性流动条件。 PCG 包括两个预处理选项：修改的不完全 Cholesky 预处理，这在标量计算机上非常有效； 多项式预处理需要较少的计算机存储空间，并且经过计算机特定的修改，在向量计算机上效率最高。 使用水头变化和残差标准来确定求解器的收敛性。 使用皮卡德迭代来解决非线性问题。 PCG 包在 USGS 水资源调查报告90-    4048（Hill，1997）中进行了描述。
PCG 求解器采用两层方法一次性求解，即内部迭代和外部迭代。 外部迭代用于改变求解方法中的预条件参数矩阵。 外部迭代是在预先调节的矩阵集中更新流动系统的水文地质参数（即透射率、饱和厚度、储存性）的地方。 内部迭代持续进行，直到执行了用户定义的最大内部迭代次数，或者满足最终的收敛标准。 外部迭代继续进行，直到更新后第一次内部迭代满足最终收敛标准。
以下是 PCG 方法的求解器参数说明：

Acceleration index  加速指数                

Linearization Method Index (INEWT)：线性化方法索引 (INEWT)

Backtracking Factor (BFACT)：回溯因子 (BFACT)

Residual Reduction Factor(RESRED)：残留减少因子(RESRED)

Aquifer Geometry-based Head Dampening Index (IDMPBOT)：基于含水层几何形状的水头阻尼指数（IDMPBOT）

Inactive cell inclusion index (IBONO)：失活细胞包含指数（IBONO）
Head change criterion (HCLDSE)：水头变化标准(HCLDSE)：
Convergence output control (MUTPCG)：汇聚输出控制（MUTPCG）
Maxt. number of orthogonolizations (MIORTH)：最大 正交化次数(MIORTH)

Residual criterion (RCLOSE)：残差准则（RCLOSE）

• 最大外部迭代次数（MXITER）：[默认值 = 25] 此参数提供要执行的外部迭代次数的上限。 仅当事先未达到收敛解时才会使用最大迭代次数。 对于大多数问题来说，二十五次迭代就足够了。 然而，如果达到最大外部迭代次数并且未达到适当的质量平衡误差，则应增加该值。

•  最大内部迭代次数（ITER1）：[默认值 = 10] 此参数提供要执行的内部迭代次数的上限。 仅当事先未达到“外部”迭代中当前矩阵集的收敛解时，才会使用此迭代次数。 对于大多数问题来说，十次内部迭代应该足够了。 超过十次迭代通常不会改进解决方案，因为当返回到外部迭代时，解决方案会再次更新。

• 收敛的水头变化准则（HCOLSE）：[默认值 = 0.01] 每次外部迭代完成后，求解器会检查每个单元的解的最大变化。 如果解的最大变化低于设定的收敛容差（此处以工作单位英尺或米为单位设置），则解已收敛并且求解器停止，否则将开始新的外部迭代。

精确到 0.01 [ft.] 的解决方案 或 m] 通常足以解决大多数问题，除非建模域的最大水头差小于一英尺或米。 如果未实现适当的质量平衡，并且内部和外部迭代的数量在上面声明的最大值范围内，则可以将该值减小一个数量级，例如 0.001。

•    Residual Criterion for Convergence（RESRED）：[默认值 = 0.01] 水头变化准则用于判断求解器整体收敛性，而残差准则则用于判断求解器内部迭代的收敛性。 如果所有节点残差的最大绝对值小于此处指定的容差（单位为 L3/T），则求解器将继续进行下一次外迭代。
 

如果您注意到所有外部迭代只执行了几次内部迭代，并且未实现适当的质量平衡，则“残差准则”值可能会减少一个或多个数量级。
注意：残差标准与单位相关。 如果长度单位是英尺或米并且时间单位是秒，则默认值 0.01 通常就足够了。 如果您的时间单位不是秒，则应将默认值乘以时间单位中的秒数（例如，如果您的时间单位是天，则使用残差标准 0.01 * 86400）。

·    阻尼系数：[默认值 = 1] 该系数允许用户减少（抑制）每次连续外部迭代期间计算的水头变化。 对于大多数“适当的”和物理上真实的地下水流问题，阻尼系数为 1 是合适的。 该参数可用于使非收敛（振荡或发散）求解过程更加稳定，从而获得解。 这是通过将阻尼系数减小到 0 到 1 之间的值（很少< 0.6）来完成的。 该参数类似于其他求解器中使用的“加速度参数”。

·    打印输出间隔：[默认=10] 打印输出间隔是将解的最大水头变化（和残差）写入列表（.LST）文件之前的迭代次数。

·    如果“预处理方法”设置为“Cholesky”，则可以设置“松弛”参数。 尽管默认值为1，但在某些情况下，0.97-0.99 的值可能会减少收敛所需的迭代次数。

10.10.1.4.2 GMG

GMG 求解器基于预条件共轭梯度算法，由 USGS 开发，用于求解基于有限差分的流动模型。 与 AMG 不同，GMG 中的预处理基于称为几何多重网格的求解器方法。 经证明，与使用相当内存量的其他求解器相比，GMG 求解器可以大大缩短模型运行时间。 有关 GMG 求解器的详细信息，包括与 AMG 求解器的比较，可以在 GMG 线性方程求解器包 PDF 文档中找到。下面使用 GMG 线性方程求解器包 PDF 文档的摘录描述了几何多重网格求解器的求解器参数：

•    · 最大限度。 外迭代 (MXITER)：外迭代的最大次数。 对于线性问题，MXITER 可以设置为 1。对于非线性问题，MXITER 需要更大，但很少超过 100。只有在事先没有达到收敛解时才会使用最大迭代次数。
•    · 最大限度。 内部迭代 (IITER)：每个线性解的 PCG 迭代的最大次数。 通常，值 100就足够了。 对于非线性问题指定较小的数字通常很有用，以防止内部迭代次数过多。 仅当事先未达到“外部”迭代中当前矩阵集的收敛解时，才会使用此迭代次数。

·    自适应阻尼控制（IADAMP）： IADAMP 是控制自适应阻尼的标志。 如果 IADAMP = 0，则分配给 DAMP 的值将用作恒定阻尼参数。 如果 IADAMP 0，则 DAMP 的值用于第一次非线性迭代。 阻尼参数根据水头变化自适应变化，并使用 Cooley 方法进行后续迭代。
·    水头变化准则 (HCLOSE)：每次外部迭代完成后，求解器都会检查每个单元解中的最大变化 ≠。 如果解的最大变化低于设定的收敛容差（此处以英尺或米的工作单位设
置），则解已收敛并且求解器停止，否则将开始新的外部迭代。 精确至 0.01 [ft.] 的解决方案 或 m] 通常足以解决大多数问题，除非整个建模域的最大水头变化小于 1 英尺或米。 如果未达到适当的质量平衡并且内部和外部迭代的次数在最大值范围内，则该值可以减少一个数量级。
·    残差准则（RCLOSE）：RCLOSE 是内迭代的残差收敛准则。 PCG 算法计算残差的 l2范数并将其与 RCLOSE 进行比较通常，RCLOSE 设置为与 HCLOSE 相同的值（见下
文）。 如果 RCLOSE 设置得太高，则由于未以足够的精度求解线性方程，可能需要额外的外部迭代。 另一方面，对于非线性问题，过于严格的 RCLOSE 设置可能会强制获得不必要的精确线性解。 这可以通过 IITER 参数或阻尼来缓解。
·    弛豫参数（RELAX）：RELAX参数可用于改善ILU预调节系统的光谱条件数。 RELAX的值应约为 1。 然而，松弛参数可能会导致因式分解失败。 如果发生这种情况，GMG求解器将报告装配错误，并且应尝试使用小于 1 的 RELAX 值。 仅当 ISC = 4 时，此项目才是只读的。
·    估计上限 (NPBOL)：IOUTGMG 是控制 GMG 求解器输出的标志。 IOUTGMG 的可能值及其含义如下： 如果 IOUTGMG = 0，则仅打印求解器输入。 如果 IOUTGMG = 1，则对于每个线性求解，PCG 迭代次数、阻尼参数值、残差的 l2-范数、头部变化的最大范数及其位置（列、行、层） ）被打印。 在时间/压力周期结束时，将打印所有时间/压力周期的 GMG 调用总数、PCG 迭代以及 PCG 迭代的运行总数。 如果 IOUTGMG = 2，则打印 PCG 迭代的收敛历史记录，显示残差的 l2 范数和每次迭代的收敛因子。 IOUTGMG = 3 与 IOUTGMG = 1 相同，只是输出发送到终端而不是 MF2K LIST 输出文件。 IOUTGMG = 4 与 IOUTGMG = 2 相同，只是输出发送到终端而不是 MF2K LIST 输出文件。
·    多重网格预处理器 (ISM)：ISM 是一个标志，用于控制多重网格预处理器中使用的平滑器类型。 ISM 的可能值及其含义如下： 如果 ISM = 0，则在多重网格预处理器中实现ILU(0) 平滑。 这种平滑需要在每个多重网格级别上有一个额外的向量来存储 ILU 分解中的主元。 如果 ISM = 1，则在多重网格预处理器中实施对称高斯赛德尔 (SGS) 平滑。 该平滑器不需要额外的存储空间； 如果在使用 ISM = 0 时超出或接近超出可用内存，用户可能需要使用此选项。使用 SGS 平滑不如 ILU 平滑稳健； 减少残差可能需要额外的迭代。 在极端情况下，求解器可能无法收敛，因为残差无法充分减少。
·    多重网格预处理器 (ISC) 中的半粗化控制：控制多重网格预处理器中的半粗化的标志。 ISC 的可能值及其含义如下： 如果ISC = 0，则行、列和层都被粗化。 如果 ISC = 1，则行和列将被粗化，但层不会。 如果 ISC = 2，则列和层会被粗化，但行不会。 如果 ISC = 3，则行和层会被粗化，但列不会。 如果 ISC = 4，则没有粗化。 通常，ISC 的值应为 0或 1。如果水力传导率垂直变化较大，则应使用值 1。 如果未实施粗化 (ISC = 4)，则GMG 求解器与 Hill (1990) 中描述的 PCG2 ILU(0) 求解器相当，并且使用最少的内存。
·    阻尼因子 (DAMP)：该因子允许用户减少（抑制）每次连续外部迭代期间计算的水头变化。 对于大多数“适当的”和物理上真实的地下水流问题，阻尼因子为 1 是合适的。该参数可用于使非收敛（振荡或发散）求解过程更加稳定，从而获得解。 这是通过将阻尼系数减小到 0 到 1 之间的值（很少< 0.6）来完成的。 该参数类似于其他求解器中使用的“加速度参数”

10.10.1.4.3 WHS

WHS 求解器使用双共轭梯度稳定 (Bi-CGSTAB) 加速例程，通过 Stone 不完全分解实现地下水流偏微分方程的预处理。 与所有迭代求解器一样，该求解器通过近似解迭代地求解一大组偏微分方程。 由于地下水流的矩阵方程最初是“病态的”，因此需要对这些矩阵进行有效的预处理才能获得有效的解决方案。
WHS 求解器采用两层方法一次性求解。 外部迭代用于改变求解方法中的因式分解参数矩阵。 外部迭代是在因式分解的矩阵集中更新流动系统的水文地质参数（即透射率、饱和厚度、储存性）的地方。不同级别的分解允许对这些矩阵进行不同的初始化，以提高求解效率和模型稳定性。 内部迭代用于迭代求解外部迭代中创建的矩阵。
WHS 方法的求解器参数描述如下：

·    外部（非线性）迭代的最大数量：[默认值 = 50] 此参数提供要执行的外部迭代数量的上限。 仅当事先未达到收敛解时才会使用最大迭代次数。 对于大多数问题来说，五十次迭代应该足够了。 然而，如果达到最大外部迭代次数并且未达到适当的质量平衡误差，则应增加该值。

·    最大内部迭代次数：[默认值 = 25] 此参数提供要执行的内部迭代次数的上限。 仅当事先未达到“外部”迭代中当前矩阵集的收敛解时，才会使用此迭代次数。 对于大多数问题来说，二十五次内部迭代应该足够了。 然而，如果所有外部迭代都使用最大内部迭代次数并且未达到适当的质量平衡误差，则可以增加该值。

·    收敛的头变化准则：[默认值 = 0.01] 每次外部迭代完成后，求解器都会检查每个单元的解的最大变化。 如果解的最大变化低于设定的收敛容差（此处以英尺或米的工作单位设置），则解已收敛并且求解器停止，否则将开始新的外部迭代。 精确至 0.01 [ft.] 的解决方案 或 m] 通常足以解决大多数问题，除非整个建模域的最大水头变化小于 1 英尺或米。 如果未达到适当的质量平衡并且内部和外部迭代的次数在最大值范围内，则该值可以减少一个数量级。

残差收敛准则：[默认= 0.01] 水头变化准则用于判断求解器整体收敛性，而残差准则则用于判断求解器内迭代的收敛性。 如果连续内部迭代的变化小于此处指定的容差（以英尺或米为工作单位），则求解器将继续进行下一个外部迭代。 0.001 的残差收敛标准应该适用于大多数问题。 但是，如果您注意到每次外部迭代仅执行几次内部迭代并且未实现适当的质量平衡，则可以将该参数值减小一个或多个数量级。

·    外部迭代的阻尼因子：[默认 = 1] 该因子允许用户减少（抑制）在每次连续的外部迭代期间计算的水头变化。 对于大多数“适当的”和物理上真实的地下水流问题，阻尼因子为 1 是合适的。 该参数可用于使非收敛（振荡或发散）求解过程更加稳定，从而获得解。 这是通过将阻尼系数减小到 0 到 1 之间的值（很少< 0.6）来完成的。 该参数类似于其他求解器中使用的“加速度参数”。
·    相对残差准则：[默认值 = 0] 该参数提供了另一种检查内迭代收敛性的方法。 该方法将最近内迭代的残差与初始内迭代的残差进行比较。 一旦最近的内迭代残差低于初始内迭代残差乘以相对残差标准，则当前外迭代完成，并且将开始新的外迭代。
·    因式分解级别：[默认 = 0] WHS 求解器提供两个可用的因式分解“级别”：0 和 1。级别 0 需要更多的外部迭代，但需要更少的内存。 级别 1 需要更少的外部迭代，但需要更多的内存。 虽然分解级别为 1时求解器的收敛需要较少的迭代，但运行求解器所需的内存会随着分解级别的增加而增加。 此外，每次迭代的工作量随着级别 1 分解而增加，使得总求解时间可能不小于使用级别 0 分解的求解时间。

10.10.1.4.4 SIP

强隐式过程，也称为 SIP，是一种通过迭代求解大型联立线性方程组的方法。 SIP 求解器的优点是它非常稳定并且通常会收敛到一个解，但通常非常慢。 它不如 PCG 方法快，但计算最终解决方案所需的内存更少。 因为每个方程涉及最多七个未知的水头值，并且由于整个网格中一组未知值从一个方程到下一个方程，所以必须在每个时间步同时求解整个网格的方程。 该软件包在 Visual MODFLOW Flex 介质随附的 MODFLOW 手册第 12 章中进行了描述，位于 Manual 文件夹中。
SIP 方法的求解器参数描述如下：

·    最大迭代次数：[默认值 = 200] 这是要执行的迭代次数的上限。 仅当事先未达到收敛解时才会考虑最大迭代次数。 对于大多数问题来说，200 次迭代应该足够了。 然而，如果达到最大迭代次数并且未达到适当的质量平衡误差，则应增加该值。
·    迭代参数个数：[默认=5] 描述地下水流系统的有限差分方程可表示为矩阵形式为[A]{h}={q}。 其中 [A] 是系数矩阵，{h} 是磁头数组，{q} 是通量数组。 迭代参数的数量表示将用于将初始系数矩阵[A]变换为可分别分解为两个下三角矩阵[L]和上三角矩阵[U]的相似矩阵的参数的数量。 默认值 5 通常就足够了。
·    加速因子：[默认= 1] 加速因子控制迭代之间水头变化的幅度。 加速因子必须为正。 大于 1的值将导致迭代之间较大的水头变化； 可能会更快地接近解决方案，但也可能更容易超出解决方案。 小于 1 的值将导致较小的头部变化，需要更多迭代才能达到解决方案。

收敛的水头变化准则：[默认值 = 0.01] 每次迭代完成后，求解器都会检查每个单元的解的最大变化。 如果解的最大变化低于设定的收敛容差（此处以英尺或米的工作单位设置），则解已收敛并且求解器停止，否则将开始新的迭代。 精确至 0.01 [ft.] 的解决方案 或 m] 通常足以解决大多数问题，除非整个建模域的最大水头变化小于一英尺或米。 如果没有达到适当的质量平衡并且没有达到最大迭代次数，则该值可以减少一个数量级。

打印输出间隔：[默认=10]打印输出间隔是将解的最大水头变化（和残差）写入列表（.LST）文件之前的迭代次数。

用户种子值：[默认 = 0.01] 有两个选项：用户可以输入种子，或者种子将在模拟开始时根据问题参数计算。 迭代参数“种子”用作确定 w 值序列的基础。 w 乘以方程右侧的每一项； 并且必须在连续迭代中循环使用一系列值以获得令人满意的收敛率。 系数矩阵对角线越强，种子的选择就越不重要。

10.10.1.4.5 SOR

切片连续过度松弛是一种迭代求解大型线性方程组的方法。 它在 SOR 包中通过将有限差分网格划分为垂直切片，并将节点方程分组为离散集合，每个集合对应于一个切片来实现。 在每次迭代中，依次处理这些方程组，从而为每个切片产生一组新的估计水头值。 当处理每个切片的方程时，它们首先用连续迭代之间计算水头的变化来表示。 然后，通过高斯消去法直接求解与切片相对应的方程组，将相邻切片的项视为已知量。 然后，为切片计算的水头变化值分别乘以加速度变量 T。结果被视为该切片迭代中水头变化的最终值。 按顺序对每个切片重复此过程，直到处理完三维数组中的所有切片，从而完成域迭代。 然后重复整个序列，直到连续迭代中计算的水头值之间的差异小于网格中所有节点处的所选标准。 SOR 包在 VMOD Flex 介质随附的 MODFLOW 参考手册的第 13 章中进行了详细描述，位于 Manual 文件夹中。

SOR 方法的求解器参数描述如下：

·    最大迭代次数：[默认值 = 50] 此参数提供要执行的迭代次数的上限。 仅当事先未达到收敛解时才会使用最大迭代次数。 50 次迭代对于大多数问题来说应该足够了。 然而，如果达到最大外部迭代次数并且未达到适当的质量平衡误差，则应增加该值。
·    加速因子：[默认值 = 1] 加速因子控制迭代之间水头变化的幅度。 加速因子必须为正。 大于1 的值将导致迭代之间较大的水头变化； 可能会更快地接近解决方案，但也可能更容易超出解决方案。 小于 1 的值将导致较小的水头变化，因此需要更多迭代才能达到解决方案。
·    收敛的水头变化准则：[默认值 = 0.01] 每次迭代完成后，求解器都会检查每个单元的解的最大变化。 如果解的最大变化低于设定的收敛容差（此处以英尺或米的工作单位设置），则解已收敛并且求解器停止，否则开始新的迭代。 精确至 0.01 [ft.] 的解决方案 或 m] 通常足以解决大多数问题，除非整个模型域的最大水头变化小于 1 英尺或米。 如果没有达到适当的质量平衡并且迭代次数小于最大值，则该值可以减少一个数量级。
·    打印输出间隔：[默认=10] 打印输出间隔是将解的最大水头变化（和残差）写入列表
（.LST）文件之前的迭代次数。

10.10.1.4.6 SAMG

Visual MODFLOW Flex 支持由弗劳恩霍夫算法和科学计算研究所 (FhG-SCAI) 开发的系统求解器代数多重网格方法 (SAMG) 包。

请注意：SAMG 求解器仅与 MODFLOW-2000、-2005、-LGR 和 -USG 流量引擎兼容。 SAMG 求解器在 Visual MODFLOW PRO 中以单核模式提供，在 Visual MODFLOW Premium 中以多核模式提供。

代数多重网格 (AMG) 包求解器可从弗劳恩霍夫算法和科学计算研究所 (FhG-SCAI) 获取，仅供研究用途。
SAMG 求解器包是一个完整的多级框架，旨在克服以前基于 AMG 的求解器的高内存要求，同时保持可扩展性和快速执行时间。 使用 Visual MODFLOW 生成的多个模型对 SAMG 求解器与 PCG2 求解器进行的测试表明，求解时间快了 2.4 到 11.3 倍。 对于大型网格（超过约 40,000 个单元）和/或高度可变的水力传导率场问题，SAMG 软件包比 MODFLOW 提供的其他求解器具有一些明显的优势。 与提到的其他迭代求解器相比，多重网格方法的优点是（1）多重网格求解器的有效性不依赖于初始水头分布，以及（2）对于许多感兴趣的问题，收敛速度与大小近似线性缩放 与收敛速度非线性增加的其他求解器不同（Demmel，1997）。
求解器设置窗口包含许多用户定义的求解器设置，这些设置可以影响 AMG 求解器的速度和有效性。

·    最大限度。 迭代次数 (MXITER)：[默认值 = 50] MXITER 是调用 AMG 例程以获得解的最大次数。 MXITER 永远不会小于 2，很少超过 50。当问题是线性时，MXITER 通常等于 2（所有层都受到限制，并且没有非线性边界条件；例如，蒸散、排水和河流包会产生非线性边界） 状况）。 对于非线性问题，MXITER一般为50以下； 然而，对于更严重的非线性问题，需要接近 50 甚至更大的值。
·    最大限度。 循环 (MXCYC)：[默认值 = 50] 对于每次调用求解器，AMG 都会循环执行一个或多个粗化和细化序列。 每次调用求解器时，求解器的最大 MXCYC 周期数受到限制。 对于大多数问题，每次迭代的收敛时间小于 50 个周期，因此 MXCYC 通常可以小于 50。然而，对于高度非线性问题，通过将求解器限制为较小的周期数可以实现更好的性能，并且 增加最大迭代次数 (MXITER)。 这可以防止求解器在这些高度非线性问题的早期迭代中不必要地找到非常准确的解决方案。
·    内部迭代的剩余收敛准则（RCLOSE）。 通常 RCLOSE 设置为与 HCLOSE 相同的值。如果 RCLOSE 设置得太高，则由于未以足够的精度求解线性方程，可能需要额外的外部迭代。 同样，对于非线性问题，过于严格的 RCLOSE 设置可能会强制获得不必要的精确线性解。 这可以通过 MXCYC 参数或阻尼来缓解。
·    注意：在新的SAMG封装中，RCLOSE和HCLOSE取代了BCLOSE
·    阻尼系数 (DAMP)：[默认值 = 1] 阻尼系数可用于限制从一次迭代到下一次迭代的水头变化，这通常在非常非线性的问题中很有用。 DAMP 使解变化缓慢，从而避免中间解的非线性效应引起的虚假偏差。 DAMP 值小于 1.0 会限制头部变化（松弛不足），而值大于 1.0会加速头部变化（过度松弛）。 对于线性问题，不需要阻尼，并且 DAMP 应设置为等于1.0。 对于非线性问题，可能需要限制水头变化 (DAMP < 1.0) 以实现收敛，并且 DAMP 值通常在 0.5 到 1.0 之间就足够了。
对于某些非线性问题，为每次迭代施加固定的 DAMP 值可能会阻碍收敛。 解决这种情况的一种方法是根据头部解决方案的进展情况来调整阻尼量。 AMG 套件提供两种自适应阻尼策略； (1) 经 Huyakorn 修改的 Cooley 法，以及 (2) 相对减少残差法。 这些方法在美国地质公开报告 01-177 中有详细描述。 DAMP 值为 -1 将使用第一种方法，DAMP值为 -2 将使用第二种方法。
·    最大限度。 阻尼因子 (DUP)：[默认值 = 1] 使用自适应阻尼策略时 DAMP 的上限。

·    最小。 阻尼系数 (DLOW)：[默认值 = 0.2] 使用自适应阻尼策略时 DAMP 的下限
·    水头变化收敛准则 (HCLOSE)，与之前的求解器所述类似
·    执行共轭梯度迭代 (ICG)：[默认 = 选中] 在某些情况下，由于少数误差分量在 AMG循环期间未减少，AMG 可能表现不佳。 共轭梯度求解器的几次迭代通常可以减少这些误差分量，从而有助于收敛（Cleary 等，2000）。 在这些情况下，参数 ICG 可以设置为 1，以便在每个多重网格循环结束时执行共轭梯度迭代。 激活此选项可以减少某些问题的执行时间，但也会增加求解器使用的内存量。

打印标志 (IOUTAMG) 框允许您在各种打印选项之间进行选择。

·    控制参数 [默认 = 2]
·    1 - 未使用设置阶段的重用
·    2 - 将使用设置阶段的重用（推荐）
·    3 - 将使用设置阶段的重用，并且将使用SSC

10.10.1.4.7 NWT

NWT 求解器参数如下图所示。 有关MODFLOW-NWT的详细信息可以找到：http://water.usgs.gov/ogw/modflow-nwt/

请注意：MODFLOW-NWT 包括上游加权 (UPW) 单元间电导包，作为 BCF、LPF 和 HUF 包提供的内部流量包的替代内部流量包。 当选择 MODFLOW-NWT 作为流程引擎时，将在转换过程中生成 UPW 属性包并在模型运行中使用。

Visual MODFLOW Flex 提供了最常用参数的选项（水头变化准则、通量变化准则、迭代次数等）以及选择求解器选项（自动求解器选项）。 MODFLOW-NWT 可以灵活地指定每个参数的值，或从预定义的求解器之一中进行选择
“计划”。 当您选择预定义方案之一时，求解器参数是固定的且无法修改（所有字段将不可见）。 如果选择“SPECIFIED”选项，则完整的求解器参数集将在 .NWT 文件中定义（使用默认值），并且您将能够在翻译后在 .NWT 文件中更改这些参数。 当您选择一个预定义方案（简单、中等或复杂）时，您将仅看到一组缩写的参数。

·    SPECIFIED 表示第 1 项和第 2 项列出的可选求解器输入值将由用户在 NWT 输入文件中指定。
·    SIMPLE 表示将定义默认求解器输入值，该值适用于近线性模型。 这将用于不包含非线性应力包的模型，以及受约束或由厚度足以将地下水位包含在单个层内的单个无约束层组成的模型。 （有关将用于此选项的求解器输入值，请参阅包文档的表 2。）
·    MODERATE 表示将定义适用于中等非线性模型的默认求解器输入值。 这将用于包含非线性应力包的模型以及包含一个或多个无约束层的模型。 这当“SIMPLE”选项不能成功收敛时，应使用“MODERATE”选项。 （有关将用于此选项的求解器输入值，请参阅包文档的表 2。）

COMPLEX 表示将定义默认求解器输入值，该值适用于高度非线性模型。 这将用于包含非线性应力包的模型，以及由代表复杂地质和地表水/地下水相互作用的一个或多个无约束层组成的模型。 当“MODERATE”选项不能成功收敛时，应使用“COMPLEX”选项。 （有关将用于此选项的求解器输入值，请参阅包文档的表 2。）
HEADTOL（长度单位）—是求解非线性问题的外部迭代之间的最大水头变化。 FLUXTOL（每次长度立方单位）— 是求解非线性问题的外部迭代之间的最大均方根通量差。

MAXITEROUT 是解决外部（非线性）问题所允许的最大迭代次数。
THICKFACT 是电池厚度（长度）的一部分，用于将存储和电导系数平滑调整为零。 LINMETH 是一个标志，用于确定将使用哪个矩阵求解器。

·    1 = GMRES will be used
·    2 = xMD will be used

IBOTAV 是一个标志，指示如果单元被脱水单元包围，是否将对相对于单元底部高度的地下水头进行修正。 勾选=将进行更正； 未选中 = 不会进行更正。 此输入变量是特定于问题的，并且应测试两个选项（选择 IBOTAV/未选择）。

请注意：Visual MODFLOW Flex 目前不支持对变量 PHIRAMP 的 .WEL 包进行修改

10.10.1.5 Recharge and Evapotranspiration

Recharge

充电区域分布可以应用于任何用户指定的模型层。 如果将补给分配给顶部网格层，并且顶层中的某些单元在模拟过程中变干，或者如果顶层中的某些单元被指定为无流动单元，则MODFLOW 程序允许补给 应用于模型中最上面最活跃（湿）层的网格单元。 充值设置显示在下面的充值选项窗口中，这些设置如下所述。

充值仅适用于顶部网格层：如果第 1 层中的任何网格单元是干燥的，或者如果它们被指定为无流量单元，则分配给这些网格单元的充值值将不会向下传递到底层活动（ 湿）网格单元。 在这种情况下，不活跃或干电池就像充电的不可渗透屏障。

·    充值应用于指定层：允许用户将充值值分配到任意指定模型层
·    充值应用于最上面的活动层：如果第 1 层中的任何网格单元是干的，或者如果它们被指定

为无流量单元，则分配给这些网格单元的充值值将向下传递到最上面的活动层（湿） ）网格单元位于同一垂直列的网格单元中。

注意：恒定的水头边界条件总是会拦截补给并防止更深的渗透。

Evapotranspiration

蒸发蒸腾分布可以应用于任何用户指定的模型层。 如果分配到顶部网格层，并且顶层中的某些单元在模拟过程中变干，或者如果顶层中的某些单元被指定为无流动单元，则 MODFLOW 程序允许将蒸散量应用于 模型中最上面最活跃（湿）层的网格单元。 蒸发蒸腾设置显示在下面的蒸发蒸腾选项窗口中，这些设置如下所述。

·    蒸散量仅应用于顶部网格层：如果第 1 层中的任何网格单元是干燥的，或者如果它们被指定为无流动单元格，则分配给这些网格单元的蒸散量值将不会被向下传递到底层活动（ 湿）网格单元。 在这种情况下，不活跃或干燥的细胞就像蒸散作用的不可渗透的屏障。
·    蒸散量应用于指定图层：它允许用户将蒸散量值分配给任何指定的模型图层。
·    蒸散量应用于最上面的活动层：如果第 1 层中的任何网格单元是干燥的，或者如果它们被指定为无流动单元，则分配给这些网格单元的蒸散量值将向下传递到最上面的活动层（湿） ）网格单元位于同一垂直列的网格单元中。

10.10.1.6 Lakes

10.10.1.7 Layer Types

10.10.1.8 Cell Rewetting

10.10.1.9 Initial Heads

10.10.1.10 Anisotropy

10.10.1.11 Output Control

10.10.1.12 Advanced Settings