#!/bin/sh
# $Id: distance_contour_map,v 1.1 2002-04-25 13:33:34-06 braup Exp braup $
# Makes a contour map of distances from a specific point on the earth.
#proj=Q0/20c
proj=G-40/40/15c
region=-180/180/-90/90
inc=1
grdfile=dist.grd
out=dist.ps

pngout=boulder_distance_ortho.png
#pngout=anchorage_distance.png
#pngout=south_pole.png
#pngout=equator.png

origin="-105.3 40"     # Boulder
#origin="-150 61.16"    # Anchorage
#origin="0 -90"          # South Pole
#origin="-75.867 0.12"    # somewhere near equator

place_name='Boulder, Colorado'
#place_name='Anchorage, Alaska'
#place_name='South Pole'
#place_name='Putumayo'

echo "Origin = $origin" >&2
bndpen=0.5p/100/100/100
gmtset PAPER_MEDIA letter BASEMAP_TYPE plain FRAME_PEN 1p HEADER_FONT_SIZE 20p
gmtset D_FORMAT %lg DEGREE_FORMAT 1
grdmath -R$region -I$inc -F $origin GDIST PI MUL 180 DIV 6378 MUL = $grdfile
#pscoast -R$region -J$proj -A3000 -Dc -G200 -B30:".Contours of distance from Boulder, Colorado (km)": -U -K -V > $out
pscoast -R$region -J$proj -A3000 -Dc -N1/$bndpen -N2/$bndpen -G200 -W$bndpen -B30 -U"B. Raup" -K -V > $out
grdcontour $grdfile -J$proj -C500 -A1000 -W0.5p/200/150/100 -K -O -V >> $out
#echo "9.8 5.6 13 0 15 CM Contours of distance from $place_name (km)" |
#   pstext -R0/20/0/10 -JX20c -N -W -O -K -V >> $out
echo "$origin" | psxy -J$proj -R$region -Sc0.2c -G0/0/255 -O -V >> $out

convert -density 80 -rotate 90 -trim $out $pngout