diff --git a/gmtsph-regional.sh b/gmtsph-regional.sh index b52e07a..9f6a6fc 100755 --- a/gmtsph-regional.sh +++ b/gmtsph-regional.sh @@ -107,12 +107,12 @@ else if [[ $plotgrad == 1 ]]; then gradfile=usergradient.nc gmt grdgradient ${data} -G${gradfile} -Nt -A0/45 - gmt grdimage ${data} -R${range} -C${cptfile} -I${gradfile} -Bxag+l"${labels[0]}" -Byag+l"${labels[1]}" -JL${plot_para[0]}/${plot_para[1]}/${plot_para[2]}/${plot_para[3]}/2.5i -K -P > $psfile + gmt grdimage ${data} -R${range} -C${cptfile} -I${gradfile} -Bxag+l"${labels[0]}" -Byag+l"${labels[1]}" -JL${plot_para[0]}/${plot_para[1]}/${plot_para[2]}/${plot_para[3]}/2i -K -P > $psfile else - gmt grdimage ${data} -R${range} -C${cptfile} -Bxag+l"${labels[0]}" -Byag+l"${labels[1]}" -JL${plot_para[0]}/${plot_para[1]}/${plot_para[2]}/${plot_para[3]}/2.5i -K -P > $psfile + gmt grdimage ${data} -R${range} -C${cptfile} -Bxag+l"${labels[0]}" -Byag+l"${labels[1]}" -JL${plot_para[0]}/${plot_para[1]}/${plot_para[2]}/${plot_para[3]}/2i -K -P > $psfile fi #-C${cptfile}+Uk 使用km(色标单位除1000) - gmt psscale -Dx0.3i/-0.3i+w1.85i/0.08i+h -C${cptfile} -Bxa -By+l${unit} -O >> $psfile + gmt psscale -Dx0.06i/-0.3i+w1.85i/0.08i+h -C${cptfile} -Bxa -By+l${unit} -O >> $psfile gmt psconvert $psfile -A -TEG -E300 # 删除临时文件 使用linux终端rm命令 rm $psfile gmt.history gmt.conf