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

globk主要作用一个生成速度场估计，第二是生成一个比较好的时间序列在不用的参考框架下，第三就是合并大型的gps网络站，当然可以得到坐标站点解

定义参考框架针对全球的地质问题与itrf08相结合，需要igs参考站 100个左右，对于区域的问题，用一个板块相符合的站或者不用站点

对于globk分析时，gamit的解不能改变，local站时间序列的形成与gamit的结算有关。

globk是一个平滑卡曼滤波器可以包含随机游走过程的噪音在估计中，合并单天解，合并许多年的数据生成位置，速度，偏移，震后形变的参数等 glred等于globk？ globk激发glorg

Glorg 应用约束，实现稳定框架，将站点吻合到所选的框架下，计算卡曼滤波的协方差矩阵

对于小的网络只选择pos_organd rate_org: xtran ytran ztran, NOT xrot, yrot, zrot

在使用globk的重要的命令为：sh_glred

利用sh_glred选项

－H 激发htoglb 将h文件转化h二进制文件  htoglb . /dev/null ../[0-3][0-9][0-9]/h*a.*

－G激发 glred生成各自的org解

－T 激发sh_plot_pos生成pos文件 进行时间序列画图

-ncomb合并解的天数

-expt 工程目录名字

-local仅仅合并自己的h文件

-archive接收h文件的机构名字

-net h文件的网络名字

-prt为画ptr而不是org文件

-detrend为删除趋势临界值为32天

-cmd生成两个cmd文件

-globk_cmdfile 自己的编辑的cmd文件

-glorg_cmdfile 自己编辑的glorgcmd文件

-mb 为不删除mb文件

－psbase 用命令 sh_baseline

-s为开始到截止日期

－d单天解

-opt 可选项为H为生成二进制h文件

R为删除globk中h文件中的gl文件

F为寻找mit和sopac的网络h文件

E为画图

在globk的cmd中 org_cmd 选项

opt选项输出

·       ERAS-- erase file before writing (normally files appended)‏

·       NOPR-- Do not write output ( e.g., for globk when invoking glorg )‏

·      BLEN -- Baselinelengths

·      BRAT -- baselinerates when velocities estimated

·      RNRP -- generatesreports on differences in parameter estimates after renames.

·       FIXA-- makes apriori coordinates and velocities consistent when equates are used inglorg (can sometimes fail in complicated rename scenarios--best if apr_file isprovided with consistent values)‏

·       VSUM-- Lat/long summary of velocity (needed to plot velocities)‏

·      PSUM -- Lat/longposition summary

·      GDLF --Include listof hfiles and chi**2 increments from run

·      CMDS -- Echos globkcommand file into output file

·      以上命令需要细细参考globk的cookbook

下面介绍一下，如何生成正确的坐标点位解：

1、在完成gamit解后，得到单天解的h文件，在gsoln中将tables中的globkcomb的cmd以及glorg的cmd文件拷贝过来，然后直接运行

sh_glred-s 1999 235 1999 250 -expt emed -opt H G E

2、查看时间序列好坏。然后查看时间序列，删除粗差大的单天解

然后运行在gsoln目录下：htoglb . /dev/null ../[0-3][0-9][0-9]/h*a.*

      ls ../glbf/h*.glx >expt.glx.gdl

glred6 glred_20150811.prt glred_20150811.log 201407_NSFBay.gdl

sh_plotcrd -f *.org-expt emed00

如果满意时间序列后，没有明显粗差后，满意所在的时间序列后

3、运行globk，删除其他的prt以及org文件，然后运行

globk6 globk_nmtbx_00100.prt globk_nmtbx_00100.log L000409_nmtbx.gdl globk_fixNaNd

grep'2 for' sina.org 得到误差比例因子

4、sed ‘s/$/0.034 /g' sina.gdl 得到新的gdl后，将误差因子编辑到gdl文件后，再次运行globk命令

globk6 sina.prt sina.log sina.gdl globk_combcmd

5、得到结算点的坐标

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

以上文件类似（L000409_nmtbx.gdl）都是随意，自己结算时需要使用自己的文件，提醒其它gamit／glbok的错误使用方式，本文为准。