分步处理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在下面网站下载
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 *.17ocp 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