解完基线后的平差，我一般使用globk分布处理。

写了个批处理的globk脚本，根据输入的参数glred/globk，可以分别得到时间序列和速度并画图。

使用前需新建glred和globk文件夹，并将提前准备好的两个cmd和apr文件放到glred/globk文件夹内。

将脚本放到工程目录执行即可。

#!/bin/bash
#############
##提前准备好两个cmd,apr文件放到globk/glred文件里，之后执行脚本
if [ $# -eq 0 ]
then
echo '##########################################'
echo '#                                        #'
echo '# Usage： bash sh_globk globk/glred      #'
echo '#                                        #'
echo '##########################################'
exit
fi
###########
run=$1 
if [ ${run} == 'glred' -o ${run} == 'globk' ];then
  if [ -e ./${run}/globk_简.cmd  ];then
  mkdir hfiles 
  cp ./tables/svnav.dat hfiles
  find ./ -name "hexpt*" | xargs -i cp {} ./hfiles/
  htoglb ./glbf/ ./tables/svs_myexp.svs ./hfiles/h*
  cd ./${run} 
  ls ../glbf/h*.glx >> expt.gdl
  ${run} 6 ${run}.prt ${run}.log expt.gdl globk_简.cmd
####  glred for TimeSeries #################
    if [ ${run} == 'glred' ];then
        echo -e "\033[42;37m "++++++++   TimeSeries  +++++++"\033[0m"
	echo -e "\033[42;37m "++++++++   TimeSeries  +++++++"\033[0m"
        mkdir plot 
	cp ./${run}.org plot 
	cd plot 
	mkdir pos_xyz
	tssum ./pos_xyz cors -R glred.org
	xyzsave 0 ./pos_xyz glred.org 
	sh_plot_pos -f ${run}.org
	gmt psconvert -Tt -A0.5c *.ps
	rm *.ps
#####   globk for Vel   ###################	
    elif  [ ${run} == 'globk' ];then
	echo -e "\033[42;37m "++++++++  globk  +++++++"\033[0m"
        echo -e "\033[42;37m "++++++++  globk  +++++++"\033[0m"
	grep 'Unc.' globk.org >> result.xyz
	mkdir vel
	cp globk.org vel
	cd vel
	sh_org2vel -file globk.org
	sh_plotvel -file globk.vel
    fi
######################################
  else 
       echo -e "\033[42;37m "+++++ do not exist ${run} or CMD files     +++++"\033[0m"
   fi
else 
        echo .
	echo -e "\033[42;37m "+++++ please input glred/globk   +++++"\033[0m"
	echo .
fi