pro reflect_alpha,alphas,help=help

if keyword_set(help) then begin
	print,'PURPOSE:'
	print,'	Reflects flare data over y=x line, in effect changing the '
	print,'	old alpha data (atan(y/x)).  New values are positive, and '
	print,'	thus can be compared to new values of alpha.'
	print,'HELP:'
	print,'	pro reflect_alpha'
	print,'	inputs: none'
	print,'	outputs:alphas'
	print,'AUTHOR:'
	print,'	Rick Pernak, Goddard Space Flight Center'
return
endif

openr,1,'/home/pernak/RHESSI/paper/spreadsheet/mem.csv'
mem_arr = strarr(25)
readf,1,mem_arr
close,1

openr,1,'/home/pernak/RHESSI/paper/spreadsheet/clean.csv'
clean_arr = strarr(25)
readf,1,clean_arr
close,1

openr,1,'/home/pernak/RHESSI/paper/spreadsheet/pixon.csv'
pixon_arr = strarr(25)
readf,1,pixon_arr
close,1

malpha = dblarr(25)
calpha = dblarr(25)
palpha = dblarr(25)

;do MEM
for ctr=0,24 do begin
	split = str_sep(mem_arr(ctr),9B)
	x1 = double(split[5])
	y1 = double(split[6])
	x2 = double(split[7])
	y2 = double(split[8])
	dx = x1-x2
	dy = y1-y2
	if ((dy lt 0) and (dx gt 0)) or ((dy gt 0) and (dx lt 0)) then begin
		arg=-dy/dx
		malpha[ctr] = atan(arg)
	endif else begin
		arg=dy/dx
		malpha[ctr] = atan(arg)
	endelse
endfor

;do Clean
for ctr=0,24 do begin
	split = str_sep(clean_arr(ctr),9B)
	x1 = double(split[5])
	y1 = double(split[6])
	x2 = double(split[7])
	y2 = double(split[8])
	dx = x1-x2
	dy = y1-y2
	if ((dy lt 0) and (dx gt 0)) or ((dy gt 0) and (dx lt 0)) then begin
		arg=-dy/dx
		calpha[ctr] = atan(arg)
	endif else begin
		arg=dy/dx
		calpha[ctr] = atan(arg)
	endelse
endfor

;do Pixon
for ctr=0,24 do begin
	split = str_sep(pixon_arr(ctr),9B)
	x1 = double(split[5])
	y1 = double(split[6])
	x2 = double(split[7])
	y2 = double(split[8])
	dx = x1-x2
	dy = y1-y2
	if ((dy lt 0) and (dx gt 0)) or ((dy gt 0) and (dx lt 0)) then begin
		arg=-dy/dx
		palpha[ctr] = atan(arg)
	endif else begin
		arg=dy/dx
		palpha[ctr] = atan(arg)
	endelse
endfor

alphas = {mem:malpha,clean:calpha,pixon:palpha}

return
end