#!/bin/bash
#
#       $Id: spotter_2.sh 12395 2013-10-28 17:12:37Z fwobbe $

ps=spotter_2.ps

# Example 2 - Using gmt hotspotter
#
# We will use gmt hotspotter to create a CVA image for the Pacific.
# It will look similar to the ones we have published but we will
# here use only seamounts with a VGG amplitude of at least 100 Eotvos.

POLES=${src}/../../src/spotter/WK97.d # Rotation poles to use

gmt hotspotter ${src}/../../src/spotter/seamounts.d -I10m -R130/260/-66/60 -E${POLES} -Gspotter_2.nc -T -N145

# Make a suitable color table

gmt makecpt -Chot -T0/3000/300 -Z > t.cpt

gmt grdimage spotter_2.nc -JM6i -P -K -Ct.cpt > $ps
gmt pscoast -R -J -O -Gdarkgreen -A500 -Dl -W0.25p -B20 -BWSne >> $ps