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

ps=sph_2.ps

gmt makecpt -Crainbow -T-9000/9000/1000 -Z > tt.cpt
gmt sphinterpolate lun2.txt -Rg -I1 -Q0 -Gtt.nc
gmt grdimage tt.nc -JH0/4.5i -B30g30 -B+t"-Q0" -Ctt.cpt -X0.8i -Y5.5i -K --MAP_TITLE_OFFSET=0i --FONT_TITLE=18p > $ps
gmt sphinterpolate lun2.txt -Rg -I1 -Q1 -Gtt.nc
gmt grdimage tt.nc -J -B30g30 -B+t"-Q1" -Ctt.cpt -X4.9i -O -K  --MAP_TITLE_OFFSET=0i --FONT_TITLE=18p >> $ps
gmt sphinterpolate lun2.txt -Rg -I1 -Q2 -Gtt.nc
gmt grdimage tt.nc -J -B30g30 -B+t"-Q2" -Ctt.cpt -X-4.9i -Y-5i -O -K  --MAP_TITLE_OFFSET=0i --FONT_TITLE=18p >> $ps
gmt sphinterpolate lun2.txt -Rg -I1 -Q3 -Gtt.nc
gmt grdimage tt.nc -J -B30g30 -B+t"-Q3" -Ctt.cpt -X4.9i -O -K  --MAP_TITLE_OFFSET=0i --FONT_TITLE=18p >> $ps
gmt psxy -Rg -J -O -K lun2.txt -Sc0.02i -G0 -B30g30 -X-2.45i -Y2.5i >> $ps
gmt psxy -Rg -J -O -T >> $ps