GAMIT基线解算_分部处理

分步处理2017年001天（chan shao bifs）

数据准备工作：

1）mkdir demo

2）cd demo

3）mkdir rinex brdc igs

4）cd brdc

两种下载方式：

第一种方式：采用gamit自带的下载脚本

sh_get_nav -archive sopac -yr 2017 -doy 001 -ndays 1

第二种方式：采用wget下载

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

解压:uncompress *

注：采用自带脚本下载无需解压，后面下载精密星历和观测文件类似

5）cd igs

第一种方式：采用gamit自带的下载脚本

sh_get_orbits -archive sopac -yr 2017 -doy 001 -ndays 1

第二种方式：采用wget下载

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

解压:uncompress *

6）cd rinex

第一种方式：采用gamit自带的下载脚本

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

第二种方式：采用wget下载

wget --http-user=anonymous --http-passwd=***@126.com -nd http://garner.ucsd.edu/pub/rinex/2017/001/chan0010.17d.Z

wget --http-user=anonymous --http-passwd=***@126.com -nd http://garner.ucsd.edu/pub/rinex/2017/001/shao0010.17d.Z

wget --http-user=anonymous --http-passwd=***@126.com -nd http://garner.ucsd.edu/pub/rinex/2017/001/bjfs0010.17d.Z

解压：sh_crx2rnx -f <filename>  #此方法只能一个个站点解压，有点扎心）

7）链接表文件

sh_setup -yr 2017 001

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

Use otl.list = Y         ; Ocean tidal loading list file from OSO

Use otl.grid = N         ; Ocean tidal loading grid file, GAMIT-format converted from OSO）

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 # 提取观测文件中的概略坐标

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

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

10）在001目录下链接文件

ln -s ../rinex/*.17o ./

ln -s ../brdc/brdc0010.17n ./

ln -s ../igs/igs19300.sp3 ./

cp ../tables/* ./

下面开始正式解算

1）sh_makexp -expt demo -orbt igsf -yr 2017 -doy 001 -sess 99 -srin -nav brdc0010.17n -apr lfile. -sinfo 30 00 00 2880

2)  sh_sp3fit -f igs19300.sp3 -o igsf -d 2017 001 -r BERNE -t -u

3)  sh_check_sess -sess 001 -type gfile -file gigsf7.001

4)  makej brdc0010.17n jbrdc7.001

5)  sh_check_sess -sess 001 -type jfile -file jbrdc7.001

6)  makex demo.makex.batch

7)  fixdrv ddemo7.001

8)  csh bdemo7.bat

运行以上步骤即可完成基线分步处理解算任务

若想生成z文件，可在解算前sestbl. 文件中设置
Output met = Y

若后面进行BLOBK平差，则需把sestbl.更改为

Choice of Experiment = RELAX.   ; BASELINE/RELAX./ORBIT #此方式用于后面的时间序列分析

Choice of Experiment = BASELINE  ; BASELINE/RELAX./ORBIT#采用默认可保证点位精度更高