#!/bin/bash

# Make an equidistant azimuthal projection

# Tallinn, Estonia:

  LAT=59.434
  LON=24.728

# Base map:

  pscoast \
    -K \
    -R0/360/0/90 \
    -JE$LON/$LAT/7i \
    -B15g15 \
    -Dl \
    -A1000 \
    -G200 \
    -W0.25p \
    -P > equidistant.eps

# Rings:

  grdmath \
    -Rg \
    -I1 \
    $LON \
    $LAT \
    SDIST \
    111.13 MUL \
    = dist.grd

  grdcontour \
    dist.grd \
    -A1000+v+ukm+kblack \
    -Glz-/z+ \
    -S8 \
    -C500 \
    -O \
    -J \
    -Wathin,red \
    -Wcthinnest,red,- >> equidistant.eps

  ps2pdf equidistant.eps
  rm equidistant.eps
  rm dist.grd