本文转载http://www.gnss.help/2016/10/02/grids-in-gamit/

原作者：山科研究生jiangyingming@live.com

0、前言

在使用 GAMIT/GLOBK 程序进行数据处理时，你肯定已经注意到了其 tables 文件夹中存在许多 *.grid 文件。如果你检查过 GAMIT/GLOBK 安装目录的 tables 文件夹，你会发现更多这样的文件。它们多数体积都很大，一般为数百 MB。这样的文件称为格网模型文件，一般在 MIT 的 FTP 服务器上提供下载。

本文将试图逐个地介绍其中的格网模型文件，可以算是介绍其格网模型文件的一份文档。

1、atl.grid

全球的大气潮格网模型文件，由澳大利亚国家大学的 Paul Tregoning 创建。

该文件包含在 GAMIT/GLOBK 程序的安装包中且不需要更新。也就是说，不管你处理哪一年的观测数据，都可以引入这个文件。

要在 GAMIT 数据处理中应用大气潮格网模型，你首先需要保证在进行数据处理的工程文件夹的 tables 目录中，其下的 atl.grid 已经正确链接到 GAMIT/GLOBK 安装目录中的 tables/atl.grid 文件。然后在 sestbl. 中修改设置为：

Use atl.list = N 　　 ; Atmospheric tides, list file, not yet available
Use atl.grid = Y 　　 ; Atmospheric tides, grid file

2、atml.grid

全球的无潮汐大气负载参数格网模型文件，该文件由澳大利亚国家大学的 Paul Tregoning 根据卢森堡大学的 Tonie van Dam 提供的数据创建。

该模型一般每年对应一个文件，GAMIT 使用的该文件的命名方式一般为：atmdisp_cm.<year>。如用于 2015 年的文件名为 atmdisp_cm.2015。

要在数据处理中应用无潮汐大气负载参数格网模型，首先从 MIT 的 FTP 服务器下载适用于当年的模型文件放入 GAMIT/GLOBK 安装目录的 tables 文件夹内，然后确保在进行数据处理的工程文件夹的 tables 目录中，其下的 atml.grid 正确链接到当年的模型文件。最后在 sestbl. 中设置：

Use atml.list = N 　　 ; Atmospheric (non-tidal) loading list file from LU
Use atml.grid = Y 　　 ; Atmospheric (non-tidal) loading grid file from LU, converted to GAMIT format

3、gpt.grid

全球气压和温度模型文件，该模型由 Boehm 和 Schuh 提出。在估计对流层参数时，GAMIT 从该模型中内插获取观测所在地区的气压和温度。

该模型一般不会更新，并且已经被包含在 GAMIT/GLOBK 程序的安装包中。

要在数据处理中使用全球气压与温度模型，首先确保你进行数据处理的工程文件夹的 tables 目录中，其下的 gpt.grid 已经正确链接到 GAMIT/GLOBK 安装目录中的 tables/gpt.grid 文件。然后在 sestbl. 中设置：

Met obs source = GPT 50 　　 ; hierarchical list with humidity value at the end; e.g. RNX UFL GPT 50 ; default GTP 50

这样，在进行数据处理时，GAMIT 就将从该模型中获取气压和温度参数。但内插出的精度毕竟不如实测数据，因此如果你有气象观测数据的话，最好还是将其引入。所以我推荐如下的设置：

Met obs source = RNX UFL GPT 50 　　 ; hierarchical list with humidity value at the end; e.g. RNX UFL GPT 50 ; default GTP 50

如此设置后，需要测站所在地区的气压和温度数据时，GAMIT 程序将首先检查有没有对应的 RINEX 格式气象数据，如果没有的话再检查有没有对应的 U-文件，最后才会尝试使用全球气压与温度模型提供的数据。

4、map.grid

全球的大气映射函数模型文件，该文件由澳大利亚国家大学的 Paul Tregoning 根据维也纳科技大学的 Johannes Boehm 提供的数据创建。

GAMIT 程序目前支持 GMF、NMFW 和 VMF1 三种映射函数，但只有 VMF1 需要引入模型文件。该模型逐年更新，每年生成一个文件。你可以根据需要在 MIT 的 FTP 服务器上下载，其命名方式为：vmf1grd.<year>。

要在数据处理时应用 VMF1 映射函数模型，需要首先下载适用于当年的模型文件放入 GAMIT/GLOBK 安装目录的 tables 文件夹内。然后确保在进行数据处理的工程文件夹的 tables 目录中，其下的 map.grid 正确链接到当年的模型文件。最后在 sestbl. 中设置：

DMap = VMF1 　　 ; GMF(default)/VMF1/NMFH; GMF now invokes GPT2 if gpt.grid is available (default)
WMap = VMF1 　　 ; GMF(default)/VMF1/NMFW; GMF now invokes GPT2 if gpt.grid is available (default)
Use map.list = N 　　 ; VMF1 list file with mapping functions, ZHD, ZWD, P, Pw, T, Ht
Use map.grid = Y 　　 ; VMF1 grid file with mapping functions and ZHD

5、met.grid

虽然各版本的 GAMIT 一直保留了这个模型文件链接，但该模型并没有什么用。有 sestbl. 中对该项配置的说明为证：

Use met.list = N 　　 ; not yet supported
Use met.grid = N 　　 ; not yet supported

6、otl.grid

全球的海潮模型文件，GAMIT 现在支持多种海潮模型。MIT 的 FTP 服务器上提供的海潮模型文件有：

• otl_CSR4.grid

• otl_FES2004.grid

• otl_FES99.grid

• otl_GOT00.grid

• otl_NAO99b.grid

要在 GAMIT 数据处理中应用全球海潮模型文件，首先下载对应的海潮模型文件放入 GAMIT/GLOBK 安装目录的 tables 文件夹内。然后确保在进行数据处理的工程文件夹的 tables 目录中，其下的 otl.grid 正确链接到对应的模型文件。最后在 sestbl. 中设置：

Tides applied = 31 　 ; Binary coded: 1 earth  2 freq-dep  4 pole  8 ocean  16 remove mean for pole tide
　　　　　　　　　　; 32 atmosphere ;  default = 31
Use otl.list = N 　　　 ; Ocean tidal loading list file from OSO
Use otl.grid = Y 　　 ; Ocean tidal loading grid file, GAMIT-format converted from OSO

这里的 “Tides applied”，以二进制的形式表示解算中应用的选项。如果你知道 UNIX/Linux 系统中的权限位表示方法，这一项应该很容易理解。否则，可以简易地以数字的对应关系计算：

• 1 ：earth;

• 2 ：freq-dep;

• 4 ：pole;

• 8 ：ocean;

• 16：remove mean for pole tide;

• 32：atmosphere

对于你要在解算中应用的项，只需将其对应的数字的和作为配置即可。例如这里的 31：

31 = 1 + 2 + 4 + 8 + 16

因此数据处理中应用到的选项为：

• earth;

• freq-dep;

• pole;

• ocean;

• remove mean for pole tide;

7、补充

文行至此，已经将 GAMIT/GLOBK 程序中所有的 .grid 文件都已做过简要的介绍。这里以一个小技巧作为最后的补充。

或许你还发现，在 tables 文件夹内还存在一些与 .grid 同名的 .list 文件（除了 gpt.grid）。这些 .list 文件是从模型文件中提取出的所有 IGS 站点对应的参数值。如果你要处理的全都是 IGS 站的观测数据，完全可以使用这些 .list 文件，只需在对应的配置项中将使用 .list 文件的选项设置为 “Y”，并将使用 .grid 文件的选项设置为 “N”。