采用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

結(jié)果圖：