Skip to content

/ pds-tools

Permalink
master
Switch branches/tags

pds-tools/solar.py /
Jump to
Code definitions
Code navigation index up-to-date

Go to file
  • Go to file
  • Copy path
  • Copy permalink
 
 
Cannot retrieve contributors at this time
3102 lines (3059 sloc) 73.3 KB
Raw Blame
################################################################################
# Models for the solar flux density at 1 AU, to be used for image
# calibration.
#
# Mark Showalter, PDS Rings Node, December 2011
################################################################################
import numpy as np
import tabulation as tab
# Class constants
AU = 149597870.7 # km
C = 299792.458 # km/sec
# Converts from W/m^2 to erg/s/cm^2
TO_CGS = 1.e7 / 1.e4
# Converts from flux per micron to flux per Angstrom
TO_PER_ANGSTROM = 1.e-4
# Converts from flux per micron to flux per nanometer
TO_PER_NM = 1.e-3
################################################################################
# Models for the solar spectrum at 1 AU in absolute units of W/m^2/micron
################################################################################
# From Colina, L., R. C. Bohlin, and F. Castelli, 1996. The 0.12-2.5 micron
# absolute flux distribution of the Sun for comparison with solar analog stars.
# Astron. J. 112, 307-315.
#
# Column 1 is wavelength in microns
# Column 2 is solar flux at 1 AU in W/m^2/Hz
COLINA_ARRAY = np.array([
0.119500, 1.34269e-18,
0.120500, 1.06999e-17,
0.121500, 9.52491e-17,
0.122500, 2.28803e-17,
0.123500, 1.01461e-18,
0.124500, 6.70275e-19,
0.125500, 5.64521e-19,
0.126500, 7.06755e-19,
0.127500, 4.16554e-19,
0.128500, 3.42433e-19,
0.129500, 8.30623e-19,
0.130500, 2.49545e-18,
0.131500, 8.22167e-19,
0.132500, 8.66106e-19,
0.133500, 3.65061e-18,
0.134500, 9.71820e-19,
0.135500, 8.17521e-19,
0.136500, 6.78593e-19,
0.137500, 6.47147e-19,
0.138500, 9.14775e-19,
0.139500, 1.83748e-18,
0.140500, 1.42836e-18,
0.141500, 9.79902e-19,
0.142500, 1.10336e-18,
0.143500, 1.21897e-18,
0.144500, 1.23414e-18,
0.145500, 1.35663e-18,
0.146500, 1.75199e-18,
0.147500, 2.12249e-18,
0.148500, 2.07099e-18,
0.149500, 2.04349e-18,
0.150500, 2.25837e-18,
0.151500, 2.43692e-18,
0.152500, 3.13014e-18,
0.153500, 3.53133e-18,
0.154500, 5.73032e-18,
0.155500, 6.01256e-18,
0.156500, 5.15069e-18,
0.157500, 4.85219e-18,
0.158500, 4.78339e-18,
0.159500, 5.05043e-18,
0.160500, 5.65788e-18,
0.161500, 6.59066e-18,
0.162500, 7.57332e-18,
0.163500, 8.75741e-18,
0.164500, 9.22951e-18,
0.165500, 1.34151e-17,
0.166500, 1.25718e-17,
0.167500, 1.23288e-17,
0.168500, 1.45822e-17,
0.169500, 1.90190e-17,
0.170500, 2.16453e-17,
0.171500, 2.28768e-17,
0.172500, 2.43741e-17,
0.173500, 2.59527e-17,
0.174500, 3.11314e-17,
0.175500, 3.74844e-17,
0.176500, 4.29370e-17,
0.177500, 5.00440e-17,
0.178500, 5.56467e-17,
0.179500, 5.76053e-17,
0.180500, 6.77727e-17,
0.181500, 8.38527e-17,
0.182500, 8.46909e-17,
0.183500, 8.54608e-17,
0.184500, 7.92253e-17,
0.185500, 8.95792e-17,
0.186500, 1.05919e-16,
0.187500, 1.19293e-16,
0.188500, 1.31747e-16,
0.189500, 1.44541e-16,
0.190500, 1.53733e-16,
0.191500, 1.72064e-16,
0.192500, 1.72116e-16,
0.193500, 1.60183e-16,
0.194500, 2.22012e-16,
0.195500, 2.33932e-16,
0.196500, 2.60995e-16,
0.197500, 2.71977e-16,
0.198500, 2.74069e-16,
0.199500, 2.97805e-16,
0.200500, 3.30848e-16,
0.201500, 3.66079e-16,
0.202500, 3.79881e-16,
0.203500, 4.28110e-16,
0.204500, 4.80783e-16,
0.205500, 5.00059e-16,
0.206500, 5.30275e-16,
0.207500, 6.12630e-16,
0.208500, 7.31097e-16,
0.209500, 1.03990e-15,
0.210500, 1.37120e-15,
0.211500, 1.66168e-15,
0.212500, 1.56699e-15,
0.213500, 1.66689e-15,
0.214500, 2.03697e-15,
0.215500, 1.85841e-15,
0.216500, 1.67279e-15,
0.217500, 1.84036e-15,
0.218500, 2.34110e-15,
0.219500, 2.49932e-15,
0.220500, 2.54485e-15,
0.221500, 2.14450e-15,
0.222500, 2.70889e-15,
0.223500, 3.46657e-15,
0.224500, 3.27573e-15,
0.225500, 2.93781e-15,
0.226500, 2.23724e-15,
0.227500, 2.27021e-15,
0.228500, 2.98963e-15,
0.229500, 2.76912e-15,
0.230500, 3.05319e-15,
0.231500, 3.01122e-15,
0.232500, 3.17036e-15,
0.233500, 2.69786e-15,
0.234500, 2.36627e-15,
0.235500, 3.13942e-15,
0.236500, 2.99997e-15,
0.237500, 3.01701e-15,
0.238500, 2.63206e-15,
0.239500, 2.79599e-15,
0.240500, 2.52204e-15,
0.241500, 3.29617e-15,
0.242500, 4.57020e-15,
0.243500, 4.34437e-15,
0.244500, 4.02143e-15,
0.245500, 3.29116e-15,
0.246500, 3.38153e-15,
0.247500, 3.75057e-15,
0.248500, 3.11131e-15,
0.249500, 3.86988e-15,
0.250500, 4.17228e-15,
0.251500, 3.22581e-15,
0.252500, 2.92139e-15,
0.253500, 3.65574e-15,
0.254500, 4.25814e-15,
0.255500, 5.69637e-15,
0.256500, 7.36479e-15,
0.257500, 9.15304e-15,
0.258500, 9.20300e-15,
0.259500, 7.75243e-15,
0.260500, 6.41940e-15,
0.261500, 6.75103e-15,
0.262500, 7.89616e-15,
0.263500, 1.24909e-14,
0.264500, 1.92332e-14,
0.265500, 1.96370e-14,
0.266500, 1.95139e-14,
0.267500, 1.97859e-14,
0.268500, 1.93795e-14,
0.269500, 1.91195e-14,
0.270500, 2.15520e-14,
0.271500, 1.82195e-14,
0.272500, 1.61400e-14,
0.273500, 1.62983e-14,
0.274500, 1.10451e-14,
0.275500, 1.46209e-14,
0.276500, 2.04461e-14,
0.277500, 1.98793e-14,
0.278500, 1.36119e-14,
0.279500, 7.50850e-15,
0.280500, 8.33037e-15,
0.281500, 1.82368e-14,
0.282500, 2.58056e-14,
0.283500, 2.77411e-14,
0.284500, 2.10194e-14,
0.285500, 1.46680e-14,
0.286500, 2.90368e-14,
0.287500, 3.05115e-14,
0.288500, 2.93994e-14,
0.289500, 4.35028e-14,
0.290500, 5.57442e-14,
0.291500, 5.42579e-14,
0.292500, 4.89752e-14,
0.293500, 5.06817e-14,
0.294500, 4.84786e-14,
0.295500, 5.29174e-14,
0.296500, 4.90460e-14,
0.297500, 4.92694e-14,
0.298500, 4.55690e-14,
0.299500, 4.77355e-14,
0.300500, 4.17886e-14,
0.301500, 4.55973e-14,
0.302500, 4.92160e-14,
0.303500, 6.29952e-14,
0.304500, 6.15369e-14,
0.305500, 6.10752e-14,
0.306500, 5.75185e-14,
0.307500, 6.40597e-14,
0.308500, 6.41336e-14,
0.309500, 5.24908e-14,
0.310500, 6.63080e-14,
0.311500, 7.78755e-14,
0.312500, 7.05428e-14,
0.313500, 7.56877e-14,
0.314500, 7.22200e-14,
0.315500, 6.93743e-14,
0.316500, 6.97244e-14,
0.317500, 8.58353e-14,
0.318500, 7.43056e-14,
0.319500, 7.99990e-14,
0.320500, 9.09963e-14,
0.321500, 7.96163e-14,
0.322500, 7.88488e-14,
0.323500, 7.62016e-14,
0.324500, 8.82206e-14,
0.325500, 1.02289e-13,
0.326500, 1.14863e-13,
0.327500, 1.12472e-13,
0.328500, 1.09015e-13,
0.329500, 1.26718e-13,
0.330500, 1.22099e-13,
0.331500, 1.16758e-13,
0.332500, 1.16214e-13,
0.333500, 1.12689e-13,
0.334500, 1.17996e-13,
0.335500, 1.16750e-13,
0.336500, 1.00873e-13,
0.337500, 1.05583e-13,
0.338500, 1.16227e-13,
0.339500, 1.21617e-13,
0.340500, 1.28307e-13,
0.341500, 1.18606e-13,
0.342500, 1.28263e-13,
0.343500, 1.23574e-13,
0.344500, 1.04127e-13,
0.345500, 1.23301e-13,
0.346500, 1.22391e-13,
0.347500, 1.19825e-13,
0.348500, 1.20039e-13,
0.349500, 1.19796e-13,
0.350500, 1.41755e-13,
0.351500, 1.33257e-13,
0.352500, 1.24247e-13,
0.353500, 1.42071e-13,
0.354500, 1.53410e-13,
0.355500, 1.45359e-13,
0.356500, 1.27124e-13,
0.357500, 1.12005e-13,
0.358500, 9.83348e-14,
0.359500, 1.41312e-13,
0.360500, 1.38583e-13,
0.361500, 1.26451e-13,
0.362500, 1.44377e-13,
0.363500, 1.44965e-13,
0.364500, 1.48160e-13,
0.365500, 1.69598e-13,
0.366500, 1.82285e-13,
0.367500, 1.75615e-13,
0.368500, 1.64757e-13,
0.369500, 1.84702e-13,
0.370500, 1.72668e-13,
0.371500, 1.78653e-13,
0.372500, 1.61434e-13,
0.373500, 1.41720e-13,
0.374500, 1.39414e-13,
0.375500, 1.65017e-13,
0.376500, 1.69657e-13,
0.377500, 2.01179e-13,
0.378500, 2.09086e-13,
0.379500, 1.69239e-13,
0.380500, 1.91331e-13,
0.381500, 1.72726e-13,
0.382500, 1.27217e-13,
0.383500, 1.10944e-13,
0.384500, 1.55798e-13,
0.385500, 1.64857e-13,
0.386500, 1.65317e-13,
0.387500, 1.65536e-13,
0.388500, 1.63269e-13,
0.389500, 1.95500e-13,
0.390500, 2.08640e-13,
0.391500, 2.25647e-13,
0.392500, 1.72778e-13,
0.393500, 9.93242e-14,
0.394500, 1.75534e-13,
0.395500, 2.25800e-13,
0.396500, 1.47468e-13,
0.397500, 1.60744e-13,
0.398500, 2.65634e-13,
0.399500, 2.91693e-13,
0.400500, 2.96474e-13,
0.401500, 3.09929e-13,
0.402500, 3.13796e-13,
0.403500, 3.02228e-13,
0.404500, 3.01731e-13,
0.405500, 3.02091e-13,
0.406500, 2.89469e-13,
0.407500, 2.97414e-13,
0.408500, 3.18617e-13,
0.409500, 3.15998e-13,
0.410500, 2.68910e-13,
0.411500, 3.27357e-13,
0.412500, 3.23895e-13,
0.413500, 3.19299e-13,
0.414500, 3.17381e-13,
0.415500, 3.18365e-13,
0.416500, 3.39962e-13,
0.417500, 3.08676e-13,
0.418500, 3.13687e-13,
0.419500, 3.18362e-13,
0.420500, 3.30576e-13,
0.421500, 3.39690e-13,
0.422500, 3.00393e-13,
0.423500, 3.26366e-13,
0.424500, 3.38807e-13,
0.425500, 3.26382e-13,
0.426500, 3.28497e-13,
0.427500, 3.05024e-13,
0.428500, 3.09960e-13,
0.429500, 2.89486e-13,
0.430500, 2.23780e-13,
0.431500, 3.33873e-13,
0.432500, 3.27484e-13,
0.433500, 3.45948e-13,
0.434500, 3.35327e-13,
0.435500, 3.47538e-13,
0.436500, 3.90983e-13,
0.437500, 3.67797e-13,
0.438500, 3.20515e-13,
0.439500, 3.75061e-13,
0.440500, 3.53505e-13,
0.441500, 4.00406e-13,
0.442500, 4.12402e-13,
0.443500, 3.99450e-13,
0.444500, 4.14671e-13,
0.445500, 3.84529e-13,
0.446500, 4.01065e-13,
0.447500, 4.42385e-13,
0.448500, 4.22167e-13,
0.449500, 4.35629e-13,
0.450500, 4.62764e-13,
0.451500, 4.57250e-13,
0.452500, 4.22779e-13,
0.453500, 4.30978e-13,
0.454500, 4.34853e-13,
0.455500, 4.48877e-13,
0.456500, 4.60358e-13,
0.457500, 4.67485e-13,
0.458500, 4.40757e-13,
0.459500, 4.51195e-13,
0.460500, 4.60136e-13,
0.461500, 4.65526e-13,
0.462500, 4.78667e-13,
0.463500, 4.66151e-13,
0.464500, 4.53513e-13,
0.465500, 4.70642e-13,
0.466500, 4.44727e-13,
0.467500, 4.68434e-13,
0.468500, 4.65549e-13,
0.469500, 4.66603e-13,
0.470500, 4.42051e-13,
0.471500, 4.77191e-13,
0.472500, 4.84666e-13,
0.473500, 4.74825e-13,
0.474500, 4.91167e-13,
0.475500, 4.84843e-13,
0.476500, 4.72429e-13,
0.477500, 5.03203e-13,
0.478500, 4.89279e-13,
0.479500, 5.07671e-13,
0.480500, 4.99747e-13,
0.481500, 5.15359e-13,
0.482500, 5.00952e-13,
0.483500, 5.02038e-13,
0.484500, 4.91664e-13,
0.485500, 4.58932e-13,
0.486500, 4.09091e-13,
0.487500, 4.62720e-13,
0.488500, 4.85889e-13,
0.489500, 4.99576e-13,
0.490500, 5.13617e-13,
0.491500, 4.87262e-13,
0.492500, 4.89247e-13,
0.493500, 4.89168e-13,
0.494500, 5.35261e-13,
0.495500, 5.03041e-13,
0.496500, 5.28875e-13,
0.497500, 5.31270e-13,
0.498500, 4.93330e-13,
0.499500, 5.22843e-13,
0.500500, 4.94904e-13,
0.501500, 4.84875e-13,
0.502500, 5.08781e-13,
0.503500, 5.21567e-13,
0.504500, 5.06087e-13,
0.505500, 5.41715e-13,
0.506500, 5.35150e-13,
0.507500, 5.22235e-13,
0.508500, 5.27862e-13,
0.509500, 5.29114e-13,
0.510500, 5.39765e-13,
0.511500, 5.55762e-13,
0.512500, 5.21707e-13,
0.513500, 5.22067e-13,
0.514500, 5.27753e-13,
0.515500, 5.37138e-13,
0.516500, 4.73555e-13,
0.517500, 4.91587e-13,
0.518500, 4.72950e-13,
0.519500, 5.24888e-13,
0.520500, 5.27773e-13,
0.521500, 5.51446e-13,
0.522500, 5.29520e-13,
0.523500, 5.52194e-13,
0.524500, 5.72987e-13,
0.525500, 5.66970e-13,
0.526500, 4.93540e-13,
0.527500, 5.41179e-13,
0.528500, 5.63682e-13,
0.529500, 5.72064e-13,
0.530500, 5.84380e-13,
0.531500, 5.89882e-13,
0.532500, 5.34037e-13,
0.533500, 5.82250e-13,
0.534500, 5.64731e-13,
0.535500, 6.07009e-13,
0.536500, 5.72935e-13,
0.537500, 5.78445e-13,
0.538500, 5.87369e-13,
0.539500, 5.67317e-13,
0.540500, 5.49893e-13,
0.541500, 5.86775e-13,
0.542500, 5.71458e-13,
0.543500, 5.90491e-13,
0.544500, 5.92666e-13,
0.545500, 6.01791e-13,
0.546500, 5.97028e-13,
0.547500, 5.84584e-13,
0.548500, 5.96298e-13,
0.549500, 6.08727e-13,
0.550500, 6.00333e-13,
0.551500, 6.05421e-13,
0.552500, 5.99521e-13,
0.553500, 6.13395e-13,
0.554500, 6.20834e-13,
0.555500, 6.22748e-13,
0.556500, 6.00018e-13,
0.557500, 6.10421e-13,
0.558500, 5.93086e-13,
0.559500, 6.02187e-13,
0.560500, 6.16008e-13,
0.561500, 6.11852e-13,
0.562500, 6.22762e-13,
0.563500, 6.28685e-13,
0.564500, 6.28551e-13,
0.565500, 6.11779e-13,
0.566500, 6.24500e-13,
0.567500, 6.46526e-13,
0.568500, 6.22402e-13,
0.569500, 6.41800e-13,
0.570500, 6.12630e-13,
0.571500, 6.33493e-13,
0.572500, 6.59708e-13,
0.573500, 6.56431e-13,
0.574500, 6.55570e-13,
0.575500, 6.44852e-13,
0.576500, 6.52737e-13,
0.577500, 6.58896e-13,
0.578500, 6.34903e-13,
0.579500, 6.53135e-13,
0.580500, 6.58966e-13,
0.581500, 6.66620e-13,
0.582500, 6.76115e-13,
0.583500, 6.72659e-13,
0.584500, 6.76054e-13,
0.585500, 6.50361e-13,
0.586500, 6.69739e-13,
0.587500, 6.78617e-13,
0.588500, 6.44549e-13,
0.589500, 5.95857e-13,
0.590500, 6.72616e-13,
0.591500, 6.65244e-13,
0.592500, 6.75317e-13,
0.593500, 6.73114e-13,
0.594500, 6.66759e-13,
0.595500, 6.72390e-13,
0.596500, 6.83334e-13,
0.597500, 6.76157e-13,
0.598500, 6.69680e-13,
0.599500, 6.78403e-13,
0.600500, 6.69572e-13,
0.601500, 6.73724e-13,
0.602500, 6.63640e-13,
0.603500, 6.92510e-13,
0.604500, 6.90542e-13,
0.605500, 6.87771e-13,
0.606500, 6.88484e-13,
0.607500, 6.89973e-13,
0.608500, 6.86353e-13,
0.609500, 6.89005e-13,
0.610500, 6.75053e-13,
0.611500, 6.94327e-13,
0.612500, 6.80280e-13,
0.613500, 6.73717e-13,
0.614500, 6.87935e-13,
0.615500, 6.90176e-13,
0.616500, 6.50480e-13,
0.617500, 6.92241e-13,
0.618500, 7.01385e-13,
0.619500, 6.96732e-13,
0.620500, 7.10014e-13,
0.621500, 6.94271e-13,
0.622500, 7.05964e-13,
0.623500, 6.88847e-13,
0.624500, 6.86920e-13,
0.625500, 6.79159e-13,
0.626500, 7.08402e-13,
0.627500, 7.10665e-13,
0.628500, 7.12932e-13,
0.629500, 7.06794e-13,
0.631000, 6.94112e-13,
0.633000, 7.03621e-13,
0.635000, 7.10214e-13,
0.637000, 7.13833e-13,
0.639000, 7.17023e-13,
0.641000, 7.05388e-13,
0.643000, 7.12867e-13,
0.645000, 7.19957e-13,
0.647000, 7.13769e-13,
0.649000, 6.98078e-13,
0.651000, 7.23971e-13,
0.653000, 7.25259e-13,
0.655000, 6.99209e-13,
0.657000, 6.35702e-13,
0.659000, 7.15609e-13,
0.661000, 7.30158e-13,
0.663000, 7.27121e-13,
0.665000, 7.33860e-13,
0.667000, 7.26480e-13,
0.669000, 7.36067e-13,
0.671000, 7.26143e-13,
0.673000, 7.32881e-13,
0.675000, 7.31925e-13,
0.677000, 7.35297e-13,
0.679000, 7.34756e-13,
0.681000, 7.36138e-13,
0.683000, 7.34034e-13,
0.685000, 7.26390e-13,
0.687000, 7.36647e-13,
0.689000, 7.37920e-13,
0.691000, 7.35624e-13,
0.693000, 7.39889e-13,
0.695000, 7.38015e-13,
0.697000, 7.31960e-13,
0.699000, 7.40832e-13,
0.701000, 7.24743e-13,
0.703000, 7.29933e-13,
0.705000, 7.48331e-13,
0.707000, 7.44628e-13,
0.709000, 7.40313e-13,
0.711000, 7.45032e-13,
0.713000, 7.42756e-13,
0.715000, 7.43132e-13,
0.717000, 7.40204e-13,
0.719000, 7.30078e-13,
0.721000, 7.35800e-13,
0.723000, 7.49317e-13,
0.725000, 7.54583e-13,
0.727000, 7.56509e-13,
0.729000, 7.45452e-13,
0.731000, 7.53517e-13,
0.733000, 7.53085e-13,
0.735000, 7.52042e-13,
0.737000, 7.54987e-13,
0.739000, 7.42286e-13,
0.741000, 7.34657e-13,
0.743000, 7.55029e-13,
0.745000, 7.54977e-13,
0.747000, 7.61405e-13,
0.749000, 7.57749e-13,
0.751000, 7.57014e-13,
0.753000, 7.59246e-13,
0.755000, 7.60866e-13,
0.757000, 7.60646e-13,
0.759000, 7.59780e-13,
0.761000, 7.61946e-13,
0.763000, 7.68427e-13,
0.765000, 7.60042e-13,
0.767000, 7.41550e-13,
0.769000, 7.56717e-13,
0.771000, 7.61289e-13,
0.773000, 7.67780e-13,
0.775000, 7.59012e-13,
0.777000, 7.68059e-13,
0.779000, 7.66223e-13,
0.781000, 7.70163e-13,
0.783000, 7.66956e-13,
0.785000, 7.73495e-13,
0.787000, 7.75470e-13,
0.789000, 7.77435e-13,
0.791000, 7.70759e-13,
0.793000, 7.64653e-13,
0.795000, 7.62479e-13,
0.797000, 7.78452e-13,
0.799000, 7.70849e-13,
0.801000, 7.79478e-13,
0.803000, 7.74481e-13,
0.805000, 7.68718e-13,
0.807000, 7.75997e-13,
0.809000, 7.62488e-13,
0.811000, 7.79522e-13,
0.813000, 7.84073e-13,
0.815000, 7.82297e-13,
0.817000, 7.84016e-13,
0.819000, 7.58680e-13,
0.821000, 7.74548e-13,
0.823000, 7.73295e-13,
0.825000, 7.78503e-13,
0.827000, 7.83007e-13,
0.829000, 7.83882e-13,
0.831000, 7.84006e-13,
0.833000, 7.62015e-13,
0.835000, 7.79735e-13,
0.837000, 7.82731e-13,
0.839000, 7.79007e-13,
0.841000, 7.84977e-13,
0.843000, 7.75897e-13,
0.845000, 7.83371e-13,
0.847000, 7.80994e-13,
0.849000, 7.43387e-13,
0.851000, 7.71483e-13,
0.853000, 7.51952e-13,
0.855000, 6.81020e-13,
0.857000, 7.88634e-13,
0.859000, 7.81358e-13,
0.861000, 7.85001e-13,
0.863000, 7.90233e-13,
0.865000, 7.70877e-13,
0.867000, 7.02664e-13,
0.869000, 7.75619e-13,
0.870000, 7.70980e-13,
0.872000, 7.68881e-13,
0.874000, 7.74844e-13,
0.876000, 7.80022e-13,
0.878000, 7.63958e-13,
0.880000, 7.68264e-13,
0.882000, 7.61855e-13,
0.884000, 7.76922e-13,
0.886000, 7.67948e-13,
0.888000, 7.73092e-13,
0.890000, 7.73217e-13,
0.892000, 7.69098e-13,
0.894000, 7.75094e-13,
0.896000, 7.82825e-13,
0.898000, 7.64933e-13,
0.900000, 7.64047e-13,
0.902000, 7.69173e-13,
0.904000, 7.63917e-13,
0.906000, 7.67301e-13,
0.908000, 7.47073e-13,
0.910000, 7.60033e-13,
0.912000, 7.70437e-13,
0.914000, 7.66729e-13,
0.916000, 7.65637e-13,
0.918000, 7.69878e-13,
0.920000, 7.49886e-13,
0.922000, 7.46837e-13,
0.924000, 7.55515e-13,
0.926000, 7.63339e-13,
0.928000, 7.72122e-13,
0.930000, 7.70865e-13,
0.932000, 7.81558e-13,
0.934000, 7.72883e-13,
0.936000, 7.62253e-13,
0.938000, 7.72982e-13,
0.940000, 7.69719e-13,
0.942000, 7.58875e-13,
0.944000, 7.77229e-13,
0.946000, 7.75778e-13,
0.948000, 7.63805e-13,
0.950000, 7.68946e-13,
0.952000, 7.65456e-13,
0.954000, 7.58053e-13,
0.956000, 7.74811e-13,
0.958000, 7.75135e-13,
0.960000, 7.70552e-13,
0.962000, 7.69734e-13,
0.964000, 7.70868e-13,
0.966000, 7.64102e-13,
0.968000, 7.71772e-13,
0.970000, 7.70793e-13,
0.972000, 7.78405e-13,
0.974000, 7.63606e-13,
0.976000, 7.73159e-13,
0.978000, 7.71948e-13,
0.980000, 7.71297e-13,
0.982000, 7.77247e-13,
0.984000, 7.74315e-13,
0.986000, 7.73697e-13,
0.988000, 7.74058e-13,
0.990000, 7.65718e-13,
0.992000, 7.75843e-13,
0.994000, 7.77882e-13,
0.996000, 7.76726e-13,
0.998000, 7.77503e-13,
1.00000, 7.77085e-13,
1.00200, 7.77741e-13,
1.00400, 7.63526e-13,
1.00600, 7.40869e-13,
1.00800, 7.74896e-13,
1.01000, 7.76769e-13,
1.01200, 7.75583e-13,
1.01400, 7.71375e-13,
1.01600, 7.68066e-13,
1.01800, 7.68347e-13,
1.02000, 7.67327e-13,
1.02200, 7.66495e-13,
1.02400, 7.76764e-13,
1.02600, 7.76822e-13,
1.02800, 7.78847e-13,
1.03000, 7.70221e-13,
1.03200, 7.75646e-13,
1.03400, 7.58176e-13,
1.03600, 7.73752e-13,
1.03800, 7.65345e-13,
1.04000, 7.67055e-13,
1.04200, 7.69905e-13,
1.04400, 7.75093e-13,
1.04600, 7.60280e-13,
1.04800, 7.68206e-13,
1.05000, 7.72487e-13,
1.05200, 7.71776e-13,
1.05400, 7.67566e-13,
1.05600, 7.73787e-13,
1.05800, 7.56614e-13,
1.06000, 7.61280e-13,
1.06200, 7.71343e-13,
1.06400, 7.65913e-13,
1.06600, 7.60727e-13,
1.06800, 7.56467e-13,
1.07000, 7.47025e-13,
1.07200, 7.61030e-13,
1.07400, 7.52664e-13,
1.07600, 7.54572e-13,
1.07800, 7.52379e-13,
1.08000, 7.67966e-13,
1.08200, 7.35189e-13,
1.08400, 7.48009e-13,
1.08600, 7.64466e-13,
1.08800, 7.29987e-13,
1.09000, 7.59994e-13,
1.09200, 7.54488e-13,
1.09400, 7.10099e-13,
1.09600, 7.31668e-13,
1.09800, 7.32810e-13,
1.10000, 7.53167e-13,
1.10200, 7.39340e-13,
1.10400, 7.50516e-13,
1.10600, 7.52031e-13,
1.10800, 7.51041e-13,
1.11000, 7.53541e-13,
1.11200, 7.45236e-13,
1.11400, 7.48728e-13,
1.11600, 7.52343e-13,
1.11800, 7.51935e-13,
1.12000, 7.39024e-13,
1.12200, 7.53811e-13,
1.12400, 7.51976e-13,
1.12600, 7.35279e-13,
1.12800, 7.48094e-13,
1.13000, 7.39450e-13,
1.13200, 7.46645e-13,
1.13400, 7.44080e-13,
1.13600, 7.52130e-13,
1.13800, 7.38265e-13,
1.14000, 7.34401e-13,
1.14200, 7.40252e-13,
1.14400, 7.41946e-13,
1.14600, 7.54398e-13,
1.14800, 7.51080e-13,
1.15000, 7.51305e-13,
1.15200, 7.49650e-13,
1.15400, 7.55196e-13,
1.15600, 7.50879e-13,
1.15800, 7.52035e-13,
1.16000, 7.37603e-13,
1.16200, 7.30451e-13,
1.16400, 7.40477e-13,
1.16600, 7.42897e-13,
1.16800, 7.48293e-13,
1.17000, 7.38459e-13,
1.17200, 7.51396e-13,
1.17400, 7.52896e-13,
1.17600, 7.34453e-13,
1.17800, 7.32650e-13,
1.18000, 7.53040e-13,
1.18200, 7.44843e-13,
1.18400, 7.22246e-13,
1.18600, 7.47360e-13,
1.18800, 7.24513e-13,
1.19000, 7.41124e-13,
1.19200, 7.49974e-13,
1.19400, 7.51226e-13,
1.19600, 7.42256e-13,
1.19800, 7.10901e-13,
1.20000, 7.32947e-13,
1.20200, 7.49770e-13,
1.20400, 7.21925e-13,
1.20600, 7.46362e-13,
1.20800, 7.14375e-13,
1.21000, 7.32725e-13,
1.21200, 7.43239e-13,
1.21400, 7.46823e-13,
1.21600, 7.49254e-13,
1.21800, 7.46965e-13,
1.22000, 7.45715e-13,
1.22200, 7.45989e-13,
1.22400, 7.44534e-13,
1.22600, 7.45692e-13,
1.22800, 7.36164e-13,
1.23000, 7.45193e-13,
1.23200, 7.41753e-13,
1.23400, 7.43028e-13,
1.23600, 7.44193e-13,
1.23800, 7.46781e-13,
1.24000, 7.38000e-13,
1.24200, 7.38378e-13,
1.24400, 7.37365e-13,
1.24600, 7.41591e-13,
1.24800, 7.42417e-13,
1.25000, 7.44604e-13,
1.25200, 7.35462e-13,
1.25400, 7.40681e-13,
1.25600, 7.38137e-13,
1.25800, 7.37776e-13,
1.26000, 7.40341e-13,
1.26200, 7.37112e-13,
1.26400, 7.37111e-13,
1.26600, 7.36791e-13,
1.26800, 7.27618e-13,
1.27000, 7.40958e-13,
1.27200, 7.40272e-13,
1.27400, 7.35379e-13,
1.27600, 7.40361e-13,
1.27800, 7.38387e-13,
1.28000, 7.28681e-13,
1.28200, 6.54730e-13,
1.28400, 7.18978e-13,
1.28600, 7.29205e-13,
1.28800, 7.32614e-13,
1.29000, 7.27365e-13,
1.29200, 7.34588e-13,
1.29400, 7.31495e-13,
1.29600, 7.37546e-13,
1.29800, 7.31092e-13,
1.30000, 7.36586e-13,
1.30200, 7.29567e-13,
1.30400, 7.26690e-13,
1.30600, 7.34370e-13,
1.30800, 7.33400e-13,
1.31000, 7.32505e-13,
1.31200, 7.16767e-13,
1.31400, 7.32099e-13,
1.31600, 7.11538e-13,
1.31800, 7.28914e-13,
1.32000, 7.33202e-13,
1.32200, 7.29788e-13,
1.32400, 7.32798e-13,
1.32600, 7.28077e-13,
1.32800, 7.28023e-13,
1.33000, 7.15809e-13,
1.33200, 7.18618e-13,
1.33400, 7.32919e-13,
1.33600, 7.29491e-13,
1.33800, 7.27149e-13,
1.34000, 7.19411e-13,
1.34200, 7.27608e-13,
1.34400, 7.29677e-13,
1.34600, 7.26449e-13,
1.34800, 7.30018e-13,
1.35000, 7.20116e-13,
1.35200, 7.28502e-13,
1.35400, 7.27629e-13,
1.35600, 7.09621e-13,
1.35800, 7.25615e-13,
1.36000, 7.17495e-13,
1.36200, 7.24479e-13,
1.36400, 7.17437e-13,
1.36600, 7.24926e-13,
1.36800, 7.15838e-13,
1.37000, 7.14058e-13,
1.37200, 7.21997e-13,
1.37400, 7.20119e-13,
1.37600, 7.12581e-13,
1.37800, 7.22109e-13,
1.38000, 7.22719e-13,
1.38200, 7.16327e-13,
1.38400, 7.23254e-13,
1.38600, 7.16321e-13,
1.38800, 7.22771e-13,
1.39000, 7.17933e-13,
1.39200, 7.22498e-13,
1.39400, 7.19660e-13,
1.39600, 7.18153e-13,
1.39800, 7.23402e-13,
1.40000, 7.09560e-13,
1.40200, 7.14030e-13,
1.40400, 7.11328e-13,
1.40600, 7.12893e-13,
1.40800, 7.16188e-13,
1.41000, 7.17706e-13,
1.41200, 7.01446e-13,
1.41400, 7.05653e-13,
1.41600, 7.13299e-13,
1.41800, 7.14400e-13,
1.42000, 7.09634e-13,
1.42200, 6.94749e-13,
1.42400, 7.05456e-13,
1.42600, 6.95005e-13,
1.42800, 6.95410e-13,
1.43000, 7.00601e-13,
1.43200, 7.05738e-13,
1.43400, 7.05736e-13,
1.43600, 7.09899e-13,
1.43800, 7.13036e-13,
1.44000, 6.83420e-13,
1.44200, 6.97987e-13,
1.44400, 6.93385e-13,
1.44600, 7.08350e-13,
1.44800, 7.04322e-13,
1.45000, 7.00715e-13,
1.45200, 6.93039e-13,
1.45400, 6.94802e-13,
1.45600, 6.93813e-13,
1.45800, 7.07452e-13,
1.46000, 7.04405e-13,
1.46200, 6.98032e-13,
1.46400, 6.98972e-13,
1.46600, 6.96551e-13,
1.46800, 7.00362e-13,
1.47000, 6.97267e-13,
1.47200, 6.95181e-13,
1.47400, 6.86028e-13,
1.47600, 6.89137e-13,
1.47800, 6.85748e-13,
1.48000, 7.01278e-13,
1.48200, 6.97518e-13,
1.48400, 6.89143e-13,
1.48600, 6.97572e-13,
1.48800, 6.45399e-13,
1.49000, 7.01932e-13,
1.49200, 6.98762e-13,
1.49400, 7.05557e-13,
1.49600, 6.86136e-13,
1.49800, 6.97511e-13,
1.50000, 6.98205e-13,
1.50200, 6.71491e-13,
1.50400, 6.63009e-13,
1.50600, 6.58236e-13,
1.50800, 6.82318e-13,
1.51000, 7.02941e-13,
1.51200, 6.91403e-13,
1.51400, 6.87744e-13,
1.51600, 6.92874e-13,
1.51800, 6.96682e-13,
1.52000, 6.95945e-13,
1.52200, 6.78982e-13,
1.52400, 6.85717e-13,
1.52600, 6.97169e-13,
1.52800, 7.00251e-13,
1.53000, 6.72957e-13,
1.53200, 6.96530e-13,
1.53400, 6.85148e-13,
1.53600, 6.96404e-13,
1.53800, 6.90052e-13,
1.54000, 6.77828e-13,
1.54200, 6.86078e-13,
1.54400, 6.97590e-13,
1.54600, 6.98047e-13,
1.54800, 6.80357e-13,
1.55000, 6.80521e-13,
1.55200, 6.91451e-13,
1.55400, 6.71547e-13,
1.55600, 6.78804e-13,
1.55800, 6.83176e-13,
1.56000, 6.82619e-13,
1.56200, 6.88693e-13,
1.56400, 6.82395e-13,
1.56600, 6.74025e-13,
1.56800, 6.72402e-13,
1.57000, 6.77539e-13,
1.57200, 6.82788e-13,
1.57400, 6.69961e-13,
1.57600, 6.48635e-13,
1.57800, 6.42191e-13,
1.58000, 6.82538e-13,
1.58200, 6.59314e-13,
1.58400, 6.70182e-13,
1.58600, 6.73788e-13,
1.58800, 6.45933e-13,
1.59000, 6.16510e-13,
1.59200, 6.58976e-13,
1.59400, 6.81467e-13,
1.59600, 6.55355e-13,
1.59800, 6.83954e-13,
1.60000, 6.65312e-13,
1.60200, 6.67991e-13,
1.60400, 6.68640e-13,
1.60600, 6.66424e-13,
1.60800, 6.74217e-13,
1.61000, 6.43323e-13,
1.61200, 6.59802e-13,
1.61400, 6.68733e-13,
1.61600, 6.37973e-13,
1.61800, 6.71160e-13,
1.62000, 6.46681e-13,
1.62200, 6.61522e-13,
1.62400, 6.71745e-13,
1.62600, 6.70677e-13,
1.62800, 6.78164e-13,
1.63000, 6.68887e-13,
1.63200, 6.67216e-13,
1.63400, 6.60924e-13,
1.63600, 6.54317e-13,
1.63800, 6.36776e-13,
1.64000, 6.36171e-13,
1.64200, 6.36295e-13,
1.64400, 6.36740e-13,
1.64600, 6.55151e-13,
1.64800, 6.57653e-13,
1.65000, 6.45953e-13,
1.65200, 6.49024e-13,
1.65400, 6.56352e-13,
1.65600, 6.57130e-13,
1.65800, 6.58138e-13,
1.66000, 6.52600e-13,
1.66200, 6.56297e-13,
1.66400, 6.53466e-13,
1.66600, 6.44878e-13,
1.66800, 6.36959e-13,
1.67000, 6.54164e-13,
1.67200, 6.28981e-13,
1.67400, 6.44757e-13,
1.67600, 6.25197e-13,
1.67800, 6.36140e-13,
1.68000, 6.08532e-13,
1.68200, 5.99633e-13,
1.68400, 6.29073e-13,
1.68600, 6.32920e-13,
1.68800, 6.33499e-13,
1.69000, 6.29329e-13,
1.69200, 6.41040e-13,
1.69400, 6.38429e-13,
1.69600, 6.37710e-13,
1.69800, 6.36531e-13,
1.70000, 6.27541e-13,
1.70200, 6.24004e-13,
1.70400, 6.26401e-13,
1.70600, 6.28900e-13,
1.70800, 6.19666e-13,
1.71000, 6.25262e-13,
1.71200, 6.01648e-13,
1.71400, 6.30154e-13,
1.71600, 6.24652e-13,
1.71800, 6.23657e-13,
1.72000, 6.15724e-13,
1.72200, 6.16649e-13,
1.72400, 5.97561e-13,
1.72600, 6.11763e-13,
1.72800, 6.10277e-13,
1.73000, 6.15037e-13,
1.73200, 6.04261e-13,
1.73400, 5.75066e-13,
1.73600, 5.61680e-13,
1.73800, 5.82457e-13,
1.74000, 5.99830e-13,
1.74200, 5.95712e-13,
1.74400, 6.06618e-13,
1.74600, 5.95533e-13,
1.74800, 5.94371e-13,
1.75000, 5.97109e-13,
1.75200, 5.95446e-13,
1.75400, 6.03188e-13,
1.75600, 6.01572e-13,
1.75800, 6.01717e-13,
1.76000, 6.00923e-13,
1.76200, 5.85140e-13,
1.76400, 5.95785e-13,
1.76600, 6.00909e-13,
1.76800, 5.95708e-13,
1.77000, 5.92095e-13,
1.77200, 5.95963e-13,
1.77400, 5.97032e-13,
1.77600, 5.86732e-13,
1.77800, 5.83875e-13,
1.78000, 5.86063e-13,
1.78200, 5.85395e-13,
1.78400, 5.85874e-13,
1.78600, 5.83927e-13,
1.78800, 5.91676e-13,
1.79000, 5.87261e-13,
1.79200, 5.84003e-13,
1.79400, 5.80302e-13,
1.79600, 5.81186e-13,
1.79800, 5.81319e-13,
1.80000, 5.79641e-13,
1.80200, 5.70133e-13,
1.80400, 5.78704e-13,
1.80600, 5.76185e-13,
1.80800, 5.76595e-13,
1.81000, 5.76716e-13,
1.81200, 5.71652e-13,
1.81400, 5.66361e-13,
1.81600, 5.52284e-13,
1.81800, 5.13590e-13,
1.82000, 5.50944e-13,
1.82200, 5.47623e-13,
1.82400, 5.65915e-13,
1.82600, 5.66850e-13,
1.82800, 5.60324e-13,
1.83000, 5.66457e-13,
1.83200, 5.62532e-13,
1.83400, 5.60541e-13,
1.83600, 5.61536e-13,
1.83800, 5.55460e-13,
1.84000, 5.57606e-13,
1.84200, 5.47818e-13,
1.84400, 5.53545e-13,
1.84600, 5.58365e-13,
1.84800, 5.36435e-13,
1.85000, 5.58834e-13,
1.85200, 5.55610e-13,
1.85400, 5.52216e-13,
1.85600, 5.44047e-13,
1.85800, 5.51977e-13,
1.86000, 5.43388e-13,
1.86200, 5.40151e-13,
1.86400, 5.43588e-13,
1.86600, 5.35598e-13,
1.86800, 5.36669e-13,
1.87000, 5.42668e-13,
1.87200, 5.26123e-13,
1.87400, 5.12178e-13,
1.87600, 4.46945e-13,
1.87800, 5.23286e-13,
1.88000, 5.29141e-13,
1.88200, 5.37328e-13,
1.88400, 5.35457e-13,
1.88600, 5.25875e-13,
1.88800, 5.38056e-13,
1.89000, 5.24510e-13,
1.89200, 5.23547e-13,
1.89400, 5.18181e-13,
1.89600, 5.08546e-13,
1.89800, 5.27070e-13,
1.90000, 5.24053e-13,
1.90200, 5.36481e-13,
1.90400, 5.25945e-13,
1.90600, 5.29258e-13,
1.90800, 5.25797e-13,
1.91000, 5.30504e-13,
1.91200, 5.19116e-13,
1.91400, 5.29884e-13,
1.91600, 5.25962e-13,
1.91800, 5.24187e-13,
1.92000, 5.15858e-13,
1.92200, 5.26272e-13,
1.92400, 5.19278e-13,
1.92600, 5.23667e-13,
1.92800, 5.11027e-13,
1.93000, 5.10823e-13,
1.93200, 5.05342e-13,
1.93400, 5.14171e-13,
1.93600, 5.13959e-13,
1.93800, 5.16829e-13,
1.94000, 5.00300e-13,
1.94200, 5.07327e-13,
1.94400, 4.64664e-13,
1.94600, 4.59172e-13,
1.94800, 5.00963e-13,
1.95000, 5.01190e-13,
1.95200, 4.87167e-13,
1.95400, 5.12530e-13,
1.95600, 5.10517e-13,
1.95800, 5.10310e-13,
1.96000, 5.07188e-13,
1.96200, 5.03740e-13,
1.96400, 4.89288e-13,
1.96600, 4.98149e-13,
1.96800, 5.08054e-13,
1.97000, 5.04290e-13,
1.97200, 4.89168e-13,
1.97400, 4.98210e-13,
1.97600, 5.02646e-13,
1.97800, 4.88209e-13,
1.98000, 4.93516e-13,
1.98200, 4.94868e-13,
1.98400, 4.95495e-13,
1.98600, 4.81433e-13,
1.98800, 5.00503e-13,
1.99000, 4.99770e-13,
1.99200, 4.92971e-13,
1.99400, 4.84237e-13,
1.99600, 4.92170e-13,
1.99800, 4.87685e-13,
2.00000, 4.94915e-13,
2.00200, 4.90296e-13,
2.00400, 4.87558e-13,
2.00600, 4.86984e-13,
2.00800, 4.92660e-13,
2.01000, 4.91268e-13,
2.01200, 4.89777e-13,
2.01400, 4.87242e-13,
2.01600, 4.88962e-13,
2.01800, 4.88320e-13,
2.02000, 4.82207e-13,
2.02200, 4.73967e-13,
2.02400, 4.82121e-13,
2.02600, 4.83177e-13,
2.02800, 4.80140e-13,
2.03000, 4.69698e-13,
2.03200, 4.82197e-13,
2.03400, 4.75758e-13,
2.03600, 4.75972e-13,
2.03800, 4.63999e-13,
2.04000, 4.78615e-13,
2.04200, 4.75565e-13,
2.04400, 4.78015e-13,
2.04600, 4.77576e-13,
2.04800, 4.74283e-13,
2.05000, 4.74405e-13,
2.05200, 4.68979e-13,
2.05400, 4.71633e-13,
2.05600, 4.72506e-13,
2.05800, 4.70368e-13,
2.06000, 4.61923e-13,
2.06200, 4.71428e-13,
2.06400, 4.53109e-13,
2.06600, 4.69738e-13,
2.06800, 4.67044e-13,
2.07000, 4.58812e-13,
2.07200, 4.63702e-13,
2.07400, 4.66966e-13,
2.07600, 4.66063e-13,
2.07800, 4.65679e-13,
2.08000, 4.55948e-13,
2.08200, 4.62139e-13,
2.08400, 4.54139e-13,
2.08600, 4.61937e-13,
2.08800, 4.60898e-13,
2.09000, 4.61020e-13,
2.09200, 4.46627e-13,
2.09400, 4.51786e-13,
2.09600, 4.50793e-13,
2.09800, 4.54576e-13,
2.10000, 4.55058e-13,
2.10200, 4.54163e-13,
2.10400, 4.55588e-13,
2.10600, 4.56493e-13,
2.10800, 4.54663e-13,
2.11000, 4.46787e-13,
2.11200, 4.52948e-13,
2.11400, 4.51579e-13,
2.11600, 4.43895e-13,
2.11800, 4.49794e-13,
2.12000, 4.45070e-13,
2.12200, 4.47130e-13,
2.12400, 4.40436e-13,
2.12600, 4.43971e-13,
2.12800, 4.47359e-13,
2.13000, 4.45428e-13,
2.13200, 4.46850e-13,
2.13400, 4.45599e-13,
2.13600, 4.36500e-13,
2.13800, 4.44875e-13,
2.14000, 4.42963e-13,
2.14200, 4.42987e-13,
2.14400, 4.40377e-13,
2.14600, 4.37239e-13,
2.14800, 4.41022e-13,
2.15000, 4.40132e-13,
2.15200, 4.37460e-13,
2.15400, 4.37087e-13,
2.15600, 4.37303e-13,
2.15800, 4.36489e-13,
2.16000, 4.33322e-13,
2.16200, 4.29252e-13,
2.16400, 4.20974e-13,
2.16600, 3.78746e-13,
2.16800, 4.14123e-13,
2.17000, 4.25559e-13,
2.17200, 4.28047e-13,
2.17400, 4.29030e-13,
2.17600, 4.25064e-13,
2.17800, 4.18623e-13,
2.18000, 4.30053e-13,
2.18200, 4.21488e-13,
2.18400, 4.27992e-13,
2.18600, 4.28594e-13,
2.18800, 4.14800e-13,
2.19000, 4.23539e-13,
2.19200, 4.26734e-13,
2.19400, 4.25392e-13,
2.19600, 4.24615e-13,
2.19800, 4.25316e-13,
2.20000, 4.23993e-13,
2.20200, 4.22401e-13,
2.20400, 4.23343e-13,
2.20600, 4.01944e-13,
2.20800, 4.12401e-13,
2.21000, 4.20265e-13,
2.21200, 4.20399e-13,
2.21400, 4.19323e-13,
2.21600, 4.15336e-13,
2.21800, 4.15640e-13,
2.22000, 4.18398e-13,
2.22200, 4.16883e-13,
2.22400, 4.15840e-13,
2.22600, 4.07600e-13,
2.22800, 4.13185e-13,
2.23000, 4.13555e-13,
2.23200, 4.13935e-13,
2.23400, 4.14172e-13,
2.23600, 4.12628e-13,
2.23800, 4.04914e-13,
2.24000, 4.08110e-13,
2.24200, 4.08163e-13,
2.24400, 4.09332e-13,
2.24600, 4.09809e-13,
2.24800, 4.02980e-13,
2.25000, 4.07687e-13,
2.25200, 4.06088e-13,
2.25400, 3.99432e-13,
2.25600, 4.03032e-13,
2.25800, 4.05353e-13,
2.26000, 4.05240e-13,
2.26200, 3.96348e-13,
2.26400, 4.02010e-13,
2.26600, 3.99723e-13,
2.26800, 4.02059e-13,
2.27000, 4.01889e-13,
2.27200, 3.99818e-13,
2.27400, 3.98066e-13,
2.27600, 4.00711e-13,
2.27800, 4.00697e-13,
2.28000, 3.99013e-13,
2.28200, 3.82156e-13,
2.28400, 3.95939e-13,
2.28600, 3.94788e-13,
2.28800, 3.96657e-13,
2.29000, 3.96396e-13,
2.29200, 3.95144e-13,
2.29400, 3.81872e-13,
2.29600, 3.85539e-13,
2.29800, 3.85464e-13,
2.30000, 3.88166e-13,
2.30200, 3.89663e-13,
2.30400, 3.89708e-13,
2.30600, 3.89602e-13,
2.30800, 3.89369e-13,
2.31000, 3.87762e-13,
2.31200, 3.87987e-13,
2.31400, 3.88549e-13,
2.31600, 3.83863e-13,
2.31800, 3.80911e-13,
2.32000, 3.85783e-13,
2.32200, 3.78113e-13,
2.32400, 3.70837e-13,
2.32600, 3.71319e-13,
2.32800, 3.78927e-13,
2.33000, 3.80441e-13,
2.33200, 3.76653e-13,
2.33400, 3.79304e-13,
2.33600, 3.74665e-13,
2.33800, 3.73766e-13,
2.34000, 3.79136e-13,
2.34200, 3.79169e-13,
2.34400, 3.78885e-13,
2.34600, 3.77304e-13,
2.34800, 3.77992e-13,
2.35000, 3.78001e-13,
2.35200, 3.63148e-13,
2.35400, 3.62496e-13,
2.35600, 3.69338e-13,
2.35800, 3.64649e-13,
2.36000, 3.69969e-13,
2.36200, 3.71970e-13,
2.36400, 3.70799e-13,
2.36600, 3.70861e-13,
2.36800, 3.70827e-13,
2.37000, 3.66784e-13,
2.37200, 3.69144e-13,
2.37400, 3.56081e-13,
2.37600, 3.67741e-13,
2.37800, 3.70410e-13,
2.38000, 3.69127e-13,
2.38200, 3.65147e-13,
2.38400, 3.46836e-13,
2.38600, 3.47862e-13,
2.38800, 3.60854e-13,
2.39000, 3.62207e-13,
2.39200, 3.60859e-13,
2.39400, 3.61537e-13,
2.39600, 3.60933e-13,
2.39800, 3.63447e-13,
2.40000, 3.58613e-13,
2.40200, 3.62200e-13,
2.40400, 3.61443e-13,
2.40600, 3.60373e-13,
2.40800, 3.60272e-13,
2.41000, 3.60351e-13,
2.41200, 3.60728e-13,
2.41400, 3.45870e-13,
2.41600, 3.36563e-13,
2.41800, 3.54201e-13,
2.42000, 3.52558e-13,
2.42200, 3.54307e-13,
2.42400, 3.52192e-13,
2.42600, 3.52508e-13,
2.42800, 3.52698e-13,
2.43000, 3.54692e-13,
2.43200, 3.53261e-13,
2.43400, 3.50906e-13,
2.43600, 3.47889e-13,
2.43800, 3.48242e-13,
2.44000, 3.52491e-13,
2.44200, 3.51656e-13,
2.44400, 3.51277e-13,
2.44600, 3.43920e-13,
2.44800, 3.41251e-13,
2.45000, 3.40821e-13,
2.45200, 3.40944e-13,
2.45400, 3.45357e-13,
2.45600, 3.39709e-13,
2.45800, 3.36053e-13,
2.46000, 3.46337e-13,
2.46200, 3.45695e-13,
2.46400, 3.46047e-13,
2.46600, 3.44790e-13,
2.46800, 3.44971e-13,
2.47000, 3.34644e-13,
2.47200, 3.44618e-13,
2.47400, 3.38220e-13,
2.47600, 3.43582e-13,
2.47800, 3.41573e-13,
2.48000, 3.27900e-13,
2.48200, 3.28511e-13,
2.48400, 3.39852e-13,
2.48600, 3.27966e-13,
2.48800, 3.38553e-13,
2.49000, 3.39015e-13,
2.49200, 3.38429e-13,
2.49400, 3.35157e-13,
2.49600, 3.36748e-13,
2.49800, 3.35981e-13,
2.50000, 3.37610e-13
])
COLINA_ARRAY = COLINA_ARRAY.reshape(COLINA_ARRAY.size//2,2)
COLINA_WAVELENGTH_MICRON = COLINA_ARRAY[:,0]
COLINA_FLUX_PER_HZ = COLINA_ARRAY[:,1] * np.pi # Column is F, not pi*F
# w = wavelength in microns
# f = frequency in Hz
#
# flux_w dw = flux_f df
#
# flux_w = flux_f |df/dw|
#
# f = C/w
#
# |df/dw| = C/w^2
SPEED_OF_LIGHT_MICRON_PER_SEC = C * 1.e9
COLINA_FLUX_PER_MICRON = COLINA_FLUX_PER_HZ * (SPEED_OF_LIGHT_MICRON_PER_SEC /
COLINA_WAVELENGTH_MICRON**2)
# W/m^2/micron
# Returns flux density in W/m^2/micron given wavelength in microns
FLUX_DENSITY_COLINA = tab.Tabulation(COLINA_WAVELENGTH_MICRON,
COLINA_FLUX_PER_MICRON)
# From Bohlin, Dickinson, & Calzetti 2001, Astron. J., for HST STIS
# To read the file instead of using the hard-wired table...
# HST_SUN_REFERENCE = os.environ["HST_SUN_REFERENCE"]
# f = pyfits.open(HST_SUN_REFERENCE)
# wavelength_angstrom = f[1].data.WAVELENGTH # Angstrom
# flux_cgs_per_angstrom = f[1].data.FLUX # erg/s/cm**2/A
# f.close()
STIS_ARRAY = np.array([
1195.000, 8.8615008E-03,
1205.000, 6.9449998E-02,
1215.000, 0.6080999,
1225.000, 0.1437000,
1235.000, 6.2695001E-03,
1245.000, 4.0755002E-03,
1255.000, 3.3780001E-03,
1265.000, 4.1625001E-03,
1275.000, 2.4150000E-03,
1285.000, 1.9545001E-03,
1295.000, 4.6680002E-03,
1305.000, 1.3810000E-02,
1315.000, 4.4809999E-03,
1325.000, 4.6494999E-03,
1335.000, 1.9305000E-02,
1345.000, 5.0630001E-03,
1355.000, 4.1964999E-03,
1365.000, 3.4324999E-03,
1375.000, 3.2259999E-03,
1385.000, 4.4944999E-03,
1395.000, 8.8989995E-03,
1405.000, 6.8195001E-03,
1415.000, 4.6124998E-03,
1425.000, 5.1210001E-03,
1435.000, 5.5789999E-03,
1445.000, 5.5705002E-03,
1455.000, 6.0394998E-03,
1465.000, 7.6934998E-03,
1475.000, 9.1944998E-03,
1485.000, 8.8510001E-03,
1495.000, 8.6170007E-03,
1505.000, 9.3970001E-03,
1515.000, 1.0006500E-02,
1525.000, 1.2685000E-02,
1535.000, 1.4125000E-02,
1545.000, 2.2624999E-02,
1555.000, 2.3435000E-02,
1565.000, 1.9819999E-02,
1575.000, 1.8435000E-02,
1585.000, 1.7945001E-02,
1595.000, 1.8710000E-02,
1605.000, 2.0700000E-02,
1615.000, 2.3815000E-02,
1625.000, 2.7030000E-02,
1635.000, 3.0874999E-02,
1645.000, 3.2145001E-02,
1655.000, 4.6160001E-02,
1665.000, 4.2739999E-02,
1675.000, 4.1414998E-02,
1685.000, 4.8404999E-02,
1695.000, 6.2390000E-02,
1705.000, 7.0175000E-02,
1715.000, 7.3305003E-02,
1725.000, 7.7200003E-02,
1735.000, 8.1254996E-02,
1745.000, 9.6354999E-02,
1755.000, 0.1147000,
1765.000, 0.1299000,
1775.000, 0.1497000,
1785.000, 0.1646000,
1795.000, 0.1685000,
1805.000, 0.1960500,
1815.000, 0.2399000,
1825.000, 0.2396500,
1835.000, 0.2392000,
1845.000, 0.2193500,
1855.000, 0.2453500,
1865.000, 0.2870000,
1875.000, 0.3198000,
1885.000, 0.3494500,
1895.000, 0.3793500,
1905.000, 0.3992500,
1915.000, 0.4422000,
1925.000, 0.4377500,
1935.000, 0.4032000,
1945.000, 0.5531000,
1955.000, 0.5768500,
1965.000, 0.6370500,
1975.000, 0.6571500,
1985.000, 0.6555500,
1995.000, 0.7052000,
2005.649, 0.7756500,
2015.651, 0.8497500,
2025.653, 0.8731000,
2035.654, 0.9743000,
2045.656, 1.083500,
2055.658, 1.116000,
2065.660, 1.172000,
2075.662, 1.341000,
2085.664, 1.585000,
2095.666, 2.233000,
2105.668, 2.916500,
2115.670, 3.501000,
2125.672, 3.270500,
2135.674, 3.446500,
2145.676, 4.172500,
2155.678, 3.771500,
2165.680, 3.363500,
2175.682, 3.666500,
2185.684, 4.621500,
2195.686, 4.889000,
2205.688, 4.933000,
2215.690, 4.119500,
2225.692, 5.157001,
2235.694, 6.540500,
2245.697, 6.125500,
2255.698, 5.445000,
2265.701, 4.110000,
2275.703, 4.134000,
2285.705, 5.396500,
2295.707, 4.955000,
2305.709, 5.416000,
2315.712, 5.295500,
2325.714, 5.527500,
2335.716, 4.663500,
2345.718, 4.055500,
2355.720, 5.335000,
2365.723, 5.055000,
2375.725, 5.041000,
2385.727, 4.361000,
2395.729, 4.594000,
2405.732, 4.109500,
2415.734, 5.326500,
2425.736, 7.324500,
2435.739, 6.905500,
2445.741, 6.340000,
2455.743, 5.146500,
2465.746, 5.245000,
2475.748, 5.770500,
2485.750, 4.748500,
2495.752, 5.859000,
2505.755, 6.266500,
2515.757, 4.806500,
2525.760, 4.318500,
2535.762, 5.361500,
2545.764, 6.196000,
2555.767, 8.224000,
2565.769, 10.55000,
2575.771, 13.01000,
2585.774, 12.98000,
2595.776, 10.85000,
2605.778, 8.915500,
2615.781, 9.304500,
2625.783, 10.80000,
2635.785, 16.95500,
2645.788, 25.91000,
2655.790, 26.25500,
2665.793, 25.89500,
2675.795, 26.06000,
2685.797, 25.33500,
2695.800, 24.81000,
2705.802, 27.76000,
2715.805, 23.29500,
2725.807, 20.48500,
2735.810, 20.53500,
2745.812, 13.81500,
2755.814, 18.15500,
2765.817, 25.20500,
2775.819, 24.33000,
2785.822, 16.54000,
2795.824, 9.058500,
2805.826, 9.978500,
2815.829, 21.69000,
2825.832, 30.47500,
2835.834, 32.53000,
2845.836, 24.47500,
2855.839, 16.96000,
2865.841, 33.34000,
2875.844, 34.79000,
2885.846, 33.29000,
2895.849, 48.92000,
2905.851, 62.25500,
2915.854, 60.18000,
2925.856, 53.95000,
2935.859, 55.45000,
2945.861, 52.68000,
2955.864, 57.11500,
2965.866, 52.58000,
2975.868, 52.46500,
2985.871, 48.20000,
2995.874, 50.15500,
3005.876, 43.61500,
3015.878, 47.27500,
3025.881, 50.69000,
3035.883, 64.45500,
3045.886, 62.55000,
3055.888, 61.67500,
3065.891, 57.70500,
3075.893, 63.85000,
3085.896, 63.51000,
3095.898, 51.64500,
3105.901, 64.82000,
3115.904, 75.64000,
3125.906, 68.08000,
3135.908, 72.58000,
3145.911, 68.81500,
3155.914, 65.68500,
3165.916, 65.60001,
3175.918, 80.25000,
3185.921, 69.03500,
3195.924, 73.86000,
3205.926, 83.49000,
3215.929, 72.59500,
3225.931, 71.45000,
3235.934, 68.62500,
3245.936, 78.95999,
3255.939, 90.99000,
3265.941, 101.5500,
3275.944, 98.83000,
3285.947, 95.21001,
3295.949, 110.0000,
3305.951, 105.3500,
3315.954, 100.1350,
3325.957, 99.07000,
3335.959, 95.49000,
3345.962, 99.39000,
3355.964, 97.75500,
3365.967, 83.96001,
3375.969, 87.36000,
3385.972, 95.60001,
3395.974, 99.44500,
3405.977, 104.3000,
3415.979, 95.85000,
3425.982, 103.0500,
3435.985, 98.70500,
3445.987, 82.69000,
3455.990, 97.35000,
3465.992, 96.07500,
3475.995, 93.52000,
3485.998, 93.14999,
3496.000, 92.43001,
3506.003, 108.7500,
3516.005, 101.6500,
3526.008, 94.24001,
3536.010, 107.1500,
3546.013, 115.0500,
3556.016, 108.4000,
3566.018, 94.27000,
3576.021, 82.59500,
3586.023, 72.11000,
3596.026, 103.0500,
3606.028, 100.5000,
3616.031, 91.19499,
3626.034, 103.5500,
3636.036, 103.4000,
3646.039, 105.1000,
3656.041, 119.6500,
3666.044, 127.9000,
3676.047, 122.5500,
3686.049, 114.3500,
3696.052, 127.5000,
3706.054, 118.5500,
3716.057, 122.0000,
3726.060, 109.6500,
3736.062, 95.74500,
3746.065, 93.68500,
3756.067, 110.3000,
3766.070, 112.8000,
3776.073, 133.0500,
3786.075, 137.5500,
3796.078, 110.7500,
3806.080, 124.5500,
3816.083, 111.8500,
3826.085, 81.95000,
3836.088, 71.09500,
3846.091, 99.32000,
3856.093, 104.5500,
3866.096, 104.3000,
3876.098, 103.9000,
3886.101, 101.9500,
3896.104, 121.4500,
3906.106, 128.9500,
3916.109, 138.7500,
3926.111, 105.7000,
3936.114, 60.45500,
3946.117, 106.3000,
3956.119, 136.0500,
3966.122, 88.40500,
3976.125, 95.88000,
3986.127, 157.6500,
3996.130, 172.2500,
4006.132, 174.2000,
4016.135, 181.2000,
4026.137, 182.5500,
4036.140, 174.9500,
4046.143, 173.8000,
4056.146, 173.1500,
4066.148, 165.1000,
4076.151, 168.8000,
4086.153, 179.9500,
4096.156, 177.6000,
4106.159, 150.4000,
4116.161, 182.2000,
4126.164, 179.4000,
4136.167, 176.0000,
4146.169, 174.1000,
4156.171, 173.8000,
4166.174, 184.7000,
4176.177, 166.9000,
4186.180, 168.8000,
4196.182, 170.5000,
4206.185, 176.2000,
4216.188, 180.2000,
4226.190, 158.6000,
4236.192, 171.5000,
4246.195, 177.2000,
4256.198, 169.9000,
4266.201, 170.2000,
4276.203, 157.3000,
4286.206, 159.1000,
4296.208, 147.9000,
4306.211, 113.8000,
4316.213, 169.0000,
4326.216, 165.0000,
4336.219, 173.5000,
4346.222, 167.4000,
4356.224, 172.7000,
4366.227, 193.4000,
4376.229, 181.1000,
4386.232, 157.1000,
4396.235, 183.0000,
4406.237, 171.7000,
4416.240, 193.6000,
4426.243, 198.5000,
4436.245, 191.4000,
4446.248, 197.8000,
4456.250, 182.6000,
4466.253, 189.6000,
4476.256, 208.2000,
4486.258, 197.8000,
4496.261, 203.2000,
4506.264, 214.9000,
4516.266, 211.4000,
4526.269, 194.6000,
4536.271, 197.5000,
4546.274, 198.4000,
4556.277, 203.9000,
4566.279, 208.2000,
4576.282, 210.5000,
4586.285, 197.6000,
4596.288, 201.4000,
4606.290, 204.5000,
4616.293, 206.0000,
4626.295, 210.9000,
4636.298, 204.5000,
4646.301, 198.1000,
4656.303, 204.7000,
4666.306, 192.6000,
4676.309, 202.0000,
4686.311, 199.9000,
4696.314, 199.5000,
4706.316, 188.2000,
4716.319, 202.3000,
4726.322, 204.6000,
4736.325, 199.6000,
4746.327, 205.6000,
4756.330, 202.1000,
4766.333, 196.1000,
4776.335, 208.0000,
4786.338, 201.4000,
4796.340, 208.1000,
4806.343, 204.0000,
4816.346, 209.5000,
4826.348, 202.8000,
4836.351, 202.4000,
4846.354, 197.4000,
4856.356, 183.5000,
4866.359, 162.9000,
4876.362, 183.5000,
4886.364, 191.9000,
4896.367, 196.5000,
4906.370, 201.2000,
4916.372, 190.1000,
4926.375, 190.1000,
4936.377, 189.3000,
4946.380, 206.3000,
4956.383, 193.1000,
4966.386, 202.2000,
4976.388, 202.3000,
4986.391, 187.1000,
4996.394, 197.5000,
5006.396, 186.2000,
5016.399, 181.7000,
5026.401, 189.9000,
5036.404, 193.9000,
5046.407, 187.4000,
5056.410, 199.8000,
5066.412, 196.6000,
5076.415, 191.1000,
5086.417, 192.4000,
5096.420, 192.1000,
5106.423, 195.2000,
5116.425, 200.2000,
5126.428, 187.2000,
5136.431, 186.6000,
5146.434, 187.9000,
5156.436, 190.5000,
5166.438, 167.3000,
5176.441, 173.0000,
5186.444, 165.8000,
5196.447, 183.3000,
5206.449, 183.6000,
5216.452, 191.1000,
5226.455, 182.8000,
5236.458, 189.9000,
5246.460, 196.3000,
5256.463, 193.5000,
5266.465, 167.8000,
5276.468, 183.3000,
5286.471, 190.2000,
5296.473, 192.3000,
5306.476, 195.7000,
5316.479, 196.8000,
5326.481, 177.5000,
5336.484, 192.8000,
5346.487, 186.3000,
5356.489, 199.5000,
5366.492, 187.6000,
5376.495, 188.7000,
5386.498, 190.9000,
5396.500, 183.7000,
5406.502, 177.4000,
5416.505, 188.6000,
5426.508, 183.0000,
5436.511, 188.4000,
5446.513, 188.4000,
5456.516, 190.6000,
5466.519, 188.4000,
5476.521, 183.8000,
5486.524, 186.8000,
5496.527, 190.0000,
5506.529, 186.7000,
5516.532, 187.6000,
5526.535, 185.1000,
5536.538, 188.7000,
5546.540, 190.3000,
5556.543, 190.2000,
5566.545, 182.6000,
5576.548, 185.1000,
5586.551, 179.2000,
5596.553, 181.3000,
5606.556, 184.8000,
5616.559, 182.9000,
5626.562, 185.5000,
5636.564, 186.6000,
5646.567, 185.9000,
5656.569, 180.3000,
5666.572, 183.4000,
5676.575, 189.2000,
5686.578, 181.5000,
5696.580, 186.5000,
5706.583, 177.4000,
5716.585, 182.8000,
5726.588, 189.7000,
5736.591, 188.1000,
5746.594, 187.2000,
5756.596, 183.5000,
5766.599, 185.1000,
5776.602, 186.2000,
5786.604, 178.8000,
5796.607, 183.3000,
5806.609, 184.3000,
5816.612, 185.8000,
5826.615, 187.8000,
5836.618, 186.2000,
5846.620, 186.5000,
5856.623, 178.8000,
5866.625, 183.5000,
5876.628, 185.3000,
5886.631, 175.4000,
5896.634, 161.6000,
5906.636, 181.8000,
5916.639, 179.2000,
5926.642, 181.3000,
5936.645, 180.1000,
5946.647, 177.8000,
5956.650, 178.7000,
5966.652, 181.0000,
5976.655, 178.5000,
5986.658, 176.2000,
5996.661, 177.9000,
6006.663, 175.0000,
6016.666, 175.5000,
6026.668, 172.3000,
6036.671, 179.2000,
6046.674, 178.1000,
6056.677, 176.8000,
6066.679, 176.4000,
6076.682, 176.2000,
6086.685, 174.7000,
6096.688, 174.8000,
6106.690, 170.7000,
6116.693, 175.0000,
6126.695, 170.9000,
6136.698, 168.7000,
6146.701, 171.7000,
6156.704, 171.7000,
6166.706, 161.3000,
6176.709, 171.1000,
6186.711, 172.8000,
6196.714, 171.1000,
6206.717, 173.8000,
6216.719, 169.4000,
6226.722, 171.7000,
6236.725, 167.0000,
6246.728, 166.0000,
6256.730, 163.6000,
6266.733, 170.1000,
6276.735, 170.1000,
6286.738, 170.1000,
6296.741, 168.1000,
6311.745, 164.3000,
6331.750, 165.5000,
6351.756, 166.0000,
6371.761, 165.8000,
6391.767, 165.5000,
6411.772, 161.8000,
6431.777, 162.5000,
6451.783, 163.1000,
6471.788, 160.7000,
6491.793, 156.2000,
6511.799, 161.0000,
6531.804, 160.3000,
6551.810, 153.6000,
6571.815, 138.8000,
6591.820, 155.3000,
6611.826, 157.5000,
6631.831, 155.9000,
6651.836, 156.4000,
6671.842, 153.9000,
6691.847, 155.0000,
6711.853, 152.0000,
6731.858, 152.5000,
6751.863, 151.4000,
6771.869, 151.2000,
6791.874, 150.2000,
6811.879, 149.6000,
6831.885, 148.3000,
6851.890, 145.9000,
6871.896, 147.1000,
6891.901, 146.5000,
6911.906, 145.2000,
6931.912, 145.2000,
6951.917, 144.0000,
6971.922, 142.0000,
6991.928, 142.9000,
7011.933, 139.0000,
7031.938, 139.2000,
7051.944, 141.9000,
7071.949, 140.4000,
7091.955, 138.8000,
7111.960, 138.9000,
7131.965, 137.7000,
7151.971, 137.0000,
7171.977, 135.7000,
7191.982, 133.1000,
7211.987, 133.4000,
7231.993, 135.1000,
7251.998, 135.3000,
7272.003, 134.9000,
7292.009, 132.2000,
7312.014, 132.9000,
7332.020, 132.1000,
7352.025, 131.2000,
7372.030, 131.0000,
7392.036, 128.1000,
7412.041, 126.1000,
7432.046, 128.9000,
7452.052, 128.2000,
7472.057, 128.6000,
7492.063, 127.3000,
7512.068, 126.5000,
7532.073, 126.2000,
7552.079, 125.8000,
7572.084, 125.1000,
7592.089, 124.3000,
7612.095, 124.0000,
7632.101, 124.4000,
7652.106, 122.4000,
7672.111, 118.8000,
7692.117, 120.6000,
7712.122, 120.7000,
7732.127, 121.1000,
7752.133, 119.1000,
7772.138, 119.9000,
7792.144, 119.0000,
7812.149, 119.0000,
7832.154, 117.9000,
7852.160, 118.3000,
7872.165, 118.0000,
7892.170, 117.7000,
7912.176, 116.1000,
7932.181, 114.6000,
7952.187, 113.7000,
7972.192, 115.5000,
7992.198, 113.8000,
8012.203, 114.5000,
8032.208, 113.2000,
8052.214, 111.8000,
8072.219, 112.3000,
8092.225, 109.8000,
8112.230, 111.7000,
8132.235, 111.8000,
8152.241, 111.0000,
8172.246, 110.7000,
8192.252, 106.6000,
8212.257, 108.3000,
8232.263, 107.6000,
8252.268, 107.8000,
8272.273, 107.9000,
8292.278, 107.5000,
8312.284, 107.0000,
8332.289, 103.5000,
8352.295, 105.4000,
8372.300, 105.3000,
8392.306, 104.3000,
8412.312, 104.6000,
8432.316, 102.9000,
8452.322, 103.4000,
8472.327, 102.6000,
8492.333, 97.20000,
8512.338, 100.4000,
8532.344, 97.39999,
8552.349, 87.80000,
8572.354, 101.2000,
8592.359, 99.80000,
8612.365, 99.80000,
8632.370, 100.0000,
8652.376, 97.10000,
8672.382, 88.10001,
8692.387, 96.80000,
8702.390, 96.00000,
8722.395, 95.30000,
8742.400, 95.60000,
8762.405, 95.80000,
8782.411, 93.40000,
8802.417, 93.50000,
8822.422, 92.30000,
8842.428, 93.70000,
8862.433, 92.20000,
8882.438, 92.40000,
8902.443, 92.00000,
8922.449, 91.10000,
8942.454, 91.40000,
8962.460, 91.90000,
8982.465, 89.40000,
9002.471, 88.90000,
9022.476, 89.10000,
9042.481, 88.10000,
9062.487, 88.10000,
9082.492, 85.40000,
9102.498, 86.50000,
9122.503, 87.30000,
9142.509, 86.50000,
9162.514, 86.00000,
9182.520, 86.10000,
9202.524, 83.50000,
9222.530, 82.80000,
9242.535, 83.40000,
9262.541, 83.90000,
9282.547, 84.50000,
9302.552, 84.00000,
9322.558, 84.80000,
9342.563, 83.50000,
9362.568, 82.00000,
9382.573, 82.80000,
9402.579, 82.10000,
9422.584, 80.60000,
9442.590, 82.20000,
9462.595, 81.70000,
9482.601, 80.10000,
9502.606, 80.30000,
9522.611, 79.60000,
9542.617, 78.50000,
9562.622, 79.90000,
9582.628, 79.60000,
9602.633, 78.80000,
9622.639, 78.38935,
9642.644, 78.17941,
9662.649, 77.17265,
9682.654, 77.62556,
9702.660, 77.20778,
9722.666, 77.64973,
9742.671, 75.86095,
9762.677, 76.49548,
9782.682, 76.06364,
9802.688, 75.68953,
9822.692, 75.96313,
9842.698, 75.36926,
9862.703, 75.00387,
9882.709, 74.73538,
9902.715, 73.63179,
9922.720, 74.30486,
9942.726, 74.20058,
9962.730, 73.79315,
9982.736, 73.57117,
10002.74, 73.23773,
10022.75, 73.00728,
10042.75, 71.38766,
10062.76, 68.99416,
10082.76, 71.87687,
10102.77, 71.76549,
10122.77, 71.37299,
10142.78, 70.70595,
10162.79, 70.12578,
10182.79, 69.87607,
10202.80, 69.50996,
10222.80, 69.16306,
10242.81, 69.81610,
10262.81, 69.54941,
10282.82, 69.45961,
10302.82, 68.42381,
10322.83, 68.63899,
10342.83, 66.83368,
10362.84, 67.94361,
10382.84, 66.94673,
10402.85, 66.83846,
10422.86, 66.82946,
10442.86, 67.02232,
10462.87, 65.49031,
10482.87, 65.92071,
10502.88, 66.03579,
10522.88, 65.72432,
10542.89, 65.11803,
10562.89, 65.39734,
10582.90, 63.70441,
10602.90, 63.85567,
10622.91, 64.45622,
10642.92, 63.76211,
10662.92, 63.09301,
10682.93, 62.50493,
10702.93, 61.49424,
10722.94, 62.41352,
10742.94, 61.49778,
10762.95, 61.42463,
10782.95, 61.01913,
10802.96, 62.05273,
10822.96, 59.18496,
10842.97, 59.99493,
10862.97, 61.08934,
10882.98, 58.11977,
10902.99, 60.28704,
10922.99, 59.63124,
10943.00, 55.91789,
10963.00, 57.40627,
10983.01, 57.28666,
11003.01, 58.66410,
11023.02, 57.37829,
11043.02, 58.03478,
11063.03, 57.94183,
11083.03, 57.65679,
11103.04, 57.64046,
11123.04, 56.80031,
11143.05, 56.86176,
11163.06, 56.93167,
11183.06, 56.69742,
11203.07, 55.52505,
11223.07, 56.43431,
11243.08, 56.09680,
11263.08, 54.65652,
11283.09, 55.41208,
11303.09, 54.57809,
11323.10, 54.91460,
11343.10, 54.53310,
11363.11, 54.92915,
11383.12, 53.72720,
11403.12, 53.25866,
11423.13, 53.49507,
11443.13, 53.43018,
11463.14, 54.13747,
11483.14, 53.71170,
11503.15, 53.54111,
11523.15, 53.23780,
11543.16, 53.44590,
11563.17, 52.95668,
11583.17, 52.85518,
11603.18, 51.66225,
11623.18, 50.98537,
11643.19, 51.50771,
11663.19, 51.49890,
11683.20, 51.69552,
11703.20, 50.84182,
11723.21, 51.55614,
11743.21, 51.48318,
11763.22, 50.05137,
11783.22, 49.75911,
11803.23, 50.97070,
11823.24, 50.24544,
11843.24, 48.55664,
11863.25, 50.07570,
11883.25, 48.38157,
11903.26, 49.32462,
11923.26, 49.74626,
11943.27, 49.66248,
11963.27, 48.90552,
11983.28, 46.68336,
12003.28, 47.97074,
12023.29, 48.90864,
12043.29, 46.93597,
12063.30, 48.36390,
12083.31, 46.13804,
12103.31, 47.16687,
12123.32, 47.68590,
12143.32, 47.75809,
12163.33, 47.75607,
12183.33, 47.45395,
12203.34, 47.21936,
12223.34, 47.08216,
12243.35, 46.83695,
12263.35, 46.75682,
12283.36, 46.00920,
12303.37, 46.42212,
12323.37, 46.05795,
12343.38, 45.98770,
12363.38, 45.91085,
12383.39, 45.92178,
12403.39, 45.23552,
12423.40, 45.11307,
12443.40, 44.90644,
12463.41, 45.01888,
12483.41, 44.92472,
12503.42, 44.91299,
12523.42, 44.21993,
12543.43, 44.39183,
12563.44, 44.09855,
12583.44, 43.93694,
12603.45, 43.94984,
12623.45, 43.61958,
12643.46, 43.48156,
12663.46, 43.32548,
12683.47, 42.65122,
12703.47, 43.29649,
12723.48, 43.12050,
12743.49, 42.70112,
12763.49, 42.85569,
12783.50, 42.60778,
12803.50, 41.91643,
12823.51, 37.54504,
12843.51, 41.10097,
12863.52, 41.55607,
12883.52, 41.62078,
12903.53, 41.19453,
12923.53, 41.47491,
12943.54, 41.17268,
12963.54, 41.38526,
12983.55, 40.89680,
13003.56, 41.07743,
13023.56, 40.56111,
13043.57, 40.27733,
13063.57, 40.57840,
13083.58, 40.40097,
13103.58, 40.22859,
13123.59, 39.24432,
13143.59, 39.96184,
13163.60, 38.72156,
13183.60, 39.54686,
13203.61, 39.65905,
13223.62, 39.35505,
13243.62, 39.39804,
13263.63, 39.02625,
13283.63, 38.90590,
13303.64, 38.13821,
13323.64, 38.17302,
13343.65, 38.81599,
13363.65, 38.51889,
13383.66, 38.28053,
13403.66, 37.76020,
13423.67, 38.07667,
13443.67, 38.07141,
13463.68, 37.79042,
13483.69, 37.86345,
13503.69, 37.23928,
13523.70, 37.56161,
13543.70, 37.40583,
13563.71, 36.37254,
13583.71, 37.08289,
13603.72, 36.56012,
13623.72, 36.80767,
13643.73, 36.34309,
13663.74, 36.61498,
13683.74, 36.05030,
13703.75, 35.85577,
13723.75, 36.14877,
13743.76, 35.94990,
13763.76, 35.47024,
13783.77, 35.84025,
13803.77, 35.76661,
13823.78, 35.34777,
13843.78, 35.58649,
13863.79, 35.14372,
13883.79, 35.35804,
13903.80, 35.02038,
13923.81, 35.14186,
13943.81, 34.90345,
13963.82, 34.73063,
13983.82, 34.88447,
14003.83, 34.11927,
14023.83, 34.23630,
14043.84, 34.00965,
14063.84, 33.98758,
14083.85, 34.04776,
14103.85, 34.02321,
14123.86, 33.15823,
14143.87, 33.26282,
14163.87, 33.52831,
14183.88, 33.48542,
14203.88, 33.16840,
14223.89, 32.38137,
14243.89, 32.78812,
14263.90, 32.21186,
14283.90, 32.14039,
14303.91, 32.28980,
14323.91, 32.43576,
14343.92, 32.34525,
14363.92, 32.44547,
14383.93, 32.49826,
14403.94, 31.06198,
14423.94, 31.63613,
14443.95, 31.34056,
14463.95, 31.92846,
14483.96, 31.65925,
14503.96, 31.41031,
14523.97, 30.98067,
14543.97, 30.97409,
14563.98, 30.84512,
14583.99, 31.36523,
14603.99, 31.14463,
14624.00, 30.77849,
14644.00, 30.73577,
14664.01, 30.54582,
14684.01, 30.62929,
14704.02, 30.41102,
14724.02, 30.23772,
14744.03, 29.75865,
14764.03, 29.81258,
14784.04, 29.58573,
14804.04, 30.17402,
14824.05, 29.93131,
14844.06, 29.49228,
14864.06, 29.77270,
14884.07, 27.47194,
14904.07, 29.79812,
14924.08, 29.58411,
14944.08, 29.79185,
14964.09, 28.89441,
14984.09, 29.29504,
15004.10, 29.24605,
15024.10, 28.05221,
15044.11, 27.62426,
15064.12, 27.35258,
15084.12, 28.27812,
15104.13, 29.05570,
15124.13, 28.50323,
15144.14, 28.27753,
15164.14, 28.41337,
15184.15, 28.49429,
15204.15, 28.38928,
15224.16, 27.62458,
15244.16, 27.82543,
15264.17, 28.21600,
15284.18, 28.26660,
15304.18, 27.09387,
15324.19, 27.96978,
15344.19, 27.44102,
15364.20, 27.81924,
15384.20, 27.49386,
15404.21, 26.93670,
15424.21, 27.19389,
15444.22, 27.57859,
15464.22, 27.52530,
15484.23, 26.75849,
15504.24, 26.69592,
15524.24, 27.05480,
15544.25, 26.20844,
15564.25, 26.42358,
15584.26, 26.52553,
15604.26, 26.43602,
15624.27, 26.60296,
15644.27, 26.29231,
15664.28, 25.90353,
15684.28, 25.77528,
15704.29, 25.90608,
15724.29, 26.04037,
15744.30, 25.48631,
15764.31, 24.61243,
15784.31, 24.30618,
15804.32, 25.76792,
15824.32, 24.82823,
15844.33, 25.17383,
15864.33, 25.24549,
15884.34, 24.14087,
15904.34, 22.98330,
15924.35, 24.50474,
15944.35, 25.27755,
15964.36, 24.24809,
15984.37, 25.24293,
16004.37, 24.49355,
16024.38, 24.53082,
16044.38, 24.49345,
16064.39, 24.35154,
16084.39, 24.57504,
16104.40, 23.39075,
16124.40, 23.93041,
16144.41, 24.19427,
16164.42, 23.02428,
16184.42, 24.16216,
16204.43, 23.22344,
16224.43, 23.69786,
16244.44, 24.00487,
16264.44, 23.90775,
16284.45, 24.11529,
16304.45, 23.72708,
16324.46, 23.60983,
16344.46, 23.32997,
16364.47, 23.04031,
16384.47, 22.36790,
16404.48, 22.29219,
16424.49, 22.24225,
16444.49, 22.20370,
16464.50, 22.79022,
16484.50, 22.82176,
16504.51, 22.36143,
16524.51, 22.41337,
16544.52, 22.61166,
16564.52, 22.58381,
16584.53, 22.56391,
16604.54, 22.32018,
16624.54, 22.39261,
16644.54, 22.24247,
16664.55, 21.89748,
16684.56, 21.57675,
16704.56, 22.10652,
16724.57, 21.20468,
16744.57, 21.68461,
16764.58, 20.97663,
16784.58, 21.29293,
16804.59, 20.32038,
16824.59, 19.97563,
16844.60, 20.90660,
16864.61, 20.98457,
16884.61, 20.95403,
16904.62, 20.76685,
16924.62, 21.10332,
16944.63, 20.96777,
16964.63, 20.89481,
16984.64, 20.80707,
17004.64, 20.46495,
17024.65, 20.30181,
17044.65, 20.33200,
17064.66, 20.36527,
17084.66, 20.01927,
17104.67, 20.15284,
17124.68, 19.34646,
17144.68, 20.21582,
17164.69, 19.99265,
17184.69, 19.91436,
17204.70, 19.61534,
17224.70, 19.59920,
17244.71, 18.94848,
17264.71, 19.35390,
17284.72, 19.26221,
17304.72, 19.36758,
17324.73, 18.98432,
17344.74, 18.02546,
17364.74, 17.56532,
17384.75, 18.17318,
17404.75, 18.67224,
17424.76, 18.50149,
17444.76, 18.79701,
17464.77, 18.41127,
17484.77, 18.33333,
17504.78, 18.37572,
17524.79, 18.28273,
17544.79, 18.47822,
17564.79, 18.38675,
17584.80, 18.34938,
17604.81, 18.28353,
17624.81, 17.76292,
17644.82, 18.04510,
17664.82, 18.15909,
17684.83, 17.96122,
17704.83, 17.81196,
17724.84, 17.88788,
17744.84, 17.87955,
17764.85, 17.53155,
17784.86, 17.40696,
17804.86, 17.43295,
17824.87, 17.37401,
17844.87, 17.34926,
17864.88, 17.25292,
17884.88, 17.44278,
17904.89, 17.27395,
17924.89, 17.13979,
17944.90, 16.99322,
17964.90, 16.98123,
17984.91, 16.94735,
18004.92, 16.86091,
18024.92, 16.54754,
18044.93, 16.75907,
18064.93, 16.64920,
18084.94, 16.62420,
18104.94, 16.59097,
18124.95, 16.40900,
18144.95, 16.22128,
18164.96, 15.78327,
18184.96, 14.64521,
18204.97, 15.67585,
18224.97, 15.54719,
18244.98, 16.03128,
18264.99, 16.02261,
18284.99, 15.80350,
18305.00, 15.94158,
18325.00, 15.79657,
18345.01, 15.70635,
18365.01, 15.69997,
18385.02, 15.49632,
18405.02, 15.52237,
18425.03, 15.21681,
18445.04, 15.34255,
18465.04, 15.44265,
18485.05, 14.80401,
18505.05, 15.38884,
18525.06, 15.26703,
18545.06, 15.14106,
18565.07, 14.88495,
18585.07, 15.06940,
18605.08, 14.80303,
18625.08, 14.68327,
18645.09, 14.74499,
18665.10, 14.49713,
18685.10, 14.49505,
18705.11, 14.62574,
18725.11, 14.14955,
18745.12, 13.74512,
18765.12, 11.96891,
18785.13, 13.98346,
18805.13, 14.10985,
18825.14, 14.29771,
18845.14, 14.21771,
18865.15, 13.93368,
18885.15, 14.22623,
18905.16, 13.83874,
18925.17, 13.78415,
18945.17, 13.61408,
18965.18, 13.33276,
18985.18, 13.78930,
19005.19, 13.68154,
19025.19, 13.97655,
19045.20, 13.67330,
19065.20, 13.73057,
19085.21, 13.61219,
19105.21, 13.70531,
19125.22, 13.38305,
19145.23, 13.63212,
19165.23, 13.50299,
19185.24, 13.42937,
19205.24, 13.18848,
19225.25, 13.42672,
19245.25, 13.22075,
19265.26, 13.30484,
19285.26, 12.95677,
19305.27, 12.92476,
19325.28, 12.75961,
19345.28, 12.95572,
19365.29, 12.92362,
19385.29, 12.96899,
19405.30, 12.52835,
19425.30, 12.67815,
19445.31, 11.58812,
19465.31, 11.42763,
19485.32, 12.44211,
19505.32, 12.42224,
19525.33, 12.04993,
19545.33, 12.65135,
19565.34, 12.57589,
19585.35, 12.54513,
19605.35, 12.44296,
19625.36, 12.33317,
19645.36, 11.95497,
19665.37, 12.14672,
19685.37, 12.36308,
19705.38, 12.24656,
19725.38, 11.85526,
19745.39, 12.04993,
19765.39, 12.13262,
19785.40, 11.76034,
19805.40, 11.86418,
19825.41, 11.87267,
19845.42, 11.86376,
19865.42, 11.50386,
19885.43, 11.93550,
19905.43, 11.89407,
19925.44, 11.70873,
19945.44, 11.47821,
19965.45, 11.64289,
19985.46, 11.51370,
20005.46, 11.66103,
20025.46, 11.52914,
20045.47, 11.44188,
20065.48, 11.40564,
20085.48, 11.51559,
20105.49, 11.46023,
20125.49, 11.40274,
20145.50, 11.32120,
20165.50, 11.33863,
20185.51, 11.30131,
20205.51, 11.13775,
20225.52, 10.92579,
20245.53, 11.09180,
20265.53, 11.09415,
20285.54, 11.00269,
20305.54, 10.74220,
20325.55, 11.00636,
20345.55, 10.83804,
20365.56, 10.82162,
20385.56, 10.52871,
20405.57, 10.83909,
20425.57, 10.74892,
20445.58, 10.78317,
20465.58, 10.75222,
20485.59, 10.65722,
20505.60, 10.63919,
20525.60, 10.49701,
20545.61, 10.53586,
20565.61, 10.53484,
20585.62, 10.46679,
20605.62, 10.25893,
20625.63, 10.44973,
20645.63, 10.02420,
20665.64, 10.37199,
20685.64, 10.29256,
20705.65, 10.09162,
20725.66, 10.17950,
20745.66, 10.23139,
20765.67, 10.19194,
20785.67, 10.16394,
20805.68, 9.932425,
20825.68, 10.04796,
20845.69, 9.855086,
20865.69, 10.00508,
20885.70, 9.963475,
20905.71, 9.947034,
20925.71, 9.618074,
20945.71, 9.710592,
20965.72, 9.670780,
20985.73, 9.733356,
21005.73, 9.725116,
21025.74, 9.687522,
21045.74, 9.699467,
21065.75, 9.700283,
21085.75, 9.643069,
21105.76, 9.458072,
21125.76, 9.570342,
21145.77, 9.523367,
21165.78, 9.343638,
21185.78, 9.449934,
21205.79, 9.333052,
21225.79, 9.358583,
21245.80, 9.201120,
21265.80, 9.257522,
21285.81, 9.310640,
21305.81, 9.253057,
21325.82, 9.265182,
21345.82, 9.221927,
21365.83, 9.016720,
21385.84, 9.172522,
21405.84, 9.116041,
21425.85, 9.099525,
21445.85, 9.029051,
21465.86, 8.947999,
21485.86, 9.008626,
21505.87, 8.973722,
21525.87, 8.902678,
21545.88, 8.878581,
21565.88, 8.866479,
21585.89, 8.833586,
21605.89, 8.753266,
21625.90, 8.655013,
21645.91, 8.472415,
21665.91, 7.608472,
21685.92, 8.303815,
21705.92, 8.517407,
21725.93, 8.551423,
21745.93, 8.555308,
21765.94, 8.460640,
21785.94, 8.317131,
21805.95, 8.528557,
21825.96, 8.343378,
21845.96, 8.456617,
21865.97, 8.453033,
21885.97, 8.166020,
21905.98, 8.322844,
21925.98, 8.370333,
21945.99, 8.328810,
21965.99, 8.298452,
21986.00, 8.297029,
22006.00, 8.256200,
22026.01, 8.210254,
22046.02, 8.213644,
22066.02, 7.784324,
22086.03, 7.972390,
22106.03, 8.109716,
22126.04, 8.097629,
22146.04, 8.062320,
22166.05, 7.971262,
22186.05, 7.962715,
22206.06, 8.001108,
22226.06, 7.957787,
22246.07, 7.923604,
22266.07, 7.752663,
22286.08, 7.844775,
22306.09, 7.837720,
22326.09, 7.830868,
22346.10, 7.821339,
22366.10, 7.778253,
22386.11, 7.619207,
22406.11, 7.665628,
22426.12, 7.652956,
22446.12, 7.661192,
22466.13, 7.656466,
22486.13, 7.515487,
22506.14, 7.589761,
22526.15, 7.546570,
22546.15, 7.409714,
22566.16, 7.463243,
22586.16, 7.492931,
22606.17, 7.477594,
22626.17, 7.300589,
22646.18, 7.391807,
22666.18, 7.336788,
22686.19, 7.366651,
22706.20, 7.350574,
22726.20, 7.299822,
22746.21, 7.255056,
22766.21, 7.290438,
22786.22, 7.277381,
22806.22, 7.234087,
22826.23, 6.916335,
22846.23, 7.153233,
22866.24, 7.119960,
22886.24, 7.141175,
22906.25, 7.124009,
22926.25, 7.089127,
22946.26, 6.839079,
22966.27, 6.892733,
22986.27, 6.879397,
23006.28, 6.915575,
23026.28, 6.930194,
23046.29, 6.918959,
23066.29, 6.905087,
23086.30, 6.889012,
23106.30, 6.848699,
23126.31, 6.840815,
23146.31, 6.838899,
23166.32, 6.744756,
23186.33, 6.681328,
23206.33, 6.755139,
23226.34, 6.609423,
23246.34, 6.471088,
23266.35, 6.468368,
23286.35, 6.589565,
23306.36, 6.604533,
23326.36, 6.527568,
23346.37, 6.562235,
23366.38, 6.470885,
23386.38, 6.444325,
23406.38, 6.525742,
23426.39, 6.515162,
23446.40, 6.499184,
23466.40, 6.461030,
23486.41, 6.461789,
23506.41, 6.450948,
23526.42, 6.186943,
23546.42, 6.165344,
23566.43, 6.271040,
23586.43, 6.180926,
23606.44, 6.260482,
23626.45, 6.283697,
23646.45, 6.253309,
23666.46, 6.243799,
23686.46, 6.232675,
23706.47, 6.154332,
23726.47, 6.183488,
23746.48, 5.954613,
23766.48, 6.139264,
23786.49, 6.173412,
23806.49, 6.141697,
23826.50, 6.065282,
23846.51, 5.751458,
23866.51, 5.758802,
23886.52, 5.963880,
23906.52, 5.976233,
23926.53, 5.944046,
23946.53, 5.945261,
23966.54, 5.925430,
23986.54, 5.956750,
24006.55, 5.867729,
24026.55, 5.916553,
24046.56, 5.894377,
24066.56, 5.867162,
24086.57, 5.855767,
24106.58, 5.847337,
24126.58, 5.843753,
24146.59, 5.593778,
24166.59, 5.434240,
24186.60, 5.709578,
24206.60, 5.673697,
24226.61, 5.692431,
24246.61, 5.649125,
24266.62, 5.644866,
24286.63, 5.638613,
24306.63, 5.661162,
24326.64, 5.629043,
24346.64, 5.582340,
24366.65, 5.525257,
24386.65, 5.521800,
24406.66, 5.580003,
24426.66, 5.557671,
24446.67, 5.542601,
24466.67, 5.417651,
24486.68, 5.366827,
24506.69, 5.351323,
24526.69, 5.344524,
24546.70, 5.404882,
24566.70, 5.307821,
24586.71, 5.242172,
24606.71, 5.393799,
24626.72, 5.375060,
24646.72, 5.371800,
24666.73, 5.343614,
24686.73, 5.337766,
24706.74, 5.169591,
24726.74, 5.315050,
24746.75, 5.207941,
24766.76, 5.281967,
24786.76, 5.242604,
24806.77, 5.024642,
24826.77, 5.025889,
24846.78, 5.191024,
24866.78, 5.001424,
24886.79, 5.154565,
24906.79, 5.153316,
24926.80, 5.136147,
24946.80, 5.078345,
24966.81, 5.094283,
24986.82, 5.074543,
25006.82, 5.090985,
25056.83, 5.053445,
25156.86, 4.941578,
25256.89, 4.894886,
25356.92, 4.822793,
25456.94, 4.742045,
25556.97, 4.667750,
25657.00, 4.591386,
25757.03, 4.565696,
25857.05, 4.443347,
25957.08, 4.423011,
26057.11, 4.367141,
26157.13, 4.284614,
26257.16, 4.086294,
26357.19, 4.168105,
26457.22, 4.105946,
26557.24, 4.057827,
26657.27, 3.974363,
26757.30, 3.944365,
26857.32, 3.913914,
26957.35, 3.858138
])
STIS_ARRAY = STIS_ARRAY.reshape(STIS_ARRAY.size//2,2)
STIS_WAVELENGTH_MICRON = STIS_ARRAY[:,0] * 1.e-4
STIS_FLUX_PER_MICRON = STIS_ARRAY[:,1] * 1.e-7 * 1.e4 * 1.e4 # W/m^2/micron
# Returns flux density in W/m^2/micron given wavelength in microns
FLUX_DENSITY_STIS = tab.Tabulation(STIS_WAVELENGTH_MICRON,
STIS_FLUX_PER_MICRON)
# Combined function
SOLAR_FLUX_MODELS = {"COLINA": FLUX_DENSITY_COLINA,
"STIS": FLUX_DENSITY_STIS}
def flux_density(model="STIS"):
"""Returns the solar flux density at 1 AU in units of Watts per square meter
per micron of wavelength.
Input:
model name of the model, either "STIS" or "COLINA". "STIS" is
the default.
Return: a Tabulation of flux density (W/m^2/micron) vs.
wavelength (micron) at 1 AU.
"""
return SOLAR_FLUX_MODELS[model]
################################################################################
# Means over filter bandpasses
################################################################################
def bandpass_flux_density(bandpass, sun_range=1., model="STIS"):
"""Returns the solar flux density averaged over a filter bandpass.
Input:
bandpass a Tabulation of a filter bandpass.
sun_range the distance from the Sun to the target body in units of
AU.
model the ame of the model for solar flux density, either
"STIS" or "COLINA"; "STIS" is the default.
Return: the mean solar flux density within the filter bandpass,
in units of W/m^2/micron.
"""
# Multiply the bandpass and the solar spectrum
product = bandpass * flux_density(model)
# Resample the bandpass at the same wavelengths for a more reliable
# normalization
bandpass = bandpass.resample(product.x)
# Return the ratio of integrals
return product.integral() / (bandpass.integral() * sun_range**2)
def mean_flux_density(center, width, sun_range=1., model="STIS"):
"""Returns the solar flux density averaged over the bandpass of a "boxcar"
filter, given its center and full width.
Input:
center the center of the bandpass in microns.
width the full width of the bandpass in microns.
sun_range the distance from the Sun to the target body in units of
AU.
model the ame of the model for solar flux density, either
"STIS" or "COLINA"; "STIS" is the default.
Return: the mean solar flux density within the filter bandpass,
in units of W/m^2/micron.
"""
# Create a boxcar filter Tabulation
bandpass = tab.Tabulation((center - width/2., center + width/2.), (1.,1.))
# Return the mean over the filter
return bandpass_flux_density(bandpass, sun_range, model)
def bandpass_f(bandpass, sun_range=1., model="STIS"):
"""Returns the solar F averaged over a filter bandpass. F is defined such
that pi*F is the solar flux density.
Input:
bandpass a tuple (wavelength,throughput) of 1-D arrays,
describing the throughput of a filter as a function
of wavelength in microns.
sun_range the distance from the Sun to the target body in units of
AU.
model the ame of the model for solar flux density, either
"STIS" or "COLINA"; "STIS" is the default.
Return: the mean solar flux density within the filter bandpass,
in units of W/m^2/micron.
"""
return bandpass_flux_density(bandpass, sun_range, model) / np.pi
def mean_f(center, width, sun_range=1., model="STIS"):
"""Returns the solar F averaged over the bandpass of a "boxcar" filter,
given its center and full width. F is defined such that pi*F is the solar
flux density.
Input:
center the center of the bandpass in microns.
width the full width of the bandpass in microns.
sun_range the distance from the Sun to the target body in units of
AU.
model the ame of the model for solar flux density, either
"STIS" or "COLINA"; "STIS" is the default.
Return: the mean solar flux density within the filter bandpass,
in units of W/m^2/micron.
"""
return mean_flux_density(center, width, sun_range, model) / np.pi
################################################################################