#!/bin/csh -f

set J=M4.0
set R=-131/-113.5/39/52
set B=a2/a1
set PS=hypo_section.ps
set ddirseis=../mapdata/events/hist_seis

# plot map showing the location of the fault
psbasemap -J$J -R$R -B$B -X2.0 -Y6 -G200 -K -P -V >! $PS 
pscoast -J$J -R$R -S255 -Na/2ta -Dl -W1 -P -V -O -K  >> $PS 

# Physiographic province boundaries from Lowry, Ribe, and Smith, 2000
set pfile=../mapdata/Yellowstone-SRP/phyprv.gmt.new
psxy $pfile -: -R$R -J$J -W3/0 -M -O -V -N -K >> $PS

#Plot cross-section 
psxy -V -J$J -R$R -W6/0/0/255 -: -K -O << FAULT  >> $PS 
39.5 -126
39.5 -113.5
FAULT

#Plot stations
#awk '{print $2,$1}' jh_oneweek.txt|psxy  -R -P -W0.01 -Jm -Sc0.01 -V -O -K >> $PS
#awk '{print $1,$2,$3,$4,$5,$6,$8}' stlist.txt|pstext  -R -P -Jm -V -O -K >> $PS
#awk '{print $1,$2}' stlist.txt|psxy -R -P -W0.1 -G0/0/250  -St0.1 -Jm -V -O -K>> $PS

# projection for first cross section
#project -Q -C-120.3500/35.8200 -A137 xyz.dat -Fxyzpqrs -L-600/600 -W-50/50  > xyzpqrs
project -Q -C-119.75/39.5 -E-113.5/39.5 $ddirseis/master_ANSS.dat -Fxyzpqrs -L-536.42/536.42 -W-25/25  > xyzpqrs
awk '{ print $1,$2,$4/40,$3,$4}' xyzpqrs >! plot2.xy
#project -Q -C-120.82/43.47 -E-116.14/41 ../mapdata/events/moment_master.dat -Fxyzpqrs -L-467/467 -W-25/25  > xyzpqrs
#awk '{ print $1,$2,$3,$4,$5,$6,$7,$8,$9}' xyzpqrs >! plot2.xy

#project -Q -C-119.75/39.5 -E-113.5/39.5 ../mapdata/events/momentBSL.dat -Fxyzpqrs -L-536.42/536.42 -W-25/25  > xyzpqrs1
#awk '{ print $1,$2,$3,$4,$5,$6,$7,$8,$9}' xyzpqrs1 >! plot2_1.xy
#
# plot epicenters  selected for first cross section (overwrite)
psxy -K -O -J$J -R$R -Sc -W0.5p/0/0/0 -G200/0/0 plot2.xy >> $PS
#psmeca  plot2.xy -R$R -J$J -G255/0/0 -Sa0.5 -C -W2 -V -K -O -H0 >> $PS
#psmeca  plot2_1.xy -R$R -J$J -G0/0/255 -Sa0.5 -C -W2 -V -K -O -H0 >> $PS
#

# draw cross section 1
psbasemap -K -O -X-1. -Y-3.5 -JX6/-2 -R-536.42/536.42/0/150 -B50/50 >> $PS
awk '{ print $5,$3,$4/40}' xyzpqrs > plot4.xy
#awk '{print $10,$3,$3,$4,$5,$6,$7,$8,$9}' xyzpqrs > plot4.xy
#awk '{print $10,$3,$3,$4,$5,$6,$7,$8,$9}' xyzpqrs1 > plot4_1.xy

psxy <<END -JX -R -Sc.025 -G155/155/155 -W0.5p/0/0/0 -O -K -V >> $PS
-450 90
END
psxy <<END -JX -R -Sc.05 -G155/155/155 -W0.5p/0/0/0  -O -K -V >> $PS
-400 90 
END
psxy <<END -JX -R -Sc.075 -G155/155/155 -W0.5p/0/0/0 -O -K -V >> $PS
-350 90
END
psxy <<END -JX -R -Sc.1 -G155/155/155 -W0.5p/0/0/0  -O -K -V >> $PS
-300 90
END
psxy <<END -JX -R -Sc.125 -G155/155/155 -W0.5p/0/0/0 -O -K -V >> $PS
-250 90
END
psxy <<END -JX -R -Sc.15 -G155/155/155 -W0.5p/0/0/0  -O -K -V >> $PS
-200 90
END
psxy <<END -JX -R -Sc.175 -G155/155/155 -W0.5p/0/0/0 -O -K -V >> $PS
-150 90
END
psxy <<END -JX -R -Sc.2 -G155/155/155 -W0.5p/0/0/0  -O -K -V >> $PS
-100 90
END

psxy -O -JX -R -Sc -W0.5p/0/0/0 -G200/0/0 plot4.xy >> $PS
#psmeca  plot4.xy -R -JX -G255/0/0 -Sa0.5 -C -W2 -V -O -K -H0 >> $PS

#project -Q -C-119.75/39.5 -E-113.5/39.5 ../mapdata/events/momentBSL.dat -Fxyzpqrs -L-536.42/536.42 -W-25/25  >! xyzpqrs
#awk '{ print $1,$2,$3,$4,$5,$6,$7,$8,$9}' xyzpqrs1 >! plot2.xy
#psmeca  plot2.xy -R$R -J$J -G0/0/255 -Sa0.5 -C -W2 -V -K -O -H0 >> $PS
#awk '{print $10,$3,$3,$4,$5,$6,$7,$8,$9}' xyzpqrs >! plot4.xy
#psmeca  plot4.xy -R -JX -G0/0/255 -Sa0.5 -C -W2 -V -O -H0 >> $PS

echo "done"
ggv  $PS&