GMT 绘图模板

Batch scripts for GMT v6.4.0

1×1

图1. GMT绘图模板 (1×1)

gmt begin

REM 设置全局字体
gmt set FONT 10p,4
REM 设置底图边框的类型
gmt set MAP_FRAME_TYPE plain 
REM 设置底图边框的画笔属性
gmt set MAP_FRAME_PEN 0.5p,black
REM 设置标注的开始位置与刻度尾端间的距离
gmt set MAP_ANNOT_OFFSET_PRIMARY 5p
REM 设置刻度的画笔属性
gmt set MAP_TICK_PEN 0.5p,black
REM 设置刻度的长度属性
gmt set MAP_TICK_LENGTH_PRIMARY 2p/2p
REM 设置网格线的画笔属性
gmt set MAP_GRID_PEN 0.5p,black
REM 设置可输出中文字体
gmt set PS_CHAR_ENCODING Standard+

REM 绘图
gmt figure Test_1 png E+600

gmt makecpt -Cgoogle/turbo -T0/60/3 -D
gmt grdimage TEC.nc -JX6.36c/3.28c -R-180/180/-90/90 -Blbrt -C -E1000
gmt coast -BWSrt -Bxa90 -Bya45 -A100000 -W0.2p
echo 0 90 10p,4 0 2022-01-01 | gmt text -D0p/8p -N -F+f+a
gmt colorbar -C -Bxa20 -DJRM+w3.28/0.25+o0.2c/0c+v
echo 180 90 TECU| gmt text -D10p/8p -N

gmt end

REM 等值线
REM gmt grdcontour TEC.nc -W0.1p,gray -C3
REM 单个点
REM echo -57.9 -63.3 | gmt plot -Sa6p -W0.2p,black -Glightred
REM echo -57.9 -63.3 OHI3 | gmt text -D12.5p/0p -F+f8p,4,black -N 
REM 多个点
REM gmt plot P1.txt -Sc3p -W0.1p,black -Glightblue -C
REM gmt text P2.txt -D0p/-6p -F+f8p,4,white -N
REM 区域边框
REM echo  30 -45 >  area1.txt
REM echo 160 -45 >> area1.txt
REM echo 160  45 >> area1.txt
REM echo  30  45 >> area1.txt
REM gmt plot area1.txt -W0.5p,red -A -L
REM del area1.txt

2×1

图2. GMT绘图模板 (2×1)

gmt begin

REM 设置全局字体
gmt set FONT 10p,4
REM 设置底图边框的类型
gmt set MAP_FRAME_TYPE plain 
REM 设置底图边框的画笔属性
gmt set MAP_FRAME_PEN 0.5p,black
REM 设置标注的开始位置与刻度尾端间的距离
gmt set MAP_ANNOT_OFFSET_PRIMARY 5p
REM 设置刻度的画笔属性
gmt set MAP_TICK_PEN 0.5p,black
REM 设置刻度的长度属性
gmt set MAP_TICK_LENGTH_PRIMARY 2p/2p
REM 设置网格线的画笔属性
gmt set MAP_GRID_PEN 0.5p,black
REM 设置可输出中文字体
gmt set PS_CHAR_ENCODING Standard+

REM 绘图
gmt figure Test_2 png E+600
gmt subplot begin 2x1 -JX6.36c/3.28c -R-180/180/-90/90 -Fs6.36c/3.28c -M0c/0.1c
gmt makecpt -Cgoogle/turbo -T0/60/3 -D

REM 01
gmt subplot set 0,0
gmt grdimage TEC.nc -JX6.36c/3.28c -R-180/180/-90/90 -Blbrt -C -E1000
gmt coast -BWbrt -Bya45 -A100000 -W0.2p
echo 0 90 10p,4 0 2022-01-01 | gmt text -D0p/8p -N -F+f+a
gmt colorbar -C -Bxa20 -DJRM+w3.28/0.25+o0.2c/0c+v
echo 180 90 TECU| gmt text -D10p/8p -N

REM 02
gmt subplot set 1,0
gmt grdimage TEC.nc -JX6.36c/3.28c -R-180/180/-90/90 -Blbrt -C -E1000
gmt coast -BWSrt -Bxa90 -Bya45 -A100000 -W0.2p
echo 0 90 10p,4 0 2022-01-01 | gmt text -D0p/8p -N -F+f+a
gmt colorbar -C -Bxa20 -DJRM+w3.28/0.25+o0.2c/0c+v
echo 180 90 TECU| gmt text -D10p/8p -N

gmt subplot end
gmt end

REM 等值线
REM gmt grdcontour TEC.nc -W0.1p,gray -C3
REM 单个点
REM echo -57.9 -63.3 | gmt plot -Sa6p -W0.2p,black -Glightred
REM echo -57.9 -63.3 OHI3 | gmt text -D12.5p/0p -F+f8p,4,black -N 
REM 多个点
REM gmt plot P1.txt -Sc3p -W0.1p,black -Glightblue -C
REM gmt text P2.txt -D0p/-6p -F+f8p,4,white -N
REM 区域边框
REM echo  30 -45 >  area1.txt
REM echo 160 -45 >> area1.txt
REM echo 160  45 >> area1.txt
REM echo  30  45 >> area1.txt
REM gmt plot area1.txt -W0.5p,red -A -L
REM del area1.txt

1×2

图3. GMT绘图模板 (1×2)

gmt begin

REM 设置全局字体
gmt set FONT 12p,4
REM 设置底图边框的类型
gmt set MAP_FRAME_TYPE plain 
REM 设置底图边框的画笔属性
gmt set MAP_FRAME_PEN 0.5p,black
REM 设置标注的开始位置与刻度尾端间的距离
gmt set MAP_ANNOT_OFFSET_PRIMARY 5p
REM 设置刻度的画笔属性
gmt set MAP_TICK_PEN 0.5p,black
REM 设置刻度的长度属性
gmt set MAP_TICK_LENGTH_PRIMARY 2p/2p
REM 设置网格线的画笔属性
gmt set MAP_GRID_PEN 0.5p,black
REM 设置可输出中文字体
gmt set PS_CHAR_ENCODING Standard+

REM 绘图
gmt figure Test_3 png E+600
gmt subplot begin 1x2 -JX6.36c/3.28c -R-180/180/-90/90 -Fs6.36c/3.28c -M-0.22c/0.05c
gmt makecpt -Cgoogle/turbo -T0/60/3 -D

REM 01
gmt subplot set 0,0
gmt grdimage TEC.nc -JX6.36c/3.28c -R-180/180/-90/90 -Blbrt -C -E1000
gmt coast -BWSrt -Bxa90 -Bya45 -A100000 -W0.2p
echo 0 90 12p,4 0 2022-01-01 | gmt text -D0p/10p -N -F+f+a

REM 02
gmt subplot set 0,1
gmt grdimage TEC.nc -JX6.36c/3.28c -R-180/180/-90/90 -Blbrt -C -E1000
gmt coast -BlSrt -Bxcmark_lon.txt -A100000 -W0.2p
echo 0 90 12p,4 0 2022-01-01 | gmt text -D0p/10p -N -F+f+a
gmt colorbar -C -Bxa20 -DJRM+w3.28/0.25+o0.2c/0c+v
echo 180 90 TECU| gmt text -D10p/10p -N

gmt subplot end
gmt end

REM 等值线
REM gmt grdcontour TEC.nc -W0.1p,gray -C3
REM 单个点
REM echo -57.9 -63.3 | gmt plot -Sa6p -W0.2p,black -Glightred
REM echo -57.9 -63.3 OHI3 | gmt text -D12.5p/0p -F+f8p,4,black -N 
REM 多个点
REM gmt plot P1.txt -Sc3p -W0.1p,black -Glightblue -C
REM gmt text P2.txt -D0p/-6p -F+f8p,4,white -N
REM 区域边框
REM echo  30 -45 >  area1.txt
REM echo 160 -45 >> area1.txt
REM echo 160  45 >> area1.txt
REM echo  30  45 >> area1.txt
REM gmt plot area1.txt -W0.5p,red -A -L
REM del area1.txt

2×2

图4. GMT绘图模板 (2×2)

gmt begin

REM 设置全局字体
gmt set FONT 12p,4
REM 设置底图边框的类型
gmt set MAP_FRAME_TYPE plain 
REM 设置底图边框的画笔属性
gmt set MAP_FRAME_PEN 0.5p,black
REM 设置标注的开始位置与刻度尾端间的距离
gmt set MAP_ANNOT_OFFSET_PRIMARY 5p
REM 设置刻度的画笔属性
gmt set MAP_TICK_PEN 0.5p,black
REM 设置刻度的长度属性
gmt set MAP_TICK_LENGTH_PRIMARY 2p/2p
REM 设置网格线的画笔属性
gmt set MAP_GRID_PEN 0.5p,black
REM 设置可输出中文字体
gmt set PS_CHAR_ENCODING Standard+

REM 绘图
gmt figure Test_4 png E+600
gmt subplot begin 2x2 -JX6.36c/3.28c -R-180/180/-90/90 -Fs6.36c/3.28c -M-0.22c/0.05c
gmt makecpt -Cgoogle/turbo -T0/60/3 -D

REM 01
gmt subplot set 0,0
gmt grdimage TEC.nc -JX6.36c/3.28c -R-180/180/-90/90 -Blbrt -C -E1000
gmt coast -BWbrt -Bya45 -A100000 -W0.2p
echo 0 90 12p,4 0 2022-01-01 | gmt text -D0p/10p -N -F+f+a

REM 02
gmt subplot set 0,1
gmt grdimage TEC.nc -JX6.36c/3.28c -R-180/180/-90/90 -Blbrt -C -E1000
gmt coast -Blbrt -A100000 -W0.2p
echo 0 90 12p,4 0 2022-01-01 | gmt text -D0p/10p -N -F+f+a
gmt colorbar -C -Bxa20 -DJRM+w3.28/0.25+o0.2c/0c+v
echo 180 90 TECU| gmt text -D10p/10p -N

REM 03
gmt subplot set 1,0
gmt grdimage TEC.nc -JX6.36c/3.28c -R-180/180/-90/90 -Blbrt -C -E1000
gmt coast -BWSrt -Bxa90 -Bya45 -A100000 -W0.2p
echo 0 90 12p,4 0 2022-01-01 | gmt text -D0p/10p -N -F+f+a

REM 04
gmt subplot set 1,1
gmt grdimage TEC.nc -JX6.36c/3.28c -R-180/180/-90/90 -Blbrt -C -E1000
gmt coast -BlSrt -Bxcmark_lon.txt -A100000 -W0.2p
echo 0 90 12p,4 0 2022-01-01 | gmt text -D0p/10p -N -F+f+a
gmt colorbar -C -Bxa20 -DJRM+w3.28/0.25+o0.2c/0c+v
echo 180 90 TECU| gmt text -D10p/10p -N

gmt subplot end
gmt end

REM 等值线
REM gmt grdcontour TEC.nc -W0.1p,gray -C3
REM 单个点
REM echo -57.9 -63.3 | gmt plot -Sa6p -W0.2p,black -Glightred
REM echo -57.9 -63.3 OHI3 | gmt text -D12.5p/0p -F+f8p,4,black -N 
REM 多个点
REM gmt plot P1.txt -Sc3p -W0.1p,black -Glightblue -C
REM gmt text P2.txt -D0p/-6p -F+f8p,4,white -N
REM 区域边框
REM echo  30 -45 >  area1.txt
REM echo 160 -45 >> area1.txt
REM echo 160  45 >> area1.txt
REM echo  30  45 >> area1.txt
REM gmt plot area1.txt -W0.5p,red -A -L
REM del area1.txt