#!/bin/bash
file=${1-tmp}
#
# convert data in $file.dat to gridded representation using GMT sphinterpolate
#
reg=-R-45/45/-45/45		# region
proj=-JH0/7			# projection for plotting
inc=-I0.25			# resolution of grid file


# grid using triangulation
sphinterpolate -Q2 -T $file.data $reg $inc -G$file.sphi.grd > /dev/null

# reuse colormap from _surface!

s=sphi

grdimage -Y0 -K $file.$s.grd -C$file.cpt $reg $proj \
    -Ba45f5WeSn:".$s":  -V  > $file.$s.ps
psxy $file.data -Sc0.02 -G0 -O $reg $proj  >> $file.$s.ps


