pro plotbobs
  restore,  '/project/ena/spike/retrieve/posterdata.sav'

;Join_input - input ion distribution
;Join - retrieved test ion distribution (noise)
;Jion2 - retrieved test ion distribution (no noise)
;Join_data - ion distribution from the real data

;Pim0 - test image (with noise)
;Pim - real data image
;Simim - 2d array with noise free simulated data

  !p.thick = 2
  !p.charthick = 2

  colors, datacolors = datacolors

;  openps, /la, file = '~/agu2001_fall/bobsposter/testion.ps'
;  plotion, jion_input, /log, min = 4, max = 7, elmax = 12, elcl = 10, /col, /noer
;  xyouts, .5, .95, 'TEST INPUT ION DISTRIBUTION', CHARS = 2, /norm, align = .5
;  closeps, /qui

;  openps, /la, file = '~/agu2001_fall/bobsposter/inv_noise.ps'
;  plotion, jion, /log, min = 4, max = 7, elmax = 12, elcl = 10, /col, /noer
;  xyouts, .5, .95, 'INVERSION FROM SIMULATED ENA IMAGE!C(With Poisson noise)', $
;    CHARS = 2, /norm, align = .5
;  closeps, /qui

;  openps, /la, file = '~/agu2001_fall/bobsposter/inv_nonoise.ps'
;  plotion, jion2, /log, min = 4, max = 7, elmax = 12, elcl = 10, /col, /noer
;  xyouts, .5, .95, 'INVERSION FROM SIMULATED ENA IMAGE!C(No Poisson noise)', $
;    CHARS = 2, /norm, align = .5
;  closeps, /qui

;  openps, /la, file = '~/agu2001_fall/bobsposter/inv_real.ps'
;  plotion, jion_data, /log, min = 4, max = 7, elmax = 12, elcl = 10, /col, /noer
;  xyouts, .5, .95, 'INVERSION FROM HENA IMAGE', CHARS = 2, /norm, align = .5
;  closeps, /qui

  openps, /la, file = '~/agu2001_fall/bobsposter/sim_ena_noise.ps'
  bobskymap, pim0, /log, min = 1., max = 1e3
  xyouts, .5, .95, 'SIMULATED TEST ENA IMAGE!C(With Poisson noise)', $
    CHARS = 2, /norm, align = .5
  closeps, /qui
  
  openps, /la, file = '~/agu2001_fall/bobsposter/real_ena.ps'	
  bobskymap, pim, /log, min = 1., max = 1e3
  xyouts, .5, .95, 'OBSERVED HENA IMAGE', CHARS = 2, /norm, align = .5
  closeps, /qui

  openps, /la, file = '~/agu2001_fall/bobsposter/sim_test.ps'
  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.
  min = 1.
  max = 1.e3
  xmargin = 10
  ymargin = 10
  skymap, reverse(simim, 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, $
    top = datacolors-1, /sphere, /colorbar,  $
    bartitle = 'ARB. UNITS (LOG)',  $
    min = min, max = max, /noerase
  xyouts, .6, .02, pim.time, /normal
  xyouts, .5, .95, 'SIMULATED TEST ENA IMAGE!C(No Poisson noise)', $
    CHARS = 2, /norm, align = .5
  closeps, /qui

  !p.thick = 1
  !p.charthick = 1

  return
end