德克萨斯大学埃尔帕索地质研究所GMT实例之五

huolongguo

本例是绘制重力拟合剖面图，搞重磁的朋友经常干的活。

bat代码：

echo 0.0 14 13 0 2 1 Calculated >tmp
echo 0.0 13.5 13 0 2 1 Observed >>tmp
pstext tmp -R-3/28/0/18 -Jx1 -X3.0 -K -V> dptal.ps

rem Plot symbols for observed and calculated
echo -0.4 14.2| psxy -R -Jx -L -St0.2 -G0/255/0 -K -O -V>> dptal.ps

echo -0.4 13.7 | psxy -R -Jx -L -Sc0.2 -G255/0/0 -K -O -V>> dptal.ps

rem # Plot observed and calculated gravity
psbasemap -Jx0.01/0.015 -R0/2346.7/-350/0 -X1.5 -Y10.0 -Bf20a200:.""::"":/f10a100:"Bouguer Gravity (mGals)":WSne -K -O >> dptal.ps
psxy deep.obs -R -Jx -St0.12 -G0/255/0 -O -K >> dptal.ps
psxy deep.calc -R -Jx -Sc0.12 -G255/0/0 -O -K >> dptal.ps

psbasemap -Jx0.01/0.05 -R0/2346.7/-100/0 -Y-6.75 -Bf20a200:"Distance (km)":/f4.0a20:"Depth (km)":WSne -K -O >> dptal.ps

rem #Plot bodies

psxy dpuc.tal -R -Jx -W6 -G0/245/245 -L -M -O -K >> dptal.ps
psxy dplcfill.tal -R -Jx -W6 -G0/150/245 -L -M -O -K >> dptal.ps
psxy dplcWY.tal -R -Jx -W6 -G194/153/230 -L -M -O -K >> dptal.ps
psxy dpmanCO.tal -R -Jx -W6 -G205/0/137 -L -M -O -K >> dptal.ps
psxy dpmanSlab.tal -R -Jx -W6 -G205/100/137 -L -M -O -K >> dptal.ps
psxy dpmanProt.tal -R -Jx -W6 -G205/0/0 -L -M -O -K >> dptal.ps
psxy dpmanArchean.tal -R -Jx -W6 -G105/0/137 -L -M -O -K >> dptal.ps
psxy dpbatholith.tal -R -Jx -W6 -G255/0/255 -L -M -O -K >> dptal.ps
psxy dpintrusion.tal -R -Jx -W6 -G255/0/255 -L -M -O -K >> dptal.ps
psxy dp245basin.tal -R -Jx -W6 -G255/255/0 -L -M -O -K >> dptal.ps
psxy dp252basin.tal -R -Jx -W6 -G255/255/0 -L -M -O -K >> dptal.ps
psxy dp260basin.tal -R -Jx -W6 -G255/255/0 -L -M -O >> dptal.ps

ps2raster dptal.ps -A -P -Tg

del .gmt*
gsview32 dptal.ps

效果图：