GAMIT基线解算_批处理

分步处理2017年001,002天（shao hkkt hkks，后面两个站点为香港CORS站）

数据准备工作：

1）mkdir demo

2）cd demo

3）mkdir rinex brdc igs

4）cd brdc

采用wget下载

wget -c ftp://nfs.kasi.re.kr/gps/data/daily/2017/001/17n/brdc0010.17n.Z

wget -c ftp://nfs.kasi.re.kr/gps/data/daily/2017/002/17n/brdc0010.17n.Z

解压:uncompress *

5）cd igs

采用wget下载

wget -nd ftp://cddis.gsfc.nasa.gov/pub/gps/products/1930/igs19300.sp3.Z

wget -nd ftp://cddis.gsfc.nasa.gov/pub/gps/products/1930/igs19301.sp3.Z

解压:uncompress *

6）cd rinex

shao站采用gamit自带的下载脚本

sh_get_rinex -archive sopac cddis  -yr 2017 -doy 001 -ndays 1 -sites shao

sh_get_rinex -archive sopac cddis  -yr 2017 -doy 002 -ndays 1 -sites shao

hkkt hkks在下面网站下载

https://www.geodetic.gov.hk/tc/rinex/downv.aspx

7）链接表文件

sh_setup -yr 2017

（在 sestbl. 文件中将海潮模型设置为：

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）

删除otl.grid链接：rm -f tables/otl.grid

重新将海潮模型格网文件链接到otl.grid下：ln -s ~/gg/tables/otl_NAO99b.grid* otl.grid

8）更新station.info

在rinex目录下：

sed -n '2,6p' ../tables/station.info > station.info
sh_upd_stnfo -files *.17o

cp station.info ../tables

9）更新L文件

在table表目录下执行

grep "APPROX POSITION XYZ"$ ../rinex/*o > lfile.rnx # 提取观测文件中的概略坐标

#上面这一步莫似对hkkt和hkks站用不上，暂时不知道啥原因，目前先采用手动添加的方式，打开o文件，找到“APPROX POSITION XYZ”这一行，手动将hkkt和hkks站点近似坐标添加在由上一条命令生成的lfile.rnx内，完成此步骤后继续看下面

rx2apr lfile.rnx 2017 001 # 生成 lfile.rnx.apr 文件

gapr_to_l lfile.rnx.apr lfile. 2017 001 # 生成 lfile. 文件

(可以仅采用001天的观测文件的近似坐标生成L文件)

完成以上数据准备工作后就可以采用批处理命令了，一条解决

sh_gamit -expt text -d 2017 001 -orbit IGSF -pres ELEV -noftp -dopt D ao x c > sh_gamit.log

*****

此外，在fig目录内可查看天空图与相位残差图

如：xdg-open HKKS.002_elev_res.png

        xdg-open HKKS.002.png