GAMIT 10.35中文文档初稿(待完善)

版本 10.3 麻省理工学院 地球、大气、行星科学系 2009 年 6 月 1 日
题记这份gamit翻译文档建议仅作为英文原版的参考，由于翻译水平有限，如有误译， 请告之，我们翻译小组会及时改正。联系人：zzh_my@21867a440066f5335a8121cc （最后修改日期：2010-3-17） 。 翻译组员： 耿长江(武汉大学) 何丽娜(同济大学) 王晶(武汉大学) 焦龙(河北省地震局) 赵忠海(中国地震局地震研究所)
李健(解放军信息工程大学) 周命端(武汉大学) 统稿：赵忠海 校正：
A
学习 GAMIT 的路上，不再孤单
目录索引题记 .................................................................................................................................................A 目录索引.......................................................................................................................................... B 1．文档结构与命名规则.................................................................................................................1 1.1 介绍.....................................................................................................................................1 1.2 测站有关文件....................................................................................................................1 1.3 工程或时段文件.................................................................................................................3 控制文件...........................................................................................................................4 用于观测模型的文件.......................................................................................................6 清洗数据时使用的文件.................................................................................................10 估计参数时用到的文件.................................................................................................10 1.4 全局文件..........................................................................................................................12 1.5 文件汇总..........................................................................................................................16 2. 创建数据输入文件....................................................................................................................18 2.1 介绍和文件组织..............................................................................................................18 2.2 准备坐标文件(L-) ...........................................................................................................20 2.3 创建台站信息文件..........................................................................................................22 2.4 创建方案文件..................................................................................................................2

5 2.5 使用makexp .....................................................................................................................26 2.6 通过外部星历来创建G-和T-文件。 ..............................................................................28 2.7 创建卫星钟(J-)文件 ........................................................................................................31 2.8 运行 makex .....................................................................................................................33 2.9 X-file的描述 .....................................................................................................................36 2.10 生成站点钟差(K-和I-)文件 ..........................................................................................38 3. 批处理........................................................................................................................................41 3.1 引言 .................................................................................................................................41 3.2 运行fixdrv.........................................................................................................................42 3.3 执行批处理.......................................................................................................................58 3.4 评定结果..........................................................................................................................69 4 数据的编辑.................................................................................................................................85 4.1 介绍..................................................................................................................................85 4.2 使用autcln进行自动编辑 .................................................................................................86 4.3 扫描残差来发现周跳......................................................................................................97 4.4 用cview进行人机交互编辑数据 ...................................................................................102 4.5 策略编辑.........................................................................................................................112 5. 分步运行各个模块..................................................................................................................115 5.1 介绍.................................................................................................................................115 5.2 运行arc ..........................................................................................................................115 5.3 运行model......................................................................................................................118 5.4 运行 cfmrg ......................................................

...............................................................120 5.5 运行solve .......................................................................................................................121 5.6 运行 autcln....................................................................................................................125学无止境 B
学习 GAMIT 的路上，不再孤单
Autcln 命令文件..........................................................................................................126 （1）基于伪距数据的测站与卫星的钟差估计。 .....................................................132 （2）基于相位数据的测站与卫星的钟差估计。 .....................................................132 (3) 默认紧接着的步骤是加入偏差标记到所有数据中的跳跃的间隙中。.............132 （4）默认紧接着的步骤是扫描所有连续单差数据和形成双差并且用三差检查是否 在双差中有周跳。.......................................................................................................132 （5）剔除数据中的周跳，要使用尽可能多的观测值。 .........................................133 （6）数据整理。.........................................................................................................133 （7）标准点（可选） 。...............................................................................................134 钟统计...........................................................................................................................134 扫描总结.......................................................................................................................134 单程（one-way）残差统计 .........................................................................................135 数据编辑与偏差标记重命名.......................................................................................135 高度角统计...................................................................................................................135 偏差标记报告...............................................................................................................135 编辑报告.......................................................................................................................135 DD_Trials行..................................................................................................................135 偏差标记行（总是输出）...........................................................................................135 6. 自动批处理..............................................................................................................................136 6.1 概述................................................................................................................................136 6.2 控制文件.....................

...................................................................................................138 6.3 使用sh_gamit .................................................................................................................141 6.4 使用sh_glred..................................................................................................................143 7 对流层延迟模型.......................................................................................................................146 7.1 大气延迟描述................................................................................................................146 7.2 大气延迟算法................................................................................................................147 7.3 估计天顶延迟参数........................................................................................................149 7.4 估计梯度........................................................................................................................150 7.5 提取可降水估值............................................................................................................151 8. 功能程序和辅助表..................................................................................................................153 8.1 利用GAMIT解算结果进行绘图和统计.......................................................................153 8.2 从NGS ARGOS文件创建RINEX或FICA文件 ............................................................156 8.3 从FICA文件创建RINEX文件 ......................................................................................158 8.4 从C文件创建X文件 ......................................................................................................159 8.5 从X文件创建RINEX文件.............................................................................................159 8.6 创建和维持基准，时间，航天器，星历表 ................................................................159 航天器参数(svnav.dat) .................................................................................................160 TAI-UTC(leap.sec)........................................................................................................161 TAI-UT1 表：..............................................................................................................162 极移表：.......................................................................................................................162 月亮表：.......................................................................................................................163 太阳表：.......................................................................................................................163 章动表：...................................................................................................

....................163 附录 ..............................................................................................................................................164学无止境 C
学习 GAMIT 的路上，不再孤单
A.1. 天线技术规范...............................................................................................................164 A.1.1 介绍...................................................................................................................164 A.1.2 TI 4100 天线..................................................................................................164 A.1.3 天宝（Trimble）天线 ....................................................................................164 A.1.4 Rogue天线 ......................................................................................................164 A.1.5 阿斯泰克（Ashtech）天线...........................................................................164 A.1.6 Macrometer 天线 .........................................................................................164 A.1.7 SERCEL 天线..................................................................................................164 A.1.8 莱卡（Leica）天线 ........................................................................................164 A.2. 数据交换格式说明.......................................................................................................164 A.2.1 RINEX..............................................................................................................164 A.2.2 FICA .................................................................................................................164 A.2.3 导航文件...........................................................................................................165 A.2.4 SP3 轨道文件..................................................................................................165 A.3 由于SA政策模拟卫星钟变化.......................................................................................166 参考文献.......................................................................................................................................167 9.部分资料翻译............................................................................................................................170 9.1 如何创建21867a440066f5335a8121cc......................................................................................................170 9.2 sh_gamit自动处理脚本命令参数 ..................................................................................170 9.3 gamit安装后检验安装正确性........................................................................................173 9.4 关于安装gamit用到的编译器gfortran..........................................................................174 9.5 sh_sp3fit使用帮助 ........................................

..................................................................174 10 常见问题.................................................................................................................................176 11 常用网址.................................................................................................................................177 12 学习资料.................................................................................................................................178 13 本文档勘误及更新发布网址.................................................................................................178 14 后记.........................................................................................................................................178
学无止境
D
学习 GAMIT 的路上，不再孤单
1． 文档结构与命名规则1.1 介绍一个工程（译者注：TGO 中为项目）的各种文件命名都要遵守明确的规则， 从而使所有程序模块顺利读取这些文件。这种统一的文件命名规则，方便数据文 件管理、易于交互处理、发现问题并解决它。四种类型文件有： 1）与测站有关文件 2）与时段或工程有关的文件 3）全局性的文件 类型 1、类型 2 文件通过文件名中第一个特定字符区分每一种文件，类型 3 有一 个唯一的文件名。 类型 1 文件名中其它字符由 4 字符的测站名与观测年纪日组成。 类型 2 文件名中其余字符由 4 字符的工程名（计算者指定）与年纪日组成。类型 3 文件有一个特定的文件名，基本上是通过链接的方式存在于工程中的。上述命 名规则使得该软件能够按既定处理顺序处理海量数据。 下面三节描述了每种文件类型的内容与格式， 以及这些文件如何被创建及被那个 程序使用。1.5 节有所有类型文件的列表，是按字母排列，还给出了一个表格， 展示了 GAMIT 模块如何读写这些文件。
1.2 测站有关文件RINEX obs 文件：观测数据文件中包含了 L1、L2 载波相位、伪距、信号振幅、 初始测站坐标、天线偏移、起止时刻和由接收机通道跟踪到的卫星编号。 命名 例子 ：sitedayn.yyo ： vndn0020.87o. 数 据 来 自 VNDN 测 站 (1987 年 第 2 天 的 Vandenberg0 测站全天数据). 类型 创建 ：ASCII 码(译者注：能用记事本能打开的文件) ：由随机软件、AIUB 或 UNAVCO(teqc)读取从接收机下载的原 始文件（二进制文件）转化而成。 输入学无止境
：makex 和可选的 cview1
学习 GAMIT 的路上，不再孤单
RINEX 气象文件：在测站采集的气象数据。 文件名 例子 ：sitedayn.yym ： vndn0020.87m. 数 据 来 自 VNDN 站 ( 在 1987 年 第 2 天 的 Vandenberg0 站). 类型 创建 输入 ：ASCII ：程序读取下载自气象传感器的原始(二进制)文件。 ：model 和 sh_met_util
X-file
：GAMIT

观测文件，与 RINEX 文件相似，对于所有测站文件来 说，除了在开始和结束时段有不同(带有空历元的观察获得不到 此文件)。
文件名 ：xsitey.day 例子 类型 创建 输入 ：xvndn7.002。1987 年第二天的测站 VNDN 数据。 ：ASCII ：makex 和 可选实用程序 ctox ：makek 和 fixdrv;可选 bxtot 和 cview
C-文件
： 数据分析主要文件， model 从 X-文件创建， 由 用于输入 autcln,cview 和 solve；包含观测文件(O’s)、近似残差(O-C’s,观测-计算值)、偏导 数和附属信息。
文件名 例子 注释
：csitey.day ：cvndn7.002 ：直接和 X-文件相对应，区别是二进制的并带有偏导数。假如周 跳需要用 cview 手动修复，为了进一步处理可用 ctox 把干净的 c文件转化为 x-文件，尽管现在这很少见了。
类型 产生于 输入
：二进制 ：model,autcln,cview ：autlcn,scandd,cview,solve
学无止境
2
学习 GAMIT 的路上，不再孤单
K-文件
：利用近似（名义）坐标、广播星历和伪距通过 makex 或 makek 计算的接收机钟数据。 fixdrv 用于估计在时段期间钟行为的线性 被 或三次多项式模型的系数。
文件名 例子 注释 类型 产生于 输入
：ksitey.day ：kvndn7.002 ：与测站有关的 K-文件的参数、格式将在 2.9 节说明。 ：ASCII ：makex,makek ：fixdrv 和实用程序 calck 和 plotk.
P-文件 文件名 例子 注释 类型 产生于
：打印输出的文件，提供 model 的运行运行记录。 ：psitey.day ：pvndn7.002 ：直接和 X-文件和 C-文件相对应。 ：ASCII ：model
Z-文件
：记录 model 操作的文件，提供给外部程序一个完整记录，包括大 气值和处理过程中使用的模型。
文件名 例子 注释 创建 输入
：zsitey.day ：zvndn7.002 ：仅仅在 sestbl. 文件中的 Output net file = Yes 才记录。 ：model ：sh_metutil
1.3 工程或时段文件在一个工程(实验)中的特定时段(天)或几天来说，这些文件是相同的。它们以文 件或链接的形式总是存在每天目录下。唯一（版本）的文件经常存储在/tables 目学无止境 3
学习 GAMIT 的路上，不再孤单
录下提供给工程并通过链接使所有年积日文件夹共享这个文件。
控制文件处理控制文件 ：被 sh_gamit 使用；这个文件中包括了目录文件名和一些处理 控制。参见第 10 章。 文件名 ：process.defaults 类型 创建 输入 ：ASCII ：用户从模板创建 ：sh_gamit
站处理控制文件
：由 sh_gamit 使用；文件中包括了使用测站及控制情况。
文件名 ：sites.defaults 类型 创建 输入 ：ASCII ：用户从模板创建 ：sh_gamit
时段控制表 ：输入 fixdrv 的控制文件，详细说明了分析计算的类型、先验误差 和卫星约束。参见 3.2 节 文件名 ：sestbl.
注释 ：该文件是硬链接（译者注：把/tables 文件下的该文件链接到

了年 积日文件夹下） ，链接可以用不同的类型（译者注：软链接） 。 类型 ：ASCII 创建 ：用户来自模板 输入 ：fixdrv
站控制表
：输入 fixdrv 控制文件，详细说明了每个站的先验坐标约束和可选 的钟、大气模型。参见 3.2 节
文件名 ：sittbl. 注释 ：文件是硬链接，同时段控制表 sestbl.。学无止境 4
学习 GAMIT 的路上，不再孤单
类型 ：ASCII 创建 ：用户来自模板 输入 ：fixdrv
时段信息或方案文件：要处理的卫星和时间(2.5 节) 文件名 :21867a440066f5335a8121cc
类型 ：ASCII 创建 ：用户或 makexp / sh_makexp 输入 ：makex，fixdrv
D-文件 ：Fixdrv 文件——定义在每个实验中的时段号，每个时段的接收机编号， 坐标(L)文件，星历(T)文件，站钟(I-)、卫星钟(J-)和每时段数据文件(X-或 C-)，还有各时段处理顺序。 文件名 ：dxxxxy.day 其中 xxxx 是工程(项目)文件名。文件名其它字母是
任意的，但第六个字母默认是年。其余 3 个字母是工程开始的年积 日，但这不是必须的。 例子 ：dcalf7.002 注释 ：D-文件主要输入 fixdrv。因此，它的文件名定义了所有随后工程 文件名。 所以用户能在任何情况下识别出工程名。 在执行 fixdrv 前， 用户用 makexp 程序或手动创建 D-文件。参见 3.2 节。 类型 ：ASCII 创建 ：makexp 或用户 输入 ：fixdrv,grdtab,autcln B-文件 ：主要批处理文件——控制数据批处理(自动)模式。 文件名 ： bxxxxx.bat 其中 xxxxx 是 D-文件名除首字母外的 5 个开始字符
例子 ：bcalf7.bat 注释 ：主 B-文件包含从属 B-文件顺序，从属 B-文件包含着模块 fixdrv 指 定的文件执行顺序。它的文件名与 D-文件相对应。
学无止境
5
学习 GAMIT 的路上，不再孤单
类型 ：ASCII 输入 ：fixdrv
B-文件
：从属批处理文件——控制程序模块的执行。 ： bxxxxx.nnn 其中 xxxxx 是 D-文件名除首字母外的开始 5 个字符，
文件名
nnn 是批文件处理顺序号 例子 ：bcalf7.001……bcalf7.015 注释 :每个从属批处理文件包含程序模块的执行输入流。例如， dbcalf7.bat 的第一行可能是 arc < bcalf7.001。这就是说程序模块 arc 会接受来自 bcalf7.001 的指令。 类型 输出 输入 ：ASCII ：fixdrv :arc,model,cfmrg,solve
用于观测模型的文件站信息文件 ：接收机、天线和每个时段的观测时间(参见 2.3 节) 文件名 ：21867a440066f5335a8121cc
注释 ：根据站日志准备(替换掉不再使用的 hi.raw 和 sited.文件。) 类型 ：ASCII 创建者 ：用户，可选用 make_stnfo 或 conv_stnfo,或者 SOPAC 创建. 输入 ：makex,fixdrv,model
导航(或者 E-)文件 ：RINEX 或 FICA 格式的广播星历数据。被 makej、makex 和 makek 用于创建卫星和接收机钟文件，也可被 bctot 创建初始 G-文件和(或)T-文件(通过 IGS 产品来做这个很少见了)

文件名 例子 注释学无止境
：brdcdayn.yyn(或者 sitedayn.yyn esitey.day) ：vndn0020.87n 或 evndn7.002 ：在附录 2 中描述了这些文件的参数和格式。6
学习 GAMIT 的路上，不再孤单
类型 输出自 输入
：ASCII ：RINEX 转换器或者实用工具 ficachop ：makej,makex,makek,bctot
G-文件 文件名 例子 注释
：在 T-文件中所有卫星的轨道初始坐标文件 ：gxxxxy.day ：gigsf5.065 ：G-文件包含着特定 UTC 历元的每个卫星初始坐标和无重力参 数。 G-文件的初始参数在数值积分卫星轨道时用于起始点,并生成 T-文件。g-文件的文件名是任意的，但通常指明轨道来源信息和 并包含初始参数的年、 天信息， 例如上面例子(2005 年， 65 天)； 第 在任何情况下它都应和 T-文件相匹配。与 solve 创建的 G-文件文 件名是相同的，除了第 6 个字符被增加了个字母。g-文件的格式 在 2.10 节有描述。
类型
：ASCII
输出自 ：orbfit / sh_sp3fit.bctot,solve,ttog 输入 I-文件 ：arc ：记录着每个站、每个时段的测站钟偏、速率和加速度，由 model 可选使用。 文件名 ：ixxxxx.xxx 例子 注释 类型 ：icalf7.002 ：I-文件的文件名由 D-文件指定 ：ASCII
输出自 ：fixdrv 输入 J-文件 ：model ：卫星传送的卫星钟参数，由接收机记录 文件名 ：jxxxxy.day 例子 注释 ：jcalf7.002 ：该文件用于 model 计算每个历元的接收机改正，当所有测站是
学无止境
7
学习 GAMIT 的路上，不再孤单
非同步观测时并且存在大的卫星钟漂(例如：在美国 SA 政策下)的 情况下，用于相位模型改正。 类型 ：ASCII
创建自 ：makej 输入 ：model
L-文件
：站坐标文件——包含着在一个工程中所使用测站的可获得的、 最好的坐标列表。 文件名 ：在工作目录中 lxxxxx.day,其中 xxxx 是工程的文件名，来自于 D文件名。 例子 注释 ：lsv4x7.002 ：L-文件可以是球坐标格式(老 GAMIT 风格)，仅有坐标。也可是 笛卡尔坐标格式(与 GLOBK apr 文件相同)，有坐标和可选的速度 值。L-文件由 solve 更新。 类型 ：ASCII
创建自 ：globk/glorg,gapr_to_l,sh_rx2apr,solve,tform 输入 :makex,makek,fixdrv,model,grdtab,tform
T-文件
：所有卫星在某一时段或时段序列的星历表（T）文件——包含 相等时间间隔(对于 arc 默认 15 分钟)的卫星状态矢量，用于后面 model 模块插值。它的文件名应与 G-文件匹配。 文件名 ：txxxxy.day 例子 ：tigsf4.289,该文件与利用 IGS 最终 sp3 文件生成的 g-文件和 2004 年第 289 天的初始坐标有关。 类型 ：二进制
输出自 ：arc,bctot,ngstot,sh_bctot,sh_sp3fit 输入 ：fixdrv,model,ttongs,ttoasc
A-文件
：ASCII 版本的 T-文件，可选生成，用于分析检查或信息输出。
学无止境
8
学习 GAMIT 的路上，不再孤单
文件名 ：Axxxxy.day 类型 ：

ASCII
输出自 ：ttoasc 输入 ：无
U-文件
：在解算中用到的每个测站的海洋、大气负载和气象资料，由全 球格网(例如：otl.grid,atlm.grid,met.grid)内插生成或对于具体测站 从文件列表中读取其值(例如：IGS 网，otl.list,atml.list,met.list)。
文件名 ：uxxxxy.day 类型 ：ASCII
输出自 ：grdtab 输入 ：model
Y-文件
：卫星偏航（Yaw）文件，有两种类型。ASCII 类型(第 6 个字符 是年的最后一位)是在观测时段中对于每个被观测的卫星给出日 食(月食)时间、偏航速率，由 arc 模块读取全球 svnav.dat 文件计 算偏航信息，然后计算生成。二进制类型(第六个字符是 t)在时段 的观测历元中给出偏离正常航道的角度，利用 T-文件和 ASCII 类 型 y-文件作为输入文件，由 yawtab 计算生成。这两种文件类型将 在 5.2 节中讨论。 文件名 ：yxxxxy.doy 例子 类型 ：ypgga5.267 ：ASCII 输入于 ：model
输出自 ：arc 输入于 ：yawtab 文件名 ：yxxxxy.doy 例子 类型 ：ypggat.267 ：二进制
输出自 ：yawtab
学无止境
9
学习 GAMIT 的路上，不再孤单
清洗数据时使用的文件Autcln(详细)输出文件 ： 编辑处理的完整记录； 假如解算成功该文件可以忽略并删
除。参见 4.2 和 5.6 节 文件名 ：autcln.out
类型 ：ASCII 创建 ：autcln
Autcln 总结文件 ：编辑后的总结文件,评价结果时很有用。参见 4.2 节和 5.6 节 文件名 ：autcln.pref.sum 或者 autcln.post.sum
类型 ：ASCII 创建 ：autcln V-文件 ： scanm 的打印输出文件(scanrms 取两者之一， 包含 rms 汇总和每一个
双差组合的 jumps 值 文件名 ：vxxxx1.day.sort,vxxxxa.day.sort,vxxxx1.day.worst,vxxxxa.day.worst
类型 ：ASCII
估计参数时用到的文件M-文件 ：合并文件 —— 为进行 solve 最小二乘分析而建立的数据和参数。 ：mxxxx1.day,mxxxxa.day
文件名
例子 ：mcalfa.002 注释 :M-文件名由 D-文件名衍生而来。M-文件最初由 cfmrg（用 fixdrv 写的
批处理文件）创建，其建立了 Solve,autcln,或 cview 最初的运行设置。在估计调整参数后，solve 用相同的文件名写了一个新 M-文件，并包含调整参数。通常 fixdrv 生成处理 顺序，来自初始解的 M-文件加上.autcln 后缀来重命名，并可用于 cview 查看残差，用于 autcln 进行在编辑。
类型
：二进制
输出自 ：cfmrg 输入学无止境
：solve,autcln,scandd10
学习 GAMIT 的路上，不再孤单
N-文件
：由 autcln 进行 postfit 编辑后，包含站信息、与高度有关值的噪声 文件在 solve 中用于重新加权相位观测值。在 10.2 版本的 LC_AUTCLN 模 式下，它也包含由 autcln 调整的双差模糊值。
格式 例子 注释 文件输入 类型
：nxxxx1.day nxxxxa.day ：ncalf1.002,ncalfa.002 ：该件包含误差模型结构：solve 批处理文件部分。solve 读取并重

写批
：ASCII
输出自 ： Shell 脚本 sh_sigelv 输入 Q-文件 ：solve ：打印输出 solve 运行文件——包含了分析记录。
文件名 ：qxxxx1.day,qxxxxa.day 例子 注释 ：qcalfp002,qcalfa.002 ： “p” 版本的 Q-文件由最初(先验)解生成， “a“版本的 Q-文件由最后(后 验)解生成。 类型 ：ASCII
输出自 :solve
O-文件
：solve 运行的解输出文件，也是 Q-文件的简略形式,主要用于绘图、统
计和输入 sh_met_util. 文件名 ：oxxxx1.day oxxxxa.day
例子 ：ocalf1.002,ocalfa.002 类型 ：ASCII 输出自：solve 输入 ：网平差，统计和绘图程序
H-文件
： 在松弛约束下生成的协方差矩阵和解的参数平差值， 用于输入 GLOBK ：hxxxxy.day
文件名
例子 ：hcalf7.002 类型 ：ASCII学无止境 11
学习 GAMIT 的路上，不再孤单
输出自：solve 输入 ：htoglb,htosnx,htoh
1.4 全局文件在一定意义上讲，这些文件是全局的，他们可以在许多相隔一段时间的工程项 目(至少一年)中仍然有效。这些文件的文件名必须是下面明确指出的。除 ftp_info 文件外，所有文件都可在 gg/tables 中找到，必须把这些文件链接(或拷贝)到每天目 录下。
ftp_info 类型 创建 输入
：一个信息表，包含了从网络下载文件要用到的网络地址和协议 ：ASCII ：MIT,使用者按需要修改 ： sh_get_hfiles、 sh_get_nav、 sh_get_orbits、 sh_get_rinex、 sh_get_stinfo、 sh_update_eop
rcvant.dat
：一个信息表，包含了天线和接收机的 GAMIT 6 字符编码，该编码 与 RINEX 和 SINEX 文件中全名(20 字符)对应。
注释 类型 创建 输入
：参见 4.3 节和附录 7 ：ASCII ：由 MIT、SIO 或者用户根据 IGS 标准创建 ：model,sh_upd_stnfo/mstinf2
guess_rcvant.dat ：由 sh_gamit 可选使用，从 RINEX 头文件中不精确的接收机和天 线类型的 20 字符名称来确定其 GAMIT 编码 类型 创建 输入 ：ASCII ：MIT,可按使用者需要修改。 ：sh_upd_stnfo/mstinf2
antmod.dat 格式学无止境
：天线相位中心偏移，作为计算高度角和方位角的函数变量 ：有两种版本可以，一种是 GAMIT 特有的，但很快就要被废除的，12
学习 GAMIT 的路上，不再孤单
另一种是 IGS ANTEX 注释 类型 创建 输入 ：参见 2.3，3.2 节 ：ASCII ：IGS/MIT ：model
svnav.dat
：此表给出了导航卫星编号、区组号(I 或 II)、飞行器质量和每个卫星 的偏航参数(用 PRN 号列出)(译者注：PRN:伪随机号，pseudo random number)
注释 类型 创建 输入
：参见 8.6 节 ：ASCII ：MIT ：makex,arc,model
svs_exclude.dat：该表给出处理过程中某天（年+年积日）应该排除的问题卫星， 类型 创建 输入 ：ASCII ：MIT,SOPAC,用户 ：sh_sp3fit/ orbfit
gdetic.dat 注释
：大地水准面参数表（注：主要为椭球参数） ：该文件的格式在 7.8 节有描述。此文件可以

在任何时候增补来包含 另外的水准面。
类型 创建 输入
：ASCII ：MIT,用户 ：tform
ut1. 格式 注释 类型学无止境
：UT1 表以表格的形式包含了 TAI-UT1 值 ：ut1.——在工作目录，例如在 tables 目录中的 ut1.iers ：要定期更新。参见 8.6 节 ：ASCII13
学习 GAMIT 的路上，不再孤单
创建 输入
：MIT,SOPAC,或者用户从 IERS、USNO 等用 sh_update_eop 创建 ：arc,model,sh_sp3fit/orbfit,ngstot,bctot,ttongs
pole.
：极表以表格的形式包含极移值，用于在 model、arc 和 bctot 中内插 值
格式 注释 类型 创建 输入
：pole.——在工作目录，例如在 tables 目录中的 pole.iers ：要定期更新。参见 8.6 节 ：ASCII ：MIT,SOPAC,或者用户从 IERS、USNO 等用 sh_update_eop 创建 ：arc,model,sh_sp3fit/orbfit,ngstot,bctot,ttongs
leap.sec 格式 注释 类型 创建 输入
：自从 1982 年 1 月 1 日在 TAI-UTC 的跳(跳秒)表 ：leap.esc ：参见 8.6 节 ：ASCII ：SIO,MIT,或者用户从诸如 IERS、USNO 公告中创建 ：fixdrv,model,arc,bctot,ngstot,ttongs
nutabl.
：章动表——以表格形式包含了章动参数，这些参数用于在惯性系统 和地固系统之间转换。
格式 注释 类型 创建 输入
：nutabl.——在工作目录；例如在 tables 目录中的 nutabl.91 ：参见 8.6 节 ：ASCII ：MIT ：arc,model,sh_sp3fit/orbfit,ngstot,bctot,ttongs
luntab. 格式 注释
：月亮星历表 ：luntab.——在工作目录；例如在 tables 目录中的 luntab.95.J2000 ：惯性参考系统(B1950 或 J2000)必须与输入控制匹配（注：覆盖工程 时间长度） 。用 21867a440066f5335a8121cc 或 sh_links.tables 链接。参见 8.6 节。
学无止境
14
学习 GAMIT 的路上，不再孤单
类型 创建 输入
：ASCII ：MIT ：arc,model
soltab. 格式 注释
：太阳星历表 ：soltab.——在工作目录；例如在 tables 目录中的 soltab.95.J2000 ：惯性参考系统(B1950 或 J2000)必须与输入控制匹配（注：覆盖工程 时间长度） 。用 21867a440066f5335a8121cc 或 sh_links.tables 链接。参见 8.6 节。
类型 创建 输入
：ACSII ：MIT ：arc,model
Otl.gird,otl.list :来自全球格网或各测站的海洋潮汐参数 类型 创建 输入 ：二进制(格网)或 ASCII(表列) ：在 Onsala 空间天文台的 Hans-Georg Scherneck ：grdlab
atl.grid,atl.list ：来自全球格网或个测站的大气潮汐参数 类型 创建 输入 ：二进制(网格)或 ASCII(表列) ：澳大利亚国家大学的 Paul Tregoning ：grdtab
atml.grid,atml.list 类型 创建
:来自全球格网或各测站的无潮汐的大气负载参数
：二进制(格网)或 ASCII(表列) ：在澳大利亚国家大学的 Paul Tregoning,其中数据由卢森堡大学的 Tonie van Dam 提供。
输入
：grdtab
map.grid,map.list 类型学无止境
:基于全球模型的大气映射函数系数。目前仅由 VMF1 提供。
：二进制(格网)或 ASCII(表列)15
学习 GAM

IT 的路上，不再孤单
创建
：澳大利亚大学的 Paul Tregoning,其中数据由维也纳科技大学的 Johannes Boehm 提供。
输入
：grdtab
1.5 文件汇总A - 文件：二进制的 T- 文件(星历表) B - 文件：数据批处理模式的控制文件 C - 文件：观测值-计算值，偏导值 D - 文件：时段和接收机驱动文件 E - 文件：广播星历，RINEX 格式导航文件或 FICA BLK 9 格式文件 G - 文件：轨道初始参数和无重力参数值 H - 文件：输入 GLOBK 的平差值和完整的方差-协方差矩阵 I - 文件：接收机钟多项式参数 J - 文件：卫星钟多项式系数 K - 文件：在观测期间，来自伪距的接收机钟偏差值 L - 文件：站坐标 M - 文件：控制合并数据(C)文件，用于解算和编辑项目 N - 文件：重写数据权重，用于 solve,由 autcln.sum.postfit 创建 O - 文件：用于后处理分析的分析记录(Q-文件的简化形式) P - 文件：model 运行记录 Q - 文件：(运行 solve 的)分析记录 S - 文件：不再使用 T - 文件：星历表 U - 文件：用于 model 的加载和气象数据 V - 文件：SCANRMS 的编辑输入 W - 文件：RINEX 格式气象数据文件 X - 文件：输入观测值 Y - 文件：卫星原始参数 Z - 文件：输出气象数据学无止境 16
学习 GAMIT 的路上，不再孤单
学无止境
17
学习 GAMIT 的路上，不再孤单
2. 创建数据输入文件
2.1 介绍和文件组织 在解算 GPS 数据中第一步同时也是最繁琐的一步是收集测站信息和接受机输 出数据，通过这样一个过程使得数据可以被程序有效的处理。在这一过程中，你必 须尝试性的确定解算的天数，包括哪些台站，多大的时间跨度，采样频率。总而言 之，这阶段你要详细查看站点记录并设计解算策略。这时你可能发现缺少某些测站 的日志、数据文件，或者先验坐标，你必须耐心的通过电子邮件等手段解决这些问 题。 主要的 GAMIT 模块(由 fixdrv, arc, model 开始)需要七种类型输入文件： （1）ASCII 格式的 X-files 文件，包含原始相位、伪距数据(每一测站的每一时 段需要一个文件)。 （2）L-file 格式的台站坐标。 （3）每个测站接受机和天线信息(21867a440066f5335a8121cc 文件) （4）卫星列表和解算方案(21867a440066f5335a8121cc) （5）卫星轨道初始参数，在 G-file 中(或者 T-file 文件的星历列表) （6）卫星和测站的时钟值(I-，J- and K-files) （7）解算控制文件(sestbl. 和 sittbl.) （8） “标准”表格形式的月亮、太阳星历表，地球转动文件，大地水准面文件， 卫星和仪器的信息文件(下面有介绍)。 X-file 是组织结构的关键，因为对于同一时段的 X-file 文件有相同的开始、结 束时间、选择的卫星和采样间隔。这一强加的硬性规定有很多优点。最主要的一个 是，建立 X-files 的过程(makex 程序)就像

一个过滤器，能够捕捉到大多数问题。例 如：丢失的或者无效的数据、不匹配时间标记、工作质量低的接收机时钟，如果采 用了有问题的文件，当发现的时候，已经浪费了大量时间。
学无止境
18
学习 GAMIT 的路上，不再孤单
第一步是为数据处理创建工作目录。建议的组织形式是(sh_gamit 命令采用的) 采用工程项目，或者顺序编号的观测网目录如(例如： emed98, scec99, scign99)。 在这一目录下面，你应该有“day”目录(e.g. /312 for day 312).，对应每一天或者一时 段。平行于 day 目录的目录是存放 RINEX 文件的目录(例如： /rinex )，存放卫星轨 道文件的目录,包括 /brdc 存放导航文件 （译者注： 广播星历） /igs 存放 IGS SP3 和 文件（译者注：IGS 精密星历） ，存放工程相关的 GAMIT 表目录( /tables )。GLOBK 处理的目录/gsoln 和 /glbf,也可以在这一层。应用 sh_gamit 批处理时，如果它们以 前不存在，这些目录大多数被自动建立。 你开始预处理通过在 day 目录内建立链接文件，链接数据文件和表格，这些数 据和表格都是批处理需要的。这些通常在两步内完成。第一步，建立链接文件，链 接 GAMIT 的全球文件：大地测量数据(gdetic.dat),月亮和太阳星历(luntab. 和 soltab.), 章动表 (nutabl.), 地球旋转 (ut1. 和 pole.), 海潮改正模型 (stations.oct and grid.oct 译者注： 这两个已经不在试用， 现在用的是 otl.grid)， 跳秒值 (leap.sec), 和 卫星数目编号信息 ，接收机和天线特征 (svnav.dat, antmod.dat, rcvant.dat). 这些文件在/tables 目录里面, 执行脚本命令 links.tables(译者注：现在已经没有这 个命令了，这是以前的版本。), 该命令将会为全球文件建立连接文件，链接到 ~/gg/tables. 第二步，在每天的目录下执行 links.day(如在 example 文件夹中建立 037 文件夹，在 037 目录中运行 links.day 2000 037 jiao)，该命令将会把../tables 中的全球 文件和项目必须文件链接到当天目录下。Links.day 也会在 day 目录中建立../tables 中 6 个 gamit 所需要的控制文件链接：21867a440066f5335a8121cc、 21867a440066f5335a8121cc、 sestbl.、 sittbl.、 autcln.cmd 和 lfile. 。因此，你就可以保证/tables 中文件的唯一复本，避免处理不同 时段数据而采用了不同的输入控制。如果你只是对一天中的数据做一个试验，你可 以执行 21867a440066f5335a8121cc，来链接~/gg/tables，实现复制所有必须文件到该 day 目录中。如果 你使用 sh_gamit 命令，这两个标准链接任务会被自动完成。 我们将在 2.5 章讲述 RINEX 文件的链接，在 2.2-2.4 节讲诉如何准备 L-file、 21867a440066f5335a8121cc、21867a440066f5335a8121cc 文件，在第三章讲诉 fixdrv 所需的控制文件(sestbl. 和 sittbl.)。导航文件仅仅是一个 RINEX 格式的“n”文件，命名为 aaaannnn.YYn 或者 eaaaY.DDD，需要

从/rinex 或者/brdc 文件夹被链接或复制到 day 目录中，或者直接 从 IGS 数据中心下载。学无止境 19
学习 GAMIT 的路上，不再孤单
2.2 准备坐标文件(L-) L-fiel 包括试验（工程）中需要的所有的测站的坐标文件。对于 10.1 版本，有 两种格式的 L-file 被支持： （1）观测历元的球坐标(传统的 GAMIT 格式)， （2）包含 笛卡尔坐标、指定历元速率格式(GLOBK apr-file 格式)。后者更受欢迎，因为它容 许测站在观测期间中运动，并且可以更好的辨别到由于地震引起位移的测站。(参考 GLOBK 3.1 节定义地震和重命名测站)* nafd_plate_scec.apr : itrf00_noam + updates from vel_020123a * North American stations for stabilization * GOLD_GPS -2353614.1450 -4641385.3890 3676976.4750 -0.00216 0.00649 0.00457 1997.0000 MOJ1_GPS -2356424.5422 -4646613.6634 3668462.2248 -0.00216 0.00649 0.00457 1991.0600 MOJ1_GLA -2356424.5553 -4646613.5858 3668462.2288 -0.00216 0.00649 0.00457 1998.9610 GOLD_GHT -2353614.1949 -4641385.3488 3676976.4648 -0.00216 0.00649 0.00457 2000.6690
在这个例子中，有多行输入，前面 4 个字符 GOLD 和 MOJ1 表示测站代码， 后四个字符代表测站表示：地震之前的观测（_GPS） ，Landers 地震之后的观测 (_GLA)， Hector Mine 地震之后的观测(_GHT)。这些代码以 GLOBK 中的 eq_file 文件中给定的地震数据（位置、日期、影响半径）为基础由 GLOBK 自动产生的。 决定用哪一个(地震之前的还是地震之后的)，GAMIT 可以读取 eq_file 中数据日期， 与要处理的数据日期做对比来确定。 每一行后面的日期， GLOBK 解的历元时间， 是 有当时的位置和速率。 一个老版本指定历元的 GAMIT L-file 是采用 GLOBK 的 apr 文件， 利用 gapr_to_l 程序产生的。 gapr_to_l 其中： 指输入的GLOBK a priori 文件的文件名
输出的 l-file 的文件名 是一个文件的文件名，这个文件给出了 GAMIT 中 4 个代码测站 的全称，你可以经常省略这一文件，通过双引号“”来代替，此时 GAMIT 测站全 称由 21867a440066f5335a8121cc 得到，不是由 l-file 得到的。 < date >是一个日期，l-file 文件坐标参考日期。这一日期的指定，可以采用下面 三种方法之一：学无止境 20
学习 GAMIT 的路上，不再孤单
1 如果给定一个值，被认定为十进制的年 2 如果给定两个值，他们被认为是年和这一年的第几天（年积日） 3 如果有三个值，他们被指定为年，月，日。 老版本的 GAMIT l-file 有下面的格式：Epoch 1993.5479: From file itrf00.aprGOLD GOLD GPS N35 14 36.99253 W116 53 21.29202 6371978.7987 MOJ1 MOJ1 GPS N35 09 01.05171 W116 53 26.91387 6371920.4323 Updated from lscec3.282
这里只有前面的 4 个代码有意义， 后面的 12 个代码是站点的文件名， 只是用来 显示。

对于这种格式的 l-file，必须确保输入到 gapr_to_l 中的 apr 文件中，对于处 理历元的一个测站只能有一个测站位置。 The GLOBK 形式的 L-file 是自由格式，没有空格的第一行是一个说明。而旧 格式的 L-file 有一个严格的格式。 (A4,1X,A12,A1,I2,1X,I2,1X,F8.5,1X,A1,I3,1X,I2,1X,F8.5,F14.4) 没有提供清楚的头文件行或注释。因为在处理过程中，这一行的前 4 个字母与 测站名相匹配时，这一行才被读取，所以注释可以放在这个文件中，没有危害，只 要注释的前 4 个字母和测站的文件名不一样就可以了。 在每天的目录中，l-file 文件(old 版或者新版)文件名除了第一字母以外通常与 D-file 相同(例如 lprojy.ddd)（译者注： lprojy.ddd 是这一天中观测的所有测站坐标 组成的 lfile） ，该文件是链接到../tables 里 l-file 的链接，../tables 里面的 l-file 文件包 含所有的(工程所用到的)测站坐标。 利用脚本 sh_rx2apr，可以获得观测网络中新测站的坐标，这个脚本调用 svpos
和 svdiff 程序，利用 RINEX 文件，进行伪距迭代求解，如下所示sh_rx2apr -site -nav