PRO plot_mdi, map, just_data=just_data

loadct,0
IF KEYWORD_SET(just_data) THEN BEGIN
    plot_map,map,dmin=-500,dmax=500
    GOTO, finish
ENDIF

plot_map,map,dmin=-500,dmax=500,grid=10,/limb

nar=get_nar2(map.time,err=err,/unique)
IF err NE '' THEN GOTO, finish

nar = drot_nar2( nar, map.time, err=err )
IF err NE '' THEN GOTO, finish

;solve the multiple NOAA AR problem
ars=nar.noaa
indx=fltarr(N_ELEMENTS(ars))

FOR i=0,N_ELEMENTS(ars)-1 DO BEGIN & $ 
    inx=min(where(ars eq ars(i))) & $
    indx(inx)=1 & $
ENDFOR
nar=nar(where(indx eq 1))

loadct,3
oplot_nar2,nar,color=128
loadct,0

finish:

END