GMT5绘制折线图--以陆地水储量为例

采用GMT5绘制折线图，具体代码如下：

set PS="hydro.ps"

REM set parameters for plot

set R=98/108/26/35

set ticks=f60ma120m

set J=M12c

REM gmt plot data_year.dat -W1p,black -O -K  >>%ps%

REM gmt plot data_re.dat -W1p,red -O -K>>%ps%

REM psbasemap -R%R% -J%J% -Ba2f1 -K  -V -U -P -Xc  -Yc>%ps%

REM pscoast  -J%J% -R%R%  -W1 -I1 -O -K -P -Xc  -Yc>>%ps%

psbasemap -R1979/2017/0.7/1 -JX14c/5c -Bxa5f5+l"Year"  -BWenS -Bya0.1+l"Sea ice concentration" -K -Xc -Yc>%ps%【此处采用线性坐标系】

REM psxy data_re.dat -W1p,black -O -K -Xc -Yc>>%PS%

gmt psxy data_re.dat   -R1978/2018/0.7/1 -JX14c/5c  -W0.5,red -Xc -Yc  -O -K >> %ps%

gmt psxy data_year.dat   -R1978/2018/0.7/1 -JX14c/5c  -W0.5,red -Xc -Yc  -O -K >> %ps%

gmt psxy data_year.dat   -R1978/2018/0.7/1 -JX14c/5c -Sc0.1c  -W0.5,red -Xc -Yc  -O -K >> %ps%

psbasemap -R1979/2017/0.7/1 -JX14c/5c -Bxa5f5+l"Year"  -BWenS -Bya0.1+l"Sea ice concentration" -O -K -Xc -Y6.5c>>%ps%

gmt psxy data_re.dat   -R1978/2018/0.7/1 -JX14c/5c  -W0.5,red -X0c -Y0c  -O -K >> %ps%

gmt psxy data_year.dat   -R1978/2018/0.7/1 -JX14c/5c  -W0.5,red -X0c -Y0c  -O -K >> %ps%

gmt psxy data_year.dat   -R1978/2018/0.7/1 -JX14c/5c -Sc0.1c  -W0.5,red -X0c -Y0c  -O -K >> %ps%

psbasemap -R2002/2022.5/-10/10 -JX14c/5c -Bxa2f5+l"Year"  -BWenS -Bya3+l"TWSC/mm" -O -K -Xc -Y-13c>>%ps%

gmt psxy oo.txt  -R2002/2022.5/-10/10 -JX14c/5c  -W0.5,blue -X0c -Y0c  -O -K >> %ps%

gmt psxy oo.txt  -R2002/2022.5/-10/10 -JX14c/5c -Sc0.08c  -W0.5,red -X0c -Y0c  -O >> %ps%

psconvert %PS% -A  -Tg  -E300

结果图：