Thanks! So this is some sped up proximity check? I’ve also noticed strange behavior for sphinterpolate where numerous datapoints ar flagged as “identical” when they are numerically distinct.
The distance in radias
echo 43.782 41.121 202.5 47.25 | gawk -f ~/awk/dist_on_sphere.awk
1.56435210975521666e+00
is quite different from zero, but I presume it’s too close to pi/2 for something to be messed up?
Where is that check? I would assume the gmtmap_geodesic_dist_degree works fine.