subroutine xskdc0 (ear,ne,param,ifl,photar,photer)
implicit none
integer ni,ngp,nc,ng,nr,np,nb
parameter (ni=2,ngp=5,nc=2,ng=100,nr=100,np=5,nb=100)
integer ne, ifl
real ear(0:ne), param(6), photar(ne), photer(ne)
c This model convolves the input spectrum with
c the line profile from a Keplerian Accretion disk around
c both Schwarzschild & Kerr (arbitrary spin) black holes. In this case, the
c accretion disk is assumed to extend down to the marginally stable orbit for
c a particular black hole spin. The value of the spin parameter is determined
c by inversion of Rms(a) (see e.g. Bardeen, Press & Teulkolsky, 1974). The
c model can be used to fit to the value of Rms & hence determine a, which is
c displayed by setting chatter>=12 in XSPEC.
c
c It requires the use of an extensive set of transfer functions, which must be
c downloaded from the XSPEC website. These transfer functions should be placed
c in $LHEASOFT/../spectral/xspec/manager/XSKDLINE/
c
c The model accepts the following parameters:
c 1: Rms The location of the innermost stable orbit in units of GMc^{-2}
c This parameter can lie in the range 1.0 -> 6.0 & is used to
c determine the black hole spin, by inversion of Rms=Rms(a)
c 2: Rout The location of the outer edge of the disc in units of GMc^{-2}
c 3: Incl The inclination of the observer relative to the spin axis
c of the black hole in Degrees, in the range 1 (face-on) -> 89
c (edge-on), measured in degrees.
c 4: Index Index of radial power law used to generate the line, taken in
c the form r^{-index}
c 5: F(r) Specifies the type of radial emissivity, current options:
c 0: No dependence, F(r) = constant
c 1: Power law depence (see above)
c 2: Zero Torque IBC (Page & Thorne, 1974)
c 6: G(m) Specifies the type of angular emissivity, current options:
c 0: Optically thick, constant, G(m) = cons.
c 1: Optically thin, limb brightened, G(m) \propto 1/m
c 2: Optically thick, limb darkened, G(m) \propto (1+(2.06*m))/2.06
c 3: G(m) \propto m*log(1+m^{-1})
c For more info, see Beckwith & Done (2004)
real eparam(8)
integer i
eparam(1)=0.0d0
eparam(2)=param(1)
eparam(3)=param(2)
eparam(4)=param(3)
eparam(5)=param(4)
eparam(6)=param(5)
eparam(7)=param(6)
eparam(8)=0.0d0
call xskdconv (ear, ne, eparam, ifl, photar, photer)
param(1)=eparam(2)
param(2)=eparam(3)
param(3)=eparam(4)
param(4)=eparam(5)
param(5)=eparam(6)
param(6)=eparam(7)
return
end