# Nature论文中的三维图原来是这么画的

## 绘图过程简介

GMT绘图复杂组合图，建议采用先拆分后组合的策略！

### Figure1(a,b)

Figure 1 (a)Figure 1 (b)

## 图1的不完整绘图代码示例

function preset(){    path_data=../../data/raw    bathy=${path_data}/bathy/bathy_Longqi_cart.nc basalt=${path_data}/sites/site_basalt_cart_zhao_bathy.txt    # 数据范围    zmin=-6    zmax=-1.5    zrange=$zmin/$zmax    xmin=-18.35    xmax=33.35    ymin=-29.9826    ymax=11.653     range=$xmin/$xmax/$ymin/$ymax     range_xyz=${range}/$zrange    xmin_eq=0    xmax_eq=20    ymin_eq=-10    range_xyz_eq=$xmin_eq/$xmax_eq/$ymin_eq/$ymax_eq/$zmin_eq/$zmax_eq    # 视角    angle_h=-40     angle_v=45    angle_view=${angle_h}/$angle_v    # 图片大小    width_x_fig=8}gmt begin Figure1 pdf     # (a)    move_x=0    move_y=12    # gmt basemap -JM$figwidth_inset -R$range_inset -Bwsen  --MAP_FRAME_TYPE=plain    # gmt basemap     gmt coast -Ggray -JL$lon0_inst/$lat0_inst/$lat_min_inset/$l    gmt grdcontour $bathy_AUV -C0.02 -A+f3p,Helvetica+ggreen -W0.2p -GD500e -p # 采样点 symbol_sample=c gmt psxy$path_data/sites/M-active.txt -S$symbol_sample$si    # legend    # (d)    move_x=-1    move_y=-4    gmt basemap -JX$width_x_fig_eq/$width_y_fig_eq -JZ$width_z gmt grdview$bathy -JZ -p -C$bathy_cpt -N${zmin_eq}+g$color_profile -Qi -I0+a45 -Wf0.5p,black awk '{print$1/1000.0, $2/1000.0}'$earthquake  | gmt grdtrack -G$bathy | gmt psxyz -JZ -p -Sc$size_eq -Ggreen -Wblack    # 东西方向    gmt basemap -JX$width_x_fig_eq/$width_z_fig_eq -JZ$width_y_# boundary between title and legend# symbol space legend type length (-) thickness,color space LegendName# N 3H 10p Helvetica Legend of figure (c)N 2G${gap_rows}cS        ${left_symbol}c a 0.3c$color_active 0.2p,$color_active${symbol_text}c  Active hydrothermal fieldS        1c   c      ${size_symbol}c$color_basalt 0gmt endrm tmp* gmt.history box.eq 

http://wap.sciencenet.cn/blog-621290-1244180.html

## 全部精选博文导读

GMT+8, 2020-10-23 05:45