# 各緯度の幅で深さ方向の震源の分布を書く# 図の軸や軸の名称、メモリ、グリッドライン等を用意する
# -R x軸、y軸の範囲
# -JX ふつーの図を書くときの図法
# -B ｘ軸の情報/y軸の情報/図のタイトル情報/どこに表示させるか
# x軸の情報 -Ba5(目盛り文字の間隔)f1(目盛の間隔)g5(グリッドの間隔):"軸名":
# 図のタイトル情報 :."図のタイトル":
# wsnE 大文字=全ての軸情報の表示 小文字=目盛り文字を表示しない
# W(左側)S(下側)N（上側）E(右側)
psbasemap -R130/145/0/400 -JX10c/-2c -Ba5f1g5:"longitude":/a100f50g100:"Source Depth(km)"::."Lat.35-36":wSnE -P -Y1.5c -K > depth.ps
# 複数条件はかつなら&&、またはなら||
awk '$2 <= 36 && $2 >= 35 {print $1,$3,($4-4)*0.1}' earth.dat | psxy -R -JX -Sc -O -K -L -W >> depth.ps
# あとは繰り返し
psbasemap -R130/145/0/400 -JX10c/-2c -Ba5f1g5/a100f50g100:."Lat.36-37":wsnE -P -Y3c -O -K >> depth.ps
awk '$2 <= 37 && $2 >= 36 {print $1,$3,($4-4)*0.1}' earth.dat | psxy -R -JX -Sc -O -K -L -W >> depth.ps
psbasemap -R130/145/0/400 -JX10c/-2c -Ba5f1g5/a100f50g100:."Lat.37-38":wsnE -P -Y3c -O -K >> depth.ps
awk '$2 <= 38 && $2 >= 37 {print $1,$3,($4-4)*0.1}' earth.dat | psxy -R -JX -Sc -O -K -L -W >> depth.ps
