stratified_blr.f90

module stratified_blr
  ! BLR template from Finke 2016, Appendix
  use const, only: dp, h, me, c
  implicit none
  public :: n_BLR, lambda_BLR, eps_BLR, lumi_BLR, rad_BLR
  private
  integer, parameter :: n_BLR = 26
  real(dp), dimension(n_BLR), parameter :: &
    ! wavelength in cm, converted from \AA
    lambda_BLR = [937.80_dp, 949.74_dp, 977.02_dp, 990.69_dp, 1025.72_dp, &
      1033.83_dp, 1066.66_dp, 1215.67_dp, 1304.35_dp, 1306.82_dp, 1396.76_dp, &
      1402.06_dp, 1549.06_dp, 1718.55_dp, 1721.89_dp, 1908.73_dp, 2423.83_dp, &
      2798.75_dp, 3188.67_dp, 4102.89_dp, 4341.68_dp, 4687.02_dp, 4862.68_dp, &
      5539.43_dp, 5877.29_dp, 6564.61_dp]*1e-8_dp, &
    eps_BLR = (h*c/lambda_BLR)/(me*c**2), &
    ! luminosity in units of L(Hβ)
    lumi_BLR = [0.24_dp, 0.24_dp, 0.60_dp, 0.60_dp, 1.1_dp, 1.1_dp, 0.094_dp, &
      12.0_dp, 0.23_dp, 0.23_dp, 1.0_dp, 1.0_dp, 2.9_dp, 0.030_dp, 0.030_dp, &
      1.8_dp, 0.051_dp, 1.7_dp, 0.051_dp, 0.12_dp, 0.30_dp, 0.016_dp, 1.0_dp, &
      0.039_dp, 0.092_dp, 3.6_dp], &
    ! radius in units of R(Hβ)
    rad_BLR = [2.7_dp, 2.8_dp, 0.83_dp, 0.85_dp, 1.2_dp, 1.2_dp, 4.5_dp, &
      0.27_dp, 4.0_dp, 4.0_dp, 0.83_dp, 0.83_dp, 0.83_dp, 3.8_dp, 3.8_dp, &
      0.46_dp, 5.8_dp, 0.45_dp, 4.3_dp, 3.4_dp, 3.2_dp, 0.63_dp, 1.0_dp, &
      4.8_dp, 0.39_dp, 1.3_dp]
end module stratified_blr