gamit10.71解算教程（入门级别，仅供参考）

1. 数据准备
   1.1 数据下载
   本次实例拟使用aira bjnm chan cmum hkws mizu twtf ulab等IGS站数据进行解算，计算日本的两个站mizu, aira的位置，解算时间为2019年129天。
   

在主目录下创建工程文件目录
mkdir test
链接tables文件夹下中的文件
sh_setup -yr 2019 -apr igs14_comb.apr
（注：igs14_comb.apr文件中保存IGS站的先验信息，采用itrf14_comb.apr也可以，两个文件有细小的差别主要是地球参考框架的实现方式不同）
进入test文件夹下创建精密星历igs 广播星历brdc 观测值rinex文件夹 电离层格网文件夹ionex
（注：电离层格网文件是gamit10.71计算新加入的必须的文件）
cd test
mkdir igs brdc rinex ionex
首先下载rinex观测值数据，实例中均为IGS站，观测值可以从武汉大学IGS数据中心下载,将下载好的文件放入rinex文件夹下(注：计算前要确保下载的文件中接收机、天线类型、天线高、文件名等是否正确）。然后准备广播星历数据，同样可以在武汉大学IGS数据中心下载，下载后放入brdc文件下（注：此次示例仅采用GPS计算，使用brdc广播星历即可，若要计算多系统的需要在cddis下载brdm多系统广播星历）。下载精密星历放在igs文件下（网解仅采用精密轨道文件即可，无需精密钟差）。电离层格网文件需要从cddis中心下载，文件名为igsg1290.19i


1.2 表文件更新（注：表文件对于解算精度很重要，这里仅展示如何计算不考虑精度，因此省略）

2. 修改表文件
2.1 sites.defaults
在tables文件夹下找到并打开该文件
cd tables
vi sites.defaults
修改工程名（test) 与测站名(需要解算的测站）

2.2 station.info
vi station.info
首先删除文件中的多余信息
:7,$d （注：7代表开始行，美元符号代表结尾行)
键入回车即可删除

保存后退出
输入命令编辑station.info文件
sh_upd_stnfo -files . ./rinex/*.*o
生成的station.info文件如下：

（注：该文件原本存在IGS站信息，重新生成主要是为了非IGS站）

2.3 lfile.文件
该文件中主要保存各个站的初始坐标，示例中均为IGS基准站，不作更新即可，直接采用igs14_comb.apr文件

2.4 sittbl.文件
起算点坐标分量约束为5cm, 非起算点为30m(本示例中均为IGS站，不改也可以）

2.5 sestbl. 文件
建议仔细学习一下其中每行的意思，这里不做阐述，采用默认值
如果想采用各种模型，如海潮模型，大气模型等可以在MIT的ftp服务器上下载后放到tables中
ftp://everest.mit.edu/pub/GRIDS/

2.6 process.defaults
这个文件主要是用于控制计算的处理进程，大部分配置不需要更改，可以自己研究一下。可以在这个文件里面设置计算时常，开始时间等。
3.批处理
sh_gamit -d 2019 129 -orbit IGSF -expt test -noftp

运行报错：

这个问题是因为配置文件默认改正海潮模型，需要otl.list和otl.grid文件，这两个文件需要手动下载，感兴趣的话可以自己下载一下这两个文件放入tables文件下即可。作为演示，我们这里更改配置文件不考虑海潮模型
进入tables文件下
vi sestbl.
修改如下配置：

重新运行sh_gamit命令即可。

基线计算结束

4. 计算结果

解算完成后需要检查计算结果。计算完成会生成以年纪日为目录名称的文件夹，这里是129，在该文件夹下有一个sh_gamit_ddd.summary文件统计计算结果。其中postfit RMS值在3~10mm,prefit nrms值应在0.2左右。还有宽巷和窄巷模糊度固定率。本示例中的站有在日本的，日本的测站由于板块活动明显数据质量较差，同时为考虑海潮模型，对流层模型改正，部分测站如CMUM站较差，模糊度固定率也偏低。

具体基线解算结果在q文件和o文件中，o文件是简略版。

5 GLOBK平差
在test目录下执行
htoglb glbf /dev/null 129/h* a.*
进入gslon目录下，在安装目录~/gg/tables/文件夹下拷贝globk.cmd和glorg.cmd文件到该目录下，打开globk.cmd文件设置平差的起算点，修改apr_site如下。网平差需要进行两次，第一行apr_site all 10 10 10 0 0 0用于第一次网平差，需要给松约束10 10 10，第二次网平差需要给定义的基准点紧约束。

这里采用IGS站的坐标作为起算点，计算日本的两个站的坐标，可以给较强的约束。起算点个数看实际情况给出。
修改glorg.cmd文件中的stab_site行添加固定点。固定的个数同样看实际情况几个都可以。

计算多天的数据可以在该文件中固定点约束前面行加上stab_site clear命令
在gsoln目录下执行 ls . ./glbf/*.glx > test.gdl
最后执行 globk 6 globk.org globk.log test.gdl globk.cmd
在globk.org文件中即可看到最终的结果,apr文件是最终的平差结果，位置估计部分可以看到基准点使用个数和调整结果。也可以在在终端使用命令 grep ‘POS S’ globk.org查看使用的基准点等。