博文
喷流强度计算（IDL程序）
|||
;+
;NAME:
; jetpower_falcke99
;PURPOSE:
; compute the jet power
;CALLING SEQUENCE:
; jetpower=jetpower_falcke99(lumdist, flux, mass_bh, nu=nu,
; angle=angle)
;INPUT:
; lumdist ---- Luminosity distance in Unit: Mpc
; flux ---- flux in Unit: mJy
; mass_bh ---- Black Hole Mass in Unit: Solar Mass
; nu ---- frequency in Unit: GHz
; angle ---- angle in Unit: degree
;OUTPUT:
; result ---- jet power in Unit: erg/s
;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
;-
;=============================================================================
Function correct_factor, angle=angle, select=select
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*angle0
endcase
return,result
End
;============================================================================
Function jetpower_falcke99, lumdist, flux, mass_bh, nu=nu, angle=angle
cor0=correct_factor(angle=angle,select=0)
cor1=correct_factor(angle=angle,select=1)
cor2=correct_factor(angle=angle,select=2)
cor3=correct_factor(angle=angle,select=3)
cor4=correct_factor(angle=angle,select=4)
aa=0.024^(cor0/cor1)*1.56^(cor4/cor1)
bb=(lumdist*1.0D2)^(1.6/cor1)*(flux)^(0.79/cor1)
cc=(cor3)^(0.79/cor1)*(mass_bh/33)^(0.15*cor2/cor1)*(nu/8.5)^(0.15*cor2/cor1)
jetpower=1.0D+39*aa*bb/cc
return,jetpower
End
https://blog.sciencenet.cn/blog-456360-658892.html
上一篇：用x射线的积分流量计算单色流量（IDL程序）
下一篇：黑洞质量计算（IDL程序）
