PSMECA

psmeca − 地図上に発震機構を描画する

psmeca files −Jparameters −Rwest/east/south/north[r] [ −B[p|s]parameters ] [ −Cpen ] [ −Ddepmin/depmax ] [ −Efill] [ −Gfill] [ −H[i][nrec] ] [ −K ] [ −L[pen] ] [ −M ] [ −N ] [ −O ] [ −P ] [ −S<symbol><scale>[/d]] [ −Tnum_of_plane[pen] ] [ −U[/dx/dy/][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 はデータ値を files [又は標準入力] から読み込み発震機構を地図上に 描画する PostScript コードを作る。大部分のオプションは psxy と同じで あ る。 PostScript コードは標準出力に書き出される。

files 1 つ以上のファイル名の一覧。ファイルが与えられないときは、 psmeca は標準入力から読み込む。

−J

地図の図法を選択する。尺度は UNIT/degree 、1:xxxxx 、 又 は UNIT 単 位 の 幅 (オプションを大文字にする) である。 UNIT は .gmtdefaults4 で設定された MEASURE_UNIT によって cm 、インチ 、 又は m になるが、 c, i, m を scale/width の値に付けることによっ て、コマンドライン上で上書き変更することができる。地図の高さ 、 最 大 の大きさ、最小の大きさを指定するには、それぞれ h, +, - を width に付けること。

詳細は psbasemap の man ページを参照すること。

円筒図法:

−Jclon0/lat0/scale (カッシーニ図法)
−Jjlon0/scale (ミラー図法)
−Jmscale (メルカトル図法 - 投影中心はグリニッジ子午線と赤道)
−Jmlon0/lat0/scale (メルカトル図法 - 中心子午線と標準緯度 を 設 定)
−Joalon0/lat0/azimuth/scale (斜めメルカトル図法 - 投影中心と方 位角を設定)
−Joblon0/lat0/lon1/lat1/scale (斜めメルカトル図法 - 2 点を設定)
−Joclon0/lat0/lonp/latp/scale (斜めメルカトル図法 - 中心点と極 を設定)
−Jqlon0/scale (等距円筒図法 (Plate Carree))
−Jtlon0/scale (TM - 横メルカトル図法、赤道は y = 0)
−Jtlon0/lat0/scale (TM - 横メルカトル図法、原点を設定)
−Juzone/scale (UTM - ユニバーサル横メルカトル図法)
−Jylon0/lats/scale (基本円筒図法)

方位図法:

−Jalon0/lat0/scale (ランベルト正積方位図法)
−Jelon0/lat0/scale (正距方位図法)
−Jflon0/lat0/horizon/scale (心射図法)
−Jglon0/lat0/scale (正射図法)
−Jslon0/lat0/[slat/]scale (一般平射図法)

円錐図法:

−Jblon0/lat0/lat1/lat2/scale (アルベルス正積円錐図法)
−Jdlon0/lat0/lat1/lat2/scale (正距円錐図法)
−Jllon0/lat0/lat1/lat2/scale (ランベルト正角円錐図法)

様々な図法:

−Jhlon0/scale (ハンメル図法)
−Jilon0/scale (正弦曲線図法＝サンソン図法)
−Jk[f|s]lon0/scale (エッケルト IV (f) 及び VI (s) 図法)
−Jnlon0/scale (ロビンソン図法)
−Jrlon0/scale (ヴィンケル第 3 式図法)
−Jvlon0/scale (ファン・デア・グリンテン図法)
−Jwlon0/scale (モルワイデ図法)

非地理的図法:

−Jp[a]scale[/origin][r|z] (極座標 (theta,r))
−Jxx-scale[d|l|ppow|t|T][/y-scale[d|l|ppow|t|T]] (線形、対数 、 指数スケーリング)

−R

west, east, south, 及び north は対象とする領域を指定する。 小数値を含む度又は [+-]dd:mm[:ss.xxx][W|E|S|N] のフォーマットで 指 定できる。東西南北の代わりに地図座標の左下と右上を使うときは r を付けること。 2 つの略号 −Rg −Rd は全球領域(それ ぞ れ 経 度 0/360 又は -180/+180 と緯度 -90/+90 )を表す。

−S

データファイルの列の意味を選択する。

同じファイルを断面を描画するのに使うために、深さが 3 列目にある 。それでも、 3 列目に深さのない「古い形式の」 psvelomeca の入力 ファイルを -o オプションによって使うことができる。

−Sascale[c/i][/fontsize[/offset[u]]]

安芸・リチャーズ方式の発震機構。 scale はマグニチュードに比例す る「ビーチボール」の半径のスケーリングを調整する。スケー ル は( c, i, m, 又は p が付かない限り) インチ単位のマグニチュード = 5 に対する大きさである。節面と円周だけを描いてビーチボールを透 明 に するには −T オプションを使うこと。圧縮象限のカラー又はグレー スケールは −G オプションで指定できる。引張象限のカラー又はグ レ ー スケールは −E オプションで指定できる。パラメータは以下の列に あると想定される。

1,2:

イベントの経度、緯度( −: オプションで順番が交換される)

3:

イベントのキロメートル単位の深さ

4,5,6:

度単位の走向、傾斜、レイク

7:

マグニチュード

8,9:

ビーチボールを置く経度、緯度。この列の内容には −C オプ シ ョンが必要である。 0,0 を 8 及び 9 列目に使うとビーチボールを 1 及び 2 列目に与えられた経度、緯度に描画する。 −: オプションで列 (1,2) と (8,9) の順番が交換される。

10:

ビーチボールの上に表示する文字列(オプション)。

−Scscale[c/i][/fontsize[/offset[u]]]

ハ ーバード CMT 方式の発震機構。 scale はマグニチュードに比例す る「ビーチボール」の半径のスケーリングを調整する。スケー ル は( c, i, m, 又は p が付かない限り) インチ単位のマグニチュード = 5 (すなわち M0 = 4.0E23 dynes-cm )に対する大きさである。節面と 円 周 だけを描いてビーチボールを透明にするには −T オプションを使う こと。圧縮象限のカラー又はグレースケールは −G オプションで指 定 で きる。引張象限のカラー又はグレースケールは −E オプションで指 定できる。パラメータは以下の列にあると想定される。

1,2:

イベントの経度、緯度( −: オプションで順番が交換される)

3:

イベントのキロメートル単位の深さ

4,5,6:

平面 1 の走向、傾斜、レイク

7,8,9:

平面 2 の走向、傾斜、レイク

10,11:

dyne-cm 単位のモーメントの仮数と指数

12,13:

ビーチボールを置く経度、緯度。この列の内容には −C オ プ ションが必要である。 0,0 を 12 及び 13 列目に使うとビーチボール を 1 及び 2 列目に与えられた経度、緯度に描画する。 −: オプシ ョ ンで列 (1,2) と (12,13) の順番が交換される。

14:

ビーチボールの上に表示する文字列(オプション)。

−Smscale[c/i][/fontsize[/offset[u]]]

地震モーメントテンソル(ゼロ・トレースのハーバード CMT )。 scale はマグニチュードに比例する「ビーチボール」の半径のスケーリン グ を 調整する。スケールは( c, i, m, 又は p が付かない限り) インチ 単位のマグニチュード = 5 (すなわち地震モーメント テ ン ソ ル = 4.0E23 dynes-cm )に対する大きさである ( −T0 オプションは最良の ダブルカップルを透過的に上書きする)。モーメントテンソルのダブル カ ッ プ ル の 部 分 だ け を 描 画 す る に は −Sdscale[/fontsize[/offset[u]]] を使うこと。モーメントテンソ ル の 非 等 方 な 部 分 だ け( ゼ ロ ・ ト レ ー ス)を描画するには −Szscale[/fontsize[/offset[u]]] を使うこと。圧縮象限のカラー 又 は グレースケールは −G オプションで指定できる。引張象限のカラー 又はグレースケールは −E オプションで指定できる。パラメータは 以 下の列にあると想定される。

1,2:

イベントの経度、緯度( −: オプションで順番が交換される)

3:

イベントのキロメートル単位の深さ

4,5,6,7,8,9:

10* 指数 dynes-cm 単位の mrr, mtt, mff, mrt, mrf, mtf

10:

指数

11,12:

ビーチボールを置く経度、緯度。この列の内容には −C オプ ションが必要である。 0,0 を 11 及び 12 列目に使うとビーチボール を 1 及び 2 列目に与えられた経度、緯度に描画する。 −: オプショ ンで列 (1,2) と (11,12) の順番が交換される。

13:

ビーチボールの上に表示する文字列(オプション)。

−Spscale[c/i][/fontsize[/offset[u]]]

両方の平面の部分的データで与えられた発震機構。 scale はマグニチ ュ ードに比例する「ビーチボール」の半径のスケーリングを調整する 。スケールは( c, i, m, 又は p が付かない限り) インチ単位のマ グ ニチュード = 5 に対する大きさである。圧縮象限のカラー又はグレー スケールは −G オプションで指定できる。引張象限のカラー又はグ レ ー スケールは −E オプションで指定できる。パラメータは以下の列に あると想定される。

1,2:

イベントの経度、緯度( −: オプションで順番が交換される)

3:

イベントのキロメートル単位の深さ

4,5:

平面 1 の走向、傾斜

6:

平面 2 の走向

7:

正/逆断層に対して -1/+1 でなければならない

8:

マグニチュード

9,10:

ビーチボールを置く経度、緯度。この列の内容には −C オプシ ョ ンが必要である。 0,0 を 9 及び 10 列目に使うとビーチボールを 1 及び 2 列目に与えられた経度、緯度に描画する。 −: オプションで 列 (1,2) と (9,10) の順番が交換される。

11:

ビーチボールの上に表示する文字列(オプション)。

−Sxscale[c/i][/fontsize[/offset[u]]]

主軸。 scale はマグニチュードに比例する「ビーチボール」の半径の スケーリングを調整する。スケールは( c, i, m, 又は p が付かな い 限 り) インチ単位のマグニチュード = 5 (すなわち地震スカラーモー メントモーメント = 4*10e+23 dynes-cm )に対する大きさで あ る ( −T0 オプションは最良のダブルカップルを透過的に上書きする)。モー メントテンソルのダブルカップルの部分 だ け を 描 画 す る に は −Syscale[c/i][/fontsize[/offset[u]]] を使うこと。モーメントテン ソルの非等方な部分だけ(ゼロ・ ト レ ー ス) を 描 画 す る に は −Stscale[c/i][/fontsize[/offset[u]]] を使うこと。圧縮象限のカラ ー又はグレースケールは −G オプションで指定できる。引張象限の カ ラ ー又はグレースケールは −E オプションで指定できる。パラメータ は以下の列にあると想定される。

1,2:

イベントの経度、緯度( −: オプションで順番が交換される)

3:

イベントのキロメートル単位の深さ

4,5,6,7,8,9,10,11,12:

T, N, P 軸の値( 10* 指数 dynes-cm ), 方位, プランジ。

13:

指数

14,15:

ビーチボールを置く経度、緯度。この列の内容には −C オ プ ションが必要である。 0,0 を 14 及び 15 列目に使うとビーチボール を 1 及び 2 列目に与えられた経度、緯度に描画する。 −: オプシ ョ ンで列 (1,2) と (14,15) の順番が交換される。

16:

ビーチボールの上に表示する文字列(オプション)。

オプションフラグとそれに付属する引き数の間にスペースを入れないこと。

−B

地 図 の 枠 の 目 盛とケバの間隔を設定する。全ての詳細は psbasemap の man ページを参照すること。

−C[pen][pointsize]

発震機構から入力ファイルの(随意の)文字列の前の最後の 2 列に指定 さ れた経度、緯度までのオフセット。小さな円が始めの位置に描かれ 、線がビーチボールと円を結ぶ。線のスタイルや円の大きさを変え る に は pen や pointsize を指定すること [デフォルトは pen width = 1, color = 0/0/0, texture = solid; pointsize 0 ]。

−D

depmin/depmax depmin と depmax の間のイベントを描画する。

−Efill

引張象限の塗りを選択する。普通は白である。グレースケ ー ル (0−255) 又 は カラー (r/g/b) を設定すること [デフォルトは 255/255/255 ]。

−Gfill

発震機構の塗りを選択する。慣例により、発震機構のビー チ ボ ー ルの圧縮象限が塗られる。グレースケール (0−255) 又はカラー (r/g/b) を設定すること [デフォルトは 255/255/255 ]。

−H

入力ファイルにヘッダ行がある時に使う。 ヘ ッ ダ 行 の 数 は .gmtdefaults4 ファイルを編集することにより変更できる。このオプ ションが使われた場合、 GMT のデフォルトではヘッダ行は 1 行で あ る。入力データだけにヘッダ行が必要な場合は −Hi を使うこと [デフ ォルトでは入力データにヘッダ行があればそれも出力される]。

−K

さらに PostScript コードが後に続くときに使用する [デフォ ル トでは描画システムを終了する]。

−L[pen]

「 ビーチボール」の輪郭を現在のペン( −W 参照)を用いるかペン属性 を設定して描く [デフォルトは width = 1, color = 0/0/0, texture = solid ]。

−N

地図の枠からはみ出す記号をスキップしない [デフォルトでは枠 の内側の記号だけを描画する]。

−O

上書き描画モードを選択する [デフォルトでは新しい描画シス テ ムに初期化する]。

−P

縦置き描画モードを選択する [ GMT のデフォルトは横置きである 。これを変更するには gmtdefaults を参照すること]。

−T[num_of_planes]

節面と透明なバブルの輪郭を描画する。 num_of_planes が
0 ならば両方の節面が描画され 1 ならば第一節面だけが描画され 2 ならば第二節面だけが描画される。

−U

図に Unix システムのタイムスタンプを入れる。図の左下隅に 対 す るスタンプの左下隅の相対的な位置を決めることができる。随意で ラベル又は c (コマンドの文字列を入れる)を付けることができ る 。 GMT のパラメータ UNIX_TIME と UNIX_TIME_POS は外見に影響を及ぼ す。詳細は gmtdefaults の man ページを参照すること。

−V

冗長モードを選択する。標準エラー出力に経過報告を送る[デフォ ルトでは「黙って」実行する]。

−W

pen (ペン属性)はコンマで区切られた width, color, texture の リストで、それぞれはオプションである。 width (幅)は長さの単 位( ポ イ ン ト、センチメートル、インチ)又は faint, thin[ner|nest], thick[er|est], fat[ter|test], obese のいずれかで指定すること が で きる。 color (色)はグレースケール (0−255) 又はカラー (r/g/b, それぞれ範囲は 0−255; h-s-v, 範囲は 0−360, 0−1, 0−1; c/m/y/k, そ れ ぞ れ 範 囲 は 0−100%; 有効なカラー名) により指定できる。 texture (テクスチャ)はダッシュ ‘-’ と点 ‘.’ の組合せである。

−X −Y

図の原点を現在の原点に対して (x-shift,y-shift) だけ動 か す 。随意で長さの単位 (c, i, m, p) を付けることができる。値の前 に a を付けると描画後に原点を元の位置に戻し、値の前に r [デフォ ルト]を付けると現在の原点を新しい位置として再設定する。 −O を使 うとデフォルトの (x-shift,y-shift) が (0,0) になり、そうでな い ときは (r1i, r1i) 又は (r2.5c, r2.5c) になる。代わりに、 c は現 在のページサイズに基づくページの中心を使って図を中央の座 標 (x 又は y) に動かす。

−Zcptfile

カ ラーパレットファイルを与え 3 列目の z の値で圧縮部の色が決ま るようにする。

−z

ゼロ・トレース・モーメントテンソルを上書きする。

−a[size/[P_axis_symbol[T_axis_symbol]]]

P 及び T 軸を計算して記号を用いて描画する。随意で size と (別々 の) P 及び T 軸の記号を以下から指定すること。 (c) 円, (d) 菱形, (h) 六角形, (i) 逆三角形, (p) 点, (s) 正方形, (t) 三角 形, (x) 十字 [デフォルトは 0.2c/cc 又は 0.08i/cc ]。

−efill

T 軸の記号の塗りを選択する。グレースケール (0−255) 又は カラー (r/g/b) を設定すること。デフォルトは白。

−gfill

P 軸の記号の塗りを選択する。グレースケール (0−255) 又は カラー (r/g/b) を設定すること。デフォルトは黒。

-o

3 列目に深さがない psvelomeca 入力フォーマットを使う。

−p[pen]

P 軸の線を現在のペン( −W 参照)を用いるかペン属性を設定して描く 。

−t[pen]

T 軸の線を現在のペン( −W 参照)を用いるかペン属性を設定して描 く 。

−:

入出力の(経度、緯度)と(緯度、経度)を入れ替える[デフォルト は(経度、緯度)]。 i を付けると入力のみが、 o を付けると出力のみ が入れ替わる[デフォルトでは両方入れ替わる]。

−c

図のコピー枚数を指定する[デフォルトは 1 ]。

下のファイルは正断層型の CMT メカニズムを与える。
psmeca −R239/240/34/35.2 −Jm4 −Sc0.4 −H1 <<END>! test.ps
lon lat depth str dip slip st dip slip mant exp plon plat
239.384 34.556 12. 180 18 -88 0 72 -90 5.5 0 0 0
END

GMT(l), psbasemap(l), psxy(l)

Bomford, G., Geodesy, 4th ed., Oxford University Press, 1980.
Aki, K. and P. Richards, Quantitative Seismology, Freeman, 1980.
F. A. Dahlen and Jeoren Trom, Theoretical Seismology, Princeton, 1998, p.167.
Cliff Frohlich, Cliff’s Nodes Concerning Plotting Nodal Lines for P, Sh and Sv
Seismological Research Letters, Volume 67, Number 1, January-February, 1996
Thorne Lay, Terry C. Wallace, Modern Global Seismology, Academic Press, 1995, p.384.
W.H. Press, S.A. Teukolsky, W.T. Vetterling, B.P. Flannery, Numerical Recipes in C, Cambridge University press (routine jacobi)

Genevieve Patau
CNRS UMR 7580
Seismology Dept.
Institut de Physique du Globe de Paris
(patau@.ipgp.jussieu.fr)