pro bobskymap, pim, ps = ps, min = min, max = max, log = log

colors, datacolors = datacolors

  if keyword_set(ps) then begin 
    openps, /land
  endif 
  theta = 90-pim.theta
  latmin = min(theta)
  latmax = max(theta)
  lonmin = min(pim.phi)
  lonmax = max(pim.phi)
  ;; Normalize and reverse since the first column is lat=+57 deg and 
  ;; theta=90-57deg.
  obsim = reverse(pim.im/double(max(pim.im)), 2)
  if n_elements(min) eq 0 then min = min(pim.im(where(pim.im gt 0)))
  if n_elements(max) eq 0 then max = max(pim.im)
  xmargin = 10
  ymargin = 10
  skymap, reverse(pim.im, 2),  $
    latmin = latmin, latmax = latmax, lonmin = lonmin,  $
    lonmax = lonmax, sc_pos = pim.scpos, spin_axis = pim.spinax,  $
    xmargin = xmargin, ymargin = ymargin, mag = [0, 0, 1], $
    /field, /limb, /grid, /square, /terminator, log = log, $
    top = datacolors-1, /sphere, /colorbar,  $
    bartitle = 'ARB. UNITS (LOG)',  $
    min = min, max = max, /noerase
  xyouts, .6, .02, pim.time, /normal
  if keyword_set(ps) then begin 
    closeps
  endif 

  return
end