pscoupe:画震源机制剖面的GMT工具

sub

这是画震源机制剖面的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<format><scale>[/d] ] [ −s<symbol><size>[/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轴