; tabulated C0 and C1 from boreas.pro

    ngrid    = 18
    Teffgrid = dindgen(ngrid)/double(ngrid-1)*(6000.-2500.)+2500.
    C0fit = [-9.0329342,-8.8057005,-8.6187019,-8.5343736,-8.5792239,$
             -8.6915425,-8.8033701,-8.8791533,-8.9288180,-8.9604793,$
             -8.9954977,-9.0593624,-9.1566287,-9.2743908,-9.4120161,$
             -9.5781877,-9.7290674,-9.8972636]
    C1fit = [0.78081830,0.68284416,0.60471646,0.56629124,0.57113263,$
             0.59584083,0.61352883,0.61218030,0.59646729,0.57949132,$
             0.56963417,0.57219919,0.58595944,0.60671743,0.63103575,$
             0.65574884,0.67753323,0.69808401]

; perform the polynomial fits

    x = Teffgrid/1000.

    coef0 = POLY_FIT(x,C0fit,6,C0fit_attempt)
    coef1 = POLY_FIT(x,C1fit,6,C1fit_attempt)

; compare with a fit to the fit...

    C0now = 4.0872049 - 48.979961*x + 47.135345*(x^2) - 20.204336*(x^3) $
    + 4.4110828*(x^4) - 0.48112223*(x^5) + 0.020825121*(x^6)

    C1now = 24.562656 - 25.713078*x + 9.8731818*(x^2) - 1.3825986*(x^3) $
    - 0.055190235*(x^4) + 0.031616983*(x^5) - 0.0021585992*(x^6)

    end