Update gmtsph-regional.sh

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