【GAMIT/GLOBK学习】glred重复性分析

glred重复性分析

（1）在GAMIT基线解算之后，在工程目录（demo）下，运行

sh_glred -cmd

（它会将gg/tables下的两个cmd命令拷贝到gsoln目录下）

进入到gsoln目录，

cd gsoln

ls查看
分别打开vim globk.cmd和vim glorg.cmd两个命令文件，并对其中的一些内容进行编辑修改。（详情参考上一篇blog）

（2）回到工程目录，cd ..
新建一个hfile的文件夹，

mkdir hfile

需要copy到hfile文件夹的内容有两部分
（a）每日解（例：245）文件中的h文件
（b）下载SOPAC全球h文件（igs1、igs2、…、igs7）

第一部分：查找h文件:

find ./ -name "hdemo*"

复制h文件：

find ./ -name "hdemo*" | xargs -i cp {} ./hfile/

第二部分：全球h文件按用脚本批量下载（hfiles.sh）
（之后别忘记将下载好的文件拷贝到hfile文件中）

copy完之后，查看一下hfile文件，ls hfile

（3）将svnav.dat拷贝到hfile目录下

cp ./tables/svnav.dat ./hfile

运行htoglb将所有的ASCII的h文件转换成二进制h文件

htoglb ./glbf/ ./tables/svs_shaz.svs ./tables/hfile/h*

（其中，glbf是指输出到该目录下，svs_shaz.svs是生成的）

到这一步，完成了文件转化
先cd glbf，进入该目录，然后查看ls，可以看到有glx、glr两类文件。

cd ../gsoln，进入该目录，ls

（4）制作h文件列表

ls ../glbf/ *.glx > shaz.gdl

对.gdl文件进行编辑，添加权重
之后，查看文件列表

cat shaz_glx.gdl

（5）在gsoln目录下，运行glred，生成时间序列，运行，

glred 6 globk_rep.prt globk_rep.log shaz_glx.gdl globk.cmd

（6）从glred解文件.org中绘制坐标的重复性，运行，

sh_plot_pos -f globk_rep.org

之后，要想查看时间序列图，可以新建一个plot目录，将生成的图移动到plot目录下，在gsoln目录下

mkdir plot
cp *.ps ./plot

这样在一个plot目录下，方便查看。

（7）可用tssum命令产生时间序列文件，例如：

tssum ./xyz pbo.final_ifr14 -R shanzhu_rep.org

之后也可以对时间序列检查，对异常值（不能得知原因）进行剔除。

在之后的拟合等操作，可以参考GLOBK的用户手册，不再赘述。