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

在计算速度时候，我们可以直接利用我们的gps测站求解速度，但是速度具有参考性，我们得到的速度并不是在文献中常说的相对于ITRF框架或者欧亚框架的速度，所以我们此次介绍得到ITRF框架下的速度

在检查完自己的时间序列后，删除粗差后，我们来求解速度

1、首先我们得到全球解，也就是全球的h文件，

      下载全球数据利用sh_get_hfiles得到mit的H文件或者sopac的h文件

·        sh_get_hfiles-yr 2013 -doy 246 -ndays 5 -net igsall eura

·       或者直接利用sh_glred进行下载

·        sh_glred-s 2013 024 2013 025  -expt sina -opt R LA H G T >sh_gred.log 合并全球解与自己的解单天解

这样其实就是相当于你用gamit结算了全球的解，然后一起得到一个参考全球的速度，为了将自己的解与全球的解融合，这里在gamit中结算需要有几个igs站点的结算。

2、将tables的globkvel的cmd与glorg的cmd文件拷贝过来，编辑两个vel的cmd用新的apr和新的eq文件，在org_opt中加入VSUM 使用松约束加入转换旋转对地球eop参数。在glorg中加入rate_org的估计参考站为处理的几个站，globk结算速度时，gdl后面不用＋1，默认为结算速度，

3、运行globk的命令

   globk 6 sina.prt sina.log sina.gdl globk_velNaNd

4、 得到chi～比例因子 加入到gdl中，删除前面生成的prt文件与org文件重新运行

·         sed‘s/$/0.034 /g' sina.gdl

·         globk6 sina.prt sina.log sina.gdl globk_velNaNd

5、画出速度

 sh_plotvel -fcoco1213_rep.vel  -maxsigma 4 -R225/310/0/50 -page L -ps coco1213_rep

注意：

查看globk_vel.org中确信所有的chi2增量在范围0.2-1.5中，如果一些偏高，则应该查找时间序列的outlier，然后分别对每一个站点进行sig_neu,或者对于每一个h文件因子在gdl中进行调整因子，最后查看总的chi增量在chi2 increments自己的解的因子用1左右，并且参考框架的站点的rms小于2mm/yr

当然：我们也可以利用tsfit进行结算速度，这也是有效的方法之一

1、       仅仅准备globkcomb的cmd与glorg的cmd就可以，然后编辑文件，加入PBOB的命令在org_prt中，然后首先运行

glred6 sina.prt sina.log sina.gdl globk_velNaNd

2、得到pos文件

tssum . ex –R *.org

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

3、编辑tsfitNaNd文件，运行tsfit

        tsfit tsfitNaNd tsfit.sum *pos

        命令编辑如下：

eq_file ../tables/ITRF14.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、直接可以得到velfile，然后直接画图即可。实现的算法就是使用的real_sigma，也就是利用tsview以及velview的工具实现的算法，这里也就是tsfit是它们的最新替代品。

下一节会讲到 如何实现参考框架，以及转化到欧亚板块的速度。

以及如何定义自己的板块，得到自己板块区域的速度