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

博文

GAMIT批处理脚本

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

最近在解算多年的CORS数据时发现使用自带的sh_gamit批处理脚本时,如果某天报错就会停止,一点都不好用。

利用sh_gamit自己写了个循环小脚本。

脚本放到工程目录下,输入年,开始年积日,结束年积日(年积日为个位数或十位是前面不补0)



#!/bin/bash
###################
if [ $# -lt 3   ]
then
echo '#######################################'
echo '#                                     #'  
echo '# Usage:  bash sh_gamit 2017 1 365    #'
echo '#                                     #'
echo '#######################################'
exit
fi
#########
#
year=$1
doy_str=$2
doy_end=$3
mkdir sh_summary
for ((i=${doy_str};i<=${doy_end};i++));
do
{
a=$(printf "%03d" "${doy_str}")
echo -e "\033[42;37m begin process $a \033[0m"
   sh_gamit -d ${year} $a -orbit igsf > ./sh_summary/sh_gamit_$a.log
    if [ -e ./$a/qexpta.$a ];then
echo -e "\033[42;37m +++++      $a sussessful        ++++\033[0m"
   grep 'nrms' ./$a/qexpta.$a | head -n 1 
     elif [ -e ./$a/GAMIT.fatal ];then
echo -e "\033[41;37m +++++      $a ERROR             ++++\033[0m"
echo  $a error >> AAprocess_error.log          
cat ./$a/GAMIT.fatal >> AAprocess_error.log
     else  
echo -e "\033[41;37m     $a ERROR but do not exist GAMIT.fatal     \033[0m"
    fi
let doy_str++
}
done




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

上一篇:win10系统自带的ubuntu安装路径
下一篇:GLOBK批处理脚本

0

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

数据加载中...

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

GMT+8, 2019-10-18 07:05

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部