Gentle''的个人博客分享 http://blog.sciencenet.cn/u/tao164411096 九分耕耘,一份收获。

博文

GLOBK批处理脚本

已有 618 次阅读 2019-4-20 13:47 |个人分类:GAMIT/GLOBK|系统分类:科研笔记| globk批处理脚本

解完基线后的平差,我一般使用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




http://blog.sciencenet.cn/blog-3391834-1174491.html

上一篇:GAMIT批处理脚本
下一篇:bat批量质检脚本

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...
扫一扫,分享此博文

Archiver|手机版|科学网 ( 京ICP备14006957 )

GMT+8, 2019-6-17 08:44

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部