#!/bin/bash
#
#	$Id: rendering.sh 12114 2013-09-03 19:19:00Z fwobbe $

ps=rendering.ps

gmt xyz2grd -R0/360/0/90 -r -I60/30 -Gt.nc <<%
#30 45 -1
90 15 -1
30 75 0
90 45 0
150 15 0
90 75 1
150 45 1
210 15 1
150 75 0
210 45 0
270 15 0
210 75 -1
270 45 -1
330 15 -1
270 75 0
330 45 0
#
30 15 0
330 75 -1
%

render () {
gmt grdimage -R -JW0/6i -Bx60 -By30 t.nc -Ct.cpt -E100 -O -K $*
gmt pstext -R -J -F+f16p,Helvetica-Bold+jBL -O -K -N <<< "180 35 $1"
}

gmt makecpt -Cpolar -T-1/1/0.5 -Z -D > t.cpt

gmt psscale -D3i/-0.4i/6i/0.4ih -Ct.cpt -P -K -Y2i > $ps
render -nn >> $ps
render -nl -Y2i >> $ps
render -nb -Y2i >> $ps
render -nc -Y2i >> $ps
gmt psxy -R -J -T -O >> $ps