#!/bin/sh

	#GMT
	gmt=$(which GMT)

	#Definition
	output=delta
	range="-R5/10/45/50"
	proj="-JM6i"
	scale="-Lf175/-48/-47/200n"

	#Cut Topological Data
	$gmt grdcut ETOPO1_Ice_g_gmt4.grd $range -G$output.grd

	#Create Continent Shadows
	$gmt grdgradient $output.grd -A180 -Nto.3 -G${output}_shadows.grd

	#Create Ocean and Land Maps
	$gmt grdimage $output.grd -Cocean.cpt $range $proj -P -K > $output.ps
	
	$gmt pscoast -Gc -Di -R -J -K -O >> $output.ps
	$gmt grdimage -Cland.cpt -R -J -K -O $output.grd -I${output}_shadows.grd  >> $output.ps
	
	$gmt pscoast -Q -K -O >> $output.ps

	#Plot
	$gmt psbasemap -B10g5 -R -J -O $scale >> $output.ps