画震源机制剖面的GMT工具
2013-04-24 14:30阅读：
画震源机制剖面的GMT工具：pscoupe画震源机制剖面的GMT工具，也可以画其他类似东西的剖面，比如地震深度剖面。还没有具体应用，用到了会把图发上来。
语法：
pscoupe files
−Jparameters
−Rwest/east/south/north[r]
−Aparameters [
−B[p|s]parameters
] [ −Efill ] [
−Gfill ] [
−H[i][nrec] ] [
−K ] [ −L[pen] ] [
−M ] [ −N ] [ −O
] [ −P ] [
−S[/d] ] [
−s[/d] ] [
−Tn ] [
−U[just/dx/dy/][c|label]
] [ −V ] [ −Wpen ] [
−X[a|c|r][x-shift[u]]
] [
−Y[a|c|r][y-shift[u]]
] [ −Zcpt ] [
−a[size/[P_symbol/[T_symbol]]]
] [ −gfill ] [
−efill ] [
−ppen ] [ −tpen
] [ −:[i|o] ] [
−ccopies ]
GMT的震源机制工具，比如psmeca,
pscoupe一般默认是绘下半球投影。还要注意的是在利用pscoupe绘震源机制剖面的时候投影面是不同于平面上的。
-J 控制投影方式
-R 坐标范围
-A 选择剖面
-Aalon1/lat1/lon2/lat2/dip/p_width/dmin/dmax[f]
lon1,lat1,lon2,lat2分别是两个端点的经纬度
dip是投影平面的倾角
p_width是投影剖面两侧的宽度
dmin和dmax是距离水平面最小和最大的距离
加上f可以获得框格
−Ablon1/lat1/strike/p_length/dip/p_width/dmin/dmax[f]
−Acx1/y1/x2/y2/dip/p_width/dmin/dmax[f]
类似于-Aa，但用于笛卡尔坐标系下
−Adx1/y1/strike/p_length/dip/p_width/dmin/dmax[f]
类似于-Ab，但用于笛卡尔坐标系下
-S 选择数据文件
−Sascale[/fontsize[/offset[u]]]
用于AKI的震源机制约定，第1～10列分别为经度、纬度、深度、走向、倾角、滑动角、震级、
0、
0、文字说明
−Scscale
用于哈佛的震源机制约定，第1～14列分别为经度、纬度、深度、节面1的走向、倾角、滑动角、节面2的走向、倾角、滑动角、地震矩的小数和指数部分、0、0、文字说明
−Spscale[/fontsize[/offset[u]]]
当只有部分震源机制解参数的时候，可用这个。第1～11列分别为经度、纬度、深度、节面1的走向、倾角、节面2走向、+1（逆断层）/-1（正断层）、震级、0、0、文字说明
−Sm|d|zscale[/fontsize[/offset[u]]]
用于绘地震矩张量解。第1～13列分别为经度、纬度、深度、mrr, mtt, mff, mrt, mrf, mtf in
10*exponent dynes-cm、指数部分、0、0、文字说明
−Sxscale[/fontsize[/offset[u]]]
用于已知力轴参数来绘震源机制解。第1～16列分别为经度、纬度、深度、value (in 10*exponent dynes-cm),
azimuth, plunge of the T, N, and P axes、指数部分、0、0、文字说明
−ssymbol[size[/fontsize[/offset[u]]]
用其他图形来替代震源机制球。第1～4列分别为经度、纬度、深度、说明。
-B设置边框刻度
−Efill膨胀象限的颜色设定
−Gfill压缩象限的颜色
−L[pen]画笔属性
−M忽略震级大小
−T[num_of_planes]控制所绘节面的数量，0-两个节面都绘，1-节面1，2-节面2，默认是0
−a[size/[P_axis_symbol/[T_axis_symbol]]]，计算并绘P、T轴的位置，符号类型有：(c)
circle, (d) diamond, (h) hexagon,
(i) inverse triangle, (p)point,
(s) square, (t) triangle,
(x) cross. [Defaults:
0.2c/c/c or
0.08i/c/c.]
−efill 填充T轴
−gfill填充P轴
转自http://hi.baidu.com/purplepool/item/4d2a961493945e4ee75e0617
