pSolve batch 控制文件总览

编写控制文件是 pSolve 全局解算的重要步骤。这里从源代码层面整理控制文件的结构和语法规范。Leonid 曾经编写过相关说明 solve_guide_03.syn，可以结合阅读。

实际使用时，通常在已有控制文件上修改，不必从头编写。例如，可以参考 gsf2026b 的控制文件，再根据自己的解算需求调整设置。本文对应的源码版本为 psolve-20260812。文中的源码路径均相对于 psolve 目录，例如 progs/solve/batch/ctrlfl.f。

1. batch 控制文件总览

阅读控制文件时，可以把内容分为两类：预处理指令和参数配置。这两类内容不一定各占连续的一块文本，而是由不同的程序处理：CTRLS 先预处理文件、检查实验列表，再调用 CTRLFL 解析配置段。

CTRLS 还负责检查恢复状态、准备全局数据。它返回后，BATCH 保存配置并完成其余准备，随后调用 PRCES 组织各次实验的处理。

flowchart TD
    A["BATCH：初步初始化"] --> B["CTRLS"]
    B --> C["预处理控制文件"]
    C --> D{"实验数据格式"}
    D -->|GVF / GVH| E["CH_DB_GVF：解析并检查数据库"]
    D -->|旧格式| F["CHSUP 等：检查 superfile"]
    E --> G["CTRLFL：解析配置段"]
    F --> G
    G --> H["检查配置、判断恢复状态"]
    H --> I["初始化或恢复全局数据"]
    I --> J["返回 BATCH"]
    J --> K["BATCH：保存配置及其他准备"]
    K --> L["PRCES：组织实验处理"]

2. 预处理指令与 CTRLS

2.1 预处理指令

预处理指令可以省略。它们的识别范围由读取流程决定，不专属于 $SETUP 段。即使某个示例把 DEFINE 放在 $SETUP 标题下面，处理它的仍是 CTRLS，而不是 GSETUP。

这里有三点需要区分：

• CTRLS 在原始文件中识别 $ARCS，判断条件是该行从第一列开始的前五个字符为 $ARCS。因此，标题前不要缩进，也不要把“在包含文件中插入 $ARCS”当作原始文件的分界方式。
• @<file> 在这里表示直接插入文本。被插入文件中的 DEFINE、UNDEFINE、@ 或行首 ${...} 不会再次经过同一套预处理。它与 $ARCS 实验记录中由 GARC 处理的 @<file> 入口不是同一个入口。
• ARCFILE 是观测实验列表的预处理入口，不是 GARC 逐条解析的实验选项。$ARCS GVF 和 $ARCS GVH 也能引用外部列表，但处理程序和旧格式不同，不能把两条路径的细节混为一谈。

2.2 CTRLS 的处理过程

CTRLS 是 BATCH 进入实际解算前的准备程序。它整理运行参数、控制文件、实验列表和已有运行状态，为后续的 PRCES 做好准备。

1. 准备错误文件，读取运行参数。

程序先调用 RTIMER，再重新建立本次运行使用的错误文件：

PRE_SCR_DIR(1:PRE_SD_LEN)//'ERRF'//PRE_LETRS

通常可将它记作 ERRFxx，其中 xx 是本次运行的标识。随后，程序从共享变量 STRING 中依次取出运行字段：跳过前四个字段，读取进程编号 IND_PROC、进程总数 NUM_PROC，再读取控制文件名 CFNAME。

2. 预处理原始控制文件，生成 .XPND 文件。

假设原始文件为 solution.cnt，程序会把预处理结果写入 solution.cnt.XPND。在正常读取路径中，这一文件都会生成并用于后续读取，不是只有使用 ARCFILE 时才生成。在遇到原始文件中的第一个 $ARCS 标题前，程序逐行处理：

3. 解析实验列表，检查数据库并统计数量。

程序读取 $ARCS 标题行的第二个字段，也就是紧跟 $ARCS 的格式标记。若为 GVF 或 GVH，就进入 CH_DB_GVF；否则进入旧式 superfile 检查路径。

flowchart LR
    A["读取原始文件的实验列表"] --> B{"格式"}
    B -->|"GVF / GVH"| C["CH_DB_GVF：查找数据库入口"]
    B -->|"旧格式"| D["CHSUP 等：检查 superfile"]
    C --> E["写入展开记录并统计实验数"]
    D --> E

GVF/GVH 路径中的主要调用是：

CALL CH_DB_GVF ( VCAT, STRING1, N_SES, N_MIS, IUER )

CH_DB_GVF 读取 VCAT 配置，查找数据库 envelope 文件，并处理实验列表及其中的 ARCFILE <file>。写入 .XPND 的记录以 GVF 开头，包含找到的 envelope 路径及该实验剩余的配置。返回后，N_SES 赋给 NARCS；若 N_MIS > 0，程序退出。旧格式则通过 CHSUP、GET_SUPERFILE_NAME 和 INQUIRE 检查 superfile。这里检查的是实验列表及相应入口文件能否找到，尚未完整读取和验证全部观测内容。

4. 切换到 .XPND，调用 CTRLFL。

完成预处理后，程序保留原始文件名，并把当前读取对象改为展开文件：

CFNAME_ORIGINAL = CFNAME
CFNAME = EXPAND_CTLFILE
TYP = CFOPEN ( CFNAME )

随后调用 CTRLFL。从这一步开始，各配置段解析程序读取的是 .XPND；原始文件名仍用于报告及恢复检查等用途。

5. 检查部分配置组合。

CTRLFL 返回后，CTRLS 继续检查一些跨配置项的限制。例如：

• 协方差输出在相应路径下要求设置实验中间结果文件目录；针对 CGM 的输出另有条件分支。
• WEIGHTS MAKE 或 WEIGHTS APPEND 使用 MYWAY 算法时，不能同时设置 NORATE_FLAG YES。

因此，“程序已读懂这一行”与“这些设置能组合使用”是两件事。部分组合检查在配置段解析器内完成，另一些在 CTRLS 中完成。

6. 确定是否恢复中断的运行。

如果 KBIT(PRE_IP(2), INT2(8)) 为真，程序禁止尝试恢复，直接设置 RESTRT=.FALSE.；否则调用 RSTORS 读取并检查已有状态。

7. 初始化或恢复全局数据。

程序先处理与恢复有关的共享状态，再调用 DISCONT_INIT 初始化测站间断运动和分段位置模型的数据，并清除源位置、站位置和站速度替换文件的初始化标志。只有满足以下条件，才进入这里的 GLO_INIT 和 GLO_RESTORE 路径：

.NOT. TRAIN .AND. SOLTYP .NE. 'I'

GLO_INIT 初始化 GLBMEM，之后根据运行状态选择恢复方式：

新运行还会检查输出 CGM 是否与已有文件冲突。实际读取的环境变量是 PSOLVE_CGM_DIR；变量为空时，使用程序中的默认目录 CGM_DIR。

GLBSET 的调用位置取决于运行状态：部分恢复情形在全局矩阵初始化前调用 GLBSET(...,1)；新运行在末尾调用 GLBSET(...,0)。它参与设置运行状态、报告输出和进度信息，不能概括为“所有路径都在最后调用一次”。相关位置为 progs/solve/batch/ctrls.f:698 和 progs/solve/batch/ctrls.f:865。

3. 参数配置段与 CTRLFL

3.1 配置结构与书写规则

本文使用三个名称描述参数配置部分：

例如：

$SETUP
 SOLUTION  COMPLETE
 ARC_FILES /scr/psolve/sub_matrices

$SETUP 是配置段标题，SOLUTION 和 ARC_FILES 是配置关键词，后面跟取值或参数。这几行共同组成一个配置段。

书写时遵循以下规则：

• 配置段标题从第一列开始，段内配置行以空格开头。
• 第一列为 * 的行是注释。CFREAD 也会跳过空行；不要把 # 或 ! 当作所有位置都通用的注释符。
• 关键词和固定取值按本文列出的大小写书写。各解析器的大小写处理并不完全一致。
• \ 只在相应解析器支持的列表或子句中用于续行，不是任意位置都有效的通用续行符。

同一段内，相互独立的关键词通常可以调整顺序，但不能一概而论。有些关键词需要使用先前已经读取的设置。例如，$SETUP 中的 WEIGHTS MAKE 在当前 FAST_MODE 不为 NONE、又尚未读取 FAST_COV 时会报错，因此应先写 FAST_MODE、FAST_COV，再写 WEIGHTS。

建议把 $SETUP 放在第一个配置段，因为它设置基本配置并初始化后续程序使用的变量。它不必是文件第一行，前面可以有空行、注释或适用的预处理指令。CTRLFL 没有显式要求 $SETUP 必须排第一，但这不代表所有配置段可以任意换序。CTRLFL 明确检查以下顺序：

• $FLAGS 在 $CARRY 之前。
• $FLAGS 和 $CONSTRAINTS 都在 $SUPPRESSION 之前。

另外，应将 $ARCS 放在文件末尾。虽然 CTRLFL 的局部计数逻辑能在遇到下一段标题时停下，但更早的 CTRLS 已把第一个 $ARCS 后的内容交给实验列表处理，不能据此在其后继续安排普通配置段。

下面列出全部配置段。这里的“必需”只表示该配置段必须存在，不表示其中每个关键词都必须出现，也不表示必须开启相应功能。表的顺序便于查阅，不是程序强制的完整顺序。

后文用表格列出关键词，用流程图说明读取与分派关系。图中分支表示不同处理路径，不表示必须依次启用所有关键词。表格中的语法符号统一如下：

“实现异常”表示源码虽然识别该入口，但存在错误或不完整的处理；“不支持”表示程序会明确拒绝。涉及问题的写法不能仅凭解析通过就视为可靠可用。文中带 <...> 的例子用于说明结构，使用前需要替换为实际值。

3.2 SETUP（基本设置段）

$SETUP 设置解算方式、矩阵文件、实验中间结果文件、计算算法和权处理等基本选项，由 GSETUP 解析。它还初始化后续程序需要的变量，因此建议将它放在第一个配置段。GSETUP 先设置初始值，再逐项读取配置；有些关键词直接在本程序中处理，有些交给辅助解析器。读到本段末尾后，程序检查选项是否齐全、组合是否允许，并把下一段标题留给 CTRLFL。下图表示这一处理关系，不表示关键词的书写顺序。

flowchart LR
    A["读取下一个关键词或记录"] --> B{"到达本段末尾？"}
    B -->|否| C{"解析方式"}
    C -->|普通选项| D["GSETUP 直接读取"]
    C -->|复杂结构| E["调用辅助解析器"]
    D --> A
    E --> A
    B -->|是| F["检查配置并放回下一段标题"]
    F --> G["返回 CTRLFL"]

关键词列举如下：

表中的“必须出现”指解析器明确要求的关键词，不表示其余关键词可以在任意条件下省略。GSETUP 还会根据解算类型、计算方式和文件设置进行组合检查。

复杂子项列举如下：

几个常见的组合限制是：GLOBAL_ONLY 要求 TRAIN NO，SUBTRACT 要求 TRAIN YES；全局解算的 FAST_MODE 只接受 NONE / B1B3D，独立解算接受 NONE / B3D / DOT_PRODUCT。快速模式下，WEIGHTS MAKE 应写在 FAST_COV 之后；使用 MYWAY 的 MAKE / APPEND 还要求 TRAIN YES。

GSETUP 读到 CGM 时先取三个字段，等读完整个 $SETUP 段后再确定其用途：

因此，CGM 应单独成行，普通写法不要额外追加第三个参数。输入文件宜使用完整路径，避免两字符文件名被误当成旧式标签。SAVE 最终怎样命名和写出文件，由后续处理阶段决定。

SETUP.ARC_FILES 指定各次实验的中间结果文件存放在哪里，并用可选的 SAVE 控制是否在正常后向处理后保留这些文件。这里的 arc 对应一次实验（session）。全局解算的前向处理会生成这些中间结果，后向处理再使用它们，因此本文将 arc files 称为“实验中间结果文件”。例如，gsf2026b 控制文件中有：

$SETUP
  ARC_FILES /scr/psolve/sub_matrices

这表示使用 /scr/psolve/sub_matrices 存放实验中间结果。它没有指定 SAVE，也不列出待处理的实验。另一个容易混淆的入口是 $ARCS 列表处理中的 ARCFILE：它指定外部实验列表文件，由预处理或数据库列表检查程序读取，不是 GARC 的逐实验选项。

3.3 FLAGS（参数估计设置段）

$FLAGS 指定哪些参数需要估计，以及采用哪些分量和时间模型，由 GFLAGS 解析。这里的 YES / NO 通常控制参数估计，不应直接理解为是否使用相应观测；观测选择主要见 $DATA。本段必须出现的关键词有：SOURCES、STATIONS、VELOCITIES、NUTATION、UT1/PM、PRECESSION、RELATIVITY、ATMOSPHERES、CLOCKS、BASELINE_CLOCKS。必需的关键词也可以采用其允许的关闭选项。段内还有两条明确的顺序要求：SOURCES 必须先于 PROPER_MOTIONS，STATIONS 必须先于 VELOCITIES。

flowchart LR
    A["读取下一个关键词或记录"] --> B{"到达本段末尾？"}
    B -->|否| C{"解析方式"}
    C -->|普通选项| D["GFLAGS 直接读取"]
    C -->|模型或例外列表| E["调用辅助解析器"]
    D --> A
    E --> A
    B -->|是| F["检查必需项、历元和解算方式"]
    F --> G["放回下一段标题，返回 CTRLFL"]

关键词列举如下：

部分关键词的结构非常复杂。这些复杂关键词的结构性子项属于相应关键词内部，不是 $FLAGS 的独立关键词。

HEO 的 CNS 并不是任意增加的可选记录。当前解析器要求 16 个约束名全部出现一次，其名称由三部分组合：HEO_VAL / HEO_ERM、E1E2 / E3、HAR / CROSS / SHIFT / DRIFT。例如 HEO_VAL_E1E2_HAR。应写完整名称，避免子串匹配误认。

3.4 OUTPUT（输出设置段）

$OUTPUT 主要设置解算报告、残差、协方差和 SINEX 等输出。需要留意，少数关键词还会影响计算过程：NORMAL_MATRIX ZERO 会清零法方程，IONOSPHERIC_MODEL 可以启动电离层模型计算和数据库更新，不能把它们都理解为显示格式开关。CTRLFL 调用 GOUTPT 解析本段。多数关键词由 GOUTPT 直接处理，COVARIANCES、CORRELATIONS 和 SINEX 则分别交给专门的解析程序。

flowchart LR
    G["GOUTPT：读取关键词"] --> K{"关键词类型"}
    K -->|"其他关键词"| D["GOUTPT 直接解析"]
    K -->|"COVARIANCES"| C["GTCVR"]
    K -->|"CORRELATIONS"| R["CORRELATIONS_PARSER"]
    K -->|"SINEX"| S["SINEX_PARSER"]
    D --> E{"本段是否结束"}
    C --> E
    R --> E
    S --> E
    E -->|"否"| G
    E -->|"是"| F["检查必需项与组合限制，返回 CTRLFL"]

图中箭头表示解析分派和返回关系，不表示关键词的书写顺序。

关键词列举如下：

这里的“必需”指关键词必须出现，不是必须启用该功能。例如，即使使用 CORRELATIONS YES，仍需写出 COVARIANCES NO。

RESIDUALS 的 SPOOL 分支写成了 RESFILE(1:4) .EQ. 'SPOOL'。左侧只有四个字符，不可能与五个字符的 SPOOL 相等。因此，GOUTPT 会把它当作普通文件前缀，创建 SPOOL.bin 和 SPOOL.ndx。这不等于后续程序完全不识别 SPOOL，但它确实无法避开本处的文件创建操作。普通前缀对应的已有 .bin 文件也会在解析阶段被删除并重建。

COVARIANCES 的结构较为复杂：

YES 后有两种完整结构：

 COVARIANCES YES (ALL / BY_ARC / CGM) (STA / NUT / EOP / ALL / SOU)
 COVARIANCES YES <database> <version> (STA / NUT / EOP / ALL / SOU)

而CORRELATIONS 也存在复杂的结构性子项：

PARAM_INCLUDE、PARAM_EXCLUDE 从属于相应的参数块；SES_INCLUDE、SES_EXCLUDE 只出现在未实现的 CROSS_LOC 分支中。文件名由同一源码文件中的 CP_CHECK 检查；相对路径以 SAVE_DIR 为起点，不是当前工作目录。

如果要设置SINEX YES，那么必需设置14 个子句，列举如下：

四个输入文件子句都会检查文件是否存在，没有单独的 NONE 豁免分支。这里的 SINEX COVARIANCES 是 SINEX 内部子句，与 $OUTPUT 顶层的 COVARIANCES 不是同一个入口。

LISTING_OPTIONS 决定了SPOOL文件的输出内容和格式，其取值列举如下：

两项同时使用时，可以写成：

 LISTING_OPTIONS SRC_STAT POST2024 \
                 SEG_STYLE POST2026

当前 include/glbc4.i:514 将默认分段输出格式设为 POST2026。因此，即使省略 LISTING_OPTIONS，也可能触发 APRIORI_ZENDEL 的强制开启。

IONOSPHERIC_MODEL 需要设置后续子句（即使前面的开关为 NO）。当其值设置为非 NO 时，程序将按下表顺序逐项读取：

3.5 CARRY（参数全局性设置段）

$CARRY 指定已选参数是多个实验共用的全局参数，还是各次实验分别求解的局部参数。它决定参数的全局或局部属性；是否估计这些参数，由 $FLAGS 设置。CTRLFL 调用 GCARRY 解析本段，并传入 $FLAGS 中的测站分量设置。$FLAGS 必须先于 $CARRY 出现，这一顺序由 CTRLFL 明确检查。

flowchart LR
    G["GCARRY：读取关键词"] --> K{"关键词类型"}
    K -->|"SOURCES"| S["GTLST：读取开关与名称列表"]
    K -->|"STATIONS"| T["GTLST_ST_CA：读取分量与名称列表"]
    K -->|"AXIS"| A["GCARRY：读取 YES / NO"]
    K -->|"TIDES"| X["报错并停止"]
    S --> E{"本段是否结束"}
    T --> E
    A --> E
    E -->|"否"| G
    E -->|"是"| F["检查 SOURCES 和 STATIONS，返回 CTRLFL"]

关键词列举如下：

$CARRY TIDES 与 $FLAGS TIDES 的处理不同；后者保留了 NO 兼容分支。

3.6 DATA（数据设置段）

$DATA 选择参与解算的观测数据，包括源、测站、基线、观测量类型及高度角等筛选条件。它与 $FLAGS 的分工不同：$DATA 选择使用哪些数据，$FLAGS 选择估计哪些参数。CTRLFL 调用 GDATA 解析本段。SOURCES、STATIONS、TYPE 必须出现，其余关键词可选；这八个关键词均不允许重复。

flowchart LR
    G["GDATA：读取关键词"] --> K{"关键词类型"}
    K -->|"SOURCES / STATIONS"| L["GTLST：名称列表"]
    K -->|"ELEVATION"| E["GTELEV：高度角与测站分组"]
    K -->|"WVR_MASK"| W["GTWVM：掩码与测站分组"]
    K -->|"BASELINES"| B["GBASE：基线列表"]
    K -->|"TYPE / NORATE_FLAG / EVERY"| D["GDATA 直接解析"]
    L --> N{"本段是否结束"}
    E --> N
    W --> N
    B --> N
    D --> N
    N -->|"否"| G
    N -->|"是"| F["检查必需项，返回 CTRLFL"]

GBASE 文件头关于调用方和配置段的旧注释不准确，这里按实际调用关系将它归入 $DATA。关键词列举如下：

根据include/solve.i:880 中的代码表，TYPE 的取值包含以下 22 个缩写：

GRPRAT  PHSRAT  SNBRAT  GRPONL  PHSONL  SNBONL  RATONL
G_GXS   PX_GXS  PS_GXS  PX_GX   PX_GS   PS_GX   PS_GS   P_PXS
GX      GS      PX      PS      SNG_X   SNG_S   FUSED

GDATA 还单独接受四种长名称，并作如下赋值：

3.7 CALIBRATIONS（校准设置段）

$CALIBRATIONS 用于选择和控制校准项，由 CTRLFL 调用 GCALIB 解析。本段同时保留了按校准名称设置的旧式写法，以及 RESET、KEEP、ENABLE、DISABLE 等列表写法。两类写法有组合限制，不能任意混用。本段必须存在，但没有某一个始终必填的关键词；是否完整，要看段末的组合检查。

flowchart LR
    A["$CALIBRATIONS → GCALIB"] --> B{"识别入口"}
    B -->|"AVAILABLE"| C["GETAVL：加载可用校准项"]
    B -->|"USE"| D["GETUSE：读取名称对应关系"]
    B -->|"GROUP"| E["GETGRP：读取组内选择"]
    B -->|"校准名称"| F["GTCALST：读取测站选择"]
    E --> F
    B -->|"其他受支持关键词"| G["GCALIB 直接读取"]
    C --> H["继续读取；段末检查与默认设置"]
    D --> H
    F --> H
    G --> H
    B -->|"OLD_IONOSPHERE"| I["报错：不再支持"]

关键词列举如下：

3.8 MAPPING（先验模型设置段）

$MAPPING 用于指定先验参数、模型文件和相关改正设置。CTRLFL 调用 GMAP 解析本段；多数文件参数由 GMAPF 读取，较复杂的结构交给专门的解析程序。

flowchart LR
    A["$MAPPING → GMAP"] --> B{"识别入口"}
    B -->|"多数文件设置"| C["GMAPF：读取文件及从属子句"]
    B -->|"SPLINE_DISPLACEMENTS"| D["PARSE_BSP"]
    B -->|"POSITION_VARIATIONS"| E["GPOSVAR"]
    B -->|"HARMONIC_EOP / ERM"| F["分别调用 GHEO / GERM"]
    B -->|"其他受支持关键词"| G["GMAP 直接处理"]
    C --> H["继续读取；段末检查配置"]
    D --> H
    E --> H
    F --> H
    G --> H
    B -->|"独立 REF_DATE / NUTATION_MODEL"| I["报错并停止"]

图中的“其他受支持关键词”包括 FIX、METRIC_TENSOR、VTD_CONF、EXTERNAL_TRP_DELAY 和 ANTENNA_THERMAL。NUTATION_EXPANSION 也在 GMAP 内处理，但只检查名称，没有据此设置模型变量。

关键词列举如下：

这里的 <model> 不一定表示一组已在本层检查过的固定取值。例如，PLATE_MODEL 在这层只读取模型字符串，并没有像 METRIC_TENSOR 那样逐一检查名称。部分复杂子句的结构列举如下：

注意两个容易混淆的拼写：EXTERNAL_TRP_DELAY 使用 REQUIRE，ANTENNA_THERMAL … INSITU 使用 REQUIRED。EARTH_ORIENTATION 的省略行为也不完全相同：指定插值方式、但省略 zonal 选项时，程序设为 UT1S；两者都省略时，设为 NO_ZONAL。

在本配置段内存在限制：

• PLATE_MODEL 与 VELOCITIES 共用状态标志，不能作为两个可叠加的设置同时指定。
• 非空的 POSITION_VARIATIONS 不能与 SETUP.TRAIN YES 同时使用。
• 启用的 HI_FREQ_EOP 与 HARMONIC_EOP 互斥。
• POSITION_VARIATIONS、SPLINE_DISPLACEMENTS、VTD_CONF、EXTERNAL_TRP_DELAY 和 ANTENNA_THERMAL 均有重复指定检查。
• GMAPF 的普通分支读完主参数后，还会取走一个词并写入 NAME2，例如形成 SRCMAP2、STAMAP2 或 VELMAP2。

3.9 PARTIALS、CONTRIBUTIONS 与 MODE_CALIBRATIONS

这三个配置段的入口较少，因此合并介绍。它们都涉及名称选择，但语法并不相同：SET 是条目前缀，NO、NONE、IN 等是固定控制词，实际偏导项或校准项的名称则由文件或数据库提供。三者需要分开理解。

flowchart LR
    A["CTRLFL"] --> B["$PARTIALS：可选"]
    A --> C["$CONTRIBUTIONS：必需"]
    A --> D["$MODE_CALIBRATIONS：可选"]
    B --> E["GPART"]
    A -.->|"省略本段时仍初始化默认项"| E
    C --> F["GCONT"]
    D --> G["GMCAL"]
    E --> H["可选 SET；名称后可带 ON / OFF"]
    F --> I["可选 SET；固定控制词或名称列表"]
    G --> J["每条必需 SET；固定控制词或名称列表"]

图中的条目均由对应的 GPART、GCONT、GMCAL 直接解析。三者的源文件均位于 progs/solve/batch/。

3.9.1 PARTIALS（偏导数设置段）

本配置段选择要使用的偏导项，可以省略。即使省略，CTRLFL 仍会调用 GPART，加载可用偏导项及默认项。

可用名称和默认项通过 ATMAVL_N 从保存目录下的 partial_calibrations 文件读取。而开关的实际处理为：

因此，条目顺序可能影响最终结果。另一个需要注意的问题是：名称没有匹配时，当前实现不会报错，而是保留原有选择。拼写错误可能因此被忽略。相关代码在 progs/solve/batch/gpart.f:110；修正时应增加“是否匹配”的检查，并测试已知名称的 ON / OFF、默认项 OFF 和未知名称。

3.9.2 CONTRIBUTIONS（附加改正项设置段）

本段选择要加入的附加改正项，必须存在，并至少提供一个条目。每个条目前的 SET 都可以省略。

多个名称可以写成：

 SET <contribution_name_1>
 SET <contribution_name_2>

解析器也接受同一行上的多个名称，但逐条写出 SET 更便于阅读。名称区分大小写；GCONT 只将名称加入列表，不在这里核验相应改正项是否存在。

3.9.3 MODE_CALIBRATIONS（模式相关校准设置段）

本段选择模式相关校准项，可以省略。一旦出现，就必须至少提供一个条目，而且每个条目都要以 SET 开头。

例如：

 SET <mode_calibration_name_1>
 SET <mode_calibration_name_2>

DEFAULT 在 GMCAL 中没有专门的控制分支，不能照搬 $CONTRIBUTIONS 中的别名关系。

这里有两处实现问题需要区分：

• NONE 会命中前面的两字符 NO 比较，因此它能够进入关闭控制分支；单独的 NONE 比较本身无效。问题在于匹配过宽，以 NO 开头的真实名称也可能被误判。
• IN 分支设置的是 KNON，不是 KIN。独占检查仍会触发，但错误提示可能把 IN 当成 NO / NONE。 修正时应比较完整控制词，并让 IN 使用 KIN 标志。验证时需要覆盖 NO、NONE、IN 各自独立使用、重复、混用，以及以 NO 开头的名称。

这三个配置段的语法区别列举如下：

3.10 CONSTRAINTS（约束设置段）

$CONSTRAINTS 用于设置参数约束，包括约束强度、参与约束的测站或源，以及例外名单。是否估计这些参数，仍由 $FLAGS 等设置决定。本段必须存在，但没有单独必需的关键词；可以给出具体约束条目，也可以单独写 NONE。CTRLFL 调用 GCONST 解析本段。大部分关键词由 GCONST 直接读取，部分结构交给辅助程序。

flowchart LR
    A["$CONSTRAINTS：GCONST"] --> B{"当前入口"}
    B -->|STATIONS / VELOCITIES| C["GSTAVELCNST"]
    B -->|六类整体约束| D["GYNEXL"]
    B -->|NUVEL_COVAR| E["GNUVCOV"]
    B -->|其他受支持入口| F["GCONST 内部读取"]
    B -->|CLOCK_REF_REQ| X["报错并停止"]
    C --> G["继续读取；段末检查"]
    D --> G
    E --> G
    F --> G

其中，GYNEXL 先拆分整体约束的选项、开关和例外名单，再由 GCONST 检查各项选项。各入口的源码位置见下表。

关键词列举如下：

EARTH_ORIENTATION 的子句结构需要额外注意：

不过，这些部分不能任意排列。例如：

 EARTH_ORIENTATION YES RATES

这行会先进入 YES 分支的标准差读取过程，把 RATES 当成数值而报错。若要关闭速率约束，应在合适的基本设置后直接写尾部 NO，不要写旧说明书中的 RATES NO；后者同样会把 NO 当成数值读取。

STATIONS 与 VELOCITIES 的结构几乎一致。除单独的 NO 外，必须给出 XYZ、UEN 各一次：

使用时注意：

• YES 后必须有 SIGMA 和三个数值。
• NO 可以单独使用，也可以带 SIGMA 和例外名单。
• 若要在 NO 后使用 EXCEPT，前面仍须给出 SIGMA 和三个数值。
• XYZ、UEN 的顺序可交换，并支持 \ 续行。

例如，下面只展示语法结构：

 STATIONS XYZ YES SIGMA <x> <y> <z> \
          UEN NO

六类整体约束选项共同结构为：

 <keyword> <options> (YES / NO) [EXCEPT <object_list>]

GYNEXL 先收集 YES / NO 之前的选项，再读取开关和例外名单，最后把选项交回 GCONST 检查。

还有以下限制：

• 启用两类平移约束时，应显式给出 ALL 或 HORIZ。只写 YES 会留下内部标记 ??，不宜作为正常配置示例。
• 单独关闭时，可以写 <keyword> NO。
• 支持 GLOBAL、LOCAL 的三个分支可以同时设置两者；两者都未指定时，源码补设 GLOBAL。
• 开关之前的选项顺序较灵活，但 SIGMA 后的数值、RIGHT_PART 后的三个数值必须紧跟对应固定词。
• 规范写法应显式给出 YES / NO。虽然解析器在行尾没有开关时会按 YES 处理，但省略开关会使配置不易读懂。
• RIGHT_PART 使用定宽实数格式读取，数值应显式带小数点。
• GCONST 传入 M_KEYS=10；结合 GYNEXL 的循环边界，开关之前实际最多容纳 9 个词，固定词和数值都计数。

在本配置段中，NONE 只设置本配置段的控制标志，不会逐项关闭已经设置的约束，也不会清除默认值。例如，GCONST 初始化时将基线时钟约束标志设为：

KBSL_CONST = .TRUE.

因此，不能将下面的配置直接解释为“所有约束均关闭”：

$CONSTRAINTS
 NONE

若要明确关闭某一项，应使用对应的 NO，并确认该关键词没有已知的实现问题。NONE 也不应放在其他约束之后充当“全部取消”；这样写可能通过当前的独占检查，但前面设置的约束仍然保留。

3.11 SUPPRESSION（数据与参数剔除设置段）

$SUPPRESSION 用于进一步限制参与解算的参数，包括分量剔除、测站绑定和部分整体限制。它与 $FLAGS 的作用不同：$FLAGS 决定估计什么，本段则对已有参数设置进一步限制。例如，分量入口中的 YES 表示剔除相应分量，而不是开启估计。CTRLFL 调用 GSUPRS 解析本段，并传入 $FLAGS 中的测站分量和参考测站名称。$FLAGS 和 $CONSTRAINTS 都必须先于 $SUPPRESSION 出现。 本段必须存在，但没有单独必需的关键词；不使用本段设置时，应单独写 NONE。

flowchart LR
    A["$SUPPRESSION：GSUPRS"] --> B{"当前入口"}
    B -->|分量及例外| C["GVELSP / GPROPSP / GSRCSP / GSTASP"]
    B -->|整体选择及例外| D["GVLOSP / GRAOSP / GDCOSP / GORISP"]
    B -->|测站绑定组| E["GVELTIE / GSTATIE"]
    B -->|其他受支持入口| F["GSUPRS 内部读取"]
    B -->|NUTATION / TIDES| X["报错并停止"]
    C --> G["继续读取；段末检查"]
    D --> G
    E --> G
    F --> G

图中合并显示了同类入口的辅助解析程序，各关键词对应的源文件见下表。CHKCMP 是分量解析器内部使用的校验函数，没有在图中展开。段末检查存在遗漏，因此“解析通过”不一定表示配置组合合理。

关键词总表

源坐标和站坐标的分量、模式与列表结构列举如下：

分量例外列表可以写成：

 VELOCITIES NO EXCEPT U-- <station_1> <station_2> \
                      -EN <station_3>

 SOURCES NO EXCEPT R- <source_1> \
                  -D <source_2>

这里先用 NO 设置默认不剔除，再为例外对象指定要剔除的分量；例外分量不是简单地把默认开关取反。

绑定组可以写成：

 VELOCITY_TIE YES <station_1> <station_2> \
              AND <station_3> <station_4>

这两组彼此独立。绑定约束作用于估计的改正量，不能直接理解为让两站的绝对坐标或绝对速度相等。这些示例只展示语法，使用前需要替换占位符。同一测站应避免出现在多个绑定组中。VELOCITY_TIE 会明确拒绝跨组重复；STATION_TIE 则会用后一次组号覆盖前一次，不会合并两个组。

GSUPRS 的段末检查没有覆盖所有入口：

因此，更应该将 NONE 规定为独立使用。它表示保留本段初始化后的状态，不是“取消此前所有设置”的命令；也不应利用上述检查遗漏，把它与其他条目混用。

3.12 ARCS（实验列表段）

$ARCS 列出本次解算要处理的实验，以及各次实验的单独设置。这部分的记录结构和选项较多，将另文说明；这里先交代它与前面各段的连接关系。

因此，前文的 ARCFILE <file>、$SETUP 中的 ARC_FILES 和实验记录里的选项分属不同用途：ARCFILE 引入实验列表，ARC_FILES 设置实验中间结果文件，而 GARC 负责解析具体实验的选项。

3.13 CTRLFL 的处理过程

CTRLFL 按文件顺序识别配置段，调用相应解析程序，并检查配置段是否重复、缺失或违反部分顺序要求。它接手的文件已经由 CTRLS 预处理，并通过 CFOPEN 打开。

flowchart LR
    A["初始化配置及共享状态"] --> B{"CFREAD：下一条记录"}
    B -->|"配置段标题"| C{"检查标题并分派"}
    C -->|"普通配置段"| D["调用相应解析程序"]
    C -->|"$ARCS"| E["保存位置并计数"]
    D --> B
    E --> B
    B -->|"文件结束"| F["补充默认值并检查必需段"]
    F --> G["返回 CTRLS"]

图中省略了报错终止分支。CFREAD 返回下一条有效记录；在主循环中，这条记录应当是配置段标题。

1. 初始化部分配置和共享状态。

程序先设置几个初始值：

USER_PROG      = ' '
USER_PART_PROG = ' '
KUSER_PART     = .FALSE.
KGLOBONLY      = .FALSE.
USER_BUFF      = 'N'
EOPMID         = 0

随后读取、清零并写回 GLBC3 对应的共享数据区，同时重置若干约束标志和索引。关键调用为：

CALL USE_GLBFIL_3 ( 'OR' )
CALL NOUT ( JGLBC3_BLOCKS*INT2(BLOCK_BYTES), BEGMARK_GLBC3_I2 )
! 此处还重置若干标志和索引。
CALL USE_GLBFIL_3 ( 'WC' )

所以，解析控制文件会修改共享状态，并伴随状态文件的读写。这里也只初始化了一部分变量；其他默认值仍由各段的解析程序设置。

2. 通过 CFREAD 取得配置段标题。

主循环结构为：

LENGTH = CFREAD(STRING)
DO WHILE (.NOT. CFEOF(IDUM))
   ! 检查段标题并调用对应的解析程序。
   LENGTH = CFREAD(STRING)
ENDDO

每次进入循环，CTRLFL 都要求当前记录的第一列为 $，再用 SPLITSTRING 取出段名。

主循环读取段标题，段内配置由被调用的解析程序继续读取。 这些程序通常在遇到下一段标题后，通过 CFUNRD 把标题放回读取缓冲区，使 CTRLFL 的下一次 CFREAD 能取到它。因此，CTRLFL 并不是先读出全部配置行，再逐行分发。

3. 按文件中的顺序调用解析程序。

遇到 $SETUP 就调用 GSETUP，遇到 $FLAGS 就调用 GFLAGS，其他配置段也按同样方式分派。多数分支采用以下结构：

IF ( KFLAGS ) CALL FERR(...)
CALL GFLAGS(...)
KFLAGS = .TRUE.

这里的 KFLAGS 表示“是否已经处理过 $FLAGS 段”，不是某个参数的估计开关。各个段标志主要用于检查重复段，以及在文件末尾检查必需段是否齐全。

源码中各分支的排列顺序不等于实际读取顺序。实际顺序由控制文件决定，程序只对第 3.1 节列出的部分先后关系做了显式检查。

4. 对 ARCS 保存位置并计数。

$ARCS 是特殊分支：CTRLFL 此时不调用 GARC，而是执行下面的读取过程。

CALL CFPOS ( ARCREC )
ITARCS = 0
LENGTH = CFREAD(STRING)

DO WHILE ( STRING(1:1).NE.'$' .AND. .NOT.CFEOF(IDUM) )
   ITARCS = ITARCS + 1
   LENGTH = CFREAD(STRING)
ENDDO

CALL CFUNRD ( LENGTH, STRING )
KARCS = .TRUE.

它保存内部读取位置 ARCREC，统计 CFREAD 返回的实验记录数 ITARCS，并把读取到的下一段标题放回缓冲区。ARCREC 是 CFPOS 与 CFSPOS 使用的内部记录位置，不是字节偏移，也不是实验编号。后续的 PRCES 用它重新定位实验列表。ITARCS 与 CTRLS 中的 NARCS 来源也不同：前者统计展开文件中的有效记录，后者来自前面的实验列表解析和数据库检查。

5. 完成依赖整份文件的设置。

部分值要等整份文件读完才能确定。例如：

IF ( FWDOUT .LT. 0 ) FWDOUT = ITARCS

这对应前文的 OUTPUT.FORWARD LAST：先保存标志，知道实验总数后，再将它换成最后一次实验的编号。若省略 $PARTIALS，程序仍执行：

IF ( .NOT. KPART ) CALL GPART ( KPART )

此时 GPART 初始化可用偏导项及默认选择，但不读取该段配置。$MODE_CALIBRATIONS 没有对应的“缺省时补调用 GMCAL”路径。

6. 检查必需配置段，返回 CTRLS。

文件末尾检查 11 个必需段是否都已出现：$SETUP、$FLAGS、$OUTPUT、$CARRY、$DATA、$CALIBRATIONS、$CONTRIBUTIONS、$MAPPING、$CONSTRAINTS、$SUPPRESSION 和 $ARCS。缺少任意一个都会报错；只有 $PARTIALS 和 $MODE_CALIBRATIONS 可省略。这一检查只确认配置段存在。各段中的关键词、参数数量和组合是否有效，仍由对应解析器及后续检查决定。