C1A, C1B
   R(dome) = C1A + C1B*Nd
C2A, C2B, C2C
    1/T(dome) = C2A + C2B*log(Rdome) + C2C*(log(Rdome))^3 T in degK
C3A, C3B
   R(body) = C3A + C3B*Nb
C4A, C4B, C4C
    1/T(body) = C4A + C4B*log(Rbody) + C4C*(log(Rbody))^3 T in degK
C5A, C5B
    V(thermopile) = C5A + C5B*Nthermopile
C6A, C6B
     dT = alpha*V(thermopile)
    Ts = T(body) + dT
    LW = sigma*Ts**4 + C6A*dT + C6B*sigma*(Ts**4 - Tdome**4)
C7A, C7B
    LWcorr = C7A + C7B*LW

V(thermopile) is in microvolts. Alpha is 694e-6 degC per microvolt. 
Sigma is 5.670e-8.

If I recall correctly, log is to the base e.