;+
;inputs - map
;
;outputs - the data coords of the centres of all AR in the data
;   	 - the NOAA number of those ARs
;   	 - the event_list of all identified flares in any of the AR
;   	 - any error messages
;
;flare_stats,map,xy_data,noaa_num,event_list
;-

PRO flare_stats, map, xy_data, noaa_num, event_list, err=err

err=''

nar_list = get_nar2( map.time , /quiet, err = err, /unique)
IF err NE '' THEN BEGIN
    event_list=-1
    GOTO, no_events
ENDIF

nar_list = drot_nar2( nar_list, map.time, err=err )
IF err NE '' THEN BEGIN
    event_list=-1
    GOTO, no_events
ENDIF

noaa_num = nar_list.noaa

x_helcent = nar_list.x
y_helcent = nar_list.y
x_data = 511 + ( map.xc / map.dx ) + ( x_helcent / map.dx )
y_data = 511 + ( map.yc / map.dy ) + ( y_helcent / map.dy )

xy_data=fltarr(2,N_ELEMENTS(x_data))

xy_data(0,*) = x_data
xy_data(1,*) = y_data

event_list= events2struct( map.time, err=err )

IF err NE '' THEN BEGIN
    err='No flares. '
    event_list=-1
    GOTO, no_events
ENDIF


identified = where(event_list.region NE '',count)

IF count eq 0 THEN BEGIN
    err='No flares assigned to regions. '
    event_list=-1
    GOTO, no_events
ENDIF

event_list = event_list(identified)

no_events:

END