#!/bin/sh
# northbrook ---
# Author: Dan Jacobson https://www.jidanni.org/
# Copyright: https://www.gnu.org/licenses/gpl.html
# Created: 2023-09-03T10:46:36+0000
# Last-Updated: 2023-09-03T10:48:48+0000
#     Update #: 3
#
set -eu
pflknwsignal='42.1528926 -87.8482539'
pflknesignal='42.1528918 -87.8481149'
pfwlsignal='42.1057223 -87.8490362'
#echo $pfwlsignal $pflknesignal | geod -f %.6f +ellps=GRS80 -I +units=us-mi
az=0.76 #half and half
pfwa='42.1274655 -87.8486412' # Pfingsten & Walters, sort of in the middle of Northbrook
oo=42.15239\ -87.77995 k=800
#Too confusing to do at once. Must use steps, and used affine to swap axes....
set -x $pfwlsignal $oo $pfwa $pflknesignal
gmt begin zz
#gmt basemap -R-87/42/-86/43+r -JOa-87.8486412/42.1274655/0.76/15 -B -V
#gmt basemap -R$2/$1/$4/$3+r -JOa$6/$5/$az/15 -Bafg -V
#gmt basemap -R$2/$1/$4/$3+r -JOb$6/$5/$8/$7/11+v -Bafg -V
#gmt basemap -RUS.IL -JOb$6/$5/$8/$7/11 -Bafg -V
gmt basemap -RUS.IL -JOb-90/36/-90/43/11+v -Bafg -V
#gmt basemap -R$2/$4/$1/$3+r -JOb$6/$5/$8/$7/11+v -B -V
gmt end show
exit
echo $@
exit
gmt begin zzz
    gmt coast -R270/20/305/25+r -JOc280/25.5/22/69/12c -Bag -Di -A250 -Gburlywood -Wthinnest -TdjTR+f2+l -Sazure
gmt end show
exit

set $pfwa
set -- proj +proj=pipeline +step +proj=omerc \
    +lonc=$2 +lat_0=$1 +alpha=$az +gamma=0 +units=us-mi \
    +step +proj=affine +s11=800 +s22=800 \
    +step -f %.0f +proj=affine +xoff=-2800 +yoff=-1400 \
    +step +proj=affine +s11=-1 +s22=-1
#800 numbers to a mile
#P&W is at -2800, -1400
#only the SW quadrant is used, so we make the numbers positive

P=$@
echo $P
    $P -r <<EOF
$pfwa pfwa
$pflknwsignal SIGw
$pflknesignal SIGe
$pfwlsignal pfwllo
$oo OO
42.15289 -87.88740 0 4400w 0
EOF