gamit10.71解算GPS文件以及获取ZTD和PWV（附下载链接）

       下面的实例中，处理的数据为2020年第5天的数据，4个远距离的IGS站（chum、cusv、daej、irkj）、解算其他3个站点。

gamit10.71下载链接（https://cloud.189.cn/t/yM3MVfVryENf （访问码：5l39））

1、数据处理前更新表文件

 更新方法在之前的博文里面gamit10.71更新tables的几个方法（使用FileZilla和Ubuntu的gftp软件）-CSDN博客。

表文件的网址如下

         ftp://garner.ucsd.edu/pub/gamit/tables/

         ftp://everest.mit.edu//pub/GRIDS/

每次用更新（可能有的并不用每次吧）

         pole.usno        极移表

         ut1.usno         地球自转参数表

         svs_exclude.info  坏卫星信息文件

         svnav.dat        卫星列表 

         rcvant.dat       接收机天线对照表       

         antmod.dat       天线相位中心改正正表  

         dcb.dat          伪距改正参数

每年更新一次（按对应年份使用）

         luntab.2020.J2000      月亮表 

         soltab.2020.J2000      太阳表 

         nutabl.2020         章动表 

         vmf1grd.2020         映射函数模型             

         leap.sec             跳秒

         atmdisp_cm.2020

2、文件准备工作

2.1、创建工作文件夹

       在“opt”文件夹内新建“work”文件夹。

2.2、在work内新建数据文件夹

       新建“brdc、igs、rinex、met”4个文件夹，分别存放当天的广播星历、精密星历和观测值文件、气象文件。分别把对应的文件放在文件夹中（可不用这个met文件，本次演示没用。广播星历和精密星历的下载方法下次出一期）。

     （注意：如果文件采用.z压缩格式，则应当使用进行解压；如果解压出的观测值文件仍采用.d的压缩格式，则应使用crx2rnx软件将其转化为.o，如果不会就百度d转o吧）。

2.3、在work文件夹下右键，点击“在终端中打开”，输入“sh_setup -yr 2020”命令，会发现主文件夹下出现了“tables”文件夹。这个终端别退，后面要用呢。

2.4、打开上一步的“tables”文件夹，右键在终端中打开， 输入“ls -sl”命令查看链接情况（这里有的用不到，红色也无所谓的。比如第一个红色，我在第6步中sestbl.设置了atl.list为no就是用不到的）。

2.5、红色连接处理

       假如atmdisp_cm.2020没连接上，是红色的（这里是假如，如果你是红绿色盲，我建议立即放弃）。打开第2.3步的终端（上面说了要用的啊），依次输入下面三个命令。

% 移除红色错误链接

rm -f tables/atml.grid

% 链接正确文件，这个文件是第1步更新的表文件

ln -s ~/gg/tables/atmdisp_cm.2020* tables/atml.grid

% 这个是用来检查是否连接成功，在下图为绿色说明成功

ls -l tables/atml.grid

3、更改配置文件

3.1、更改sestbl.文件。

打开tables下的“sestbl.”，该文件内部有详细的说明，翻译下就知道大概意思。我的修改如下，所在行数已经标出来，如果好兄弟还是找不到，我劝你放弃读研（我主要是为了获取ZTD）。

      36行  Elevation Cutoff = 10 

      33行  Interval zen = 1 

      44行：Met obs source = GPT 50

      46行：Output met = Y 

      49行：DMap = VMF1

      50行：WMap = VMF1

      52行：Use map.grid = Y

      64行：Use otl.grid = Y 

      71行：Use atml.grid = Y

3.2、更改sites.defaults文件。

打开“tables”文件夹下的“sites.defaults”文件，它是用来控制需要参与解算的测站。把红圈处改为工作文件夹（我的是“work”），蓝圈处改为测站名字。

3.3、生成测站信息文件station.info文件。

      3.3.1、将work/tables下的“station.info”文件拷贝到“rinex”文件夹下，打开并编辑，仅保留以#或*开头的前几行，如下图红框的内容。

      3.3.2、删除完成，如下图一样。

      3.3.3、右键进入终端，打开终端并进入“rinex”文件夹，运行下面命令（20是年份的意思）：

sh_upd_stnfo -files *.20o

      3.3.4、命令完成后如下图，将生成的“station.info”文件拷贝到text/tables文件夹，覆盖原文件。

       注：输入此命令有时会报错：打不开guess-rcvant.dat和rcvant.dat文件，这时只要从tables把这两个文件复制到当前目录下即可。

3.4、测站近似坐标文件lfile.文件。

       igs站可不用生成lfile.文件。进入“work/rinex“文件夹内打开终端。解算的每个站点都要依次生成，我结算三个站点，就需要输入三个命令。

       sh_rx2apr -site dxin0050.20o -nav ../brdc/brdc0050.20n

       sh_rx2apr -site fjwy0050.20o -nav ../brdc/brdc0050.20n

       sh_rx2apr -site gsdh0050.20o -nav ../brdc/brdc0050.20n

注：以dxin站为例sh_rx2apr -site dxin0050.20o -nav ../brdc/brdc0050.20n中红色为o文件名字，蓝色为brdc名字。（igs站可以用这个生成lfile.的文件grep bjfs itrf14.apr > lfile.1 ）。

       将上面的站点生成“*.apr（或者*.lfile）”文件进行合并数据进行合并，输入下面命令：

       cat *.apr > lfile.

       cat lfile.* > lfile.

会组合生成“lfile.”文件，将其拷贝到text/tables文件夹，覆盖原文件。

4、解算文件

4.1、解算命令

      在前面章节2.3的终端下输入“sh_gamit -expt work -d 2020 005 -orbit IGSF -met -metutil Z -noftp”等待解算完成。（work：文件夹。-d：处理一天数据，多天改为-s。2020 055：处理的年份和年纪日。-noftp：不启用ftp。还有几个忘记了，看到了再补充）。

4.2、结果文件

       将存放在名称为年积日的文件夹内，此例中为“005”，可供参考的结果文件为：

       sh_gamit_2020_005.summary 解算总结

       qworka.005 解算记录

       oworka.005 解算记录的简略版

4.3、分析结果文件

      打开“qworka.005”，划到最底部。nrms是衡量结果的重要指标，通常postfit nrms值在3~10mm,prefit nrms值应在0.2左右，如果存在很大偏差，就需要查找原因了，重新解算了。

5、个人总结

        我个人做gamit解算是为了提取天顶总延迟 ZTD。我是大气水汽方向的，所以需要求ZTD还有可将水量（里面的每个站点m文件就是可以提取需要的ZTD和PWV）。教程或许会有一些瑕疵，还请海涵。如果有不对的地方，欢迎大家指正。本文出于帮助初学者利用gamit软件解算GPS基线的目的撰写，但受笔者水平限制，难免存在疏漏和错误，欢迎大家批评指正。