[Proj] Polar stereographic, different values on different platforms?

[Gerald I. Evenden]

On Wednesday 20 August 2008 3:21:01 pm Bruce Raup wrote:
> Frank and all,
>
> Has the problem described in the below been resolved yet?  A colleague
> of mine and I have just executed the same "proj" command on Windows
> and Linux, and are getting different results.  Windows version is Rel.
> 4.6.1, 21 Jul 2008.  Linux versions are Rel. 4.5.0, 22 Oct 2006 and
> Rel. 4.6.0, 21 Dec 2007
>
> Command:
> echo "105 40" | proj +proj=stere +lat_0=90 +lon_0=0 +lat_ts=70 +datum=WGS84
>
> On both Linux boxes, we get the output
> 5577808.93      1494569.40

For what it is worth, lproj reproduces the above value as well as the test in 
the earlier email.

> On Windows, my colleague gets:
> 5610189.12      1503245.64
>
> Using GMT's "mapproject" program, I get the same results as the Linux
> versions of Proj:
> echo "105 40" | mapproject -Js0/90/70/1:1 -C -F -R-180/180/30/90
> 5577808.933369  1494569.399231

Also, for what it is worth, lproj gives:
5577808.933368  1494569.399231

Must be a problem in the longitude computation.  ;-)

	...
-- 
The whole religious complexion of the modern world is due
to the absence from Jerusalem of a lunatic asylum.
-- Havelock Ellis (1859-1939) British psychologist