**4.5.1 RSF4输入**
对Recharge Seepage Face(RSF4)包的输入来自于IUNIT(8)中指定的单元。RSF4需要的附加输入变量显示在阴影框中。
请注意,新的RSF4包包含了原始MODFLOW(McDonald和Harbaugh,1988)的Recharge(RCH1)包。
瞬时充水条件的输入通过两个文件进行。第一个是在IUNIT(8)上读取的充水渗流面包输入文件,它指示了向模型施加瞬时充水的存在以及充水区域的数量。第二个文件包含了每个区域的充水时间序列(默认扩展名.RTS)。该文件是无格式的,当应用瞬时充水时,内部打开。
1. 数据: 格式: NRCHOP I10 IRCHCB I10 IUNCNF I10 MXSEEP I10 MXZRCH I10 IPNDPOR I10 **为每个应力周期输入2到6项。
2. 数据: 格式: INRECH I10 INIRCH I10 INPNDEL I10 INSEEP I10 **仅在INRECH $ 0时输入第3项a。
3a. 数据: 实用模块: RECH(NCOL, NROW) U2DREL。请注意,U2DREL是读取二维数组(即项3a的RECH)的实用模块。
**仅在INRECH > 0且MXZRCH > 0时输入第3项b。表格化的区域充水值将添加到上述第3项a中读取的默认充水值中。**
3b. 数据: 实用模块: IZNRCH (NCOL, NROW) U2DINT **仅在充水选项(NRCHOP)等于2且INIRCH $ 0时输入第4项。**
4. 数据: 实用模块: IRCH(NCOL, NROW) U2DINT **仅在使用非约束选项(IUNCNF> 0)并且INPNDEL $ 0时输入第5项(5a和5b)。**
5a. 数据: 实用模块: PNDEL(NCOL, NROW) U2DREL **仅在IPNDPOR=1时输入第5b项,如果池塘孔隙度标志打开。**
5b. 数据: 实用模块: PNDOR(NCOL, NROW) U2DREL **仅在INSEEP > 0时输入第6项。**
6. 数据: 格式: 层 I10 行 I10 列 I10 海拔 F10.0(通常,输入项6包括每个渗流面边界单元的一条记录。如果INSEEP为负或零,则不读取第6项。)** 输入说明中使用的字段的解释**
- NRCHOP -- 是充水选项代码(标志)。充水率在二维数组RECH中定义,每个垂直列有一个值。因此,充水应用于每个垂直列中的一个单元,选项代码确定选择哪个单元。
1 - 充水仅应用于顶部网格层。
2 - 在数组IRCH中指定了充水的垂直分布。
3 - 充水应用于每个垂直列中最高的活动单元。常数压力节点拦截了充水,防止更深的渗透。
- IRCHCB -- 是流动模拟的标志和单元数。
如果IRCHCB > 0,则是记录每当ICBCFL(请参阅输出控制)设置时,将记录每个单元的流动术语的单元号。
如果IRCHCB # 0,则不会打印或记录单元格的流动术语。
- IUNCNF -- 是非约束选项代码(标志)。
如果IUNCNF > 0,则采用了充水渗流面边界条件,充水
(RECH)将一直施加,直到水位达到PNDEL数组中描述的塘池高程。
如果IUNCNF # 0,则总是施加指定的充水(RECH),就像在受限制的情况下一样。
- MXSEEP -- 是同时激活的渗流边界单元的最大数量。
- MXZRCH -- 是用于输入充水的区域时间序列标志。如果MXZRCH = n,则充水会随时间序列变化,n是模拟中使用的充水区域的总数。
如果MXZRCH = 0,则MODFLOW的充水输入结构仅在应力周期变化时使用。
注意1:MXZRCH > 0表示充水变化的时间尺度与MODFLOW的应力周期不同(可能小得多),并且将为模拟中的每个区域提供独立的时间序列(在文件*.RTS中详细说明)。
注意2:分区索引数组标识了位于特定区域内的网格块。值为零的索引允许在该应力周期的默认应用MODFLOW的充水。
注意3:充水时间序列中的分区位置可以在应力周期之间变化。因此,分区本身可能会在应力周期基础上变化。这由标志INRECH控制,该标志还控制每个应力周期的充水率的读取。
- IPNDOR -- 是读取塘池孔隙度的标志。如果IPNDPOR = 0,则假定塘池中的整个单元区域都被浸润。
如果IPNDPOR = 1,则读取塘池孔隙度数组。塘池孔隙度是湖泊/池塘浸润的总单元面积的分数,并用于确定塘池存储的总体积,如果发生塘池。
- INRECH -- 是RECH读取标志。如果INRECH $ 0,则会读取充水率数组(RECH)。
如果INRECH < 0,则使用上一个应力周期中的充水率。
注意:如果MXZRCH > 0,则INRECH还控制区域索引的输入数组(IZNRCH)。因此,区域位置可以在每个应力周期上更改。
- INIRCH -- 是IRCH读取标志。当NRCHOP为2时,如果INIRCH $ 0,则会读取层号数组(IRCH)。
如果INIRCH < 0,则重用上一个应力周期中使用的数组(IRCH)。
注意:当NRCHOP为1或3时,将忽略INIRCH。
- INPNDEL -- 是PNDEL读取标志,当使用非约束选项(IUNCNF > 0)时。
如果INPNDEL $ 0,则会读取塘池高程数组(PNDEL)。
如果INPNDEL < 0,则使用上一个应力周期中的塘池高程。
- INSEEP -- 是渗流边界单元的标志和计数器。
如果INSEEP < 0,则将重用上一个应力周期中的渗流边界数据。
如果INSEEP $ 0,则INSEEP是当前应力周期内渗流面边界的数量。
- RECH -- 是充水率(L/T)。仅当INRECH大于或等于零时才读取。
- IZNRCH -- 是每个节点的区域值的积分器数组。
注意:将在每个应力周期中使用默认值,将其读入RECH数组时,将应用于区域值为零的单元。
- IRCH -- 是定义每个垂直列中应用充水的层号数组。仅当NRCHOP为2且INIRCH大于或等于零时才读取。
- PNDEL -- 是定义每个垂直列中塘池高程的塘池高程数组。通常,如果地表上方没有水池,塘池高程与地表高程相一致。
仅当IUNCNF大于零且INPNDEL大于或等于零时才读取。
- PNDPOR -- 是定义湖泊/池塘存在时水位超过域顶部时占用的单元区域的分数的塘池孔隙度数组。
- Layer -- 是表示渗流面边界的单元的层号。
- Row -- 是表示渗流面边界的单元的行号。
- Column -- 是表示渗流面边界的单元的列号。
- Elevation -- 是渗流面的高程。