处理过程中，为了自己的方便使用，需要编写自己的脚本（用红色标出）

第一种方法是，得到站点位置后：

1、 grep Unc. sina.org >sina.apr 得到运算点的坐标

作为apr坐标，然后输入到globkcomb的cmd的apr选项中，然后重新进行画图

2、 然后运行在gsoln目录下：

glred6 glred_20150811.prt glred_20150811.log 201407_NSFBay.gdl cmd

sh_plotcrd -f *.org-expt emed00

上述仍然需要编辑时间序列，遇到不好的天数进行重新编辑，所以我们采用：

第二种方法，当然我们介绍一种新的globk处理时间序列的算法tsfit

1、在globk的cmd中org_opt 加入PBOB 命令变为，运行

glred 6 glred_20150811.prt glred_20150811.log 201407_NSFBay.gdl

2、得到pos文件

tssum . ex –R *.org

. 是当前目录，ex是随意命名，然后得到pos文件

3、编辑tsfitNaNd文件，运行tsfit

    tsfit tsfitNaNd tsfit.sum *pos

    命令编辑如下：

eq_file ../tables/RESET.eq

*Annual signal

per 365.25

*Remove bad outliers (5-sigma)

nsigma 5

* Savethe velocity file (and an apriori coordinate file)

velfile @.vel

* Savethe edited data list

rep_edits @.rns

* Usethe realistic sigma algorithm for the velocity sigmas

real_sigma

4、将得到的eq文件  上面输出的eq文件应为tsfit.rns 输入到globkcomd的cmd文件中，重新运行glred，删除原先的org文件以及prt文件

glred6 glred.prt glred.log 201407_NSFBay.gdl globkcmd

5、得到新的org文件，进行画图

sh_plotcrd -f *.org-expt emed00

可以得到不同的图以及误差重命名的GPS站点图

也可以用

sh_plot_pos -f globk_rep.org 生成raw时间序列

备注：当然在精准的时间序列时，可以选择利用globk得到的站点坐标返回gamit进行重新再次处理，然后继续返回这样的方法进行求解时间序列可以得到更佳精确的时间序列。