;+
; NAME:
;  yl89q23
; PURPOSE: (one line)
;  radiative heating rate (Q) from Yelle and Lunine 1989 FOR 2.3 um
; DESCRIPTION:
;  radiative heating rate (Q) from Yelle and Lunine 1989 for 1 pt
; CALLING SEQUENCE:
;  yl89q23, t,n,mixrat, h, dh, p10=p10, distau=distau
; INPUTS:
; KEYWORD INPUT PARAMETERS
;  p10 - probability that a V-T transition occurs during a collision
; OUTPUTS:
;  h - heading rate (erg/cm^3/s)
;  dh - temperature derivative dH/dt (erg/cm^3/s/K)
; MODIFICATION HISTORY:
;  Written 2004 Apr, by Leslie Young, SwRI
;-

pro yl89q23, t,n,mixrat, h, dh, p10=p10, distau=distau

if not keyword_set(p10) then p10 = 1e-6
if not keyword_set(distau) then distau=30.121473

physconstants

; set up constants
wl = 3.3 ; micron
wn = 1e4/wl ; wavenumber in cm
;piF = 23. * (wl/wn) ; erg/cm^2/s/(cm^-1)
piF = 23. * (wl/wn)*(30./distau)^2 ; erg/cm^2/s/(cm^-1)
piF = 0.3 * piF * (1.025/1.2) 

Nlosch = 2.6868e19    ; Noschmidt number
b23 = 1.17e-17; band strength in cm^2 cm^-1 FOR 3.3 mu

;a10 = 4.24 ; s^-1 ; 3.3 um
a10 = 25.8

moldiam = 3.7405e-8
mun2 = 28.
much4 = 16.0426
molwgt = (mun2*much4)/(mun2+much4)

vteff = yl89_vtefficiency(t, n, a10, p10, moldiam, molwgt)
eps = vteff/(1.d + vteff)

q = (piF/4) * eps * (n*mixrat) * b

h = q
dh = 0


end