This tutorial notebook takes you through the process of loading an Elf Owl model using PICASO and then changing the chemical abundance of certain gases within an Elf Owl model and generating new spectrum using this changed chemistry. Please note that this essentially means that you are loosing the self-consistency of your model in this process. If you are not sure if it is "safe" to manipulate a certain gas abundance within the Elf Owl grid, please contact the Sagnick Mukherjee at samukher@ucsc.edu.
import numpy as np
import picaso.justdoit as jdi
import matplotlib.pyplot as plt
import xarray
teff= '450.0' #K
grav='31.0' # m/s/s
mh='-1.0'
co='0.5' # relative to solar
logkzz='9.0'
ds = xarray.load_dataset(f'/Users/sagnickmukherjee/Downloads/output_425.0_475.0/spectra_logzz_{logkzz}_teff_{teff}_grav_{grav}_mh_{mh}_co_{co}.nc')
Inspect the contents of the Elf Owl model Xarray file
ds
<xarray.Dataset> Dimensions: (wavelength: 193132, pressure: 91) Coordinates: * wavelength (wavelength) float64 15.0 15.0 15.0 15.0 ... 0.6 0.6 0.6 0.6 * pressure (pressure) float64 0.0001778 0.0002108 ... 668.0 791.9 Data variables: (12/40) flux (wavelength) float64 3.827e+08 3.806e+08 ... 7.606e+03 temperature (pressure) float64 101.9 103.2 104.8 ... 2.997e+03 3.111e+03 e- (pressure) float64 4.5e-38 4.5e-38 ... 1.737e-08 2.236e-08 H2 (pressure) float64 0.8377 0.8377 0.8377 ... 0.8322 0.8307 H (pressure) float64 4.5e-38 4.5e-38 ... 0.005537 0.007109 H+ (pressure) float64 4.5e-38 4.5e-38 ... 1.069e-18 6.708e-18 ... ... OCS (pressure) float64 4.5e-38 4.5e-38 ... 8.072e-13 7.991e-13 Li (pressure) float64 4.5e-38 4.5e-38 ... 2.044e-10 1.995e-10 LiOH (pressure) float64 4.5e-38 4.5e-38 ... 1.938e-12 1.711e-12 LiH (pressure) float64 4.5e-38 4.5e-38 ... 1.46e-10 1.506e-10 LiCl (pressure) float64 4.5e-38 4.5e-38 ... 2.163e-14 1.567e-14 LiF (pressure) float64 4.5e-38 4.5e-38 ... 8.85e-16 6.762e-16 Attributes: author: Batalha,Mukherjee contact: natasha.e.batalha@nasa.gov code: {"spectra": "PICASO", "chemistry": "visccher"} doi: ELF OWL Paper planet_params: {"logg": {"value": 31.0, "unit": "m / s2"}, "teff": {"val...
array([14.999761, 14.999511, 14.999261, ..., 0.600022, 0.600012, 0.600002])
array([1.778000e-04, 2.107689e-04, 2.498512e-04, 2.961803e-04, 3.511002e-04, 4.162036e-04, 4.933790e-04, 5.848648e-04, 6.933146e-04, 8.218739e-04, 9.742715e-04, 1.154928e-03, 1.369083e-03, 1.622947e-03, 1.923886e-03, 2.280626e-03, 2.703516e-03, 3.204821e-03, 3.799081e-03, 4.503533e-03, 5.338610e-03, 6.328533e-03, 7.502014e-03, 8.893090e-03, 1.054211e-02, 1.249690e-02, 1.481416e-02, 1.756111e-02, 2.081742e-02, 2.467753e-02, 2.925341e-02, 3.467778e-02, 4.110798e-02, 4.873051e-02, 5.776646e-02, 6.847792e-02, 8.117558e-02, 9.622772e-02, 1.140709e-01, 1.352228e-01, 1.602967e-01, 1.900201e-01, 2.252549e-01, 2.670233e-01, 3.165366e-01, 3.752311e-01, 4.448091e-01, 5.272887e-01, 6.250622e-01, 7.409657e-01, 8.783607e-01, 1.041233e+00, 1.234305e+00, 1.463179e+00, 1.734492e+00, 2.056113e+00, 2.437372e+00, 2.889327e+00, 3.425086e+00, 4.060190e+00, 4.813058e+00, 5.705530e+00, 6.763489e+00, 8.017623e+00, 9.504306e+00, 1.126666e+01, 1.335581e+01, 1.583233e+01, 1.876807e+01, 2.224818e+01, 2.637360e+01, 3.126397e+01, 3.706116e+01, 4.393330e+01, 5.207972e+01, 6.173671e+01, 7.318436e+01, 8.675472e+01, 1.028414e+02, 1.219110e+02, 1.445166e+02, 1.713138e+02, 2.030800e+02, 2.407366e+02, 2.853756e+02, 3.382920e+02, 4.010205e+02, 4.753805e+02, 5.635289e+02, 6.680224e+02, 7.918918e+02])
array([3.82700593e+08, 3.80643193e+08, 3.77353824e+08, ..., 7.66920446e+03, 7.65829284e+03, 7.60566013e+03])
array([ 101.92793327, 103.22061407, 104.8269031 , 106.80453662, 109.12564562, 111.71266816, 114.82928818, 118.61362689, 123.02148031, 127.9787556 , 133.78440036, 139.72057273, 145.46672701, 151.41083234, 157.13760399, 162.56996184, 168.24896832, 174.66178135, 181.51860015, 188.18340151, 194.26992778, 200.01912458, 205.33063221, 210.16663944, 215.12856072, 220.27598031, 225.36318537, 230.47160487, 235.68690948, 241.05093353, 246.75509577, 252.82843695, 258.90160717, 264.88524095, 270.89227767, 277.0856875 , 283.59724416, 290.55688125, 298.09520229, 306.03047196, 314.49952921, 323.8696078 , 334.53223856, 346.85374574, 361.12682171, 378.14272805, 397.89209686, 418.64867366, 440.47336763, 463.42050086, 487.5499181 , 512.92301377, 539.60249128, 567.65395505, 597.14564853, 628.14368903, 660.72367578, 694.94789407, 730.89523178, 768.62995574, 808.22508174, 849.75149355, 893.26537156, 938.8516339 , 986.54010562, 1036.41642116, 1088.51683255, 1142.90498055, 1199.63745324, 1258.7499 , 1320.32090573, 1384.36646801, 1450.96011165, 1520.14664692, 1591.97703353, 1666.51757252, 1743.78634266, 1823.85492284, 1906.73240823, 1992.44574571, 2081.07275557, 2172.45130722, 2266.65689033, 2363.64327736, 2463.18341061, 2565.37886352, 2669.99929908, 2776.84462448, 2886.05422112, 2997.36486954, 3110.61679467])
array([4.50047919e-38, 4.50046568e-38, 4.50044935e-38, 4.50042992e-38, 4.50040801e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50034577e-38, 4.50015314e-38, 4.50008048e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 1.75539429e-37, 9.42806065e-36, 1.61309766e-33, 1.51251529e-31, 1.24315417e-29, 9.59176341e-28, 6.32966742e-26, 3.54672037e-24, 2.33309090e-22, 1.04522239e-20, 2.95648129e-19, 3.97108672e-18, 3.42027816e-17, 2.10189777e-16, 9.08099237e-16, 3.21795565e-15, 1.04464749e-14, 3.19102511e-14, 9.14982745e-14, 2.51409949e-13, 6.51541102e-13, 1.59998579e-12, 3.69451125e-12, 8.09712813e-12, 1.68764556e-11, 3.36622831e-11, 6.45002917e-11, 1.19077614e-10, 2.12274978e-10, 3.65912992e-10, 6.11534714e-10, 9.90806463e-10, 1.55906541e-09, 2.38501268e-09, 3.54842201e-09, 5.12660087e-09, 7.20343174e-09, 9.87098597e-09, 1.32299096e-08, 1.73666510e-08, 2.23583606e-08])
array([0.83765961, 0.83765822, 0.83765539, 0.83765044, 0.83764226, 0.8376383 , 0.8376383 , 0.8376383 , 0.8376383 , 0.8376383 , 0.8376383 , 0.8376383 , 0.83763829, 0.83763826, 0.83763812, 0.83763706, 0.83763577, 0.83762809, 0.83760231, 0.837535 , 0.83749793, 0.83745751, 0.83745751, 0.83745751, 0.83745751, 0.83745751, 0.83745751, 0.83745751, 0.83745751, 0.83745751, 0.83745751, 0.83745751, 0.83745751, 0.83745751, 0.83745751, 0.83745751, 0.83745751, 0.83745751, 0.83745751, 0.83745751, 0.83745751, 0.83745751, 0.83745751, 0.83745751, 0.83745751, 0.83745751, 0.83745751, 0.83745751, 0.83745751, 0.83745752, 0.83748038, 0.83750079, 0.83750084, 0.83750094, 0.83750113, 0.83750164, 0.83750208, 0.83750234, 0.83750355, 0.83750525, 0.837507 , 0.83750862, 0.8375101 , 0.83751134, 0.83751246, 0.83751418, 0.83751551, 0.83751665, 0.83751739, 0.83751818, 0.83751895, 0.83751912, 0.83751986, 0.83751954, 0.83751618, 0.83750585, 0.83749125, 0.83746711, 0.83743655, 0.83739141, 0.8373035 , 0.83717727, 0.83700036, 0.83674634, 0.83641331, 0.83593824, 0.83528742, 0.83445515, 0.83341904, 0.83219591, 0.83073557])
array([4.50047919e-38, 4.50046568e-38, 4.50044935e-38, 4.50042992e-38, 4.50040801e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50034577e-38, 4.50015314e-38, 4.50008048e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 1.21019355e-37, 7.23638775e-37, 6.50942947e-36, 5.85063919e-35, 5.40223409e-34, 5.54331886e-33, 6.73691959e-32, 1.00497031e-30, 1.84070259e-29, 4.46388461e-28, 1.29190127e-26, 3.14556228e-25, 6.50868395e-24, 1.15410645e-22, 1.76812377e-21, 2.35797521e-20, 2.75634347e-19, 2.84277943e-18, 2.60286767e-17, 2.12886154e-16, 1.56226094e-15, 1.03414611e-14, 6.20931396e-14, 3.39364638e-13, 1.69605246e-12, 7.78384837e-12, 3.29283159e-11, 1.28909039e-10, 4.68251076e-10, 1.58569921e-09, 5.01248356e-09, 1.48671180e-08, 4.14397656e-08, 1.09064729e-07, 2.71493031e-07, 6.41198935e-07, 1.44139150e-06, 3.09060643e-06, 6.33913940e-06, 1.24725016e-05, 2.35788452e-05, 4.29364138e-05, 7.54496942e-05, 1.28194170e-04, 2.11054672e-04, 3.36860771e-04, 5.22348778e-04, 7.87962511e-04, 1.15696713e-03, 1.65705618e-03, 2.31634161e-03, 3.16310342e-03, 4.22846541e-03, 5.53711735e-03, 7.10942281e-03])
array([4.50047919e-38, 4.50046568e-38, 4.50044935e-38, 4.50042992e-38, 4.50040801e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50034577e-38, 4.50015314e-38, 4.50008048e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50001051e-38, 4.50004639e-38, 4.50003248e-38, 4.50001835e-38, 4.50000109e-38, 4.50003445e-38, 4.50006589e-38, 4.50003197e-38, 4.50000014e-38, 4.50004416e-38, 4.50009080e-38, 4.50009657e-38, 4.50003901e-38, 4.71584740e-38, 5.29138479e-38, 2.01808334e-36, 1.05736271e-34, 5.22794469e-33, 2.08695136e-31, 6.81304839e-30, 1.84136033e-28, 4.13627994e-27, 7.82661529e-26, 1.25813278e-24, 1.72424715e-23, 2.04842674e-22, 2.11566412e-21, 1.90659447e-20, 1.51698651e-19, 1.06880830e-18, 6.70844211e-18])
array([4.50047919e-38, 4.50046568e-38, 4.50044935e-38, 4.50042992e-38, 4.50040801e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50034577e-38, 4.50015314e-38, 4.50008048e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 3.72245505e-37, 9.11261070e-36, 1.84040694e-33, 4.08572897e-31, 5.80962671e-29, 4.91810384e-27, 1.86893936e-25, 4.31215673e-24, 6.74721711e-23, 7.06649680e-22, 5.78648643e-21, 4.16301608e-20, 2.69627524e-19, 1.57915260e-18, 8.50549802e-18, 4.17713793e-17, 1.87787894e-16, 7.68209710e-16, 2.89639833e-15, 1.00909918e-14, 3.27222223e-14, 9.94098907e-14, 2.84169162e-13, 7.66864564e-13, 1.95793088e-12, 4.74946317e-12, 1.09551759e-11, 2.41024125e-11, 5.06785786e-11, 1.01930888e-10, 1.96061247e-10, 3.61419637e-10, 6.40739226e-10, 1.09699410e-09, 1.81712115e-09, 2.91783317e-09])
array([4.50047919e-38, 4.50046568e-38, 4.50044935e-38, 4.50042992e-38, 4.50040801e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50034577e-38, 4.50015314e-38, 4.50008048e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.97059215e-38, 7.49166706e-37, 1.24958293e-34, 1.02630266e-32, 4.93690026e-31, 1.55687636e-29, 3.18079443e-28, 4.92685089e-27, 6.51008121e-26, 7.53117421e-25, 7.67256697e-24, 7.00979752e-23, 5.70280563e-22, 4.15221534e-21, 2.69328965e-20, 1.57803469e-19, 8.38361136e-19, 4.07251536e-18, 1.82231408e-17, 7.55194803e-17, 2.91027116e-16, 1.04610421e-15, 3.52535402e-15, 1.11502941e-14, 3.32338804e-14, 9.35773001e-14, 2.49201661e-13, 6.28076054e-13, 1.50169337e-12, 3.41934043e-12, 7.45223477e-12, 1.55798664e-11, 3.13153228e-11])
array([4.50047919e-38, 4.50046568e-38, 4.50044935e-38, 4.50042992e-38, 4.50040801e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50034577e-38, 4.50015314e-38, 4.50008048e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50001051e-38, 4.50004639e-38, 4.50003248e-38, 4.50001835e-38, 4.50000109e-38, 4.50003445e-38, 4.50006589e-38, 4.50003197e-38, 4.50000014e-38, 4.50004416e-38, 4.50009080e-38, 4.50009657e-38, 4.50003901e-38, 5.07230799e-38, 7.02403726e-38, 3.28081042e-36, 1.91439512e-34, 8.99161958e-33, 3.43490761e-31, 1.08063906e-29, 2.83312624e-28, 6.21327059e-27, 1.15472743e-25, 1.83373169e-24, 2.49663866e-23, 2.96192197e-22, 3.07032210e-21, 2.79049639e-20, 2.24917335e-19, 1.61218292e-18, 1.03367830e-17])
array([4.50047919e-38, 4.50046568e-38, 4.50044935e-38, 4.50042992e-38, 4.50040801e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50034577e-38, 4.50015314e-38, 4.50008048e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50001051e-38, 4.50004639e-38, 4.50003248e-38, 4.50001835e-38, 4.50000109e-38, 4.50003445e-38, 4.50006589e-38, 8.63755543e-38, 1.91563574e-37, 1.06240194e-35, 4.92404008e-34, 1.90989644e-32, 6.29543868e-31, 1.74901341e-29, 4.14751387e-28, 8.43773674e-27, 1.47223984e-25, 2.22213701e-24, 2.91849988e-23, 3.36181013e-22, 3.41998686e-21, 3.07871335e-20, 2.47391480e-19, 1.78422734e-18, 1.15742604e-17, 6.83589424e-17, 3.68100396e-16, 1.81032454e-15, 8.19336161e-15, 3.41825508e-14, 1.31960593e-13])
array([0.16238 , 0.16238 , 0.16238 , 0.16238 , 0.16238 , 0.16238 , 0.16238 , 0.16238 , 0.16238 , 0.16238 , 0.16238 , 0.16238 , 0.16238 , 0.16238 , 0.16238 , 0.16237933, 0.16237921, 0.1623752 , 0.16236872, 0.16236453, 0.16236268, 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16236 , 0.16235975, 0.16236 , 0.16236 , 0.16235855, 0.16235499, 0.16235193, 0.16235 , 0.16234971, 0.16234931, 0.16234239, 0.16233517, 0.16232382, 0.16230697, 0.16228718, 0.16225462, 0.16221418, 0.16215777, 0.16208162, 0.16199438, 0.16187954, 0.16175827])
array([8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.17225696e-05, 8.22153764e-05, 8.38513008e-05, 8.78315928e-05, 9.06945019e-05, 9.27071226e-05, 9.25398638e-05, 9.14460999e-05, 9.03818167e-05, 8.96514303e-05, 8.93523088e-05, 8.94083956e-05, 8.96268342e-05, 8.93518129e-05, 8.91846800e-05, 8.94149818e-05, 8.99966047e-05, 9.08285564e-05, 9.17734275e-05])
array([2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.28192237e-05, 2.24554580e-05, 2.19421083e-05, 2.10035268e-05, 2.02522428e-05, 1.99384732e-05, 2.00074119e-05, 1.86929195e-05, 1.73304594e-05, 1.62123200e-05, 1.54946778e-05, 1.52656376e-05, 1.55604251e-05, 1.48859990e-05, 1.38768022e-05, 1.31744570e-05, 1.28117543e-05, 1.28096556e-05, 1.31619500e-05])
array([2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 2.80503226e-07, 4.62359145e-07, 7.31779971e-07, 1.11373747e-06, 1.61296526e-06, 2.26778701e-06, 3.05329965e-06, 3.65152644e-06, 4.18647023e-06, 4.69923481e-06, 5.24665272e-06, 5.88301313e-06, 6.63932694e-06, 6.80843583e-06, 6.63710045e-06, 6.45623404e-06, 6.31414693e-06, 6.22933097e-06, 6.19348269e-06])
array([4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 4.04323366e-06, 3.99304884e-06, 3.99458027e-06, 4.03896987e-06, 4.09925212e-06, 4.09710178e-06, 4.15539995e-06, 4.26987182e-06, 4.43648347e-06, 4.65391620e-06, 4.92336970e-06, 5.09372767e-06, 5.26763145e-06, 5.50921981e-06, 5.81713255e-06, 6.19293142e-06, 6.63657914e-06])
array([4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.45691620e-06, 4.52271497e-06, 4.61190093e-06, 4.71039792e-06, 4.75633494e-06, 4.57404885e-06, 4.45320221e-06, 4.37681828e-06, 4.33064500e-06, 4.30390086e-06, 4.28809193e-06, 4.03322606e-06, 3.74290424e-06, 3.51043975e-06, 3.31984205e-06, 3.15848639e-06, 3.01520481e-06])
array([4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 7.68256220e-38, 2.61136501e-37, 1.96650946e-36, 1.62660584e-35, 1.30304871e-34, 1.04973977e-33, 8.89166826e-33, 9.39667295e-32, 8.66816802e-31, 7.82557634e-30, 7.73571003e-29, 8.85492538e-28, 1.21010847e-26, 1.96544281e-25, 4.07469275e-24, 9.74568360e-23, 1.97239780e-21, 3.41793650e-20, 5.11013410e-19, 6.64281870e-18, 7.52524779e-17, 7.52446137e-16, 6.74917360e-15, 5.35377076e-14, 3.82756553e-13, 2.47803958e-12, 1.44845890e-11, 7.77070931e-11, 3.79332157e-10, 1.70626452e-09, 6.93992559e-09, 2.62884861e-08, 6.58623080e-08, 1.37998838e-07, 1.66037019e-07, 1.66183297e-07, 1.64036792e-07, 1.61879506e-07, 1.57482200e-07, 1.52993405e-07, 1.47356801e-07, 1.39366042e-07, 1.31119725e-07, 1.22183356e-07, 1.12308607e-07, 1.02943449e-07, 9.41362954e-08, 8.59351552e-08, 7.82292227e-08, 7.09034235e-08, 6.50979237e-08, 5.98448338e-08, 5.54715763e-08, 5.21429472e-08, 4.95042022e-08, 4.76931241e-08, 4.69421025e-08, 4.68114716e-08, 4.78031420e-08, 4.91942895e-08])
array([4.94524661e-13, 8.25217343e-13, 1.58722913e-12, 3.57003465e-12, 9.09879755e-12, 2.49341423e-11, 8.10467033e-11, 3.19741611e-10, 1.44305223e-09, 6.99212207e-09, 3.87552414e-08, 1.90625587e-07, 6.90823907e-07, 1.86603654e-06, 2.22112192e-06, 2.53842017e-06, 2.53847545e-06, 2.53849451e-06, 2.53862091e-06, 2.53906838e-06, 2.53884591e-06, 2.53850336e-06, 2.53828685e-06, 2.53820000e-06, 2.53820252e-06, 2.53821970e-06, 2.53821619e-06, 2.53820000e-06, 2.53820000e-06, 2.53820000e-06, 2.53820000e-06, 2.53820000e-06, 2.53820000e-06, 2.53820000e-06, 2.53820000e-06, 2.53820000e-06, 2.53820000e-06, 2.53820000e-06, 2.53820000e-06, 2.53820000e-06, 2.53820000e-06, 2.53820000e-06, 2.53820000e-06, 2.53820000e-06, 2.53820000e-06, 2.53820000e-06, 2.53820000e-06, 2.53820000e-06, 2.53820000e-06, 2.53820000e-06, 2.53842810e-06, 2.53890619e-06, 2.53957679e-06, 2.54036433e-06, 2.54116687e-06, 2.54169427e-06, 2.54212440e-06, 2.54236594e-06, 2.54251177e-06, 2.54280732e-06, 2.54409559e-06, 2.54797149e-06, 2.57379951e-06, 2.65494886e-06, 2.71043954e-06, 2.72470666e-06, 2.72470278e-06, 2.72467379e-06, 2.72465745e-06, 2.72451841e-06, 2.72413660e-06, 2.72341763e-06, 2.71565519e-06, 2.68592059e-06, 2.61719371e-06, 2.55714442e-06, 2.52216079e-06, 2.50671234e-06, 2.51096326e-06, 2.51490042e-06, 2.51424237e-06, 2.51033638e-06, 2.50388590e-06, 2.49539570e-06, 2.48533996e-06, 2.47168077e-06, 2.45255975e-06, 2.42882395e-06, 2.40088197e-06, 2.36980352e-06, 2.33560369e-06])
array([4.50047919e-38, 4.50046568e-38, 4.50044935e-38, 4.50042992e-38, 4.50040801e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50034577e-38, 4.50015314e-38, 4.50008048e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 2.33859952e-37, 7.97496798e-36, 1.51681425e-33, 2.16713778e-31, 2.35025469e-29, 1.96742103e-27, 1.28226075e-25, 6.58048358e-24, 2.70477566e-22, 8.93487899e-21, 2.41805186e-19, 5.31537891e-18, 9.34035427e-17, 1.26133057e-15, 1.45424004e-14, 1.44369769e-13, 1.23884801e-12, 9.19136781e-12, 6.04158908e-11, 3.53858001e-10, 1.87474772e-09, 8.97019223e-09, 1.19128362e-08, 1.17844092e-08, 1.04273122e-08, 8.98228085e-09, 7.54698191e-09, 6.18388294e-09, 4.98823982e-09, 3.99541937e-09, 3.19295916e-09, 2.55539422e-09, 2.05050186e-09])
array([4.50047919e-38, 4.50046568e-38, 4.50044935e-38, 4.50042992e-38, 4.50040801e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50034577e-38, 4.50015314e-38, 4.50008048e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 1.44892038e-37, 4.07431000e-36, 5.17483194e-34, 6.80788635e-32, 6.87122289e-30, 5.38457765e-28, 3.30657377e-26, 1.61512297e-24, 6.32370352e-23, 2.00476249e-21, 5.22452369e-20, 1.12276659e-18, 2.02358144e-17, 3.05861800e-16, 3.94650305e-15, 4.36183556e-14, 4.17098702e-13, 3.47519975e-12, 2.53251319e-11, 1.42282221e-10, 3.12577626e-10, 3.25730424e-10, 2.42901928e-10, 1.78063922e-10, 1.30904015e-10, 9.73792286e-11, 7.34939941e-11, 5.62782547e-11, 4.36670016e-11, 3.41123201e-11, 2.70263344e-11, 2.17757995e-11, 1.78136615e-11, 1.47736363e-11, 1.23953894e-11])
array([4.50047919e-38, 4.50046568e-38, 4.50044935e-38, 4.50042992e-38, 4.50040801e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50034577e-38, 4.50015314e-38, 4.50008048e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.79275937e-37, 4.64826392e-35, 4.81575163e-33, 3.91657241e-31, 2.53013593e-29, 1.31079324e-27, 5.51570845e-26, 1.89960447e-24, 5.40307052e-23, 1.28080936e-21, 2.55103931e-20, 4.30110809e-19, 6.16978202e-18, 7.59927697e-17, 8.07520912e-16, 7.41047085e-15, 5.96239207e-14, 4.20044422e-13, 2.63264226e-12, 1.46622433e-11, 7.36647888e-11, 3.34778291e-10, 1.38393535e-09, 5.23722641e-09, 1.82165054e-08, 5.84305702e-08, 1.73613000e-07, 4.74815679e-07, 1.18285940e-06, 2.75650314e-06, 4.18408367e-06, 5.06158453e-06, 5.30873781e-06, 5.27809510e-06, 5.24433050e-06, 5.19887164e-06, 5.14553835e-06, 5.08772848e-06, 5.02557659e-06, 4.95941767e-06, 4.88949219e-06])
array([4.50047919e-38, 4.50046568e-38, 4.50044935e-38, 4.50042992e-38, 4.50040801e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50034577e-38, 4.50015314e-38, 4.50008048e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 1.19433935e-37, 4.45680259e-36, 4.40111012e-34, 3.42690891e-32, 2.13292772e-30, 1.07017574e-28, 4.37328413e-27, 1.46955576e-25, 4.09693822e-24, 9.55464431e-23, 1.87520132e-21, 3.12706698e-20, 4.45389793e-19, 5.42886678e-18, 5.74379354e-17, 5.28387172e-16, 4.28607814e-15, 3.06960302e-14, 1.96736790e-13, 1.13292358e-12, 5.89804646e-12, 2.79261595e-11, 1.20833896e-10, 4.79641691e-10, 1.75396927e-09, 5.87345449e-09, 1.78271973e-08, 5.03866980e-08, 9.23467214e-08, 1.34267156e-07, 1.68511716e-07, 1.99624812e-07, 2.35309599e-07, 2.75604500e-07, 3.20929971e-07, 3.71792843e-07, 4.28616812e-07, 4.91750420e-07, 5.61469006e-07])
array([4.50047919e-38, 4.50046568e-38, 4.50044935e-38, 4.50042992e-38, 4.50040801e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50034577e-38, 4.50015314e-38, 4.50008048e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 5.47336464e-37, 5.35911502e-35, 4.74944226e-33, 3.33365729e-31, 1.87944522e-29, 8.57956227e-28, 3.20346211e-26, 9.87939087e-25, 2.53889788e-23, 5.48229739e-22, 1.00124457e-20, 1.56153562e-19, 2.09057672e-18, 1.53207362e-17, 8.48514516e-17, 8.22348608e-16, 8.36365830e-15, 6.88385638e-14, 4.94366086e-13, 3.06103461e-12, 1.66107017e-11, 7.57910819e-11, 2.83878780e-10, 5.07959605e-10, 7.43753579e-10, 9.24172634e-10, 1.13465846e-09, 1.37754060e-09, 1.65734661e-09, 1.97788665e-09, 2.34330972e-09, 2.75752690e-09, 3.22417671e-09, 3.74207488e-09, 4.31472525e-09, 4.94667857e-09, 5.64116553e-09, 6.40067741e-09, 7.22659345e-09])
array([4.50047919e-38, 4.50046568e-38, 4.50044935e-38, 4.50042992e-38, 4.50040801e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50034577e-38, 4.50015314e-38, 4.50008048e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 1.12753516e-28, 3.66099393e-19, 7.58586059e-18, 1.33235153e-16, 2.00002389e-15, 2.58044139e-14, 2.89457462e-13, 2.83366490e-12, 2.43539005e-11, 1.85135155e-10, 1.25139181e-09, 7.55801994e-09, 4.07886914e-08, 1.75298599e-07, 3.15894815e-07, 3.48084282e-07, 3.48354707e-07, 3.50165116e-07, 3.52028005e-07, 3.56309249e-07, 3.60211521e-07, 3.63937980e-07, 3.65988181e-07, 3.67238269e-07, 3.67509645e-07, 3.67139888e-07, 3.66467411e-07, 3.65533664e-07, 3.64320723e-07, 3.62529415e-07, 3.60481756e-07, 3.58181693e-07, 3.55593584e-07, 3.52694226e-07, 3.49437408e-07, 3.44947045e-07, 3.39468273e-07, 3.33219931e-07, 3.26078002e-07, 3.18001566e-07, 3.08741072e-07])
array([4.50047919e-38, 4.50046568e-38, 4.50044935e-38, 4.50042992e-38, 4.50040801e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50034577e-38, 4.50015314e-38, 4.50008048e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 1.21004842e-32, 1.77003823e-26, 2.41679494e-24, 2.86984978e-22, 2.98386795e-20, 2.74666505e-18, 2.02472650e-16, 1.23687057e-14, 6.13669920e-13, 1.71613793e-11, 2.60385852e-10, 1.46323378e-09, 5.00201080e-09, 1.09829042e-08, 1.45887436e-08, 1.56217485e-08, 1.61848006e-08, 1.70422972e-08, 1.79225672e-08, 1.93222802e-08, 2.06585770e-08, 2.19334838e-08, 2.27237925e-08, 2.32744253e-08, 2.35259385e-08, 2.35903505e-08, 2.35720547e-08, 2.34823311e-08, 2.33269215e-08, 2.30541481e-08, 2.26767489e-08, 2.21846793e-08, 2.15577305e-08, 2.07842277e-08, 1.98676474e-08, 1.86667024e-08, 1.72418851e-08, 1.57051702e-08, 1.41287833e-08, 1.25899816e-08, 1.11168213e-08])
array([4.50047919e-38, 4.50046568e-38, 4.50044935e-38, 4.50042992e-38, 4.50040801e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50034577e-38, 4.50015314e-38, 4.50008048e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 3.20885381e-32, 1.11099811e-25, 1.35874425e-23, 1.45304928e-21, 1.36765332e-19, 3.64201821e-18, 6.90421145e-17, 7.27905676e-16, 6.91629926e-15, 5.77224956e-14, 4.10733956e-13, 2.33290395e-12, 8.28549489e-12, 1.91461978e-11, 2.60693473e-11, 2.81377703e-11, 2.92722600e-11, 3.10400008e-11, 3.28695039e-11, 3.58443024e-11, 3.87219166e-11, 4.15213751e-11, 4.32784031e-11, 4.45072665e-11, 4.50597671e-11, 4.51632401e-11, 4.50495989e-11, 4.47373834e-11, 4.42306899e-11, 4.33598602e-11, 4.21462485e-11, 4.05836363e-11, 3.86547837e-11, 3.63824975e-11, 3.38364600e-11, 3.07546482e-11, 2.73989973e-11, 2.40840199e-11, 2.09585145e-11, 1.81299770e-11, 1.56053969e-11])
array([4.50047919e-38, 4.50046568e-38, 4.50044935e-38, 4.50042992e-38, 4.50040801e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50034577e-38, 4.50015314e-38, 4.50008048e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 2.22838614e-32, 5.55595983e-26, 3.60410524e-24, 7.91608459e-23, 1.17595256e-21, 1.73215498e-20, 2.26666496e-19, 2.72104447e-18, 2.92547328e-17, 2.75203896e-16, 2.21566612e-15, 1.45058109e-14, 6.63468494e-14, 2.23400821e-13, 3.90467785e-13, 4.83775492e-13, 5.58255827e-13, 6.69840882e-13, 7.98467874e-13, 1.02983479e-12, 1.29933114e-12, 1.62631837e-12, 1.87597987e-12, 2.07293640e-12, 2.17120003e-12, 2.19219148e-12, 2.17615724e-12, 2.12822548e-12, 2.05269242e-12, 1.93304600e-12, 1.78356968e-12, 1.61673261e-12, 1.44130693e-12, 1.26612529e-12, 1.09904181e-12, 9.32134435e-13, 7.80719458e-13, 6.52293073e-13, 5.45396563e-13, 4.57836596e-13, 3.86088802e-13])
array([4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 4.35031114e-13, 8.49266710e-13, 1.57383202e-12, 2.77385557e-12, 4.66440726e-12, 7.51914842e-12, 1.14907810e-11, 1.69901374e-11, 2.46575984e-11, 3.56257664e-11, 4.86179821e-11, 6.41790712e-11, 7.97150654e-11, 8.75860862e-11, 9.27743873e-11, 9.70507432e-11, 1.01944401e-10, 1.08453399e-10, 1.16835042e-10, 1.14187589e-10, 1.06620442e-10, 1.00145311e-10, 9.52551233e-11, 9.19261707e-11, 8.97800594e-11])
array([2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 2.03995937e-10, 3.65889551e-10, 6.28037653e-10, 1.01752485e-09, 1.60463140e-09, 2.49838994e-09, 3.81593342e-09, 5.13548327e-09, 6.70481242e-09, 8.62773620e-09, 1.10855188e-08, 1.43507052e-08, 1.87786381e-08, 2.18999184e-08, 2.41708374e-08, 2.67275237e-08, 2.98018896e-08, 3.36153721e-08, 3.83139599e-08])
array([4.50047919e-38, 4.50046568e-38, 4.50044935e-38, 4.50042992e-38, 4.50040801e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50034577e-38, 4.50015314e-38, 4.50008048e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 1.22253294e-37, 2.38257605e-36, 1.54393302e-34, 8.16527527e-33, 3.52680039e-31, 1.25995051e-29, 3.72459263e-28, 9.22323326e-27, 1.93162283e-25, 3.43089282e-24, 5.20902662e-23, 6.80671895e-22, 7.70215937e-21, 7.58826616e-20, 6.53453380e-19, 4.95557239e-18, 3.31697839e-17, 1.97209475e-16, 1.04771369e-15, 4.97381283e-15, 2.13298859e-14, 8.31456071e-14, 2.89578038e-13, 9.13603343e-13, 2.61292350e-12, 6.61949651e-12, 1.57771009e-11, 3.65035598e-11, 8.19632438e-11, 1.49476447e-10, 2.52543906e-10, 4.09354698e-10, 6.55186174e-10, 1.05688340e-09, 1.73301089e-09, 2.38721110e-09, 2.98204435e-09, 3.69706569e-09, 4.62040565e-09, 5.87181938e-09, 7.59519156e-09])
array([4.50047919e-38, 4.50046568e-38, 4.50044935e-38, 4.50042992e-38, 4.50040801e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50034577e-38, 4.50015314e-38, 4.50008048e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 2.30961487e-31, 1.44968940e-24, 1.36959954e-23, 1.15588792e-22, 8.75990341e-22, 6.00302429e-21, 3.71902552e-20, 2.09723049e-19, 1.08214276e-18, 5.11502133e-18, 2.22417294e-17, 8.92971789e-17, 3.32101703e-16, 1.14717044e-15, 3.68807679e-15, 1.10783482e-14, 3.11229649e-14, 8.20062855e-14, 2.03469584e-13, 4.74127079e-13, 1.04636054e-12, 2.19624763e-12, 4.29862863e-12, 7.93621394e-12, 1.38015952e-11, 2.20465296e-11, 3.42925608e-11, 5.34816041e-11, 8.34699788e-11, 1.08944697e-10, 1.35371510e-10, 1.65727372e-10, 2.05349590e-10, 2.62544261e-10, 3.49106416e-10, 3.98230813e-10, 4.20386567e-10, 4.49076520e-10, 4.92142042e-10, 5.57825305e-10, 6.53959913e-10])
array([4.50047919e-38, 4.50046568e-38, 4.50044935e-38, 4.50042992e-38, 4.50040801e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50034577e-38, 4.50015314e-38, 4.50008048e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 2.11824343e-27, 3.66234842e-17, 8.08569211e-17, 1.72299114e-16, 3.54944233e-16, 7.08830106e-16, 1.37048636e-15, 2.57273135e-15, 4.69796815e-15, 8.34244859e-15, 1.44264175e-14, 2.43231293e-14, 4.00289088e-14, 6.43295760e-14, 1.01012813e-13, 1.55176542e-13, 2.33124823e-13, 3.42727680e-13, 4.93976722e-13, 6.95240870e-13, 9.61149895e-13, 1.30900126e-12, 1.71867208e-12, 2.19625357e-12, 2.72377875e-12, 3.19213509e-12, 3.74123014e-12, 4.50803632e-12, 5.56787450e-12, 5.88340057e-12, 6.04587116e-12, 6.24909816e-12, 6.66428286e-12, 7.46881577e-12, 8.86180513e-12, 9.16800705e-12, 8.91658257e-12, 8.90862018e-12, 9.25613589e-12, 1.00784407e-11, 1.14914867e-11])
array([4.50047919e-38, 4.50046568e-38, 4.50044935e-38, 4.50042992e-38, 4.50040801e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50034577e-38, 4.50015314e-38, 4.50008048e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 1.51243369e-37, 1.72751176e-35, 4.19036593e-33, 7.61380403e-31, 1.04533801e-28, 1.09810392e-26, 8.95612752e-25, 5.73582192e-23, 2.91529108e-21, 1.18540935e-19, 3.82670244e-18, 1.00782032e-16, 2.19646305e-15, 4.01586003e-14, 6.17866330e-13, 8.11915216e-12, 9.06560190e-11, 8.76268439e-10, 7.33617609e-09, 5.40114826e-08, 3.36860798e-07, 1.74635787e-06, 3.46040401e-06, 4.97030580e-06, 5.80264272e-06, 5.81387700e-06, 5.63372556e-06, 5.50920010e-06, 5.40307809e-06, 5.28261975e-06, 5.12022001e-06, 4.89281874e-06, 4.36844978e-06, 3.78138266e-06, 3.22553558e-06, 2.70001717e-06, 2.21535902e-06, 1.77975621e-06])
array([4.50047919e-38, 4.50046568e-38, 4.50044935e-38, 4.50042992e-38, 4.50040801e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50034577e-38, 4.50015314e-38, 4.50008048e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 2.30136095e-37, 2.78134277e-35, 4.13690847e-33, 5.14804205e-31, 4.95926529e-29, 3.74130547e-27, 2.23125977e-25, 1.07912033e-23, 4.21865469e-22, 1.34316650e-20, 3.53176878e-19, 7.69889887e-18, 1.41328235e-16, 2.19338728e-15, 2.91988753e-14, 3.34992707e-13, 3.33454977e-12, 2.94395960e-11, 2.32379334e-10, 1.95131724e-09, 7.42649031e-09, 1.44107254e-08, 1.98108625e-08, 2.68923197e-08, 3.60756558e-08, 4.78302008e-08, 6.27176794e-08, 8.13657484e-08, 1.04448125e-07, 1.32656612e-07, 1.66744204e-07, 2.07512691e-07, 2.55850141e-07, 3.12573656e-07, 3.78438009e-07])
array([4.50047919e-38, 4.50046568e-38, 4.50044935e-38, 4.50042992e-38, 4.50040801e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50034577e-38, 4.50015314e-38, 4.50008048e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.78904244e-38, 5.00659154e-38, 1.11625223e-37, 4.52814651e-37, 2.32115703e-36, 1.45169930e-35, 1.06349066e-34, 9.43942227e-34, 1.02079304e-32, 1.43337426e-31, 2.40107078e-30, 3.48726333e-29, 4.41759180e-28, 4.91595431e-27, 3.98122700e-26, 2.90132401e-25, 2.27302709e-24, 1.60083761e-23, 1.01806927e-22, 5.89386723e-22, 3.09319724e-21, 1.48420218e-20, 6.54716917e-20, 2.65333079e-19, 9.92463491e-19, 3.44091467e-18, 1.11789933e-17, 3.44987365e-17, 9.82100711e-17, 2.57806532e-16, 6.29306247e-16, 1.44427180e-15, 3.11826410e-15, 6.34583006e-15, 1.22184512e-14, 2.23560651e-14, 3.84003627e-14, 6.24331813e-14, 9.61638678e-14, 1.43041584e-13, 2.04578131e-13, 2.86490989e-13, 3.87592946e-13, 4.66129913e-13, 5.36866500e-13, 6.04994549e-13, 6.77867915e-13, 7.62601045e-13, 8.63396448e-13, 8.87421232e-13, 8.65548511e-13, 8.41203978e-13, 8.20860900e-13, 8.07150696e-13, 7.99056820e-13])
array([4.50047919e-38, 4.50046568e-38, 4.50044935e-38, 4.50042992e-38, 4.50040801e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50034577e-38, 4.50015314e-38, 4.50008048e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 5.28762271e-36, 7.66505080e-33, 1.99568345e-30, 3.86859862e-28, 5.66491553e-26, 6.33656105e-24, 5.50620391e-22, 3.74847929e-20, 2.01532935e-18, 7.67570674e-17, 1.91647285e-15, 2.35938332e-14, 1.29542430e-13, 5.15098846e-13, 1.32500255e-12, 2.48462305e-12, 4.26164616e-12, 7.06955899e-12, 1.13141161e-11, 1.99798166e-11, 3.45567474e-11, 6.10699665e-11, 9.66247992e-11, 1.39137663e-10, 1.75546869e-10, 1.98358972e-10, 2.14723959e-10, 2.25249590e-10, 2.31636210e-10, 2.33071570e-10, 2.33062800e-10, 2.32231187e-10, 2.30873644e-10, 2.29175251e-10, 2.27245082e-10, 2.23354117e-10, 2.18666223e-10, 2.13949042e-10, 2.09209358e-10, 2.04423792e-10, 1.99547884e-10])
array([4.50047919e-38, 4.50046568e-38, 4.50044935e-38, 4.50042992e-38, 4.50040801e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50034577e-38, 4.50015314e-38, 4.50008048e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 1.72386571e-31, 2.02808454e-24, 1.07440965e-22, 4.62008950e-21, 1.62912811e-19, 4.74934707e-18, 1.15730250e-16, 2.37227752e-15, 4.11067161e-14, 5.37984416e-13, 4.90720929e-12, 2.34031413e-11, 5.26659930e-11, 9.04464651e-11, 1.05729304e-10, 9.45833901e-11, 8.09356465e-11, 7.00064149e-11, 6.08128950e-11, 6.06524428e-11, 6.14282596e-11, 6.58740593e-11, 6.53334783e-11, 6.12217930e-11, 5.24679698e-11, 4.25157048e-11, 3.33851683e-11, 2.57638187e-11, 1.94746465e-11, 1.45708398e-11, 1.10578315e-11, 8.55840758e-12, 6.76217343e-12, 5.44933524e-12, 4.46909587e-12, 3.66915693e-12, 3.04941072e-12, 2.58083692e-12, 2.22006925e-12, 1.93765523e-12, 1.71148610e-12])
array([4.50047919e-38, 4.50046568e-38, 4.50044935e-38, 4.50042992e-38, 4.50040801e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50034577e-38, 4.50015314e-38, 4.50008048e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 7.21486057e-36, 1.34392167e-32, 2.98532788e-30, 4.99349027e-28, 6.37833753e-26, 6.28730264e-24, 4.86418029e-22, 2.97666392e-20, 1.45189307e-18, 5.06195813e-17, 1.16691890e-15, 1.33741544e-14, 6.89107458e-14, 2.59101597e-13, 6.34966326e-13, 1.14241736e-12, 1.89280166e-12, 3.05316770e-12, 4.78036764e-12, 8.30908018e-12, 1.42237501e-11, 2.50127911e-11, 3.95817860e-11, 5.72766643e-11, 7.29471671e-11, 8.35648308e-11, 9.20749398e-11, 9.86862172e-11, 1.04071269e-10, 1.07747721e-10, 1.11215819e-10, 1.14747032e-10, 1.18458183e-10, 1.22431798e-10, 1.26736438e-10, 1.30345782e-10, 1.33824637e-10, 1.37610676e-10, 1.41684426e-10, 1.46028605e-10, 1.50589754e-10])
array([4.50047919e-38, 4.50046568e-38, 4.50044935e-38, 4.50042992e-38, 4.50040801e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50034577e-38, 4.50015314e-38, 4.50008048e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 2.01875486e-27, 7.51139831e-17, 8.49711481e-16, 7.42279721e-15, 5.07226084e-14, 2.71815799e-13, 1.29589196e-12, 5.36317460e-12, 1.98773850e-11, 6.05137067e-11, 1.42850097e-10, 2.06226745e-10, 1.87095494e-10, 1.59618381e-10, 1.75609824e-10, 2.06477071e-10, 2.43942568e-10, 2.61416070e-10, 2.74704760e-10, 2.62102933e-10, 2.39699098e-10, 2.00233908e-10, 1.40742854e-10, 9.05590761e-11, 5.13465821e-11, 2.67441476e-11, 1.38731083e-11, 7.24620701e-12, 3.85727839e-12, 2.08572543e-12, 1.16225458e-12, 6.68933740e-13, 3.97185580e-13, 2.43112382e-13, 1.53370536e-13, 9.86982232e-14, 6.50576876e-14, 4.40241157e-14, 3.05072157e-14, 2.16290654e-14, 1.56659296e-14])
array([4.50047919e-38, 4.50046568e-38, 4.50044935e-38, 4.50042992e-38, 4.50040801e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50040000e-38, 4.50034577e-38, 4.50015314e-38, 4.50008048e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 4.50000000e-38, 2.41051183e-29, 1.80808629e-20, 3.66792623e-19, 6.31873561e-18, 9.31824382e-17, 1.18273333e-15, 1.30643256e-14, 1.26061206e-13, 1.06866242e-12, 7.18928108e-12, 3.49063188e-11, 8.94782372e-11, 1.09675265e-10, 1.05345589e-10, 7.30720555e-11, 4.00036615e-11, 2.15090966e-11, 1.18579191e-11, 6.70777178e-12, 4.43692941e-12, 3.04377538e-12, 2.25284136e-12, 1.57178962e-12, 1.04796179e-12, 6.41072145e-13, 3.66811284e-13, 2.09528530e-13, 1.20556937e-13, 7.05486232e-14, 4.18201989e-14, 2.54508690e-14, 1.59326846e-14, 1.02494422e-14, 6.77099255e-15, 4.59257556e-15, 3.16649335e-15, 2.22863077e-15, 1.60492862e-15, 1.17997696e-15, 8.85019163e-16, 6.76212993e-16])
PandasIndex(Float64Index([14.999760753054142, 14.999510759124895, 14.999260769362177, 14.999010783765922, 14.998760802336058, 14.99851082507252, 14.998260851975232, 14.998010883044127, 14.997760918279138, 14.997510957680193, ... 0.6000915069772242, 0.6000815055354544, 0.6000715042603735, 0.6000615031519793, 0.6000515022102686, 0.6000415014352392, 0.6000315008268879, 0.600021500385212, 0.6000115001102087, 0.600001500001875], dtype='float64', name='wavelength', length=193132))
PandasIndex(Float64Index([ 0.0001778, 0.000210768924604, 0.000249851178733, 0.0002961803388785, 0.0003511001772457, 0.0004162036377185, 0.0004933790390226, 0.0005848648452022, 0.0006933145920246, 0.0008218738524933, 0.000974271473848, 0.0011549277323696, 0.0013690825430084, 0.0016229474425423, 0.001923885754519, 0.0022806261616478, 0.0027035158802829, 0.0032048207802987, 0.0037990811552991, 0.0045035334622373, 0.0053386102629582, 0.0063285328684122, 0.0075020138751206, 0.0088930899708864, 0.0105421091651885, 0.0124969010787681, 0.0148141642365284, 0.0175611106020272, 0.0208174150531019, 0.0246775263429574, 0.0292534065758858, 0.0346777786558133, 0.0411079758995613, 0.0487305054724325, 0.0577664580080703, 0.0684779203180215, 0.0811755771909401, 0.096227722770782, 0.1140709428880285, 0.1352227781837813, 0.1602967353192546, 0.1900200816691478, 0.2252549396320263, 0.2670232923958693, 0.3165366264482683, 0.3752310706090316, 0.4448090507889792, 0.5272886686666364, 0.625062236505918, 0.7409656659107748, 0.8783607231300077, 1.0412325367183388, 1.2343051857527547, 1.4631787212276317, 1.734491594919194, 2.056113206950672, 2.4373721568791646, 2.889326818701921, 3.425086087780413, 4.060189602911503, 4.813058471845216, 5.705529573493154, 6.763488934204631, 8.017622549119784, 9.504306426090642, 11.266661667796392, 13.355805194593795, 15.832332385181749, 18.76807463890974, 22.24818283763207, 26.373597137690155, 31.263974728068654, 37.061160474028576, 43.93330111185604, 52.07972232649315, 61.73671017115252, 73.18436451835693, 86.75472332599927, 102.841393361868, 121.9109667247463, 144.51655429702564, 171.31382866514352, 203.0800418309876, 240.7365693209032, 285.3756345817088, 338.2919888018924, 401.02046502772737, 475.3804958273229, 563.5288857325677, 668.022369117856, 791.8917680000002], dtype='float64', name='pressure'))
Let's plot up the T(P) profile and the abundances of the Elf Owl Model first
fig,ax=plt.subplots(nrows=1,ncols=2,figsize=(20,6))
ax[0].set_ylim(800,1.77e-4)
ax[0].semilogy(ds['temperature'],ds['pressure'],linewidth=3)
ax[0].set_ylabel("Pressure [bars]",fontsize=20)
ax[0].set_xlabel("Temperature [K]",fontsize=20)
ax[1].set_ylim(800,1.77e-4)
ax[1].set_xlim(1e-12,1e-1)
ax[1].loglog(ds['CH4'],ds['pressure'],label="CH4",linewidth=3)
ax[1].loglog(ds['CO'],ds['pressure'],label="CO",linewidth=3)
ax[1].loglog(ds['NH3'],ds['pressure'],label="NH3",linewidth=3)
ax[1].loglog(ds['PH3'],ds['pressure'],label="PH3",linewidth=3)
ax[1].loglog(ds['CO2'],ds['pressure'],label="CO2",linewidth=3)
ax[1].loglog(ds['H2O'],ds['pressure'],label="H2O",linewidth=3)
ax[1].legend(fontsize=20)
ax[1].set_ylabel("Pressure [bars]",fontsize=20)
ax[1].set_xlabel("VMR",fontsize=20)
for i in range(2):
ax[i].minorticks_on()
ax[i].tick_params(axis='both',which='major',color="k",length =18, width=2,direction='in',labelsize=16)
ax[i].tick_params(axis='both',which='minor',color="k",length =12, width=1,direction='in',labelsize=16)
plt.show()
Let's plot up the spectrum of the Elf Owl Model
fig,ax=plt.subplots(nrows=1,ncols=1,figsize=(20,6))
ax.set_xlim(0.6,14)
ax.set_ylim(1e4,1e11)
ax.loglog(ds['wavelength'],ds['flux'],linewidth=3)
ax.set_ylabel("Pressure [bars]",fontsize=20)
ax.set_xlabel("Temperature [K]",fontsize=20)
ax.set_ylabel(r"Flux [erg/s/cm$^2$/cm]",fontsize=20)
ax.set_xlabel("Wavelength [microns]",fontsize=20)
ax.minorticks_on()
ax.tick_params(axis='both',which='major',color="k",length =18, width=2,direction='in',labelsize=16)
ax.tick_params(axis='both',which='minor',color="k",length =12, width=1,direction='in',labelsize=16)
plt.show()