;
;  fsnspp: north stereographic polar projection
;
FUNCTION fsyspp, plam, pphi
;
rpi = double( 2.*asin( 1. ))
rad = double( rpi/180.)
;
z = double( 2. * sin( rad*plam ) * tan( rpi/4. - rad*pphi/2. ))
;
return, z
END