#!/bin/csh
#
# Creates a Western US map
# This version adds the earthquakes
#
set file = ps.wus_faults
#
set cpath = /usr/contrib/data/maps
set data = /data/a1/lind/graphics/gmt/data
set black = 0/0/0
set white = 255/255/255
set grey = 130
set shade = 230
set red = 255/0/0
set green = 0/255/0
set blue = 0/0/255
set tour = 0/255/255
set yellow = 255/255/0
set pblue = 220/240/255
set lblue = 180/200/255
set pyellow = 255/255/160
#
set water = 80/80/255
set land = 215/200/150
set state = 1/155/150/0to
#
pscoast -R-128/-105/30/50 -JM5.5 -Bf1a5/f1a5::WESN -G$land -C$water -S$water -N1 -N2/$state -W1 -A150 -P -X1.5 -Y2.3 -K > $file
#
# add the faults
#
psxy $cpath/cafaults.asc -JM -R -Ms -W5/$red -: -O -K  >> $file
psxy $data/wus_faults.asc -JM -R -M$ -W5/$red -O -K >> $file
#
# or add the plate boundaries
#
psxy $data/nuvel1_plates.asc -JM -R -MS -W7/$black -: -O >> $file
#
# STOP HERE
#