#!/bin/sh
#
#   make_mecmap.plで作成したサンプル
#
#   make_mecmap.pl -start:200010010000 -end:200010102359 -s:34.5 -n:36 -w:132.5 -e:134.5 -d:0 -D:100 -m:4.0 -M:8.0 -o:mec.ps -path:data > tmp.sh
#   の出力です。
#
gmtdefaults -D > .gmtdefaults
gmtset DOTS_PR_INCH 1200
#--- topography ---
pscoast -JL133.5/35.2/35.2/35.6/16 -R132.5/134.5/34.5/36.0 -Ba1.0f0.5WeSn -Na -Df -G255/235/205 -S224/255/255 -W1 -P -X2.0  -Y12.8 -K -U/-1.0/-1.5/:"REDC" > mec.ps
psxy /jma/data/map/trench_data_jma -JL -R -M -W1 -P -O -K -: >> mec.ps
psxy /jma/data/map/border_data -: -JL -R -M -W1 -P -O -K >> mec.ps
#--- mechanism distribution ---
psmeca -R -Sa0.750000 -JL -W1 -P -G255/0/0   -O -K << END >> mec.ps
133.351700 35.271000 8.960000 130 83 -6 7.3 133.351700 35.271000
END
psmeca -R -Sa0.750000 -JL -W1 -P -G255/0/255   -O -K << END >> mec.ps
133.351000 35.276800 14.510000 314 60 6 4.6 133.351000 35.276800
END
psmeca -R -Sa0.750000 -JL -W1 -P -G255/0/255   -O -K << END >> mec.ps
133.411500 35.200500 12.490000 341 77 -13 4.1 133.411500 35.200500
END
psmeca -R -Sa0.750000 -JL -W1 -P -G255/0/0   -O -K << END >> mec.ps
133.427000 35.189800 8.950000 42 80 176 4.4 133.427000 35.189800
END
psmeca -R -Sa0.750000 -JL -W1 -P -G255/0/0   -O -K << END >> mec.ps
133.323700 35.314200 7.960000 78 70 167 4.2 133.323700 35.314200
END
psmeca -R -Sa0.750000 -JL -W1 -P -G255/0/0   -O -K << END >> mec.ps
133.311700 35.358200 7.680000 354 78 -6 4.3 133.311700 35.358200
END
psmeca -R -Sa0.750000 -JL -W1 -P -G255/0/255   -O -K << END >> mec.ps
133.423200 35.193500 10.430000 242 47 -158 4.0 133.423200 35.193500
END
psmeca -R -Sa0.750000 -JL -W1 -P -G255/0/255   -O -K << END >> mec.ps
133.344500 35.270500 13.390000 337 41 -5 4.0 133.344500 35.270500
END
psmeca -R -Sa0.750000 -JL -W1 -P -G255/0/0   -O -K << END >> mec.ps
133.426700 35.181000 5.310000 216 69 -170 4.0 133.426700 35.181000
END
psmeca -R -Sa0.750000 -JL -W1 -P -G255/0/0   -O -K << END >> mec.ps
133.292800 35.294800 6.900000 255 40 -167 4.0 133.292800 35.294800
END
psmeca -R -Sa0.750000 -JL -W1 -P -G255/0/0   -O -K << END >> mec.ps
133.363700 35.287500 6.080000 57 54 153 4.4 133.363700 35.287500
END
psmeca -R -Sa0.750000 -JL -W1 -P -G255/0/0   -O -K << END >> mec.ps
133.327500 35.307700 8.220000 76 76 173 4.0 133.327500 35.307700
END
psmeca -R -Sa0.750000 -JL -W1 -P -G255/0/0   -O -K << END >> mec.ps
133.298800 35.356700 8.410000 73 79 -175 4.3 133.298800 35.356700
END
psmeca -R -Sa0.750000 -JL -W1 -P -G255/0/0   -O -K << END >> mec.ps
133.288000 35.382300 7.710000 74 72 173 4.0 133.288000 35.382300
END
psmeca -R -Sa0.750000 -JL -W1 -P -G255/0/0   -O -K << END >> mec.ps
133.380700 35.246500 8.820000 43 77 -174 4.0 133.380700 35.246500
END
psmeca -R -Sa0.750000 -JL -W1 -P -G255/0/0   -O -K << END >> mec.ps
133.313800 35.373500 8.790000 341 70 -18 4.3 133.313800 35.373500
END
psmeca -R -Sa0.750000 -JL -W1 -P -G255/0/0   -O -K << END >> mec.ps
133.324500 35.311800 7.530000 162 84 1 4.3 133.324500 35.311800
END
psmeca -R -Sa0.750000 -JL -W1 -P -G255/0/0   -O -K << END >> mec.ps
133.297500 35.336300 6.260000 78 72 170 4.0 133.297500 35.336300
END
psmeca -R -Sa0.750000 -JL -W1 -P -G255/0/0   -O -K << END >> mec.ps
133.152800 35.136200 6.800000 74 77 -174 5.5 133.152800 35.136200
END
psmeca -R -Sa0.750000 -JL -W1 -P -G255/0/0   -O -K << END >> mec.ps
133.313300 35.365700 8.300000 164 79 -3 5.0 133.313300 35.365700
END
psmeca -R -Sa0.750000 -JL -W1 -P -G255/0/0   -O -K << END >> mec.ps
133.309800 35.369200 9.860000 72 77 175 4.4 133.309800 35.369200
END
#--- explanatory notes ---
pstext -R -JL -B -P -N -O -K <<END>> mec.ps
133.500000 36.054102 15 0 5 CB Period:2000/10/01 00:00--2000/10/10 23:59 
END
psmeca -R0/21/0/30 -Sa0.750000 -JX21/30 -W1 -P -G255/0/0   -N -O -K << END >> mec.ps
17.000000 8.172172 0.000000 327 50 10 4.0 17.000000 8.172172
END
pstext -R -JX  -P -N -O -K <<END>> mec.ps
17.900000 8.172172 14 0 5 5 4 
END
psmeca -R0/21/0/30 -Sa0.750000 -JX21/30 -W1 -P -G255/0/0   -N -O -K << END >> mec.ps
17.000000 8.997172 0.000000 327 50 10 5.0 17.000000 8.997172
END
pstext -R -JX  -P -N -O -K <<END>> mec.ps
17.900000 8.997172 14 0 5 5 5 
END
psmeca -R0/21/0/30 -Sa0.750000 -JX21/30 -W1 -P -G255/0/0   -N -O -K << END >> mec.ps
17.000000 9.972172 0.000000 327 50 10 6.0 17.000000 9.972172
END
pstext -R -JX  -P -N -O -K <<END>> mec.ps
17.900000 9.972172 14 0 5 5 6 
END
psmeca -R0/21/0/30 -Sa0.750000 -JX21/30 -W1 -P -G255/0/0   -N -O -K << END >> mec.ps
17.000000 11.097172 0.000000 327 50 10 7.0 17.000000 11.097172
END
pstext -R -JX  -P -N -O -K <<END>> mec.ps
17.900000 11.097172 14 0 5 5 7 
END
psmeca -R0/21/0/30 -Sa0.750000 -JX21/30 -W1 -P -G255/0/0   -N -O -K << END >> mec.ps
17.000000 12.372172 0.000000 327 50 10 8.0 17.000000 12.372172
END
pstext -R -JX  -P -N -O -K <<END>> mec.ps
17.900000 12.372172 14 0 5 5 8 
END
pstext -R -JX  -P -N -O -K <<END>> mec.ps
17.000000 13.422172 16 0 5 CM M 
END
psscale -Ctmp.cpt -D17.000000/3.738869/5.877738/0.200000 -O -K >> mec.ps
pstext -R0/21/0/30 -JX21/30 -P -N -O <<END>> mec.ps
17.500000 0.000000 14 0 5 CT Depth(km) 
END
# remove cptfile
rm tmp.cpt