pro pretty_mroi, epsf=epsf, mroi=mroi, fmroi=fmroi, hef=hef, pfss=pfss, foot=foot, $
  trim=trim

;
; epsf= file for eps output
; mroi= map variable mroi 
; fmroi= file containing file for mroi 
;
;

if keyword_set(epsf) then begin
	set_plot, 'ps'
	device, /encaps, /color, bits=8, /cmyk, $
	; file='~/Desktop/MROI_He_PFSS_CR2005_900.eps', $
	file=epsf, $
	/inch, xsize=6.5, ysize=3.75
	!p.font = 0
endif

if not keyword_set(mroi) then begin
  if keyword_set(fmroi) then begin
	restore, fmroi
	s1=strpos(fmroi, 'mroi', 0)
	s2=strpos(fmroi, '.', 0)
	r=execute('mroi= '+strmid(fmroi, s1, s2-s1))
	aa=size(mroi.data)
	mroi.dx=360./aa[1]
	mroi.dy=180./aa[2]
  endif
endif

loadct, 3
plot_map, mroi, dmin=1,dmax=4096, /log, smooth=3, $
   yticks=4, yminor=3, yrange=90*[-1,1], /ystyle, xticks=4, xminor=6, $
   bottom=11, $
   ; title='SOHO-SOI MDI Synoptic Map CR 2005', $
   xtitle='Carrington Longitude (degrees)', ytitle='Latitude (degrees)', $
   xtickname=['0', '90', '180', '270', '360']; , ytickname=['-90', '-45', '0', '45', '90']
plot_map_colorbar, [1,4096], /log       

if keyword_set(pfss) then begin
; Current plot neutral line
    restore, pfss
	contour, br[*,*,38], lon, lat, levels=[0], /down, c_thick=3, /over, c_color=255 
	contour, br[*,*,38], -360+lon, lat, levels=[0], /down, c_thick=3, /over, c_color=255 
endif

if keyword_set(hef) then begin
; He
	set_line_color
	fits2map, hef, h
	h.xc=180.
	lev=0.1*median(h.data[where(h.data gt 0)])
;	plot_map, h, /over, levels=lev, c_color=2, c_thick=2, smooth=9
	hx=get_map_xp(h)
	hx=hx[*,0]
	hy=get_map_yp(h)
	hy=hy[0,*]
	contour, smooth(h.data, 9), hx, hy, levels=[lev], $                                 
	  PATH_XY=xy, PATH_INFO=info, /path_data_coords, /closed
	xxx=reverse(sort(info.n))                     
    ;jj=0
    for jj=0,8 do begin
	  S = [INDGEN(info[xxx[jj]].N), 0]
	  plots, xy[0,INFO[xxx[jj]].OFFSET + S], xy[1,INFO[xxx[jj]].OFFSET + S], color=2, thick=2
	endfor
endif

if keyword_set(foot) then begin
; Footpoint
	set_line_color
	restore, /verb, foot
	dtheta1= results.dtheta[*,1] / cos(results.dtheta[*,0] * !DTOR)
	dphi1= results.dphi[*,1] / cos(results.dtheta[*,0] * !DTOR)
	if keyword_set(trim) then $
  	  yyy=where(dphi1 le trim) else $
	  yyy=where(dphi1 le 4)
	theta1= 90* sin(results.dtheta[*,0] * !DTOR)
	oplot, -360+results.dphi[yyy,0]+results.l0[yyy], theta1[yyy], psym=8, color=5, symsize=0.33
	oplot, results.dphi[yyy,0]+results.l0[yyy], theta1[yyy], psym=8, color=5, symsize=0.33
	oploterr, -360+results.dphi[yyy,0]+results.l0[yyy], theta1[yyy], $ 
	   results.dphi[yyy,1], dtheta1[yyy], /noconn, errcolor=5
	oploterr, results.dphi[yyy,0]+results.l0[yyy], theta1[yyy], $
	   results.dphi[yyy,1], dtheta1[yyy], /noconn, errcolor=5
endif

; Text
; fd= [0, 15, 30, 45, 56, 71, 85, 100, 114, 129, 145, 160, 175, 190, 205, 219, 220, 235, $
;      249, 264, 279, 294, 308, 323, 338, 354, 369, 383, 398, 413, 428]
; for i=1,27 do xyouts, results.dphi[fd[i],0]+results.l0[fd[i]], theta1[fd[i]]+12.5, $
;    strmid(results.mdidates[fd[i]], 0, strpos(results.mdidates[fd[i]], '-', 0)), $
;    color=5, chars=0.75
;for i=28,31 do xyouts, results.dphi[fd[i],0]+results.l0[fd[i]], theta1[fd[i]]+7.5, $
;   strmid(results.mdidates[fd[i]], 0, strpos(results.mdidates[fd[i]], '-', 0)), $
;   color=5, chars=1.25

if keyword_set(epsf) then begin
  device, /close
endif

return
end