;+
; Project:
;     SDAC
; Name:
;     GOES_GET_CHIANTI_EM
;
; Usage:
;     goes_get_chianti_em, fl, temperature, emission_meas, sat=goes, /photospheric, em49=em49
;
;Purpose:
;     Called by GOES_CHIANTI_TEM to derive temperature and emission measures.
;     This procedures computes the emission measure of solar plasma from the
;     temperature derived from the B4/B8 ratio together with the flux B8 in the
;     1-8 Angstrom channel using CHIANTI spectral models with coronal or
;     photospheric abundances
;     WARNING: fluxes are asssumed to be TRUE fluxes, so corrections
;     such as the (0.70,0.85) scaling of GOES 8-12 must be applied before
;     use of this routine. GOES_CHIANTI_TEM applies these corrections.
;
;     EM49 is EM in units of 10^49, used in some applications
;
;Category:
;     GOES, SPECTRA
;
;Method:
;     From the temperature the emission measure per unit B8 flux is computed
;     from a spline fit from a lookup table for 101 temperatures logT=.02 apart
;
;Inputs:
;     FL - GOES long wavelength flux in Watts/meter^2
;     TEMP - GOES temperature derived from GOES_GET_CHIANTI_TEMP in units of MK
;
;Keywords:
;     sat  - GOES satellite number, needed to get the correct response
;     photospheric - use photospheric abundances rather than the default
;             coronal abundnaces
;
;Outputs:
;     Emission_meas - Emission measure in units of cm-3 (i.e., NOT scaled)
;
;Common Blocks:
;     None.
;
;Needed Files:
;     None
;
; MODIFICATION HISTORY:
;     Stephen White, 24-Mar-2004: Initial version based on CHIANTI 4.2
;     Kim Tolbert, 9-Nov-2010, For GOES 15, use same table as GOES 14 (by ensuring gsat<13) until
;       tables for GOES 15 are generated.
;     27-Jun-2012 (online date), Kim, Richard, use CHIANTI 7.0
;     S. White, 27-Dec-2013: add GOES 15 (very similar to GOES 14), use CH 7.1
;     Note from rschwartz70@gmail.com 2-oct-2020
;       goes_get_chianti_em and goes_get_chiant_temp are created in lockstep
;       so past revisions to one whether in using transfer functions or Chianti versions
;       are locked together
;     SW, 2022 Jan: corrected GOES 13-15 for using wrong G-bar (2.5 A normalization
;                   instead of 3.5) and added 16+17 responses
;     This routine created 01/26/22 using CHIANTI version 10.0.1
;
; Contact     : Stephen White (smwhitemd at gmail.com)
;
;-
;-------------------------------------------------------------------------

pro goes_get_chianti_em, b8, temp, em, sat=sat, photospheric=photospheric, em49=em49

; interpolate tables of b8 vs temp to determine emission measure
; in units of cm^-3. Data values are responses to CHIANTI spectra for
; coronal and photospheric abundance
; EM49 is EM in units of 10^49

b8_cor=fltarr(17,101)     ; responses for 101 temps for each of 17 GOES satellites

 b8_cor[0,*]=[5.28e-05,9.84e-05,1.80e-04,3.21e-04,5.54e-04,9.30e-04,1.52e-03,2.41e-03,$
    3.74e-03,5.70e-03,8.55e-03,1.27e-02,1.85e-02,2.69e-02,3.85e-02,5.44e-02,7.58e-02,$
    1.04e-01,1.41e-01,1.88e-01,2.47e-01,3.20e-01,4.09e-01,5.16e-01,6.45e-01,7.98e-01,$
    9.79e-01,1.19e+00,1.44e+00,1.72e+00,2.04e+00,2.42e+00,2.84e+00,3.31e+00,3.83e+00,$
    4.42e+00,5.06e+00,5.76e+00,6.52e+00,7.34e+00,8.22e+00,9.15e+00,1.01e+01,1.11e+01,$
    1.22e+01,1.33e+01,1.44e+01,1.55e+01,1.67e+01,1.79e+01,1.90e+01,2.02e+01,2.14e+01,$
    2.26e+01,2.37e+01,2.48e+01,2.58e+01,2.68e+01,2.78e+01,2.87e+01,2.96e+01,3.05e+01,$
    3.14e+01,3.24e+01,3.33e+01,3.44e+01,3.54e+01,3.66e+01,3.78e+01,3.91e+01,4.04e+01,$
    4.18e+01,4.33e+01,4.49e+01,4.65e+01,4.81e+01,4.98e+01,5.16e+01,5.34e+01,5.52e+01,$
    5.70e+01,5.88e+01,6.06e+01,6.25e+01,6.42e+01,6.60e+01,6.77e+01,6.94e+01,7.10e+01,$
    7.25e+01,7.40e+01,7.53e+01,7.66e+01,7.77e+01,7.88e+01,7.98e+01,8.06e+01,8.14e+01,$
    8.20e+01,8.26e+01,8.31e+01]
 b8_cor[1,*]=[5.43e-05,1.01e-04,1.85e-04,3.30e-04,5.70e-04,9.56e-04,1.56e-03,2.47e-03,$
    3.84e-03,5.85e-03,8.78e-03,1.30e-02,1.91e-02,2.76e-02,3.95e-02,5.59e-02,7.79e-02,$
    1.07e-01,1.45e-01,1.93e-01,2.54e-01,3.28e-01,4.20e-01,5.30e-01,6.63e-01,8.20e-01,$
    1.01e+00,1.22e+00,1.47e+00,1.77e+00,2.10e+00,2.48e+00,2.91e+00,3.40e+00,3.94e+00,$
    4.54e+00,5.20e+00,5.92e+00,6.70e+00,7.54e+00,8.45e+00,9.40e+00,1.04e+01,1.15e+01,$
    1.25e+01,1.37e+01,1.48e+01,1.60e+01,1.72e+01,1.84e+01,1.96e+01,2.08e+01,2.20e+01,$
    2.32e+01,2.43e+01,2.55e+01,2.65e+01,2.76e+01,2.86e+01,2.95e+01,3.05e+01,3.14e+01,$
    3.23e+01,3.33e+01,3.43e+01,3.53e+01,3.64e+01,3.76e+01,3.88e+01,4.01e+01,4.15e+01,$
    4.30e+01,4.45e+01,4.61e+01,4.78e+01,4.95e+01,5.12e+01,5.30e+01,5.48e+01,5.67e+01,$
    5.86e+01,6.04e+01,6.23e+01,6.42e+01,6.60e+01,6.78e+01,6.96e+01,7.13e+01,7.30e+01,$
    7.45e+01,7.60e+01,7.74e+01,7.87e+01,7.99e+01,8.10e+01,8.20e+01,8.29e+01,8.36e+01,$
    8.43e+01,8.49e+01,8.54e+01]
 b8_cor[2,*]=[5.43e-05,1.01e-04,1.85e-04,3.30e-04,5.70e-04,9.56e-04,1.56e-03,2.47e-03,$
    3.84e-03,5.85e-03,8.78e-03,1.30e-02,1.91e-02,2.76e-02,3.95e-02,5.59e-02,7.79e-02,$
    1.07e-01,1.45e-01,1.93e-01,2.54e-01,3.28e-01,4.20e-01,5.30e-01,6.63e-01,8.20e-01,$
    1.01e+00,1.22e+00,1.47e+00,1.77e+00,2.10e+00,2.48e+00,2.91e+00,3.40e+00,3.94e+00,$
    4.54e+00,5.20e+00,5.92e+00,6.70e+00,7.54e+00,8.45e+00,9.40e+00,1.04e+01,1.15e+01,$
    1.25e+01,1.37e+01,1.48e+01,1.60e+01,1.72e+01,1.84e+01,1.96e+01,2.08e+01,2.20e+01,$
    2.32e+01,2.43e+01,2.55e+01,2.65e+01,2.76e+01,2.86e+01,2.95e+01,3.05e+01,3.14e+01,$
    3.23e+01,3.33e+01,3.43e+01,3.53e+01,3.64e+01,3.76e+01,3.88e+01,4.01e+01,4.15e+01,$
    4.30e+01,4.45e+01,4.61e+01,4.78e+01,4.95e+01,5.12e+01,5.30e+01,5.48e+01,5.67e+01,$
    5.86e+01,6.04e+01,6.23e+01,6.42e+01,6.60e+01,6.78e+01,6.96e+01,7.13e+01,7.30e+01,$
    7.45e+01,7.60e+01,7.74e+01,7.87e+01,7.99e+01,8.10e+01,8.20e+01,8.29e+01,8.36e+01,$
    8.43e+01,8.49e+01,8.54e+01]
 b8_cor[3,*]=[4.74e-05,8.82e-05,1.61e-04,2.88e-04,4.97e-04,8.34e-04,1.36e-03,2.16e-03,$
    3.35e-03,5.11e-03,7.67e-03,1.14e-02,1.66e-02,2.41e-02,3.45e-02,4.88e-02,6.80e-02,$
    9.34e-02,1.26e-01,1.69e-01,2.21e-01,2.87e-01,3.66e-01,4.63e-01,5.79e-01,7.16e-01,$
    8.78e-01,1.07e+00,1.29e+00,1.54e+00,1.83e+00,2.17e+00,2.54e+00,2.97e+00,3.44e+00,$
    3.96e+00,4.54e+00,5.17e+00,5.85e+00,6.59e+00,7.37e+00,8.20e+00,9.08e+00,1.00e+01,$
    1.09e+01,1.19e+01,1.29e+01,1.39e+01,1.50e+01,1.60e+01,1.71e+01,1.81e+01,1.92e+01,$
    2.02e+01,2.12e+01,2.22e+01,2.32e+01,2.41e+01,2.49e+01,2.58e+01,2.66e+01,2.74e+01,$
    2.82e+01,2.90e+01,2.99e+01,3.08e+01,3.18e+01,3.28e+01,3.39e+01,3.50e+01,3.63e+01,$
    3.75e+01,3.89e+01,4.02e+01,4.17e+01,4.32e+01,4.47e+01,4.63e+01,4.79e+01,4.95e+01,$
    5.11e+01,5.28e+01,5.44e+01,5.60e+01,5.76e+01,5.92e+01,6.07e+01,6.22e+01,6.37e+01,$
    6.50e+01,6.63e+01,6.76e+01,6.87e+01,6.97e+01,7.07e+01,7.15e+01,7.23e+01,7.30e+01,$
    7.36e+01,7.41e+01,7.45e+01]
 b8_cor[4,*]=[4.47e-05,8.31e-05,1.52e-04,2.71e-04,4.69e-04,7.86e-04,1.28e-03,2.03e-03,$
    3.16e-03,4.81e-03,7.22e-03,1.07e-02,1.57e-02,2.27e-02,3.25e-02,4.60e-02,6.41e-02,$
    8.80e-02,1.19e-01,1.59e-01,2.09e-01,2.70e-01,3.45e-01,4.36e-01,5.45e-01,6.74e-01,$
    8.27e-01,1.01e+00,1.21e+00,1.45e+00,1.73e+00,2.04e+00,2.40e+00,2.79e+00,3.24e+00,$
    3.73e+00,4.27e+00,4.87e+00,5.51e+00,6.20e+00,6.94e+00,7.73e+00,8.56e+00,9.42e+00,$
    1.03e+01,1.12e+01,1.22e+01,1.31e+01,1.41e+01,1.51e+01,1.61e+01,1.71e+01,1.81e+01,$
    1.91e+01,2.00e+01,2.09e+01,2.18e+01,2.27e+01,2.35e+01,2.43e+01,2.50e+01,2.58e+01,$
    2.66e+01,2.74e+01,2.82e+01,2.90e+01,2.99e+01,3.09e+01,3.19e+01,3.30e+01,3.42e+01,$
    3.54e+01,3.66e+01,3.79e+01,3.93e+01,4.07e+01,4.21e+01,4.36e+01,4.51e+01,4.66e+01,$
    4.82e+01,4.97e+01,5.12e+01,5.28e+01,5.43e+01,5.58e+01,5.72e+01,5.86e+01,6.00e+01,$
    6.13e+01,6.25e+01,6.37e+01,6.47e+01,6.57e+01,6.66e+01,6.74e+01,6.81e+01,6.88e+01,$
    6.93e+01,6.98e+01,7.02e+01]
 b8_cor[5,*]=[4.56e-05,8.43e-05,1.53e-04,2.72e-04,4.69e-04,7.84e-04,1.27e-03,2.02e-03,$
    3.13e-03,4.75e-03,7.12e-03,1.05e-02,1.54e-02,2.22e-02,3.18e-02,4.48e-02,6.24e-02,$
    8.56e-02,1.16e-01,1.54e-01,2.02e-01,2.62e-01,3.34e-01,4.22e-01,5.27e-01,6.51e-01,$
    7.98e-01,9.69e-01,1.17e+00,1.40e+00,1.66e+00,1.96e+00,2.30e+00,2.68e+00,3.11e+00,$
    3.58e+00,4.09e+00,4.66e+00,5.27e+00,5.93e+00,6.63e+00,7.38e+00,8.16e+00,8.98e+00,$
    9.82e+00,1.07e+01,1.16e+01,1.25e+01,1.34e+01,1.43e+01,1.53e+01,1.62e+01,1.71e+01,$
    1.81e+01,1.90e+01,1.98e+01,2.07e+01,2.14e+01,2.22e+01,2.29e+01,2.36e+01,2.43e+01,$
    2.50e+01,2.57e+01,2.65e+01,2.73e+01,2.81e+01,2.90e+01,2.99e+01,3.09e+01,3.20e+01,$
    3.31e+01,3.42e+01,3.54e+01,3.67e+01,3.80e+01,3.93e+01,4.07e+01,4.21e+01,4.35e+01,$
    4.49e+01,4.63e+01,4.77e+01,4.92e+01,5.06e+01,5.19e+01,5.33e+01,5.46e+01,5.58e+01,$
    5.70e+01,5.81e+01,5.92e+01,6.02e+01,6.11e+01,6.19e+01,6.27e+01,6.33e+01,6.39e+01,$
    6.44e+01,6.49e+01,6.53e+01]
 b8_cor[6,*]=[3.92e-05,7.42e-05,1.38e-04,2.48e-04,4.33e-04,7.33e-04,1.20e-03,1.92e-03,$
    3.00e-03,4.60e-03,6.95e-03,1.04e-02,1.53e-02,2.22e-02,3.20e-02,4.55e-02,6.38e-02,$
    8.80e-02,1.20e-01,1.60e-01,2.11e-01,2.74e-01,3.51e-01,4.45e-01,5.57e-01,6.91e-01,$
    8.49e-01,1.04e+00,1.25e+00,1.50e+00,1.79e+00,2.12e+00,2.50e+00,2.92e+00,3.39e+00,$
    3.91e+00,4.49e+00,5.13e+00,5.82e+00,6.56e+00,7.37e+00,8.22e+00,9.11e+00,1.01e+01,$
    1.10e+01,1.20e+01,1.31e+01,1.41e+01,1.52e+01,1.63e+01,1.74e+01,1.85e+01,1.97e+01,$
    2.08e+01,2.18e+01,2.29e+01,2.39e+01,2.49e+01,2.59e+01,2.68e+01,2.77e+01,2.86e+01,$
    2.95e+01,3.05e+01,3.14e+01,3.25e+01,3.35e+01,3.47e+01,3.59e+01,3.71e+01,3.85e+01,$
    3.99e+01,4.13e+01,4.29e+01,4.44e+01,4.61e+01,4.77e+01,4.94e+01,5.12e+01,5.29e+01,$
    5.47e+01,5.65e+01,5.83e+01,6.00e+01,6.18e+01,6.35e+01,6.52e+01,6.68e+01,6.84e+01,$
    6.98e+01,7.13e+01,7.26e+01,7.38e+01,7.49e+01,7.60e+01,7.69e+01,7.78e+01,7.85e+01,$
    7.92e+01,7.97e+01,8.02e+01]
 b8_cor[7,*]=[6.63e-05,1.21e-04,2.18e-04,3.82e-04,6.52e-04,1.08e-03,1.75e-03,2.75e-03,$
    4.24e-03,6.42e-03,9.57e-03,1.41e-02,2.05e-02,2.95e-02,4.19e-02,5.90e-02,8.18e-02,$
    1.12e-01,1.51e-01,2.00e-01,2.62e-01,3.39e-01,4.32e-01,5.44e-01,6.78e-01,8.37e-01,$
    1.02e+00,1.24e+00,1.49e+00,1.78e+00,2.12e+00,2.50e+00,2.92e+00,3.40e+00,3.93e+00,$
    4.52e+00,5.17e+00,5.87e+00,6.63e+00,7.45e+00,8.32e+00,9.24e+00,1.02e+01,1.12e+01,$
    1.22e+01,1.33e+01,1.44e+01,1.55e+01,1.66e+01,1.77e+01,1.89e+01,2.00e+01,2.11e+01,$
    2.22e+01,2.33e+01,2.43e+01,2.53e+01,2.62e+01,2.71e+01,2.79e+01,2.87e+01,2.95e+01,$
    3.03e+01,3.11e+01,3.20e+01,3.29e+01,3.38e+01,3.48e+01,3.59e+01,3.71e+01,3.83e+01,$
    3.95e+01,4.09e+01,4.23e+01,4.38e+01,4.53e+01,4.68e+01,4.84e+01,5.00e+01,5.17e+01,$
    5.33e+01,5.50e+01,5.67e+01,5.83e+01,6.00e+01,6.16e+01,6.31e+01,6.46e+01,6.61e+01,$
    6.75e+01,6.88e+01,7.01e+01,7.12e+01,7.23e+01,7.32e+01,7.41e+01,7.49e+01,7.56e+01,$
    7.62e+01,7.67e+01,7.72e+01]
 b8_cor[8,*]=[7.00e-05,1.27e-04,2.28e-04,3.99e-04,6.80e-04,1.13e-03,1.82e-03,2.86e-03,$
    4.40e-03,6.65e-03,9.90e-03,1.45e-02,2.11e-02,3.04e-02,4.32e-02,6.06e-02,8.40e-02,$
    1.15e-01,1.55e-01,2.05e-01,2.69e-01,3.47e-01,4.41e-01,5.56e-01,6.92e-01,8.54e-01,$
    1.04e+00,1.27e+00,1.52e+00,1.82e+00,2.16e+00,2.54e+00,2.97e+00,3.46e+00,4.00e+00,$
    4.59e+00,5.24e+00,5.96e+00,6.72e+00,7.55e+00,8.43e+00,9.35e+00,1.03e+01,1.13e+01,$
    1.24e+01,1.34e+01,1.45e+01,1.56e+01,1.68e+01,1.79e+01,1.90e+01,2.02e+01,2.13e+01,$
    2.24e+01,2.35e+01,2.45e+01,2.54e+01,2.64e+01,2.72e+01,2.80e+01,2.88e+01,2.96e+01,$
    3.04e+01,3.12e+01,3.20e+01,3.29e+01,3.38e+01,3.48e+01,3.59e+01,3.70e+01,3.82e+01,$
    3.95e+01,4.08e+01,4.22e+01,4.37e+01,4.52e+01,4.67e+01,4.83e+01,4.99e+01,5.15e+01,$
    5.32e+01,5.48e+01,5.65e+01,5.81e+01,5.98e+01,6.14e+01,6.29e+01,6.44e+01,6.59e+01,$
    6.73e+01,6.86e+01,6.98e+01,7.10e+01,7.20e+01,7.30e+01,7.38e+01,7.46e+01,7.53e+01,$
    7.59e+01,7.64e+01,7.68e+01]
 b8_cor[9,*]=[4.43e-05,8.31e-05,1.53e-04,2.75e-04,4.77e-04,8.04e-04,1.31e-03,2.09e-03,$
    3.26e-03,4.99e-03,7.51e-03,1.12e-02,1.64e-02,2.38e-02,3.42e-02,4.85e-02,6.77e-02,$
    9.33e-02,1.26e-01,1.69e-01,2.22e-01,2.88e-01,3.69e-01,4.67e-01,5.84e-01,7.24e-01,$
    8.88e-01,1.08e+00,1.31e+00,1.57e+00,1.87e+00,2.21e+00,2.59e+00,3.03e+00,3.51e+00,$
    4.05e+00,4.65e+00,5.30e+00,6.00e+00,6.77e+00,7.58e+00,8.45e+00,9.36e+00,1.03e+01,$
    1.13e+01,1.23e+01,1.34e+01,1.45e+01,1.55e+01,1.66e+01,1.78e+01,1.89e+01,2.00e+01,$
    2.11e+01,2.22e+01,2.32e+01,2.42e+01,2.52e+01,2.61e+01,2.70e+01,2.79e+01,2.88e+01,$
    2.97e+01,3.06e+01,3.16e+01,3.26e+01,3.36e+01,3.47e+01,3.59e+01,3.71e+01,3.84e+01,$
    3.98e+01,4.12e+01,4.27e+01,4.43e+01,4.59e+01,4.75e+01,4.92e+01,5.09e+01,5.26e+01,$
    5.44e+01,5.62e+01,5.79e+01,5.96e+01,6.14e+01,6.31e+01,6.47e+01,6.63e+01,6.78e+01,$
    6.93e+01,7.07e+01,7.20e+01,7.32e+01,7.43e+01,7.54e+01,7.63e+01,7.71e+01,7.78e+01,$
    7.85e+01,7.90e+01,7.95e+01]
b8_cor[10,*]=[5.53e-05,1.02e-04,1.86e-04,3.29e-04,5.67e-04,9.48e-04,1.54e-03,2.44e-03,$
    3.78e-03,5.74e-03,8.59e-03,1.27e-02,1.86e-02,2.68e-02,3.83e-02,5.41e-02,7.53e-02,$
    1.03e-01,1.40e-01,1.86e-01,2.44e-01,3.15e-01,4.03e-01,5.09e-01,6.35e-01,7.85e-01,$
    9.62e-01,1.17e+00,1.41e+00,1.69e+00,2.00e+00,2.36e+00,2.77e+00,3.23e+00,3.74e+00,$
    4.31e+00,4.93e+00,5.61e+00,6.35e+00,7.14e+00,7.99e+00,8.88e+00,9.83e+00,1.08e+01,$
    1.18e+01,1.29e+01,1.39e+01,1.50e+01,1.61e+01,1.72e+01,1.84e+01,1.95e+01,2.06e+01,$
    2.17e+01,2.28e+01,2.38e+01,2.48e+01,2.58e+01,2.67e+01,2.75e+01,2.84e+01,2.92e+01,$
    3.00e+01,3.09e+01,3.18e+01,3.27e+01,3.37e+01,3.48e+01,3.59e+01,3.71e+01,3.83e+01,$
    3.97e+01,4.10e+01,4.25e+01,4.40e+01,4.55e+01,4.71e+01,4.88e+01,5.04e+01,5.21e+01,$
    5.38e+01,5.55e+01,5.72e+01,5.89e+01,6.06e+01,6.22e+01,6.38e+01,6.54e+01,6.69e+01,$
    6.83e+01,6.97e+01,7.09e+01,7.21e+01,7.32e+01,7.42e+01,7.51e+01,7.59e+01,7.66e+01,$
    7.72e+01,7.78e+01,7.82e+01]
b8_cor[11,*]=[5.39e-05,9.99e-05,1.82e-04,3.23e-04,5.56e-04,9.31e-04,1.51e-03,2.40e-03,$
    3.71e-03,5.65e-03,8.47e-03,1.25e-02,1.83e-02,2.65e-02,3.78e-02,5.34e-02,7.44e-02,$
    1.02e-01,1.38e-01,1.84e-01,2.41e-01,3.12e-01,3.99e-01,5.04e-01,6.29e-01,7.78e-01,$
    9.53e-01,1.16e+00,1.40e+00,1.67e+00,1.99e+00,2.35e+00,2.75e+00,3.21e+00,3.72e+00,$
    4.28e+00,4.90e+00,5.58e+00,6.31e+00,7.10e+00,7.94e+00,8.84e+00,9.77e+00,1.08e+01,$
    1.18e+01,1.28e+01,1.39e+01,1.50e+01,1.61e+01,1.72e+01,1.83e+01,1.94e+01,2.06e+01,$
    2.17e+01,2.27e+01,2.38e+01,2.48e+01,2.57e+01,2.66e+01,2.75e+01,2.83e+01,2.91e+01,$
    3.00e+01,3.09e+01,3.18e+01,3.27e+01,3.37e+01,3.48e+01,3.59e+01,3.71e+01,3.84e+01,$
    3.97e+01,4.11e+01,4.25e+01,4.40e+01,4.56e+01,4.72e+01,4.88e+01,5.05e+01,5.22e+01,$
    5.39e+01,5.56e+01,5.73e+01,5.90e+01,6.07e+01,6.23e+01,6.39e+01,6.55e+01,6.70e+01,$
    6.84e+01,6.98e+01,7.10e+01,7.22e+01,7.33e+01,7.43e+01,7.52e+01,7.60e+01,7.67e+01,$
    7.73e+01,7.79e+01,7.83e+01]
b8_cor[12,*]=[3.63e-05,6.91e-05,1.29e-04,2.33e-04,4.09e-04,6.93e-04,1.14e-03,1.83e-03,$
    2.86e-03,4.39e-03,6.64e-03,9.91e-03,1.46e-02,2.14e-02,3.08e-02,4.39e-02,6.16e-02,$
    8.51e-02,1.16e-01,1.55e-01,2.04e-01,2.66e-01,3.41e-01,4.32e-01,5.42e-01,6.73e-01,$
    8.27e-01,1.01e+00,1.22e+00,1.47e+00,1.75e+00,2.07e+00,2.44e+00,2.85e+00,3.32e+00,$
    3.84e+00,4.41e+00,5.03e+00,5.71e+00,6.45e+00,7.24e+00,8.08e+00,8.97e+00,9.90e+00,$
    1.09e+01,1.19e+01,1.29e+01,1.40e+01,1.50e+01,1.61e+01,1.72e+01,1.83e+01,1.94e+01,$
    2.06e+01,2.16e+01,2.27e+01,2.37e+01,2.47e+01,2.57e+01,2.66e+01,2.76e+01,2.85e+01,$
    2.94e+01,3.04e+01,3.14e+01,3.24e+01,3.35e+01,3.46e+01,3.59e+01,3.71e+01,3.85e+01,$
    3.99e+01,4.14e+01,4.29e+01,4.45e+01,4.62e+01,4.78e+01,4.96e+01,5.13e+01,5.31e+01,$
    5.49e+01,5.67e+01,5.85e+01,6.03e+01,6.20e+01,6.38e+01,6.55e+01,6.71e+01,6.87e+01,$
    7.02e+01,7.16e+01,7.29e+01,7.42e+01,7.53e+01,7.64e+01,7.73e+01,7.82e+01,7.89e+01,$
    7.96e+01,8.01e+01,8.06e+01]
b8_cor[13,*]=[5.55e-05,1.03e-04,1.87e-04,3.31e-04,5.69e-04,9.51e-04,1.54e-03,2.44e-03,$
    3.79e-03,5.75e-03,8.61e-03,1.27e-02,1.86e-02,2.69e-02,3.84e-02,5.42e-02,7.54e-02,$
    1.03e-01,1.40e-01,1.86e-01,2.44e-01,3.16e-01,4.04e-01,5.09e-01,6.36e-01,7.86e-01,$
    9.63e-01,1.17e+00,1.41e+00,1.69e+00,2.01e+00,2.37e+00,2.78e+00,3.24e+00,3.75e+00,$
    4.31e+00,4.94e+00,5.62e+00,6.35e+00,7.15e+00,8.00e+00,8.89e+00,9.83e+00,1.08e+01,$
    1.18e+01,1.29e+01,1.39e+01,1.50e+01,1.61e+01,1.73e+01,1.84e+01,1.95e+01,2.06e+01,$
    2.17e+01,2.28e+01,2.38e+01,2.48e+01,2.58e+01,2.67e+01,2.75e+01,2.84e+01,2.92e+01,$
    3.00e+01,3.09e+01,3.18e+01,3.27e+01,3.37e+01,3.48e+01,3.59e+01,3.71e+01,3.83e+01,$
    3.96e+01,4.10e+01,4.25e+01,4.40e+01,4.55e+01,4.71e+01,4.87e+01,5.04e+01,5.21e+01,$
    5.38e+01,5.55e+01,5.72e+01,5.89e+01,6.05e+01,6.22e+01,6.38e+01,6.53e+01,6.68e+01,$
    6.83e+01,6.96e+01,7.09e+01,7.21e+01,7.31e+01,7.41e+01,7.50e+01,7.58e+01,7.65e+01,$
    7.72e+01,7.77e+01,7.82e+01]
b8_cor[14,*]=[6.42e-05,1.18e-04,2.12e-04,3.72e-04,6.37e-04,1.06e-03,1.71e-03,2.70e-03,$
    4.16e-03,6.30e-03,9.40e-03,1.38e-02,2.01e-02,2.90e-02,4.13e-02,5.81e-02,8.07e-02,$
    1.10e-01,1.49e-01,1.98e-01,2.59e-01,3.35e-01,4.27e-01,5.38e-01,6.71e-01,8.29e-01,$
    1.01e+00,1.23e+00,1.48e+00,1.77e+00,2.10e+00,2.48e+00,2.90e+00,3.38e+00,3.90e+00,$
    4.49e+00,5.13e+00,5.83e+00,6.59e+00,7.40e+00,8.27e+00,9.18e+00,1.01e+01,1.11e+01,$
    1.22e+01,1.32e+01,1.43e+01,1.54e+01,1.65e+01,1.77e+01,1.88e+01,1.99e+01,2.11e+01,$
    2.22e+01,2.32e+01,2.42e+01,2.52e+01,2.61e+01,2.70e+01,2.78e+01,2.87e+01,2.94e+01,$
    3.02e+01,3.11e+01,3.19e+01,3.28e+01,3.38e+01,3.48e+01,3.59e+01,3.71e+01,3.83e+01,$
    3.96e+01,4.09e+01,4.23e+01,4.38e+01,4.53e+01,4.69e+01,4.85e+01,5.01e+01,5.17e+01,$
    5.34e+01,5.51e+01,5.67e+01,5.84e+01,6.00e+01,6.17e+01,6.32e+01,6.48e+01,6.62e+01,$
    6.76e+01,6.90e+01,7.02e+01,7.14e+01,7.24e+01,7.34e+01,7.43e+01,7.51e+01,7.57e+01,$
    7.64e+01,7.69e+01,7.73e+01]

b8_pho=fltarr(17,101)

 b8_pho[0,*]=[3.13e-05,5.60e-05,9.85e-05,1.70e-04,2.85e-04,4.68e-04,7.50e-04,1.18e-03,$
    1.81e-03,2.73e-03,4.05e-03,5.95e-03,8.61e-03,1.23e-02,1.74e-02,2.44e-02,3.36e-02,$
    4.57e-02,6.13e-02,8.09e-02,1.05e-01,1.35e-01,1.72e-01,2.15e-01,2.67e-01,3.27e-01,$
    3.98e-01,4.81e-01,5.77e-01,6.87e-01,8.13e-01,9.56e-01,1.12e+00,1.30e+00,1.50e+00,$
    1.73e+00,1.97e+00,2.25e+00,2.54e+00,2.87e+00,3.21e+00,3.59e+00,3.99e+00,4.41e+00,$
    4.86e+00,5.33e+00,5.83e+00,6.35e+00,6.90e+00,7.46e+00,8.06e+00,8.67e+00,9.31e+00,$
    9.96e+00,1.06e+01,1.13e+01,1.20e+01,1.27e+01,1.34e+01,1.42e+01,1.49e+01,1.57e+01,$
    1.65e+01,1.73e+01,1.81e+01,1.90e+01,2.00e+01,2.09e+01,2.19e+01,2.29e+01,2.40e+01,$
    2.51e+01,2.63e+01,2.75e+01,2.87e+01,3.00e+01,3.13e+01,3.26e+01,3.39e+01,3.53e+01,$
    3.67e+01,3.81e+01,3.95e+01,4.09e+01,4.23e+01,4.37e+01,4.50e+01,4.64e+01,4.78e+01,$
    4.91e+01,5.04e+01,5.17e+01,5.29e+01,5.41e+01,5.53e+01,5.64e+01,5.75e+01,5.85e+01,$
    5.95e+01,6.05e+01,6.14e+01]
 b8_pho[1,*]=[3.22e-05,5.76e-05,1.01e-04,1.74e-04,2.93e-04,4.81e-04,7.71e-04,1.21e-03,$
    1.86e-03,2.80e-03,4.17e-03,6.11e-03,8.85e-03,1.27e-02,1.79e-02,2.51e-02,3.45e-02,$
    4.70e-02,6.30e-02,8.32e-02,1.08e-01,1.39e-01,1.76e-01,2.21e-01,2.74e-01,3.36e-01,$
    4.09e-01,4.94e-01,5.93e-01,7.06e-01,8.35e-01,9.82e-01,1.15e+00,1.34e+00,1.54e+00,$
    1.77e+00,2.03e+00,2.31e+00,2.61e+00,2.95e+00,3.30e+00,3.69e+00,4.10e+00,4.53e+00,$
    4.99e+00,5.48e+00,5.99e+00,6.53e+00,7.09e+00,7.67e+00,8.28e+00,8.91e+00,9.56e+00,$
    1.02e+01,1.09e+01,1.16e+01,1.23e+01,1.31e+01,1.38e+01,1.46e+01,1.53e+01,1.61e+01,$
    1.69e+01,1.78e+01,1.86e+01,1.96e+01,2.05e+01,2.15e+01,2.25e+01,2.36e+01,2.47e+01,$
    2.58e+01,2.70e+01,2.82e+01,2.95e+01,3.08e+01,3.21e+01,3.35e+01,3.49e+01,3.63e+01,$
    3.77e+01,3.91e+01,4.05e+01,4.20e+01,4.34e+01,4.49e+01,4.63e+01,4.77e+01,4.91e+01,$
    5.05e+01,5.18e+01,5.31e+01,5.44e+01,5.56e+01,5.68e+01,5.80e+01,5.91e+01,6.01e+01,$
    6.12e+01,6.21e+01,6.31e+01]
 b8_pho[2,*]=[3.22e-05,5.76e-05,1.01e-04,1.74e-04,2.93e-04,4.81e-04,7.71e-04,1.21e-03,$
    1.86e-03,2.80e-03,4.17e-03,6.11e-03,8.85e-03,1.27e-02,1.79e-02,2.51e-02,3.45e-02,$
    4.70e-02,6.30e-02,8.32e-02,1.08e-01,1.39e-01,1.76e-01,2.21e-01,2.74e-01,3.36e-01,$
    4.09e-01,4.94e-01,5.93e-01,7.06e-01,8.35e-01,9.82e-01,1.15e+00,1.34e+00,1.54e+00,$
    1.77e+00,2.03e+00,2.31e+00,2.61e+00,2.95e+00,3.30e+00,3.69e+00,4.10e+00,4.53e+00,$
    4.99e+00,5.48e+00,5.99e+00,6.53e+00,7.09e+00,7.67e+00,8.28e+00,8.91e+00,9.56e+00,$
    1.02e+01,1.09e+01,1.16e+01,1.23e+01,1.31e+01,1.38e+01,1.46e+01,1.53e+01,1.61e+01,$
    1.69e+01,1.78e+01,1.86e+01,1.96e+01,2.05e+01,2.15e+01,2.25e+01,2.36e+01,2.47e+01,$
    2.58e+01,2.70e+01,2.82e+01,2.95e+01,3.08e+01,3.21e+01,3.35e+01,3.49e+01,3.63e+01,$
    3.77e+01,3.91e+01,4.05e+01,4.20e+01,4.34e+01,4.49e+01,4.63e+01,4.77e+01,4.91e+01,$
    5.05e+01,5.18e+01,5.31e+01,5.44e+01,5.56e+01,5.68e+01,5.80e+01,5.91e+01,6.01e+01,$
    6.12e+01,6.21e+01,6.31e+01]
 b8_pho[3,*]=[2.81e-05,5.02e-05,8.84e-05,1.52e-04,2.56e-04,4.20e-04,6.73e-04,1.05e-03,$
    1.62e-03,2.45e-03,3.64e-03,5.33e-03,7.73e-03,1.11e-02,1.56e-02,2.19e-02,3.01e-02,$
    4.10e-02,5.50e-02,7.26e-02,9.45e-02,1.21e-01,1.54e-01,1.93e-01,2.39e-01,2.94e-01,$
    3.57e-01,4.31e-01,5.17e-01,6.16e-01,7.29e-01,8.57e-01,1.00e+00,1.17e+00,1.35e+00,$
    1.55e+00,1.77e+00,2.02e+00,2.28e+00,2.57e+00,2.88e+00,3.22e+00,3.58e+00,3.96e+00,$
    4.36e+00,4.78e+00,5.23e+00,5.70e+00,6.19e+00,6.70e+00,7.23e+00,7.78e+00,8.35e+00,$
    8.93e+00,9.53e+00,1.01e+01,1.08e+01,1.14e+01,1.20e+01,1.27e+01,1.34e+01,1.41e+01,$
    1.48e+01,1.55e+01,1.63e+01,1.71e+01,1.79e+01,1.88e+01,1.97e+01,2.06e+01,2.15e+01,$
    2.26e+01,2.36e+01,2.47e+01,2.58e+01,2.69e+01,2.80e+01,2.92e+01,3.04e+01,3.16e+01,$
    3.29e+01,3.41e+01,3.54e+01,3.66e+01,3.79e+01,3.92e+01,4.04e+01,4.16e+01,4.29e+01,$
    4.40e+01,4.52e+01,4.64e+01,4.75e+01,4.85e+01,4.96e+01,5.06e+01,5.16e+01,5.25e+01,$
    5.34e+01,5.42e+01,5.50e+01]
 b8_pho[4,*]=[2.64e-05,4.73e-05,8.33e-05,1.43e-04,2.41e-04,3.96e-04,6.34e-04,9.94e-04,$
    1.53e-03,2.30e-03,3.43e-03,5.02e-03,7.28e-03,1.04e-02,1.47e-02,2.06e-02,2.84e-02,$
    3.86e-02,5.18e-02,6.84e-02,8.91e-02,1.14e-01,1.45e-01,1.82e-01,2.25e-01,2.77e-01,$
    3.37e-01,4.06e-01,4.87e-01,5.80e-01,6.87e-01,8.08e-01,9.44e-01,1.10e+00,1.27e+00,$
    1.46e+00,1.67e+00,1.90e+00,2.15e+00,2.42e+00,2.72e+00,3.03e+00,3.37e+00,3.73e+00,$
    4.11e+00,4.51e+00,4.93e+00,5.37e+00,5.83e+00,6.31e+00,6.81e+00,7.33e+00,7.86e+00,$
    8.42e+00,8.98e+00,9.56e+00,1.01e+01,1.07e+01,1.14e+01,1.20e+01,1.26e+01,1.33e+01,$
    1.39e+01,1.46e+01,1.53e+01,1.61e+01,1.69e+01,1.77e+01,1.85e+01,1.94e+01,2.03e+01,$
    2.12e+01,2.22e+01,2.32e+01,2.43e+01,2.53e+01,2.64e+01,2.75e+01,2.87e+01,2.98e+01,$
    3.10e+01,3.22e+01,3.33e+01,3.45e+01,3.57e+01,3.69e+01,3.81e+01,3.92e+01,4.04e+01,$
    4.15e+01,4.26e+01,4.37e+01,4.47e+01,4.57e+01,4.67e+01,4.77e+01,4.86e+01,4.94e+01,$
    5.03e+01,5.11e+01,5.19e+01]
 b8_pho[5,*]=[2.71e-05,4.83e-05,8.45e-05,1.45e-04,2.43e-04,3.97e-04,6.34e-04,9.90e-04,$
    1.52e-03,2.28e-03,3.39e-03,4.96e-03,7.16e-03,1.02e-02,1.44e-02,2.01e-02,2.77e-02,$
    3.76e-02,5.04e-02,6.65e-02,8.64e-02,1.11e-01,1.40e-01,1.76e-01,2.18e-01,2.67e-01,$
    3.25e-01,3.92e-01,4.69e-01,5.58e-01,6.60e-01,7.76e-01,9.06e-01,1.05e+00,1.22e+00,$
    1.40e+00,1.60e+00,1.81e+00,2.05e+00,2.31e+00,2.59e+00,2.89e+00,3.21e+00,3.55e+00,$
    3.90e+00,4.28e+00,4.68e+00,5.09e+00,5.53e+00,5.98e+00,6.45e+00,6.94e+00,7.44e+00,$
    7.96e+00,8.49e+00,9.03e+00,9.58e+00,1.01e+01,1.07e+01,1.13e+01,1.19e+01,1.25e+01,$
    1.31e+01,1.37e+01,1.44e+01,1.51e+01,1.58e+01,1.66e+01,1.73e+01,1.82e+01,1.90e+01,$
    1.99e+01,2.08e+01,2.17e+01,2.27e+01,2.37e+01,2.47e+01,2.57e+01,2.68e+01,2.78e+01,$
    2.89e+01,3.00e+01,3.11e+01,3.22e+01,3.33e+01,3.44e+01,3.55e+01,3.65e+01,3.76e+01,$
    3.87e+01,3.97e+01,4.07e+01,4.16e+01,4.26e+01,4.35e+01,4.44e+01,4.52e+01,4.60e+01,$
    4.68e+01,4.75e+01,4.82e+01]
 b8_pho[6,*]=[2.30e-05,4.17e-05,7.43e-05,1.29e-04,2.20e-04,3.64e-04,5.89e-04,9.29e-04,$
    1.44e-03,2.19e-03,3.27e-03,4.83e-03,7.04e-03,1.01e-02,1.44e-02,2.03e-02,2.81e-02,$
    3.84e-02,5.17e-02,6.86e-02,8.97e-02,1.16e-01,1.47e-01,1.85e-01,2.30e-01,2.84e-01,$
    3.46e-01,4.19e-01,5.04e-01,6.02e-01,7.15e-01,8.43e-01,9.88e-01,1.15e+00,1.34e+00,$
    1.54e+00,1.76e+00,2.01e+00,2.28e+00,2.58e+00,2.90e+00,3.25e+00,3.62e+00,4.01e+00,$
    4.43e+00,4.87e+00,5.34e+00,5.83e+00,6.34e+00,6.88e+00,7.44e+00,8.02e+00,8.62e+00,$
    9.24e+00,9.88e+00,1.05e+01,1.12e+01,1.19e+01,1.26e+01,1.33e+01,1.40e+01,1.48e+01,$
    1.55e+01,1.63e+01,1.72e+01,1.80e+01,1.89e+01,1.99e+01,2.08e+01,2.18e+01,2.29e+01,$
    2.40e+01,2.51e+01,2.63e+01,2.74e+01,2.87e+01,2.99e+01,3.12e+01,3.25e+01,3.38e+01,$
    3.52e+01,3.65e+01,3.79e+01,3.92e+01,4.06e+01,4.20e+01,4.33e+01,4.46e+01,4.60e+01,$
    4.72e+01,4.85e+01,4.98e+01,5.10e+01,5.21e+01,5.33e+01,5.43e+01,5.54e+01,5.64e+01,$
    5.74e+01,5.83e+01,5.92e+01]
 b8_pho[7,*]=[3.99e-05,7.02e-05,1.21e-04,2.06e-04,3.42e-04,5.55e-04,8.80e-04,1.37e-03,$
    2.08e-03,3.11e-03,4.59e-03,6.68e-03,9.61e-03,1.37e-02,1.92e-02,2.66e-02,3.65e-02,$
    4.94e-02,6.59e-02,8.67e-02,1.12e-01,1.44e-01,1.82e-01,2.27e-01,2.80e-01,3.43e-01,$
    4.16e-01,5.01e-01,5.99e-01,7.11e-01,8.39e-01,9.84e-01,1.15e+00,1.33e+00,1.53e+00,$
    1.76e+00,2.01e+00,2.28e+00,2.57e+00,2.89e+00,3.23e+00,3.60e+00,3.99e+00,4.40e+00,$
    4.84e+00,5.30e+00,5.78e+00,6.28e+00,6.80e+00,7.35e+00,7.92e+00,8.50e+00,9.11e+00,$
    9.73e+00,1.04e+01,1.10e+01,1.17e+01,1.23e+01,1.30e+01,1.37e+01,1.43e+01,1.51e+01,$
    1.58e+01,1.65e+01,1.73e+01,1.81e+01,1.90e+01,1.99e+01,2.08e+01,2.17e+01,2.27e+01,$
    2.38e+01,2.48e+01,2.59e+01,2.71e+01,2.82e+01,2.94e+01,3.06e+01,3.19e+01,3.31e+01,$
    3.44e+01,3.57e+01,3.70e+01,3.82e+01,3.95e+01,4.08e+01,4.21e+01,4.34e+01,4.46e+01,$
    4.58e+01,4.70e+01,4.82e+01,4.93e+01,5.04e+01,5.15e+01,5.25e+01,5.35e+01,5.45e+01,$
    5.54e+01,5.63e+01,5.71e+01]
 b8_pho[8,*]=[4.22e-05,7.40e-05,1.28e-04,2.16e-04,3.58e-04,5.80e-04,9.17e-04,1.42e-03,$
    2.16e-03,3.23e-03,4.76e-03,6.92e-03,9.93e-03,1.41e-02,1.98e-02,2.74e-02,3.76e-02,$
    5.08e-02,6.77e-02,8.89e-02,1.15e-01,1.47e-01,1.86e-01,2.32e-01,2.86e-01,3.50e-01,$
    4.24e-01,5.11e-01,6.10e-01,7.24e-01,8.53e-01,1.00e+00,1.17e+00,1.35e+00,1.56e+00,$
    1.78e+00,2.03e+00,2.31e+00,2.60e+00,2.92e+00,3.27e+00,3.64e+00,4.03e+00,4.44e+00,$
    4.88e+00,5.34e+00,5.82e+00,6.33e+00,6.85e+00,7.40e+00,7.97e+00,8.55e+00,9.16e+00,$
    9.78e+00,1.04e+01,1.11e+01,1.17e+01,1.24e+01,1.30e+01,1.37e+01,1.44e+01,1.51e+01,$
    1.58e+01,1.66e+01,1.73e+01,1.81e+01,1.90e+01,1.99e+01,2.08e+01,2.17e+01,2.27e+01,$
    2.37e+01,2.48e+01,2.59e+01,2.70e+01,2.82e+01,2.94e+01,3.06e+01,3.18e+01,3.30e+01,$
    3.43e+01,3.56e+01,3.69e+01,3.81e+01,3.94e+01,4.07e+01,4.20e+01,4.32e+01,4.45e+01,$
    4.57e+01,4.69e+01,4.80e+01,4.92e+01,5.03e+01,5.13e+01,5.24e+01,5.33e+01,5.43e+01,$
    5.52e+01,5.61e+01,5.69e+01]
 b8_pho[9,*]=[2.61e-05,4.70e-05,8.32e-05,1.44e-04,2.44e-04,4.02e-04,6.47e-04,1.02e-03,$
    1.57e-03,2.38e-03,3.55e-03,5.22e-03,7.58e-03,1.09e-02,1.55e-02,2.17e-02,2.99e-02,$
    4.08e-02,5.48e-02,7.26e-02,9.47e-02,1.22e-01,1.55e-01,1.94e-01,2.41e-01,2.97e-01,$
    3.62e-01,4.38e-01,5.25e-01,6.27e-01,7.43e-01,8.75e-01,1.02e+00,1.19e+00,1.38e+00,$
    1.59e+00,1.82e+00,2.07e+00,2.35e+00,2.65e+00,2.98e+00,3.33e+00,3.70e+00,4.10e+00,$
    4.52e+00,4.97e+00,5.44e+00,5.93e+00,6.45e+00,6.98e+00,7.55e+00,8.13e+00,8.73e+00,$
    9.35e+00,9.99e+00,1.06e+01,1.13e+01,1.20e+01,1.27e+01,1.34e+01,1.41e+01,1.48e+01,$
    1.56e+01,1.64e+01,1.72e+01,1.81e+01,1.89e+01,1.99e+01,2.08e+01,2.18e+01,2.29e+01,$
    2.39e+01,2.50e+01,2.62e+01,2.74e+01,2.86e+01,2.98e+01,3.11e+01,3.24e+01,3.37e+01,$
    3.50e+01,3.63e+01,3.77e+01,3.90e+01,4.04e+01,4.17e+01,4.30e+01,4.44e+01,4.57e+01,$
    4.69e+01,4.82e+01,4.94e+01,5.06e+01,5.17e+01,5.29e+01,5.39e+01,5.50e+01,5.60e+01,$
    5.69e+01,5.78e+01,5.87e+01]
b8_pho[10,*]=[3.29e-05,5.86e-05,1.02e-04,1.76e-04,2.94e-04,4.80e-04,7.67e-04,1.20e-03,$
    1.83e-03,2.76e-03,4.09e-03,5.99e-03,8.65e-03,1.24e-02,1.74e-02,2.43e-02,3.35e-02,$
    4.54e-02,6.08e-02,8.02e-02,1.04e-01,1.34e-01,1.69e-01,2.12e-01,2.62e-01,3.22e-01,$
    3.91e-01,4.72e-01,5.65e-01,6.72e-01,7.95e-01,9.34e-01,1.09e+00,1.27e+00,1.46e+00,$
    1.68e+00,1.92e+00,2.18e+00,2.47e+00,2.78e+00,3.12e+00,3.47e+00,3.86e+00,4.27e+00,$
    4.70e+00,5.15e+00,5.63e+00,6.12e+00,6.64e+00,7.19e+00,7.75e+00,8.33e+00,8.94e+00,$
    9.56e+00,1.02e+01,1.08e+01,1.15e+01,1.22e+01,1.28e+01,1.35e+01,1.42e+01,1.50e+01,$
    1.57e+01,1.65e+01,1.73e+01,1.81e+01,1.90e+01,1.99e+01,2.08e+01,2.18e+01,2.28e+01,$
    2.38e+01,2.49e+01,2.60e+01,2.72e+01,2.84e+01,2.96e+01,3.08e+01,3.21e+01,3.34e+01,$
    3.47e+01,3.60e+01,3.73e+01,3.86e+01,3.99e+01,4.12e+01,4.25e+01,4.38e+01,4.51e+01,$
    4.63e+01,4.76e+01,4.87e+01,4.99e+01,5.10e+01,5.21e+01,5.32e+01,5.42e+01,5.52e+01,$
    5.61e+01,5.70e+01,5.78e+01]
b8_pho[11,*]=[3.21e-05,5.71e-05,1.00e-04,1.72e-04,2.88e-04,4.71e-04,7.52e-04,1.18e-03,$
    1.80e-03,2.72e-03,4.03e-03,5.90e-03,8.52e-03,1.22e-02,1.72e-02,2.40e-02,3.30e-02,$
    4.49e-02,6.01e-02,7.93e-02,1.03e-01,1.32e-01,1.68e-01,2.10e-01,2.60e-01,3.19e-01,$
    3.88e-01,4.68e-01,5.61e-01,6.67e-01,7.89e-01,9.27e-01,1.08e+00,1.26e+00,1.45e+00,$
    1.67e+00,1.91e+00,2.17e+00,2.46e+00,2.77e+00,3.10e+00,3.46e+00,3.84e+00,4.25e+00,$
    4.68e+00,5.13e+00,5.60e+00,6.10e+00,6.62e+00,7.16e+00,7.73e+00,8.31e+00,8.92e+00,$
    9.54e+00,1.02e+01,1.08e+01,1.15e+01,1.21e+01,1.28e+01,1.35e+01,1.42e+01,1.49e+01,$
    1.57e+01,1.65e+01,1.73e+01,1.81e+01,1.90e+01,1.99e+01,2.08e+01,2.18e+01,2.28e+01,$
    2.39e+01,2.49e+01,2.61e+01,2.72e+01,2.84e+01,2.96e+01,3.09e+01,3.21e+01,3.34e+01,$
    3.47e+01,3.60e+01,3.73e+01,3.86e+01,4.00e+01,4.13e+01,4.26e+01,4.39e+01,4.51e+01,$
    4.64e+01,4.76e+01,4.88e+01,5.00e+01,5.11e+01,5.22e+01,5.33e+01,5.43e+01,5.52e+01,$
    5.62e+01,5.71e+01,5.79e+01]
b8_pho[12,*]=[2.12e-05,3.86e-05,6.91e-05,1.21e-04,2.06e-04,3.43e-04,5.56e-04,8.80e-04,$
    1.36e-03,2.08e-03,3.12e-03,4.61e-03,6.74e-03,9.73e-03,1.39e-02,1.95e-02,2.71e-02,$
    3.71e-02,5.00e-02,6.64e-02,8.70e-02,1.12e-01,1.43e-01,1.80e-01,2.24e-01,2.76e-01,$
    3.37e-01,4.09e-01,4.92e-01,5.88e-01,6.99e-01,8.25e-01,9.68e-01,1.13e+00,1.31e+00,$
    1.51e+00,1.73e+00,1.98e+00,2.25e+00,2.54e+00,2.86e+00,3.20e+00,3.56e+00,3.96e+00,$
    4.37e+00,4.81e+00,5.27e+00,5.76e+00,6.27e+00,6.81e+00,7.37e+00,7.95e+00,8.55e+00,$
    9.17e+00,9.81e+00,1.05e+01,1.11e+01,1.18e+01,1.25e+01,1.32e+01,1.40e+01,1.47e+01,$
    1.55e+01,1.63e+01,1.71e+01,1.80e+01,1.89e+01,1.99e+01,2.08e+01,2.19e+01,2.29e+01,$
    2.40e+01,2.51e+01,2.63e+01,2.75e+01,2.87e+01,3.00e+01,3.13e+01,3.26e+01,3.39e+01,$
    3.53e+01,3.66e+01,3.80e+01,3.94e+01,4.07e+01,4.21e+01,4.35e+01,4.48e+01,4.61e+01,$
    4.74e+01,4.87e+01,5.00e+01,5.12e+01,5.24e+01,5.35e+01,5.46e+01,5.57e+01,5.67e+01,$
    5.76e+01,5.86e+01,5.95e+01]
b8_pho[13,*]=[3.31e-05,5.88e-05,1.03e-04,1.76e-04,2.95e-04,4.82e-04,7.69e-04,1.20e-03,$
    1.84e-03,2.77e-03,4.10e-03,6.00e-03,8.67e-03,1.24e-02,1.75e-02,2.44e-02,3.35e-02,$
    4.55e-02,6.09e-02,8.03e-02,1.04e-01,1.34e-01,1.70e-01,2.12e-01,2.63e-01,3.22e-01,$
    3.92e-01,4.73e-01,5.66e-01,6.73e-01,7.96e-01,9.35e-01,1.09e+00,1.27e+00,1.47e+00,$
    1.68e+00,1.92e+00,2.19e+00,2.47e+00,2.78e+00,3.12e+00,3.48e+00,3.86e+00,4.27e+00,$
    4.70e+00,5.15e+00,5.63e+00,6.13e+00,6.65e+00,7.19e+00,7.75e+00,8.34e+00,8.94e+00,$
    9.56e+00,1.02e+01,1.08e+01,1.15e+01,1.22e+01,1.28e+01,1.35e+01,1.42e+01,1.50e+01,$
    1.57e+01,1.65e+01,1.73e+01,1.81e+01,1.90e+01,1.99e+01,2.08e+01,2.18e+01,2.28e+01,$
    2.38e+01,2.49e+01,2.60e+01,2.72e+01,2.84e+01,2.96e+01,3.08e+01,3.21e+01,3.34e+01,$
    3.46e+01,3.59e+01,3.73e+01,3.86e+01,3.99e+01,4.12e+01,4.25e+01,4.38e+01,4.51e+01,$
    4.63e+01,4.75e+01,4.87e+01,4.99e+01,5.10e+01,5.21e+01,5.31e+01,5.42e+01,5.51e+01,$
    5.61e+01,5.69e+01,5.78e+01]
b8_pho[14,*]=[3.85e-05,6.79e-05,1.18e-04,2.00e-04,3.33e-04,5.41e-04,8.59e-04,1.34e-03,$
    2.04e-03,3.05e-03,4.50e-03,6.56e-03,9.44e-03,1.34e-02,1.89e-02,2.63e-02,3.60e-02,$
    4.87e-02,6.51e-02,8.56e-02,1.11e-01,1.42e-01,1.80e-01,2.24e-01,2.77e-01,3.40e-01,$
    4.12e-01,4.96e-01,5.93e-01,7.05e-01,8.32e-01,9.76e-01,1.14e+00,1.32e+00,1.52e+00,$
    1.75e+00,1.99e+00,2.26e+00,2.55e+00,2.87e+00,3.21e+00,3.58e+00,3.97e+00,4.38e+00,$
    4.82e+00,5.27e+00,5.75e+00,6.26e+00,6.78e+00,7.32e+00,7.89e+00,8.48e+00,9.08e+00,$
    9.70e+00,1.03e+01,1.10e+01,1.16e+01,1.23e+01,1.30e+01,1.36e+01,1.43e+01,1.50e+01,$
    1.58e+01,1.65e+01,1.73e+01,1.81e+01,1.90e+01,1.99e+01,2.08e+01,2.17e+01,2.27e+01,$
    2.38e+01,2.48e+01,2.59e+01,2.71e+01,2.82e+01,2.94e+01,3.07e+01,3.19e+01,3.32e+01,$
    3.44e+01,3.57e+01,3.70e+01,3.83e+01,3.96e+01,4.09e+01,4.22e+01,4.34e+01,4.47e+01,$
    4.59e+01,4.71e+01,4.83e+01,4.94e+01,5.05e+01,5.16e+01,5.27e+01,5.36e+01,5.46e+01,$
    5.55e+01,5.64e+01,5.72e+01]

; Adding GOES-R fits
b8_cor[15,*]=[6.50e-05,1.18e-04,2.10e-04,3.67e-04,6.23e-04,1.03e-03,1.66e-03,2.60e-03,$
    3.99e-03,6.03e-03,8.97e-03,1.32e-02,1.91e-02,2.75e-02,3.90e-02,5.48e-02,7.60e-02,$
    1.04e-01,1.40e-01,1.86e-01,2.44e-01,3.15e-01,4.02e-01,5.07e-01,6.32e-01,7.81e-01,$
    9.56e-01,1.16e+00,1.40e+00,1.67e+00,1.99e+00,2.35e+00,2.75e+00,3.20e+00,3.71e+00,$
    4.27e+00,4.89e+00,5.56e+00,6.30e+00,7.08e+00,7.93e+00,8.82e+00,9.76e+00,1.07e+01,$
    1.18e+01,1.28e+01,1.39e+01,1.50e+01,1.61e+01,1.73e+01,1.84e+01,1.96e+01,2.07e+01,$
    2.19e+01,2.30e+01,2.40e+01,2.51e+01,2.60e+01,2.69e+01,2.78e+01,2.87e+01,2.95e+01,$
    3.03e+01,3.12e+01,3.21e+01,3.30e+01,3.40e+01,3.50e+01,3.61e+01,3.72e+01,3.84e+01,$
    3.97e+01,4.10e+01,4.24e+01,4.38e+01,4.53e+01,4.68e+01,4.84e+01,4.99e+01,5.15e+01,$
    5.31e+01,5.47e+01,5.64e+01,5.80e+01,5.95e+01,6.11e+01,6.26e+01,6.40e+01,6.54e+01,$
    6.68e+01,6.80e+01,6.92e+01,7.03e+01,7.13e+01,7.22e+01,7.30e+01,7.37e+01,7.43e+01,$
    7.49e+01,7.53e+01,7.57e+01]
b8_cor[16,*]=[6.89e-05,1.24e-04,2.21e-04,3.85e-04,6.52e-04,1.07e-03,1.73e-03,2.70e-03,$
    4.15e-03,6.25e-03,9.29e-03,1.36e-02,1.97e-02,2.83e-02,4.02e-02,5.64e-02,7.82e-02,$
    1.07e-01,1.44e-01,1.91e-01,2.50e-01,3.23e-01,4.11e-01,5.18e-01,6.46e-01,7.98e-01,$
    9.76e-01,1.18e+00,1.43e+00,1.71e+00,2.02e+00,2.39e+00,2.80e+00,3.26e+00,3.77e+00,$
    4.34e+00,4.97e+00,5.65e+00,6.39e+00,7.18e+00,8.03e+00,8.93e+00,9.88e+00,1.09e+01,$
    1.19e+01,1.30e+01,1.40e+01,1.52e+01,1.63e+01,1.74e+01,1.86e+01,1.97e+01,2.09e+01,$
    2.20e+01,2.31e+01,2.42e+01,2.52e+01,2.62e+01,2.71e+01,2.79e+01,2.88e+01,2.96e+01,$
    3.04e+01,3.13e+01,3.21e+01,3.30e+01,3.40e+01,3.50e+01,3.61e+01,3.72e+01,3.84e+01,$
    3.97e+01,4.10e+01,4.23e+01,4.37e+01,4.52e+01,4.67e+01,4.83e+01,4.98e+01,5.14e+01,$
    5.30e+01,5.46e+01,5.62e+01,5.78e+01,5.93e+01,6.09e+01,6.24e+01,6.38e+01,6.52e+01,$
    6.65e+01,6.78e+01,6.89e+01,7.00e+01,7.10e+01,7.19e+01,7.27e+01,7.34e+01,7.40e+01,$
    7.46e+01,7.50e+01,7.54e+01]

b8_pho[15,*]=[3.93e-05,6.87e-05,1.18e-04,1.99e-04,3.29e-04,5.32e-04,8.40e-04,1.30e-03,$
    1.97e-03,2.94e-03,4.33e-03,6.28e-03,9.01e-03,1.28e-02,1.79e-02,2.49e-02,3.41e-02,$
    4.61e-02,6.14e-02,8.08e-02,1.05e-01,1.34e-01,1.70e-01,2.12e-01,2.62e-01,3.21e-01,$
    3.90e-01,4.71e-01,5.64e-01,6.70e-01,7.93e-01,9.31e-01,1.09e+00,1.26e+00,1.46e+00,$
    1.68e+00,1.92e+00,2.18e+00,2.47e+00,2.78e+00,3.12e+00,3.48e+00,3.87e+00,4.28e+00,$
    4.72e+00,5.19e+00,5.67e+00,6.18e+00,6.72e+00,7.28e+00,7.86e+00,8.46e+00,9.09e+00,$
    9.73e+00,1.04e+01,1.11e+01,1.17e+01,1.24e+01,1.31e+01,1.38e+01,1.45e+01,1.53e+01,$
    1.60e+01,1.68e+01,1.76e+01,1.84e+01,1.93e+01,2.02e+01,2.11e+01,2.21e+01,2.31e+01,$
    2.41e+01,2.51e+01,2.62e+01,2.73e+01,2.85e+01,2.97e+01,3.08e+01,3.21e+01,3.33e+01,$
    3.45e+01,3.58e+01,3.70e+01,3.83e+01,3.95e+01,4.08e+01,4.20e+01,4.33e+01,4.45e+01,$
    4.56e+01,4.68e+01,4.79e+01,4.90e+01,5.00e+01,5.10e+01,5.20e+01,5.29e+01,5.38e+01,$
    5.47e+01,5.55e+01,5.63e+01]
b8_pho[16,*]=[4.18e-05,7.28e-05,1.25e-04,2.10e-04,3.46e-04,5.57e-04,8.78e-04,1.36e-03,$
    2.05e-03,3.06e-03,4.49e-03,6.51e-03,9.33e-03,1.32e-02,1.85e-02,2.56e-02,3.51e-02,$
    4.74e-02,6.31e-02,8.29e-02,1.07e-01,1.37e-01,1.74e-01,2.17e-01,2.68e-01,3.28e-01,$
    3.98e-01,4.80e-01,5.74e-01,6.83e-01,8.07e-01,9.47e-01,1.11e+00,1.28e+00,1.48e+00,$
    1.70e+00,1.95e+00,2.21e+00,2.50e+00,2.82e+00,3.16e+00,3.52e+00,3.91e+00,4.33e+00,$
    4.77e+00,5.23e+00,5.72e+00,6.23e+00,6.77e+00,7.33e+00,7.91e+00,8.52e+00,9.14e+00,$
    9.79e+00,1.04e+01,1.11e+01,1.18e+01,1.25e+01,1.32e+01,1.39e+01,1.46e+01,1.53e+01,$
    1.61e+01,1.68e+01,1.76e+01,1.84e+01,1.93e+01,2.02e+01,2.11e+01,2.21e+01,2.30e+01,$
    2.41e+01,2.51e+01,2.62e+01,2.73e+01,2.84e+01,2.96e+01,3.08e+01,3.20e+01,3.32e+01,$
    3.44e+01,3.57e+01,3.69e+01,3.82e+01,3.94e+01,4.07e+01,4.19e+01,4.31e+01,4.43e+01,$
    4.55e+01,4.66e+01,4.77e+01,4.88e+01,4.99e+01,5.09e+01,5.18e+01,5.28e+01,5.36e+01,$
    5.45e+01,5.53e+01,5.60e+01]

if keyword_set(sat) then gsat=fix(sat-1)>0 else gsat=8-1 ; subtract 1 to get array index
if keyword_set(photospheric) then b8dat=reform(b8_pho[gsat,*]) else b8dat=reform(b8_cor[gsat,*])

; do spline fit

logtemp=findgen(101)*0.02d0     ; temp in MK as in goes_tem
b8_ftn=spl_init(logtemp,b8dat,/double)
denom=spl_interp(logtemp,b8dat,b8_ftn,alog10(temp),/double)
; assume that B8 = flux is in W/m^2, calibrate accordingly
em=1.d55*b8/denom
; per Richard should be in units of 10^49 for SolarSoft
em49=em/1.d49

end