#!/bin/bash
#
#	$Id: imgmap.sh 12713 2014-01-06 17:15:37Z pwessel $

ps=imgmap.ps
IMG=topo.15.1.img
if [ -f /tmp/gmt/$IMG ]; then # use one placed/linked from /tmp
	# PW: Run (cd /tmp/gmt; ln -s /Volumes/MacNutHD3/UH/RESOURCES/DATA/img/topo.15.1.img topo.15.1.img) before make check
	IMG=/tmp/gmt/topo.15.1.img
else
	OK=`gmt gmtwhich -C $IMG`
	if [ $OK = N ]; then
		echo "File $IMG not available for testing imgmap.sh"
		exit
	fi
fi
# Get merc grid
gmt img2grd $IMG -R180/200/-5/5 -I1m -T1 -D -S1 -Gimg.nc -M
gmt makecpt -Crainbow -T-8000/0/500 -Z > t.cpt
gmt grdimage img.nc -Jx0.25i -Ct.cpt -P -K -Xc > $ps
gmt psbasemap -R -Jm0.25i -Ba -BWSne -O -K >> $ps
# Get geo grid
gmt img2grd $IMG -R -I1m -T1 -D -S1 -Gimg.nc
gmt grdimage img.nc -Jm -Ct.cpt -O -K -Ba -BWSne -Y3.25i >> $ps
# Get resampled geo grid
gmt img2grd $IMG -R -I1m -T1 -D -S1 -Gimg.nc -E
gmt grdimage img.nc -Jm -Ct.cpt -O -K -Ba -BWSne -Y3.25i >> $ps
gmt psxy -R -J -O -T >> $ps