1.首先更新tables
2.建一个项目,mkdir test
进入项目,cd test
3.数据准备。
建立3个文件夹,分别放导航文件,观测文件,精密星历。mkdir igs brdc rinex
分别下载这三种数据。
(1)cd brdc
sh_get_nav -archive sopac -yr 2017 -doy 010 -ndays 1
(2)cd ../igs
sh_get_orbits -archive sopac -yr 2017 -doy 010 -ndays 1
(3)cd ../rinex
sh_get_rinex -archive sopac cddis -yr 2017 -doy 010 -ndays 1 -sites urum lhaz shao chao
(注:有很多其他下载数据的方法)
4.生成tables文件并链接tables表。
cd ..(回到项目目录下)
(1)sh_setup -yr 2017 (生成一个tables文件)
(2)cd tables
vim sestbl. 改两个参数的Y,N Use otl.list = Y Use otl.grid = N
5.更新station.info
在rinex 文件下
sed -n "1,6p" ../tables/station.info > ./station.info
sh_upd_stnfo -files *.17o
cp station.info ../tables
6.根据需要修改相关参数
进入tables目录
vim process.defaults,把set use_rxc的N改为Y
其他参数根据需要自行修改。
7.制作l文件
进入rinex文件
grep POSITION *.17o > lfile.rnx
rx2apr lfile.rnx 2017 010
gapr_to_l lfile.rnx.apr lfile. “” 2017 010
8.接下来进行分布解算。
(1)首先在项目下新建一个文件夹 mkdir 010
(2)进入所建的文件010,将所需要的文件链接进去
ln -s ../brdc/brdc0100.17n ./
ln -s ../igs/igs19312.sp3 ./
ln -s ../rinex/*.17o ./
cp ../tables/* ./
(3)分布解算
1)sh_makexp -expt test -orbt igsf -sp3file igs19312.sp3 -yr 2017 -doy 010 -sess 99 -srin -nav brdc0100.17n -apr lfile. -sinfo 30 00 00 2880
2)由于版本问题,首先从tables表更新的网站下载svnav.dat.gnss文件并重命名svnav.dat,将其代替原来文件svnav.dat。如果不改的话,可能出现不能识别svnav.dat版本问题。
sh_sp3fit -f igs19312.sp3 -o igsf -d 2017 010 -r BERNE -t -u
3)sh_check_sess -sess 010 -type gfile -file gigsf7.010
4)makej brdc0100.17n jbrdc7.010
5)sh_check_sess -sess 010 -type jfile -file jbrdc7.010
6)makex test.makex.batch
如果在此过程出现了cannot fund site code SHAO on L-file,可能是在将tables表内容复制到010文件时出了错。
将rinex的lfile.文件复制代替tables和010文件中的lfile.文件,再重新运行该步骤。
7)fixdrv dtest7.010
8)csh btest7.bat
出现ARC not exected时,删除解算目录下的GAMIT.fatal 和GAMIT.warning,再次重新运行此命令。