Function correct_factor, angle=angle, select=select
;+
;NAME:
; correct_factor
;PURPOSE:
; compute the correct_factor
;CALLING SEQUENCE:
; result=correct_factor(angle=angle, select=select)
;INPUT:
; angle ---- Unit:degree
; select
;OUTPUT:
; result
;METHOD:
; see Falcke, Heino; Biermann, Peter L. 1999A&A...342...49F
; The jet/disk symbiosis. III. What the radio cores in GRS
; 1915+105, NGC 4258, M 81 and SGR A* tell us about accreting
; black holes
;REVISION HISTORY:
; Original by DL.Wang,Oct-08-2008
;-
angle0=angle*!dtor
case select of
0:result=2.38-1.90*angle0+0.520*angle0*angle0
1:result=1.12-0.19*angle0+0.067*angle0*angle0
2:result=-0.155+1.79*angle0-0.634*angle0*angle0
3:result=0.33+0.60*angle0+0.045*angle0*angle0
4:result=0.68+0.50*angle0-0.177*angle0*angle0
5:result=0.09+0.80*angle0+0.103*angle0*angle0
6:result=1.19-0.29*angle0+0.101*angle0*angle
endcase
return,result
End
https://blog.sciencenet.cn/blog-456360-658883.html
上一篇:
AGN宽线区光度计算(IDL程序)下一篇:
能量不同单位表达的转化(IDL程序)