#!/bin/bash
R=70/110/15/45             # 小区域地图的范围 - R
J=M20c                      # 小区域的投影方式 - J
PS=map.ps

# 打开 GMT
gmt psxy -R$R -J$J -T -K > $PS

# 绘制小区域的地图
gmt pscoast -R$R -J$J -B1 -Df -W2p -I2p -N1/2p -A1000 -K -O >> $PS
# ...
# 此处可以任意添加其他绘图命令
# ...

# 开始绘制 insert map
Rg=72/135/16/55     # 大区域地图的范围 - R
Jg=M7c              # 大区域地图的投影方式 - J
# 绘制大区域地图的海岸线及边框
gmt pscoast -R$Rg -J$Jg -B0 -B+gwhite -Df -N1 -W -A5000 -K -O --MAP_FRAME_TYPE=plain >> $PS
# 利用 psbasemap 在大区域地图内绘制小区域所在的方框
gmt psbasemap -R$Rg -J$Jg -D$R -F+p2p,blue -K -O >> $PS

# 结束 GMT 绘图
gmt psxy -R$R -J$J -T -O >> $PS