function hxr_cnt2phot, pl_exponent=pl_exponent,time=time, edges=edges
    

;WRITE RUN_BONOMO.COM
checkvar, pl_exponent, 2.4 ;Use the value for the diffuse cosmic background

if n_elements(time) eq 0 then begin
    time= ''
    read,'Enter time of observation from 80/02/14 - 89/11/15: ',time
endif
comfile = ['$ bonomo ','time,'+time,'enbe,10.','enen,1000.', $
    'comp','1.,'+string(float(pl_exponent)),'exit',$
    '$ pur run_bonomo.com','$ exit']

openw,lu,/get,'run_bonomo.com'
for i=0,n_elements(comfile)-1 do printf,lu,comfile(i)
free_lun,lu

spawn,'@run_bonomo',result
for i=0,n_elements(result)-1 do print,result(i)


istart = (where(strpos(result,'EDGE      MEAN     FLUX ') ne -1))(0)+1
buff=result(istart:istart+15)
arr = fltarr(6,15)
max_e=0.0
reads,buff,arr,max_e
edges= [(arr(1,*))(*),max_e]
conv=(arr(5,*))(*)


width = edges(1:*) - edges

scaling = 1./(width*68.62*conv)

print, 'Scaling factor for a single power law of exponent ',pl_exponent
print, 'photon flux (ph/cm2/s/keV) = count rate (cnts/s) * scaling '

return, scaling

end