今天要画震源机制解，在网上搜了一下，还是不太懂，自己琢磨了一上午总算画出来了，具体方法如下：

command：psmeca

命令对应选项有：

psmeca files −Jparameters

−Rwest/east/south/north[r]

[ −B[p|s]parameters ] [

−C[pen][Ppointsize] ] [

−Ddepmin/depmax ] [ −Efill] [

−Gfill] [ −H[i][nrec] ] [

−K ] [ −L[pen] ] [ −M ] [ −N ] [

−O ] [ −P ] [

−S[/d]]

[ −Tnum_of_plane[pen] ] [

−U[just/dx/dy/][c|label]

] [ −V ] [ −Wpen ] [

−X[a|c|r][x-shift[u]] ] [

−Y[a|c|r][y-shift[u]] ] [

−Zcptfile] [ −z ] [

−a[size[P_symbol[T_symbol]]] ] [

−gfill ] [ −efill ] [ −o ] [

−ppen ] [ −tpen ] [

−:[i|o] ] [ −ccopies ]

普通的命令选项我就不解释了，下面说一下，psmeca对应的特殊的命令选项：

[ −Efill] [

−Gfill]：对应震源球正负初动的两种颜色E--负初动Default is

255/255/255，G--正初动Default is 0/0/0，

[ −C[pen][Ppointsize]

]：files中的最后两列是震源球位置坐标，-C是对应坐标偏移量，Defaults: pen

width = 1, color = 0/0/0, texture = solid; pointsize 0

[

−S[/d]]

画震源球是对应的输入文件有多种，对不同的文件不同。

具体如下：

−Sascale[/fontsize[/offset[u]]]：AKI

& RICHARD规定的震源机制解，对应文件格式为：

X,

Y, depth, strike, dip, rake, mag, newX, newY, event_title

35.74 141.36

19.0 17 26 91 5.4 35.74 141.36 C201106210849A

−Scscale[/fontsize[/offset[u]]]：cmt规定：对应文件格式为：

X, Y,

depth, strike1, dip1, rake1, strike2, dip2, rake2, moment, newX,

newY, event_title

35.74 141.36

19.0 17 26 91 198 64 -90

5.4 35.74 141.36 C201106210849A

−Sm|d|zscale[/fontsize[/offset[u]]矩张量，对应文件格式为：

X, Y, depth, mrr, mtt, mff, mrt, mrf, mtf, exp, newX, newY,

event_title

35.74 141.36 19.0 -0.612

0.024 0.588 0.153 0.464 0.202 24

35.74 141.36 C201106210849A

−Spscale[/fontsize[/offset[u]]]对应文件格式为：

X, Y, depth, strike1, dip1, strike2, fault, mag, newX, newY,

event_title

35.74 141.36

19.0 17 26 198 +1 5.4 35.74 141.36 C201106210849A

fault = -1/+1 for a normal/inverse fault

−Sx|y|tscale[/fontsize[/offset[u]]]Harvard

CMT，对应文件格式为：

(x) X, Y, depth, T_value, T_azim, T_plunge,

N_value, N_azim, N_plunge

P_value, P_azim, P_plunge, exp, newX, newY,

event_title

。。。。

(y) X, Y, depth, T_value, T_azim, T_plunge, N_value, N_azim,

N_plunge

P_value, P_azim, P_plunge, exp, newX, newY,

event_title

。。。。

-Tn[/]

n=1：画第一个节点平面

n=2：画第二个节点平面

n=3：两个都画

example：

psmeca cmt_gmt.lst -R96/100/21/36 -JB101/28.5/26/31/4 -D0/800

-E255/255/255 -G 255/0/0 -Sm0.8/-1 -T0 -W1/255/0/0 -O -P

>>sh.ps