最近在解算多年的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