#!/bin/sh
#
#	drft06rr_poster - Create the final regional summary poster
#	showing all dredges as well as published data
#
# Paul Wessel, 11-DEC-2001

DIR=drft06rr
prefix=drft06rr
R=243/290/-30/-10
width=32.5
gmtset PLOT_DEGREE_FORMAT ddd:mm:ssF GRIDFILE_SHORTHAND true PAPER_MEDIA archE+ HEADER_FONT 17 ANOT_FONT_SIZE +16p HEADER_OFFSET 0.35i
#--------------------------------------------------
if [ ! -f data/$prefix.grd ]; then
	grdraster 5 -R$R -Gdata/$prefix.grd -V
	grdgradient data/$prefix.grd -Nt1.2 -A90 -Gdata/$prefix.int.grd -V
fi
grep '^Dredge' waypoints.txt | grep '\-' | egrep -v ':' | tr '-' ' ' | awk '{printf "%s:%sW\t%s:%sS\t%s\t%s\n", $6, $7, $4, $5, $2, $3}' > dredges.d
# Need map height and where to center the color bar:
height=`echo 243 -10 | mapproject -R$R -JM$width | cut -f2`
Y=`gmtmath -Q $height 2.5 ADD =`
x0=`gmtmath -Q $width 0.5 MUL =`
makecpt -Crainbow -T-6000/0/500 -Z > $prefix.cpt
psxy -R$R -JM$width -K /dev/null -Y1 -X0.75 -P > $prefix.poster.ps
if [ ! -f $prefix.posterimage.ps ]; then	# We do this to save time
	grdmath data/$prefix.int.grd 0.25 ADD = t.grd
	grdimage -R -JM data/$prefix.grd -It.grd -C$prefix.cpt -Ei -O -K -V > $prefix.posterimage.ps
fi
cat $prefix.posterimage.ps >> $prefix.poster.ps
pscoast -R -JM -B5g5f1:."DRFT06RR\072 Drift Expedition Leg 6 - Easter/Salas y Gomez/Nazca Survey, Nov-Dec 2001": -G100 -O -K -W0.25p >> $prefix.poster.ps
psxy -R -JM $GMT_DATADIR/plate_boundaries_1996.d -M -O -K -W1p >> $prefix.poster.ps
psxy -R -JM 200nm_line.d -O -K -W2p/255/130/0 >> $prefix.poster.ps
psxy -R -JM -O -K -Ss0.075 -W0.25p -G150 -H Bonatti_rocks.d >> $prefix.poster.ps
psxy -R -JM -O -K -Ss0.075 -W0.25p -G255/255/50 -H misc_rocks.d >> $prefix.poster.ps
psxy -R -JM -O -K -Ss0.075 -W0.25p -G200/50/200 -H Hydros_rocks.d >> $prefix.poster.ps
psxy -R -JM -O -K -Ss0.075 -W0.25p -G0/0/150 -H Clark_and_Dymond_rocks.d >> $prefix.poster.ps
psxy -R -JM -O -K -Ss0.075 -W0.25p -G255 -H glor07mv_no_rocks.d >> $prefix.poster.ps
psxy -R -JM -O -K -Ss0.075 -W0.25p -G0/200/200 -H glor07mv_rocks.d >> $prefix.poster.ps
psxy -R -JM -O -K -Ss0.075 -W0.25p -G0/200/0 -H Oconnor_rocks.d >> $prefix.poster.ps
psxy Bonatti_dates.d -H -R -JM -O -K -St0.2i -G100/255/255 -W0.25p >> $prefix.poster.ps
psxy Oconnor_dates.d -H -R -JM -O -K -St0.2i -G255/0/255 -W0.25p >> $prefix.poster.ps
psxy glor07mv_dates.d -H -R -JM -O -K -St0.2i -G255 -W0.25p >> $prefix.poster.ps
awk '{printf "%s %s 6 0 0 LM %s Ma\n", $1, $2, $3}' Bonatti_dates.d | pstext -H -R -JM -O -K -D0.15/0 >> $prefix.poster.ps
awk '{printf "%s %s 6 0 0 LM %s Ma\n", $1, $2, $3}' Oconnor_dates.d | pstext -H -R -JM -O -K -D0.15/0 >> $prefix.poster.ps
awk '{printf "%s %s 6 0 0 LM %s Ma\n", $1, $2, $3}' glor07mv_dates.d | pstext -H -R -JM -O -K -D0.15/0 >> $prefix.poster.ps
psxy -R -JM -O -K -W0.75p waypts.d >> $prefix.poster.ps
psxy -R -JM -O -K -W0/255 waypts.d >> $prefix.poster.ps
psxy -R -JM -O -K -Sc0.04i -G0 waypts.d >> $prefix.poster.ps
#psxy -R -JM -O -K -Sc0.175i -W0.25p/255/0/0 waypts.d >> $prefix.poster.ps
#awk '{print $1, $2, 9, 0, 0, "CM", $3}' waypts.d | pstext -R -JM -O -K >> $prefix.poster.ps
echo "282.95	-12.10" > lima.d
psxy -R -JM -O -K -Sa0.25 -G0/255/0 -W0.5p lima.d >> $prefix.poster.ps
awk '{if ($NF != "D") print $1, $2}' dredges.d | psxy -R -JM -O -K -G220/180/0 -Sc0.175 -W0.25p >> $prefix.poster.ps
awk '{if ($NF == "D") print $1, $2}' dredges.d | psxy -R -JM -O -K -G255/0/0 -Sc0.175 -W0.25p >> $prefix.poster.ps
# Plot Easter and Salas y Gomez locations
echo "109:22.50W 27:07.02S" | psxy -R -JM -O -K -Sa0.25 -G255/0/0 -W0.25p >> $prefix.poster.ps
echo "105:21.96W 26:28.44S" | psxy -R -JM -O -K -Sa0.25 -G255/255/0 -W0.25p >> $prefix.poster.ps
xi=`gmtmath -Q $width 0.1 SUB =`
yi=`gmtmath -Q $height 0.1 SUB =`
psimage data/moai.ras -E600 -C$xi/$yi/TR -F1p -O -K >> $prefix.poster.ps
psscale -C$prefix.cpt -D$x0/-0.4/$width/0.1h -O -K -B/:m: >> $prefix.poster.ps
pslegend -R -JM -O -K -D243.25/-10.25/3.5/5.6/TL -F -G255 << EOF >> $prefix.poster.ps
G -0.05i
H 24 5 MAP LEGEND
D 0.25 1p
S 0.2 c 0.1i 0/0/0 0.25p 0.5 Waypoints
S 0.2 c 0.2i 255/0/0 0.25p 0.5 Datable dredge locations
S 0.2 c 0.2i 220/180/0 0.25p 0.5 Other dredge locations
D 0.25 1p
S 0.2 t 0.2i 255/0/255 0.25p 0.5 Dates from O\`Conner
S 0.2 t 0.2i 255/255/255 0.25p 0.5 Dates from Gloria07mv
S 0.2 t 0.2i 100/255/255 0.25p 0.5 Dates from Bonatti
D 0.25 1p
S 0.2 s 0.20i 0/0/150 0.25p 0.5 Clark and Dymond samples
S 0.2 s 0.20i 150/150/150 0.25p 0.5 Bonatti samples
S 0.2 s 0.20i 255/255/255 0.25p 0.5 Gloria07mv no-rock samples
S 0.2 s 0.20i 0/200/200 0.25p 0.5 Gloria07mv rock samples
S 0.2 s 0.20i 0/200/0 0.25p 0.5 OConner samples
S 0.2 s 0.20i 200/50/200 0.25p 0.5 Hydros samples
S 0.2 s 0.20i 255/255/50 0.25p 0.5 Misc samples
D 0.25 1p
S 0.2 a 0.2i 255/255/0 0.25p 0.5 Salas y Gomez Island
S 0.2 a 0.2i 255/0/0 0.25p 0.5 Easter Island
S 0.2 a 0.2i 0/255/0 0.25p 0.5 Callao, Peru
S 0.2 - 0.2i - 2p/255/100/0 0.5 Peru 200 nm boundary
EOF
psxy -R -JM -O /dev/null >> $prefix.poster.ps
#gv $prefix.poster.ps &
rm $prefix.cpt lima.d
ps2rtl -h $Y -w 34.25 -d 24 -r 300 $prefix.poster.ps poster.rtl
mv $prefix.poster.ps $prefix