diff --git a/gmtsph-regional.sh b/gmtsph-regional.sh index 83f3137..c542455 100755 --- a/gmtsph-regional.sh +++ b/gmtsph-regional.sh @@ -18,6 +18,7 @@ inverseCPT=0 coastline=0 axistick=("a" "a") #横纵坐标轴间隔 a表示自动 bartick="a" #色标轴标轴间隔 a表示自动 +polyfile='null' # 不同图幅的绘图参数 pageChoice='small' # 默认尺寸为small exSmallpage="1i,0.08i/-0.25i,0.85i/0.03i" @@ -25,7 +26,7 @@ smallPage="1.5i,0.125i/-0.3i,1.2i/0.05i" # 图大小,色标位置,色标尺寸 middlePage="2i,0.2i/-0.35i,1.5i/0.05i" largePage="2.5i,0.35i/-0.4i,1.75i/0.08i" # 从命令行获取参数 -while getopts "hi:r:u:c:a:t:v:l:G:gnpb" arg +while getopts "hi:r:u:c:a:t:v:l:G:s:gnpb" arg do case $arg in h) @@ -45,6 +46,7 @@ For further explanations, please look for GMT's manuscripts.\nAuthor: Yi Zhang ( printf "%s\t%s\n" "-G" "provide a different input grid data for the use of grdgradient. This should be used with the '-g' option at the same time" printf "%s\t%s\n" "-n" "disable the use of grd2cpt." printf "%s\t%s\n" "-p" "inverse the color pattern specified by the -c option, no use if -n option is implemented." + printf "%s\t%s\n" "-s" "plot polygons via a file." exit 0;; i) data=$OPTARG;; @@ -73,6 +75,8 @@ For further explanations, please look for GMT's manuscripts.\nAuthor: Yi Zhang ( nocpt=1;; p) inverseCPT=1;; + s) + polyfile=$OPTARG;; ?) printf "error: unknow argument\nuse -h option to see help information\n" exit 1;; @@ -155,7 +159,11 @@ else fi if [[ $coastline == 1 ]]; then - gmt pscoast -R${range} -JL -W0.25p -Dc -A5000 -K -O -P >> $psfile + gmt pscoast -R${range} -JL${plot_para[0]}/${plot_para[1]}/${plot_para[2]}/${plot_para[3]}/${pagePara[0]} -W0.25p -Dc -A5000 -K -O -P >> $psfile + fi + + if [[ $polyfile != 'null' ]]; then + gmt psxy ${polyfile} -JL${plot_para[0]}/${plot_para[1]}/${plot_para[2]}/${plot_para[3]}/${pagePara[0]} -W0.25p -R${range} -L -K -O >> $psfile fi #-C${cptfile}+Uk 使用km(色标单位除1000