博文
画出BLazars的能谱指数的IDL程序
||
Pro aox_aro_blazar
;+
;NAME:
; aox_aro_blazar
;PURPOSE:
; plot broad band energy index
;CALLING SEQUENCE:
; aox_aro_blazar
;PROCESS:
; xyouts, alog10(), textoidl(), one_arrow
;REVISION HISTORY:
; Original by DL.Wang,Aug-21-2007
;-
aox=findgen(31)*0.1
vr=4.85D9 ;4.85GHz
vo=5.9958D14 ;5000AA
vx=2.4180D17 ;1keV
arx=0.78
aro=-alog10(vx/vo)/alog10(vo/vr)*aox+alog10(vx/vr)/alog10(vo/vr)*arx
plot,aox,aro,linestyle=5,xrange=[0.0,3.0]
xtitle=textoidl('alpha_{ox}')
ytitle=textoidl('alpha_{ro}')
xyouts,0.48,0.04,xtitle,/normal,charsize=3.0
xyouts,0.05,0.50,ytitle,/normal,charsize=3.0,orientation=90.0
x1=[0.55,0.55]
y1=[0.30,0.55]
oplot,x1,y1,linestyle=1
x2=[0.55,0.85]
y2=[0.30,0.10]
oplot,x2,y2,linestyle=1
x3=[0.85,1.46]
y3=[0.10,0.10]
oplot,x3,y3,linestyle=1
x4=[1.46,1.46]
y4=[0.10,0.44]
oplot,x4,y4,linestyle=1
x5=[1.22,1.46]
y5=[0.55,0.44]
oplot,x5,y5,linestyle=1
x6=[0.55,1.20]
y6=[0.55,0.55]
oplot,x6,y6,linestyle=1
one_arrow,0.15,0.18,270,' ',charsize=2,/normal
xyouts,0.20,0.15,'Radio weak',/normal,charsize=2.0
one_arrow,0.60,0.15,0,' ',charsize=2,/normal
xyouts,0.70,0.15,'X-ray weak',/normal,charsize=2.0
k=-alog10(vx/vo)/alog10(vo/vr)
sita=atan(k)/(!pi)*180.0
one_arrow,0.25,0.70,90.0+sita,' ',charsize=2,/normal
xyouts,0.40,0.75,'LBL',/normal,charsize=2.0
one_arrow,0.25,0.70,270.0+sita,' ',charsize=2,/normal
xyouts,0.15,0.50,'HBL',/normal,charsize=2.0
End
https://blog.sciencenet.cn/blog-456360-657721.html
上一篇：在字符上面画点的IDL程序
下一篇：2MASS星等转化为SDSS AB星等的IDL程序
