#!/bin/csh

#  Draw basemap
#  ------------
pscoast -Jm0.677705740167622 -R-141.41948/-133.45142/57.01818/62.15012 \
-Ba1.0f10mwSnE \
-W3/0/0/0 -Na/1 -S153/255/255 \
-Df -Ir/1/114/187/243 \
-Lf-139.5888/57.927933/57.927933/50 \
-G230/225/200 -P -K >! Yakutat.ps

#  Station labels
#  --------------
pstext << EOT -: -R -Jm -G0 -O -K \
>> Yakutat.ps
59.6698 -138.550766 9 0.0 1 5 COMB
59.6693 -138.774066 9 0.0 1 5 FT83
59.5107 -139.560266 9 0.0 1 5 GEOS
58.4178 -135.608466 9 0.0 1 5 GUS2
59.7055 -138.856966 9 0.0 1 5 HIDD
59.7589 -138.772866 9 0.0 1 5 MIDX
59.5680 -138.902766 9 0.0 1 5 MOSE
59.5735 -138.468866 9 0.0 1 5 NOVA
59.8793 -139.046366 9 0.0 1 5 NQ1X
59.8812 -138.977166 9 0.0 1 5 NQ2X
59.8273 -138.943066 9 0.0 1 5 NQ4X
59.8236 -139.017666 9 0.0 1 5 NQ5X
59.8330 -138.953166 9 0.0 1 5 NQ6X
60.7505 -135.133566 9 0.0 1 5 WHIT
59.4963 -138.644766 9 0.0 1 5 YAKU
EOT

#  Kodiak Earthquake
#  ------------------
psxy << EOT -: -R -W5/51/51/255 -G255/255/255 \
-Jm -Sa.15 -O -K \
>> Yakutat.ps
  56.6833     -153.3000
EOT

#  Kodiak Earthquake label
#  ------------------------
pstext << EOT -: -R -Jm -G0 -O -K \
>> Yakutat.ps
57.0 -153.2 9 0.0 1 5 Kodiak_Eq_20010110
56.8 -153.2 9 0.0 1 5 M6.7
EOT

#  Kodiak Earthquake
#  ------------------
psxy << EOT -: -R -W5/51/51/255 -G255/255/255 \
-Jm -Sa.15 -O -K \
>> Yakutat.ps
  57.51     -154.67
EOT

#  Kodiak Earthquake label
#  ------------------------
pstext << EOT -: -R -Jm -G0 -O -K \
>> Yakutat.ps
57.4 -154.3 9 0.0 1 5 Kodiak_Eq_19991206
57.2 -154.3 9 0.0 1 5 M7.0
EOT

#  Velocity vectors
#  ----------------
psvelomeca << EOT -R \
-W1/255/0/63 -G70/100/190 \
-Se0.0070576610911144/.95/0 -B  \
-Jm -A0.01/0.12/0.03 -O -K \
>> Yakutat.ps
-138.6393	59.6698	-11.82	15.69	1.07 	1.06 	-0.02	COMB
-138.8626	59.6693	-52.75	20.16	196.82	197.78	-0.00	FT83
-139.6488	59.5107	-24.77	39.35	0.98 	0.94 	-0.02	GEOS
-135.6970	58.4178	121.03	88.89	143.57	153.83	0.04 	GUS2
-138.9455	59.7055	-20.16	26.99	1.05 	1.02 	-0.03	HIDD
-138.8614	59.7589	19.38	-14.08	198.81	199.18	0.00 	MIDX
-138.9913	59.5680	0.00 	0.00 	200.00	200.00	0.00 	MOSE
-138.5574	59.5735	0.00 	0.00 	200.00	200.00	0.00 	NOVA
-139.1349	59.8793	-7.78	1.50 	199.09	199.54	-0.00	NQ1X
-139.0657	59.8812	0.00 	0.00 	200.00	200.00	0.00 	NQ2X
-139.0316	59.8273	-33.02	-27.58	197.10	197.84	0.00 	NQ4X
-139.1062	59.8236	0.00 	0.00 	200.00	200.00	0.00 	NQ5X
-139.0417	59.8330	-6.55	24.06	198.46	198.58	0.00 	NQ6X
-135.2221	60.7505	28.48	39.84	134.87	144.38	0.03 	WHIT
-138.7333	59.4963	-26.78	-141.69	172.83	180.69	0.01 	YAKU
EOT

#  Station dots
#  ------------
psxy << EOT -: -R -W5/51/51/255 -G255/255/255 \
-Jm -Sc0.07 -O -K \
>> Yakutat.ps
59.6698	-138.6393
59.6693	-138.8626
59.5107	-139.6488
58.4178	-135.6970
59.7055	-138.9455
59.7589	-138.8614
59.5680	-138.9913
59.5735	-138.5574
59.8793	-139.1349
59.8812	-139.0657
59.8273	-139.0316
59.8236	-139.1062
59.8330	-139.0417
60.7505	-135.2221
59.4963	-138.7333
EOT

#  Olympia Earthquake
#  ------------------
psxy << EOT -: -R -W5/51/51/255 -G255/255/255 \
-Jm -Sa.15 -O -K \
>> Yakutat.ps
  47.167 -122.733
EOT

#  HM Earthquake
#  -------------
psxy << EOT -: -R -W5/51/51/255 -G255/255/255 \
-Jm -Sa.15 -O -K \
>> Yakutat.ps
  34.595500     -116.268167
EOT

#  HM Earthquake label
#  --   --------------
pstext << EOT -: -R -Jm -G0 -O -K \
>> Yakutat.ps
34.55 -116.268167 9 0.0 1 5 Hector_Mine_Eq
34.53 -116.268167 9 0.0 1 5 M7.1
EOT

#  Landers Earthquake
#  ------------------
psxy << EOT -: -R -W5/51/51/255 -G255/255/255 \
-Jm -Sa.15 -O -K \
>> Yakutat.ps
  34.216666667     -116.43333333333
EOT

#  Landers Earthquake label
#  ------------------------
pstext << EOT -: -R -Jm -G0 -O -K \
>> Yakutat.ps
34.18 -116.43333333333 9 0.0 1 5 Landers_Eq
34.16 -116.43333333333 9 0.0 1 5 M7.3
EOT

#  Label plot and vector scale
#  ---------------------------
psvelomeca << EOT -: -R \
-W1/255/0/0 -G70/100/190 -Se0.0070576610911144/.95/0 -B  \
-Jm -A0.03/0.12/0.03 -O -K \
>> Yakutat.ps
57.578028 -139.5888 0.0 20.0 0.0 0.0 0.0 Scale
EOT
pstext << EOT -: -R -Jm -O \
>> Yakutat.ps
57.4030755 -139.5888 9 0.0 0 5 20 mm/yr 
EOT