Gamit解算脚本

这是一个解算单天的shell脚本，对于初学者很有帮助。

首先就是需要在项目（四个字符）建立rinex brdc igs 还有以年纪日命名的目录，然后提前准备好station.info和lfile.文件，以及将tables表更新，然后将脚本放在以年纪日命名的目录中，执行即可。

 

#！/bin/bash

#AUTHOR 王会
#EMAIL    wzxwhd@126.com
#READEME
#    在时段文件下执行该脚本
#    需要准别下列文件：
    #    station.info  sh_upd_stnfo
    #    lfile.

#step 1:自定义工程
expt='wwww'    #工程名称
yr_4=2016    #工程时间
yr_2=16
yr_1=6
doy=123        #年纪日
sp3type='igsf'    #igs星历类型
sp3file='igs18951.sp3'
brdcfile='brdc1230.16n'
session=99    #处理时段，99代表所有

#step2:链接文件
links.day ${yr_4} ${doy} ${expt}

#step3:移除坏链接，然后重新链接
rm luntab.    #极移，太阳表，月亮表
rm soltab.
rm nutabl.
rm pmu.
rm ut1.
rm map.grid
rm lfile.
ln -s ../tables/lfile. ./
ln -s ../tables/luntab.2016.J2000 ./luntab.
ln -s ../tables/soltab.2016.J2000 ./soltab.
ln -s ../tables/nutabl.2016 ./nutabl.
ln -s ../tables/ut1.usno ./ut1.
ln -s ../tables/pmu.usno ./pmu.
ln -s ../tables/vmf1grd.2016 ./map.grid

#step4:链接o和sp3和n文件
rm *.${yr_2}o
rm *.${yr_2}n
rm *.sp3

ln -s ../rinex/*.${yr_2}o ./
ln -s ../brdc/$brdcfile ./
ln -s ../igs/$sp3file ./

#step5:分布处理
rm *.fatal
rm session.info

sh_makexp -expt $expt -orbt $sp3type -yr $yr_4 -doy $doy -sess 99 -srin -nav $brdcfile -apr lfile. -sinfo 30 00 00 2880
sh_sp3fit -f $sp3file -o $sp3type -d $yr_4 $doy -t -u
sh_check_sess -sess $doy -type gfile -file g${sp3type}${yr_1}.$doy
makej $brdcfile jbrdc${yr_1}.$doy
sh_check_sess -sess $doy -type jfile -file jbrdc${yr_1}.$doy
makex ${expt}.makex.batch
fixdrv d${expt}${yr_1}.$doy
csh b${expt}${yr_1}.bat