Inverse Map Projection Help

Hi Group,

I've read the previous posts on inverse map projections and the
lengthy tutorial by David Fanning, but still cannot get things to line
up quite right. So I'm posting here for help...

I have a data set (BEDMAP) with this information in the header:

ncols         1371
nrows         1371
xllcorner     -3426225.75
yllcorner     -3426225.75
cellsize      5000
NODATA_value  -9999

And this information on the website:

Polar Stereographic projection with 71=B0S as the latitude of true scale
and 0=B0E as the central meridian.

I've managed to load the data, and inverse project it approximately
such that things roughly line up. But I cannot get it accurate where
my reference for 'accurate' is the /MAP_CONTINENTS, /HIRES keywords.

My goal is to inverse project it to an equilateral lat/lon grid
(something like what comes from /CYLINDRICAL) so I can generate an
image that would work in Google Earth.

I'm pleased to post the code I've written if necessary. If anyone has
any suggestions or familiarity with this data set I would appreciate
any tips.