To keep track of bugs we report them here so that they can be reproduced easily. Additionally, as soon as they are fixed they should disappear in this notebook. To each bug report add time and your name. Add new bugs at the end of the file.
In [2]:
%pylab nbagg
import sygma as s
reload(s)
s.__file__
Populating the interactive namespace from numpy and matplotlib
Out[2]:
'/astro/adampaul/nupycee/sygma.pyc'
In [2]:
!echo $PYTHONPATH
:/Users/benoitcote/Desktop/OMEGA_SYGMA_inheritance
In [ ]:
In [4]:
s_0_1_10=sygma.sygma(imf_bdys=[0.1,10])
s_1_10=sygma.sygma(imf_bdys=[1,10])
s_2_8_10=sygma.sygma(imf_bdys=[2.8,10])
---------------------------------------------------------------------------
NameError Traceback (most recent call last)
<ipython-input-4-74a4524ab009> in <module>()
----> 1 s_0_1_10=sygma.sygma(imf_bdys=[0.1,10])
2 s_1_10=sygma.sygma(imf_bdys=[1,10])
3 s_2_8_10=sygma.sygma(imf_bdys=[2.8,10])
NameError: name 'sygma' is not defined
In [107]:
s_0_1_10.plot_sn_distr(rate_only='sn2',label2='[0.1,10]',marker2='d')
s_1_10.plot_sn_distr(rate_only='sn2',label2='[1,10]',marker2='o')
s_2_8_10.plot_sn_distr(rate_only='sn2',label2='[2.8,10]',marker2='x')
plt.xlim(1e7,6e7)
Out[107]:
(10000000.0, 60000000.0)
In [109]:
s_0_1_10.plot_sn_distr(rate_only='sn1a',label1='[0.1,10]',marker1='d')
s_1_10.plot_sn_distr(rate_only='sn1a',label1='[1,10]',marker1='o')
s_2_8_10.plot_sn_distr(rate_only='sn1a',label1='[2.8,10]',marker1='x')
plt.ylim(1e-1,1e5)
plt.xlim(1e6,1.5e10)
Out[109]:
(1000000.0, 15000000000.0)
In [110]:
print sum(s_0_1_10.history.sn1a_numbers)
print sum(s_1_10.history.sn1a_numbers)
41653255.3985
41653255.3985
In [111]:
c_2_8_10=sygma.sygma(imf_type='chabrier',imf_bdys=[2.8,10])
c_1_10=sygma.sygma(imf_type='chabrier',imf_bdys=[1,10])
c_0_1_10=sygma.sygma(imf_type='chabrier',imf_bdys=[0.1,10])
################ Start SYGMA simulation ##############
10 30
Use SNIa tables from yield_tables/sn1a_t86.txt
START SFRMDOT Run time: 2.66s
Use NuGrid yields
Total mass of the gas in stars:
1.000E+11
sum msssfac 21001804996.4
END SFRMDOT Run time: 3.03s
################ Simulation successful ##############
Run time: 4.38s
################ Start SYGMA simulation ##############
10 30
Use SNIa tables from yield_tables/sn1a_t86.txt
START SFRMDOT Run time: 2.58s
Use NuGrid yields
Total mass of the gas in stars:
1.000E+11
sum msssfac 43945020318.1
END SFRMDOT Run time: 3.03s
################ Simulation successful ##############
Run time: 4.42s
################ Start SYGMA simulation ##############
10 30
Use SNIa tables from yield_tables/sn1a_t86.txt
START SFRMDOT Run time: 2.79s
Use NuGrid yields
Total mass of the gas in stars:
4.946E+10
sum msssfac 21733178007.1
END SFRMDOT Run time: 3.27s
################ Simulation successful ##############
Run time: 4.63s
In [112]:
c_2_8_10.plot_totmasses(source='agb',marker='o',color='r',label='[2.8,10]')
c_1_10.plot_totmasses(source='agb',marker='s',color='b',label='[1,10]')
c_0_1_10.plot_totmasses(source='agb',marker='p',color='k',label='[0.1,10]')
In [113]:
s1=s.sygma()
################ Start SYGMA simulation ##############
10 30
Use SNIa tables from yield_tables/sn1a_t86.txt
START SFRMDOT Run time: 2.65s
Use NuGrid yields
Total mass of the gas in stars:
1.000E+11
sum msssfac 36877281297.2
END SFRMDOT Run time: 3.32s
################ Simulation successful ##############
Run time: 4.7s
In [114]:
s1.plot_mass_range_contributions(specie='H',prodfac=False,rebin=0.5)
[[1.25, 1.75, 2.25, 2.75, 3.25, 3.75, 4.25, 4.75, 5.25, 5.75, 6.25, 6.75, 7.25, 7.75, 8.25, 8.75, 9.25, 9.75, 10.25, 10.75, 11.25, 11.75, 12.25, 12.75, 13.25, 13.75, 14.25, 14.75, 15.25, 15.75, 16.25, 16.75, 17.25, 17.75, 18.25, 18.75, 19.25, 19.75, 20.25, 20.75, 21.25, 21.75, 22.25, 22.75, 23.25, 23.75, 24.25, 24.75, 25.25, 25.75, 26.25, 26.75, 27.25, 27.75, 28.25, 28.75, 29.25, 29.75], [1.0, 1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, 5.0, 5.5, 6.0, 6.5, 7.0, 7.5, 8.0, 8.5, 9.0, 9.5, 10.0, 10.5, 11.0, 11.5, 12.0, 12.5, 13.0, 13.5, 14.0, 14.5, 15.0, 15.5, 16.0, 16.5, 17.0, 17.5, 18.0, 18.5, 19.0, 19.5, 20.0, 20.5, 21.0, 21.5, 22.0, 22.5, 23.0, 23.5, 24.0, 24.5, 25.0, 25.5, 26.0, 26.5, 27.0, 27.5, 28.0, 28.5, 29.0, 29.5, 30.0], [2556444798.9621968, 5453087152.6765184, 3763305672.9371104, 3345537525.4494596, 2666735318.4808583, 2188452042.2160754, 1847257096.901278, 1581572830.352963, 1381321735.7575526, 1214674076.8088381, 1085191905.2978191, 936172718.10445476, 849997729.38083661, 776752474.1610961, 815074188.51327407, 752793169.6240679, 698353751.31595373, 650419876.79241407, 607936808.63175774, 570061632.77221668, 536113242.15338749, 505535746.62491274, 477871306.27747971, 452739688.3636511, 429822693.81911004, 317159079.50523001, 302225071.60204148, 288473960.29612714, 275776527.81996149, 264021307.6938931, 253111677.21813169, 242963498.43876553, 233503191.27953658, 209145755.75527218, 201446396.49401817, 194227352.65117398, 187447012.80128813, 181068328.72019631, 175058215.05138102, 169387040.4676404, 164028194.56535253, 158957717.74309888, 154153983.69979832, 132997064.11698076, 129150107.8559899, 125492791.60529315, 122012062.26762006, 118696009.93550111, 115533747.18221213, 112515303.18821818, 109631530.63268493, 106874023.60177353, 104235045.03248471, 101707462.43300822, 99284690.806062073, 96960641.857151806, 94729678.700350568, 92586575.384394348]]
In [116]:
s1.plot_mass(specie='Ni',source='sn1a')
s1.plot_mass(specie='Ni',source='massive')
read BB abundance has problems with isotopes
In [118]:
reload(s)
s1=s.sygma(iniZ=0.0)
################ Start SYGMA simulation ##############
10 30
Use SNIa tables from yield_tables/sn1a_t86.txt
START SFRMDOT Run time: 2.55s
Taking POPIII yields from yield_tables/popIII_heger10.txt
Currently only PopIII massive stars between 10Msun and 100Msun with their yield contributions included
lower IMF minimum (imf_bdys) fixed to 10Msun for now
Correct lower imf_bdy of PopIII to 10Msun: [10, 30]
TRY POPIII stars
Total mass of the gas in stars:
1.000E+11
sum msssfac 6278650311.79
END SFRMDOT Run time: 2.69s
################ Simulation successful ##############
Run time: 4.08s
In [129]:
reload(sygma)
s1=s.sygma(iolevel=1,iniZ=0.02)
################ Start SYGMA simulation ##############
Take alpha-enhanced abundance
10 30
Use SNIa tables from yield_tables/sn1a_t86.txt
Use isotopes with care
['H-1', 'H-2', 'He-3', 'He-4', 'Li-7', 'B-11', 'C-12', 'C-13', 'N-14', 'N-15', 'O-16', 'O-17', 'O-18', 'F-19', 'Ne-20', 'Ne-21', 'Ne-22', 'Na-23', 'Mg-24', 'Mg-25', 'Mg-26', 'Al-27', 'Si-28', 'Si-29', 'Si-30', 'P-31', 'Pb-206', 'Pb-207', 'S-32', 'S-33', 'S-34', 'S-36', 'Cl-35', 'Cl-37', 'Ar-36', 'Ar-38', 'Ar-40', 'K-39', 'K-40', 'K-41', 'Ca-40', 'Ca-42', 'Ca-43', 'Ca-44', 'Ca-46', 'Ca-48', 'Sc-45', 'Ti-46', 'Ti-47', 'Ti-48', 'Ti-49', 'Ti-50', 'V-50', 'V-51', 'Cr-50', 'Cr-52', 'Cr-53', 'Cr-54', 'Mn-55', 'Fe-54', 'Fe-56', 'Fe-57', 'Fe-58', 'Co-59', 'Ni-58', 'Ni-60', 'Ni-61', 'Ni-62', 'Ni-64', 'Cu-63', 'Cu-65', 'Zn-64', 'Zn-66', 'Zn-67', 'Zn-68', 'Zn-70', 'Ga-69', 'Ga-71', 'Ge-70', 'Ge-72', 'Ge-73', 'Ge-74', 'Ge-76', 'As-75', 'Se-74', 'Se-76', 'Se-77', 'Se-78', 'Se-80', 'Se-82', 'Br-79', 'Br-81', 'Kr-78', 'Kr-80', 'Kr-82', 'Kr-83', 'Kr-84', 'Kr-86', 'Rb-85', 'Rb-87', 'Sr-84', 'Sr-86', 'Sr-87', 'Sr-88', 'Y-89', 'Zr-90', 'Zr-91', 'Zr-92', 'Zr-94', 'Zr-96', 'Nb-93', 'Mo-92', 'Mo-94', 'Mo-95', 'Mo-96', 'Mo-97', 'Mo-98', 'Mo-100', 'Ru-96', 'Ru-98', 'Ru-99', 'Ru-100', 'Ru-101', 'Ru-102', 'Ru-104', 'Rh-103', 'Pd-102', 'Pd-104', 'Pd-105', 'Pd-106', 'Pd-108', 'Pd-110', 'Ag-107', 'Ag-109', 'Cd-106', 'Cd-108', 'Cd-110', 'Cd-111', 'Cd-112', 'Cd-113', 'Cd-114', 'Cd-116', 'In-113', 'In-115', 'Sn-112', 'Sn-114', 'Sn-115', 'Sn-116', 'Sn-117', 'Sn-118', 'Sn-119', 'Sn-120', 'Sn-122', 'Sn-124', 'Sb-121', 'Sb-123', 'Te-120', 'Te-122', 'Te-123', 'Te-124', 'Te-125', 'Te-126', 'Te-128', 'Te-130', 'I-127', 'Xe-124', 'Xe-126', 'Xe-128', 'Xe-129', 'Xe-130', 'Xe-131', 'Xe-132', 'Xe-134', 'Xe-136', 'Cs-133', 'Ba-130', 'Ba-132', 'Ba-134', 'Ba-135', 'Ba-136', 'Ba-137', 'Ba-138', 'La-138', 'La-139', 'Ce-136', 'Ce-138', 'Ce-140', 'Ce-142', 'Pr-141', 'Nd-142', 'Nd-143', 'Nd-144', 'Nd-145', 'Nd-146', 'Nd-148', 'Nd-150', 'Sm-144', 'Sm-147', 'Sm-148', 'Sm-149', 'Sm-150', 'Sm-152', 'Sm-154', 'Eu-151', 'Eu-153', 'Gd-152', 'Gd-154', 'Gd-155', 'Gd-156', 'Gd-157', 'Gd-158', 'Gd-160', 'Tb-159', 'Dy-156', 'Dy-158', 'Dy-160', 'Dy-161', 'Dy-162', 'Dy-163', 'Dy-164', 'Ho-165', 'Er-162', 'Er-164', 'Er-166', 'Er-167', 'Er-168', 'Er-170', 'Tm-169', 'Yb-168', 'Yb-170', 'Yb-171', 'Yb-172', 'Yb-173', 'Yb-174', 'Yb-176', 'Lu-175', 'Lu-176', 'Hf-174', 'Hf-176', 'Hf-177', 'Hf-178', 'Hf-179', 'Hf-180', 'Ta-180', 'Ta-181', 'W-180', 'W-182', 'W-183', 'W-184', 'W-186', 'Re-185', 'Re-187', 'Os-184', 'Os-186', 'Os-187', 'Os-188', 'Os-189', 'Os-190', 'Os-192', 'Ir-191', 'Ir-193', 'Pt-190', 'Pt-192', 'Pt-194', 'Pt-195', 'Pt-196', 'Pt-198', 'Au-197', 'Hg-196', 'Hg-198', 'Hg-199', 'Hg-200', 'Hg-201', 'Hg-202', 'Hg-204', 'Tl-203', 'Tl-205', 'Pb-204', 'Pb-208', 'Bi-209']
Use initial abundance of iniab2.0E-02GN93.ppn
Number of timesteps: 2.0E+02
##### Start with initial metallicity of 2.0000E-02
#################################
################## Star formation at 7.000E+06 (Z=2.0000E-02) of 1.0
Mass locked away: 1.000E+11 , new ISM mass: 0.000E+00
START SFRMDOT Run time: 2.71s
Use NuGrid yields
Stars under consideration (take into account user-selected imf ends):
[1.0, 1.65, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 15.0, 20.0, 25.0]
with (corrected) mass bdys:####
[1, 1.325, 1.825, 2.5, 3.5, 4.5, 5.5, 6.5, 8, 17.5, 22.5, 30]
Scaling factor for gas mass (mfac) 1.0
Total mass of the gas in stars:
1.000E+11
sum msssfac 36877281254.1
Ejection of mass at 1.629E+09 excluded due to timelimit
END SFRMDOT Run time: 3.59s
time and metallicity and total mass:
1.124E+07 1.2081E-01 8.9902E+09
time and metallicity and total mass:
1.851E+07 1.1290E-01 1.4334E+10
time and metallicity and total mass:
3.048E+07 1.0958E-01 1.9096E+10
time and metallicity and total mass:
5.018E+07 9.9220E-02 2.4072E+10
time and metallicity and total mass:
8.263E+07 8.5561E-02 2.9313E+10
time and metallicity and total mass:
1.361E+08 7.6909E-02 3.4362E+10
time and metallicity and total mass:
2.240E+08 7.1738E-02 3.8964E+10
time and metallicity and total mass:
3.689E+08 6.8435E-02 4.3455E+10
time and metallicity and total mass:
6.073E+08 6.5984E-02 4.8249E+10
time and metallicity and total mass:
1.000E+09 6.3524E-02 5.3241E+10
Converting isotopes into elements
################ Simulation successful ##############
Run time: 5.28s
In [130]:
s1.plot_mass_range_contributions(specie='C',prodfac=True,rebin=0.5)
[[1, 1.325], [1.325, 1.825], [1.825, 2.5], [2.5, 3.5], [3.5, 4.5], [4.5, 5.5], [5.5, 6.5], [6.5, 8], [8, 17.5], [17.5, 22.5], [22.5, 30]]
[[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], array([ 5.28646453e+09, 2.83986470e+04, 2.23021828e+06,
2.31951860e+09, 1.65491813e-01, 2.25597799e+01,
1.07618336e+08, 7.81377811e+05, 2.03673425e+07,
1.59023264e+04, 1.08992187e+08, 2.85417565e+05,
1.23703842e+05, 1.29428222e+04, 1.43910902e+07,
3.71970184e+04, 7.29858395e+06, 4.43925636e+05,
4.66994885e+06, 7.52755913e+05, 1.01893956e+06,
5.22235148e+05, 5.88466220e+06, 3.12608368e+05,
2.21018295e+05, 5.72151738e+04, 2.42484719e+01,
2.56108742e+01, 3.15871265e+06, 2.69675522e+04,
1.50951889e+05, 9.74289403e+02, 5.34657052e+04,
2.17640911e+04, 6.42847826e+05, 1.25936351e+05,
2.61661390e+02, 3.08830278e+04, 1.14716567e+02,
2.59428882e+03, 5.67743966e+05, 4.40662739e+03,
9.33646308e+02, 1.35610552e+04, 3.04365262e+01,
1.32004193e+03, 4.69971563e+02, 2.57081886e+03,
2.32581542e+03, 2.34527831e+04, 2.03043743e+03,
1.96060000e+03, 7.85957315e+00, 3.26175148e+03,
6.39413198e+03, 1.30229635e+05, 1.50665670e+04,
4.60011142e+03, 8.74112760e+04, 6.30254191e+05,
1.03668514e+07, 2.74026050e+05, 7.23561578e+04,
3.97672649e+04, 4.43353198e+05, 1.84038803e+05,
1.06530704e+04, 2.94462085e+04, 7.81377811e+03,
6.78911417e+03, 2.90283288e+03, 9.17618045e+03,
5.55951744e+03, 8.49497928e+02, 3.87254278e+03,
1.23074161e+02, 3.90402687e+02, 2.67271282e+02,
4.49993478e+02, 5.77589899e+02, 1.63087574e+02,
7.72791242e+02, 1.53069910e+02, 1.16147661e+02,
9.20480235e+00, 1.13113740e+02, 8.75830074e+01,
2.90626751e+02, 5.83314278e+02, 1.02237419e+02,
1.07618336e+02, 1.11510914e+02, 2.99728514e+00,
2.11172362e+01, 1.19067095e+02, 1.08133530e+02,
5.46907224e+02, 1.70128561e+02, 1.07217629e+02,
4.33564509e+01, 2.47980123e+00, 5.89038658e+01,
4.31618220e+01, 4.95845758e+02, 1.22902429e+02,
1.46486873e+02, 3.36536275e+01, 5.42556696e+01,
5.28875428e+01, 7.87102191e+00, 1.46200654e+01,
7.81950249e+00, 5.12904409e+00, 1.03267807e+01,
1.24276280e+01, 6.39985636e+00, 1.75223259e+01,
5.65454214e+00, 2.25025361e+00, 7.77943183e-01,
5.35687440e+00, 7.28141081e+00, 7.71073928e+00,
1.64003475e+01, 8.32324789e+00, 8.60374249e+00,
3.25831685e-01, 5.11530558e+00, 7.73936118e+00,
1.09106675e+01, 1.12827521e+01, 4.16391370e+00,
6.02204731e+00, 6.38840760e+00, 4.87659895e-01,
3.76091738e-01, 7.02953811e+00, 5.90183534e+00,
1.25478400e+01, 5.84459154e+00, 1.58851533e+01,
3.48900935e+00, 2.00410529e-01, 5.30764474e+00,
9.38798249e-01, 6.54296585e-01, 3.39283977e-01,
1.94056467e+01, 9.30784118e+00, 3.19821087e+01,
1.10137063e+01, 4.45242243e+01, 5.57440082e+00,
6.30254191e+00, 5.69919230e+00, 3.82445799e+00,
1.25077694e-01, 4.75810430e+00, 1.66808421e+00,
9.12466103e+00, 1.05729291e+01, 3.10375860e+01,
4.49363796e+01, 4.85599119e+01, 1.18838120e+01,
1.69670610e-01, 1.49864257e-01, 4.62701601e+00,
3.82789262e+01, 9.74289403e+00, 3.17760310e+01,
4.45070511e+01, 1.43166733e+01, 1.16662856e+01,
1.20040240e+01, 1.39331399e-01, 1.35438821e-01,
5.75872585e+00, 9.87455476e+00, 1.79917250e+01,
2.12031019e+01, 1.46658605e+02, 1.23360380e-02,
1.87816894e+01, 6.43992702e-02, 8.83271768e-02,
4.33507265e+01, 4.22172993e+00, 6.43420264e+00,
1.03897489e+01, 3.87311522e+00, 8.08282395e+00,
2.61890365e+00, 6.09073986e+00, 1.80146225e+00,
1.66922908e+00, 2.84215446e-01, 1.53928567e+00,
1.33835994e+00, 1.36640940e+00, 9.39370687e-01,
2.80036649e+00, 2.29547621e+00, 1.59481215e+00,
1.76826085e+00, 2.66355381e-02, 3.28407656e-01,
1.72819019e+00, 2.50441606e+00, 1.87759650e+00,
3.17645822e+00, 2.60688246e+00, 2.18900275e+00,
7.70501490e-03, 1.33893238e-02, 4.31732707e-01,
2.73968806e+00, 3.87197034e+00, 3.67848631e+00,
4.55431639e+00, 3.34761717e+00, 1.23532111e-02,
1.46887580e-01, 3.25659954e+00, 2.20331370e+00,
2.82269157e+00, 1.52898178e+00, 1.48833869e+00,
1.21585822e-02, 3.82102336e-01, 1.48891113e+00,
2.40939136e+00, 1.70872730e+00, 3.78782196e+00,
1.32462143e+00, 1.50207720e+00, 8.36904293e-02,
9.45667505e-03, 3.68192094e-01, 1.18723632e+00,
1.98635971e+00, 9.42232877e-01, 2.79578698e+00,
1.86672018e-04, 1.66464958e+00, 6.47999767e-03,
1.72017606e+00, 9.65702833e-01, 2.14835965e+00,
1.74078383e+00, 9.04451972e-01, 1.40361787e+00,
5.45876836e-03, 6.10218862e-01, 5.07637980e-01,
4.08033776e+00, 4.65792766e+00, 7.86529753e+00,
1.19467802e+01, 1.04870634e+01, 1.78199936e+01,
7.74508556e-03, 5.41526307e-01, 1.94514418e+01,
1.98635971e+01, 1.55588637e+01, 4.28870517e+00,
8.56367184e+00, 2.31436666e-02, 2.28231013e+00,
2.93717916e+00, 4.63731989e+00, 2.48838780e+00,
6.13653490e+00, 1.11339183e+00, 3.10662079e+00,
6.92077490e+00, 3.02361729e+00, 7.08678191e+01,
6.89215300e+00]), array([ 6.09122104e+09, 2.39030093e+04, 1.30153958e+06,
2.82138411e+09, 1.39588837e-01, 2.65834624e+01,
1.55261000e+08, 9.11906714e+05, 2.82967417e+07,
1.73735993e+04, 1.41878473e+08, 2.56083933e+05,
1.38720355e+05, 2.46530625e+04, 1.70538398e+07,
4.83586894e+04, 1.28022228e+07, 6.28860344e+05,
5.55434089e+06, 1.15034466e+06, 1.84670979e+06,
6.22149343e+05, 6.93207009e+06, 3.71868452e+05,
2.72940391e+05, 6.95180833e+04, 3.77158301e+01,
3.66539127e+01, 3.71078923e+06, 3.23430810e+04,
1.80170660e+05, 1.52537123e+03, 6.26491756e+04,
2.82019981e+04, 7.53211260e+05, 1.50286964e+05,
4.24372172e+02, 3.65117974e+04, 2.12620327e+02,
3.11982630e+03, 6.65968237e+05, 5.48328323e+03,
1.17758343e+03, 1.61300902e+04, 4.94245544e+01,
1.55418906e+03, 6.05569221e+02, 3.06021682e+03,
2.73690444e+03, 2.74045732e+04, 2.51346755e+03,
2.56241839e+03, 9.15854362e+00, 3.82645532e+03,
7.45315964e+03, 1.52773982e+05, 1.77052018e+04,
6.07148280e+03, 1.01849321e+05, 7.35446843e+05,
1.21429656e+07, 3.37405484e+05, 1.21666515e+05,
5.60566032e+04, 5.17141903e+05, 2.23200024e+05,
1.60787708e+04, 4.19240230e+04, 1.30864535e+04,
1.12981689e+04, 4.99377486e+03, 1.14994990e+04,
7.87161034e+03, 1.28022228e+03, 5.72408976e+03,
1.44839209e+02, 6.13069752e+02, 4.12529228e+02,
7.35841608e+02, 8.59402994e+02, 2.42740882e+02,
1.14087030e+03, 1.79460083e+02, 1.58971790e+02,
1.07257599e+01, 1.80644378e+02, 1.20600650e+02,
4.39373235e+02, 7.58343203e+02, 1.19771644e+02,
1.28495946e+02, 1.42944338e+02, 3.49129999e+00,
2.61689593e+01, 1.74288664e+02, 1.39115120e+02,
7.18471957e+02, 2.30740032e+02, 1.46852510e+02,
6.80969299e+01, 2.88888889e+00, 8.89799885e+01,
6.28070815e+01, 7.79660502e+02, 2.00974765e+02,
2.24700131e+02, 5.56223619e+01, 9.55330843e+01,
8.67298290e+01, 1.56603201e+01, 1.70262063e+01,
9.11117184e+00, 5.99252984e+00, 1.61458808e+01,
2.11159698e+01, 9.68358082e+00, 2.81941028e+01,
7.28341077e+00, 2.62044882e+00, 9.05985242e-01,
6.27676050e+00, 1.21705991e+01, 9.84938204e+00,
2.45938477e+01, 1.02362516e+01, 1.10060429e+01,
3.79566366e-01, 8.75193587e+00, 9.79806262e+00,
1.62366767e+01, 1.76104582e+01, 5.40433027e+00,
7.01497070e+00, 8.75193587e+00, 5.68066563e-01,
4.35425587e-01, 1.18311014e+01, 8.17952689e+00,
1.96829735e+01, 8.28611339e+00, 2.60742158e+01,
5.51091677e+00, 2.33384956e-01, 7.72159971e+00,
1.09349853e+00, 7.61501321e-01, 3.95159576e-01,
3.19996356e+01, 1.40615226e+01, 5.31748201e+01,
1.76420394e+01, 7.70580912e+01, 9.72700495e+00,
7.40184021e+00, 8.67298290e+00, 4.98982721e+00,
1.45707692e-01, 8.52297228e+00, 2.94849838e+00,
1.66156509e+01, 1.42904862e+01, 4.76086362e+01,
5.38064438e+01, 5.68066563e+01, 1.53642465e+01,
1.97698217e-01, 1.74525523e-01, 8.91773709e+00,
4.61085300e+01, 1.88973915e+01, 3.99501989e+01,
6.75837357e+01, 1.90987216e+01, 1.37536060e+01,
1.62603626e+01, 1.62327290e-01, 1.57748019e-01,
1.14402842e+01, 1.30272388e+01, 3.63973156e+01,
3.85053597e+01, 3.07087547e+02, 1.43654915e-02,
3.70881540e+01, 7.50447906e-02, 1.02875710e-01,
8.57034405e+01, 6.20175519e+00, 1.12113206e+01,
2.01922201e+01, 5.89778628e+00, 1.37772919e+01,
3.86316844e+00, 1.08402417e+01, 2.71716620e+00,
1.95211199e+00, 3.31010294e-01, 2.22607877e+00,
2.31924327e+00, 1.75551912e+00, 1.81552337e+00,
3.89396010e+00, 2.93823449e+00, 1.94816434e+00,
2.18344417e+00, 4.57927181e-02, 6.35966111e-01,
2.16528499e+00, 3.45142875e+00, 2.43214600e+00,
4.62269594e+00, 3.22325468e+00, 2.75269503e+00,
8.98089946e-03, 1.56405818e-02, 8.04135921e-01,
3.35944854e+00, 5.17931432e+00, 4.57137651e+00,
6.60836294e+00, 4.14108287e+00, 1.43931250e-02,
1.71130546e-01, 4.18055935e+00, 2.75072121e+00,
4.08581580e+00, 2.23594789e+00, 1.89921351e+00,
1.41641614e-02, 6.78995475e-01, 2.09817497e+00,
3.66618080e+00, 2.40885488e+00, 6.28070815e+00,
1.89487109e+00, 2.01290577e+00, 2.22252589e-01,
1.10139382e-02, 5.76751389e-01, 1.61814096e+00,
3.30891865e+00, 1.44286539e+00, 5.09641371e+00,
2.17949652e-04, 2.32398044e+00, 7.67422793e-03,
2.69426984e+00, 1.60669278e+00, 3.76368771e+00,
2.62439647e+00, 1.33391030e+00, 1.79302177e+00,
6.35966111e-03, 1.15152895e+00, 7.09392366e-01,
5.74382800e+00, 5.63724151e+00, 1.02086180e+01,
1.45036592e+01, 1.24074580e+01, 2.10172786e+01,
9.02432359e-03, 8.68087820e-01, 2.34845586e+01,
2.34648204e+01, 1.95527011e+01, 5.13983784e+00,
1.04138957e+01, 2.69545413e-02, 3.99501989e+00,
4.08581580e+00, 7.64264675e+00, 3.88606480e+00,
1.07968176e+01, 1.41325802e+00, 5.28590082e+00,
1.07928699e+01, 5.58197443e+00, 8.96510886e+01,
8.20321278e+00]), array([ 4.17793514e+09, 1.28834482e+03, 5.80565076e+05,
1.90822002e+09, 5.07278365e-02, 1.79522811e+01,
6.05257364e+07, 6.52271481e+05, 2.05637375e+07,
1.09455974e+04, 7.57164320e+07, 1.22216949e+05,
9.34158641e+04, 4.19571358e+03, 1.18562490e+07,
4.00212604e+04, 2.70034862e+06, 3.92508611e+05,
3.66631093e+06, 1.07026253e+06, 2.08007834e+06,
4.51177487e+05, 4.69548549e+06, 2.53836721e+05,
1.92283785e+05, 5.05302982e+04, 2.84455158e+01,
2.99863146e+01, 2.50281032e+06, 2.10575832e+04,
1.22533010e+05, 1.16508092e+03, 4.23917201e+04,
1.76559736e+04, 5.09056210e+05, 1.01534688e+05,
5.79182308e+02, 2.45145036e+04, 7.58349550e+01,
2.00501379e+03, 4.49597181e+05, 3.54383718e+03,
7.83831991e+02, 1.10009082e+04, 1.12893141e+02,
1.04912593e+03, 4.47819336e+02, 2.06032451e+03,
1.87365082e+03, 1.86653944e+04, 1.56766398e+03,
1.79720349e+03, 6.25801348e+00, 2.58972717e+03,
5.09253748e+03, 1.03075487e+05, 1.19471166e+04,
3.86977538e+03, 6.79531767e+04, 5.00759601e+05,
8.19783963e+06, 2.07810296e+05, 6.76371154e+04,
3.81051389e+04, 3.52210797e+05, 1.48371020e+05,
1.05801516e+04, 3.09147446e+04, 1.78061027e+04,
6.77951460e+03, 5.24266660e+03, 7.33459724e+03,
7.14100970e+03, 1.30197496e+03, 7.20422196e+03,
1.29822174e+02, 8.84181450e+02, 5.78392155e+02,
1.20439104e+03, 1.34187770e+03, 3.95471685e+02,
2.16501981e+03, 1.66781590e+02, 2.51466261e+02,
7.37410490e+00, 4.37744882e+02, 2.07612758e+02,
1.09594251e+03, 1.45901792e+03, 9.36529101e+01,
1.16508092e+02, 2.65689019e+02, 2.40206578e+00,
2.95912380e+01, 4.94043299e+02, 2.30132124e+02,
1.37605183e+03, 1.26760330e+03, 4.07719060e+02,
6.08220439e+02, 1.98723534e+00, 1.91868955e+02,
9.57270623e+01, 1.63186393e+03, 4.18188590e+02,
2.81689622e+02, 8.82601143e+01, 1.64490146e+02,
1.26819591e+02, 1.11569634e+02, 1.17001938e+01,
6.23628427e+00, 3.97842145e+00, 2.50083493e+01,
2.83862543e+01, 1.26661561e+01, 3.23567742e+01,
6.95532369e+00, 1.80115426e+00, 6.20270275e-01,
4.38930112e+00, 1.24547901e+01, 7.78300919e+00,
2.18674903e+01, 8.17808580e+00, 8.44476251e+00,
2.61145638e-01, 7.85412298e+00, 7.33657262e+00,
1.36380445e+01, 1.47501852e+01, 4.82388539e+00,
4.82586077e+00, 6.81112073e+00, 3.90730766e-01,
2.84060082e-01, 9.65962308e+00, 6.29752114e+00,
1.59591196e+01, 6.35678263e+00, 2.08007834e+01,
7.24175423e+00, 1.60341842e-01, 5.94195219e+00,
7.51633248e-01, 5.19525740e-01, 2.69639785e-01,
2.29737048e+01, 1.06631177e+01, 4.00212604e+01,
1.29703651e+01, 5.50934331e+01, 1.61507318e+01,
6.17109663e+00, 6.24813656e+00, 4.29448273e+00,
1.00250689e-01, 5.89849377e+00, 1.99316149e+00,
1.22710795e+01, 1.00724781e+01, 3.19221900e+01,
3.81446466e+01, 3.87570153e+01, 1.04418748e+01,
1.36024876e-01, 1.20004520e-01, 5.53502329e+00,
3.15271134e+01, 1.10977019e+01, 2.66479172e+01,
4.06138754e+01, 1.81162379e+01, 1.00467982e+01,
1.04359486e+01, 1.11688157e-01, 1.08488037e-01,
5.52119560e+00, 8.09709509e+00, 1.70594080e+01,
2.47713034e+01, 1.58702274e+02, 9.84135832e-03,
2.01489070e+01, 5.15970051e-02, 7.06792053e-02,
4.38139959e+01, 6.81507150e+00, 6.65901624e+00,
9.43442941e+00, 3.74927702e+00, 8.87539601e+00,
2.55021951e+00, 6.83877609e+00, 2.06625066e+00,
1.36222415e+00, 2.26773973e-01, 1.48489543e+00,
1.30691342e+00, 1.19036582e+00, 1.12399295e+00,
2.59170255e+00, 2.09588141e+00, 1.30177743e+00,
1.48904374e+00, 2.71812707e-02, 3.86384923e-01,
1.48647574e+00, 2.31909969e+00, 1.65477837e+00,
3.09542523e+00, 2.31909969e+00, 1.88471296e+00,
6.17899816e-03, 1.07599114e-02, 5.14982359e-01,
2.31514893e+00, 3.53198488e+00, 2.95319765e+00,
4.16213207e+00, 2.86628079e+00, 9.90457058e-03,
2.45737651e-01, 2.86825618e+00, 1.89617018e+00,
2.74775781e+00, 1.66228483e+00, 1.30592573e+00,
9.74653993e-03, 4.21941818e-01, 1.45151146e+00,
2.40601655e+00, 1.62099932e+00, 4.07126445e+00,
1.57161475e+00, 1.35906353e+00, 1.43709116e-01,
7.57954473e-03, 2.96504995e-01, 1.09534990e+00,
2.17292135e+00, 9.64579540e-01, 3.29691430e+00,
1.51136557e-04, 1.57220736e+00, 5.43625413e-03,
1.68579189e+00, 1.05208901e+00, 2.38033657e+00,
1.84105700e+00, 8.97811593e-01, 1.24626916e+00,
4.37547344e-03, 6.62148396e-01, 4.51572564e-01,
3.84607078e+00, 3.86977538e+00, 6.96520061e+00,
1.01179119e+01, 8.53167936e+00, 1.44657300e+01,
6.20862890e-03, 5.52712175e-01, 1.63067870e+01,
1.61843133e+01, 1.34148262e+01, 3.66236016e+00,
7.15088661e+00, 1.85290929e-02, 2.43762267e+00,
2.69837324e+00, 4.76264852e+00, 2.47910572e+00,
6.87235761e+00, 1.15816708e+00, 3.53000950e+00,
7.02051133e+00, 3.60902482e+00, 7.31681879e+01,
6.05849979e+00]), array([ 3.04297187e+09, 1.06132638e+03, 2.89093933e+05,
1.62129356e+09, 9.16141111e-02, 1.36713229e+01,
2.63213585e+07, 7.74553557e+05, 1.94624860e+07,
6.00702609e+03, 4.53892562e+07, 7.18614867e+04,
5.98497557e+04, 2.41975452e+03, 8.83297441e+06,
3.50835391e+04, 1.06445988e+06, 3.31222033e+05,
2.79461337e+06, 5.46504750e+05, 9.67553641e+05,
3.30757812e+05, 3.57334492e+06, 1.90678977e+05,
1.37061395e+05, 3.62208818e+04, 1.53425203e+01,
1.67351847e+01, 1.91259254e+06, 1.57951362e+04,
9.18230108e+04, 6.24842126e+02, 3.24722933e+04,
1.21045754e+04, 3.90178163e+05, 7.62251687e+04,
2.80854002e+02, 1.86617039e+04, 2.71105350e+01,
1.46693991e+03, 3.44104180e+05, 2.54741543e+03,
5.53700183e+02, 8.27126641e+03, 5.09483086e+01,
7.99621517e+02, 2.96869643e+02, 1.55746310e+03,
1.42515998e+03, 1.43212330e+04, 1.12991511e+03,
1.19653089e+03, 4.81513741e+00, 1.97642300e+03,
3.91918993e+03, 7.89176534e+04, 9.13239727e+03,
2.61008533e+03, 5.21901011e+04, 3.85187782e+05,
6.28555898e+06, 1.52728870e+05, 3.43523903e+04,
2.36752960e+04, 2.70873240e+05, 1.10995358e+05,
6.31921504e+03, 1.93696417e+04, 8.69370796e+03,
4.54124673e+03, 2.67623689e+03, 5.49986411e+03,
4.25110830e+03, 7.18498812e+02, 3.72189580e+03,
8.40705119e+01, 4.26155328e+02, 2.76095731e+02,
5.52307518e+02, 6.45616038e+02, 1.87893648e+02,
9.79507345e+02, 1.05784472e+02, 1.23250806e+02,
5.67742883e+00, 1.82206935e+02, 9.83453228e+01,
4.57258168e+02, 6.81245037e+02, 6.56409187e+01,
7.49601652e+01, 1.25920079e+02, 1.84876209e+00,
1.70369287e+01, 2.01820293e+02, 1.13072750e+02,
6.39465103e+02, 4.95092220e+02, 1.69556899e+02,
2.23406592e+02, 1.52960981e+00, 8.00433905e+01,
4.46232908e+01, 6.85539086e+02, 1.73038561e+02,
1.35204509e+02, 3.86696502e+01, 6.93779017e+01,
5.72965375e+01, 4.27083771e+01, 9.00589691e+00,
4.80005021e+00, 3.06386184e+00, 1.15045691e+01,
1.27312744e+01, 6.18923302e+00, 1.59576137e+01,
4.25923217e+00, 1.38686170e+00, 4.77451803e-01,
3.33078919e+00, 6.24145794e+00, 5.18767516e+00,
1.23947138e+01, 5.62520391e+00, 5.70876378e+00,
2.01007905e-01, 4.05845638e+00, 5.05769314e+00,
7.93586638e+00, 8.30376191e+00, 3.05457741e+00,
3.71377192e+00, 4.34395260e+00, 3.00699470e-01,
2.21781817e-01, 5.19115682e+00, 3.99694703e+00,
9.07436958e+00, 3.97373596e+00, 1.15405463e+01,
3.62673039e+00, 1.23482916e-01, 3.64878091e+00,
5.78768143e-01, 4.00042869e-01, 2.07623062e-01,
1.31722848e+01, 6.41438044e+00, 2.26540087e+01,
7.60162691e+00, 3.10216011e+01, 7.41941997e+00,
4.25807162e+00, 3.85071726e+00, 2.75399399e+00,
7.71420062e-02, 3.26579819e+00, 1.12039857e+00,
6.58033963e+00, 6.81825314e+00, 2.00891850e+01,
2.81550334e+01, 2.96173311e+01, 7.44843381e+00,
1.04658735e-01, 9.23684711e-02, 3.01163692e+00,
2.36869015e+01, 6.01747107e+00, 1.95553303e+01,
2.67275523e+01, 1.09289344e+01, 7.36719505e+00,
7.34862619e+00, 8.59390034e-02, 8.35018406e-02,
3.05573796e+00, 5.87356241e+00, 9.59313710e+00,
1.37641672e+01, 8.22020204e+01, 7.59930580e-03,
1.07444064e+01, 3.97141485e-02, 5.44067587e-02,
2.33735520e+01, 3.22633936e+00, 3.78224459e+00,
5.38380874e+00, 2.32923133e+00, 4.89057340e+00,
1.59924303e+00, 3.63601482e+00, 1.17680148e+00,
1.03242859e+00, 1.74663336e-01, 9.41441183e-01,
7.47048433e-01, 8.37571624e-01, 5.50334577e-01,
1.70601398e+00, 1.44024717e+00, 9.70106860e-01,
1.08604618e+00, 1.84760153e-02, 1.92303752e-01,
1.06399566e+00, 1.53425203e+00, 1.15335829e+00,
1.93232195e+00, 1.63289909e+00, 1.34508177e+00,
4.74782529e-03, 8.25269755e-03, 2.54509432e-01,
1.68628456e+00, 2.37101126e+00, 2.25959810e+00,
2.74819122e+00, 2.05301954e+00, 7.60975078e-03,
9.05464017e-02, 1.98222576e+00, 1.34740288e+00,
1.68396346e+00, 9.71035303e-01, 9.08945678e-01,
7.49369541e-03, 2.12961609e-01, 9.06856682e-01,
1.42632053e+00, 1.02975932e+00, 2.20969429e+00,
8.56372594e-01, 9.10570454e-01, 4.55401282e-02,
5.82946137e-03, 2.08087283e-01, 7.22328639e-01,
1.17448037e+00, 5.64725443e-01, 1.61665134e+00,
1.15509912e-04, 1.00817302e+00, 4.03640586e-03,
9.92157380e-01, 5.65305720e-01, 1.23831083e+00,
1.08581407e+00, 5.46388694e-01, 8.70067128e-01,
3.36676636e-03, 3.30989923e-01, 2.98610474e-01,
2.46965833e+00, 2.86308604e+00, 4.81049519e+00,
7.38808502e+00, 6.45499982e+00, 1.09753566e+01,
4.77799969e-03, 3.24722933e-01, 1.19885200e+01,
1.22322363e+01, 9.42949902e+00, 2.66231025e+00,
5.22249177e+00, 1.42632053e-02, 1.17215926e+00,
1.71065619e+00, 2.55669986e+00, 1.41935721e+00,
3.47353730e+00, 7.11999711e-01, 1.82671156e+00,
4.15246123e+00, 1.74779391e+00, 4.59463220e+01,
4.29869100e+00]), array([ 2.34052928e+09, 2.55440004e+03, 1.77572974e+05,
1.36907408e+09, 2.73137056e-02, 1.08290898e+01,
1.07010685e+07, 5.13666351e+05, 1.73355804e+07,
5.42810008e+03, 3.11317504e+07, 6.07573688e+04,
4.66524417e+04, 1.63340026e+03, 6.93573830e+06,
2.82098542e+04, 5.01617295e+05, 2.71932150e+05,
2.19895286e+06, 3.35265005e+05, 5.13365125e+05,
2.61087999e+05, 2.83303447e+06, 1.49935450e+05,
1.04374954e+05, 2.76827079e+04, 1.10248869e+01,
1.21394247e+01, 1.51893422e+06, 1.24180591e+04,
7.20458288e+04, 3.73821986e+02, 2.58151041e+04,
9.10456852e+03, 3.10263212e+05, 5.99139348e+04,
1.48354011e+02, 1.47826865e+04, 1.19737502e+01,
1.13712973e+03, 2.73438282e+05, 1.95194720e+03,
4.19231869e+02, 6.50874985e+03, 2.47231584e+01,
6.33930999e+02, 2.16732409e+02, 1.22749766e+03,
1.12733987e+03, 1.13938893e+04, 8.75815814e+02,
8.74309682e+02, 3.83009392e+00, 1.57164884e+03,
3.11694037e+03, 6.27379325e+04, 7.23545859e+03,
1.94140427e+03, 4.18780029e+04, 3.06497882e+05,
4.99960549e+06, 1.19210355e+05, 1.96399625e+04,
1.63490639e+04, 2.15828729e+05, 8.66779021e+04,
4.27892129e+03, 1.33443304e+04, 4.60349275e+03,
2.75998707e+03, 1.49483611e+03, 4.29473567e+03,
2.81646702e+03, 4.43631209e+02, 2.17033635e+03,
6.21656023e+01, 2.30588824e+02, 1.50989743e+02,
2.79989957e+02, 3.50476939e+02, 1.00684931e+02,
5.00186469e+02, 7.76411096e+01, 6.94100977e+01,
4.49580431e+00, 8.30631851e+01, 5.36183026e+01,
2.11310333e+02, 3.69604816e+02, 5.03951799e+01,
5.44617366e+01, 6.96812014e+01, 1.46772573e+00,
1.16499318e+01, 9.05938456e+01, 6.45528217e+01,
3.44151184e+02, 2.00390875e+02, 7.88460152e+01,
7.99756143e+01, 1.20942407e+00, 3.69077670e+01,
2.33224555e+01, 3.14555688e+02, 7.77917228e+01,
7.40941685e+01, 1.87588753e+01, 3.13124863e+01,
2.85336726e+01, 1.40597431e+01, 7.13680694e+00,
3.81653873e+00, 2.49189555e+00, 5.87542131e+00,
6.26927485e+00, 3.35340311e+00, 8.63013691e+00,
2.90382268e+00, 1.10022950e+00, 3.80298354e-01,
2.62066985e+00, 3.41666066e+00, 3.71562788e+00,
7.77164162e+00, 4.15315926e+00, 4.13734487e+00,
1.58821630e-01, 2.30890050e+00, 3.70734416e+00,
5.10428167e+00, 5.10126941e+00, 2.11611560e+00,
2.91888400e+00, 2.99870900e+00, 2.38119484e-01,
2.10632574e-01, 3.03410311e+00, 2.75170334e+00,
5.63444017e+00, 2.69748258e+00, 6.95531802e+00,
1.97303305e+00, 9.82751193e-02, 2.43165027e+00,
4.58993756e-01, 3.20053070e-01, 1.65599224e-01,
8.20088926e+00, 4.20512081e+00, 1.39392526e+01,
4.86028828e+00, 1.90224484e+01, 3.48518967e+00,
3.16212434e+00, 2.58376961e+00, 1.91730616e+00,
6.08326754e-02, 1.97755144e+00, 6.87775222e-01,
3.81503260e+00, 4.93107648e+00, 1.37509860e+01,
2.18765687e+01, 2.35107220e+01, 5.59076234e+00,
8.23101191e-02, 7.30172840e-02, 1.82392597e+00,
1.84049342e+01, 3.73144227e+00, 1.49935450e+01,
1.91805922e+01, 7.34164090e+00, 5.69619159e+00,
5.49211069e+00, 6.75650858e-02, 6.59911778e-02,
1.90902243e+00, 4.46643473e+00, 6.00344254e+00,
8.59248361e+00, 5.34676894e+01, 5.68941399e-03,
7.22265647e+00, 3.13953235e-02, 4.33766044e-02,
1.62662266e+01, 2.14849744e+00, 2.70651938e+00,
3.87603095e+00, 1.72678045e+00, 3.49121420e+00,
1.19059742e+00, 2.56795522e+00, 8.53223833e-01,
8.15570530e-01, 1.38865379e-01, 6.98167533e-01,
5.51771494e-01, 6.38825929e-01, 3.85946350e-01,
1.29150827e+00, 1.10851322e+00, 7.49827864e-01,
8.41174776e-01, 2.52578353e-02, 1.52269955e-01,
8.20841993e-01, 1.18532596e+00, 8.89371003e-01,
1.46772573e+00, 1.26741016e+00, 1.03621888e+00,
3.70207269e-03, 6.45076377e-03, 1.91278776e-01,
1.30280426e+00, 1.83748116e+00, 1.73657031e+00,
2.07394390e+00, 1.59273469e+00, 5.93039513e-03,
8.10299068e-02, 1.55809365e+00, 1.04600874e+00,
1.28849601e+00, 7.14132534e-01, 6.99899585e-01,
5.86111305e-03, 1.63867172e-01, 6.79416189e-01,
1.06558846e+00, 7.74904964e-01, 1.59198163e+00,
6.41311047e-01, 6.90335646e-01, 2.40152763e-02,
4.56960478e-03, 1.57993257e-01, 5.44240833e-01,
8.45693172e-01, 4.13885100e-01, 1.11303162e+00,
9.01420060e-05, 7.53819114e-01, 3.14405075e-03,
7.39284940e-01, 4.03191562e-01, 8.72050484e-01,
7.99003077e-01, 4.02589109e-01, 6.65409160e-01,
2.65003942e-03, 2.35634366e-01, 2.26597574e-01,
1.88266512e+00, 2.22079178e+00, 3.75779958e+00,
5.80162084e+00, 5.02972814e+00, 8.59248361e+00,
3.76608331e-03, 3.25776372e-01, 9.54134683e+00,
9.63924542e+00, 7.49752558e+00, 2.09126442e+00,
4.09216091e+00, 1.13035214e-02, 8.68285154e-01,
1.29828587e+00, 1.83898729e+00, 1.03621888e+00,
2.41583588e+00, 5.40400196e-01, 1.28021228e+00,
3.01904179e+00, 1.18607903e+00, 3.44151184e+01,
3.33984792e+00]), array([ 2.60730909e+09, 4.55242135e+03, 1.81412826e+05,
1.62274610e+09, 2.73829308e-02, 1.23270489e+01,
8.46302114e+06, 4.96065478e+05, 1.95238914e+07,
6.70638030e+03, 3.35610090e+07, 6.16061367e+04,
9.69281527e+04, 1.87452644e+03, 7.82265497e+06,
3.25495216e+04, 5.14257699e+05, 3.34664095e+05,
2.51853106e+06, 3.00389951e+05, 4.43817420e+05,
2.84890179e+05, 3.22147847e+06, 1.69405961e+05,
1.15920831e+05, 3.08103452e+04, 1.21232960e+01,
1.33749208e+01, 1.73117174e+06, 1.41098865e+04,
8.16466872e+04, 3.59551053e+02, 2.94204596e+04,
1.02385819e+04, 3.53802311e+05, 6.79952447e+04,
1.20869115e+02, 1.68387196e+04, 1.46410993e+01,
1.28655386e+03, 3.11741897e+05, 2.19470952e+03,
4.67831152e+02, 7.38604167e+03, 1.61619690e+01,
7.22012861e+02, 2.36717178e+02, 1.39570718e+03,
1.27927697e+03, 1.29892457e+04, 9.96933703e+02,
9.59821572e+02, 4.35885612e+00, 1.79520835e+03,
3.54748307e+03, 7.15318124e+04, 8.22288383e+03,
2.17724499e+03, 4.82093853e+04, 3.49217872e+05,
5.70216970e+06, 1.35859505e+05, 1.93055848e+04,
1.74135938e+04, 2.46104364e+05, 9.82379926e+04,
4.57279663e+03, 1.41972092e+04, 3.87639842e+03,
2.92239836e+03, 1.36659963e+03, 4.87551519e+03,
2.91584916e+03, 4.39014674e+02, 2.04043949e+03,
6.80752904e+01, 2.04407794e+02, 1.36878270e+02,
2.31041205e+02, 3.05774848e+02, 8.65222024e+01,
4.10707578e+02, 8.47029803e+01, 6.24648095e+01,
5.08945571e+00, 5.95176698e+01, 4.68704378e+01,
1.52159735e+02, 3.17345101e+02, 5.62576238e+01,
5.89282418e+01, 6.09293861e+01, 1.66640743e+00,
1.24434791e+01, 6.03254044e+01, 5.74874179e+01,
2.90420614e+02, 9.32897085e+01, 5.70580815e+01,
2.42902533e+01, 1.36587194e+00, 2.65169811e+01,
1.95020608e+01, 2.20271410e+02, 5.38708044e+01,
6.65689746e+01, 1.47356989e+01, 2.28203218e+01,
2.34606880e+01, 4.00156090e+00, 8.08462295e+00,
4.34066390e+00, 2.91075534e+00, 4.86969368e+00,
5.19860903e+00, 2.96897044e+00, 7.61890209e+00,
3.06575306e+00, 1.24871404e+00, 4.33848083e-01,
2.95951049e+00, 3.05119928e+00, 3.98191330e+00,
7.73533231e+00, 4.54878290e+00, 4.47237558e+00,
1.79448067e-01, 2.17142348e+00, 4.03285152e+00,
5.20224748e+00, 5.03051291e+00, 2.24637543e+00,
3.26877824e+00, 3.16617412e+00, 2.69463175e-01,
2.81542810e-01, 2.91730454e+00, 2.91148303e+00,
5.63012851e+00, 2.83507570e+00, 6.82499358e+00,
1.79375298e+00, 1.11991312e-01, 2.52580794e+00,
5.20588592e-01, 3.66682404e-01, 1.89053559e-01,
8.22288383e+00, 4.34648541e+00, 1.38915798e+01,
4.96356554e+00, 1.90145092e+01, 2.75502993e+00,
3.47398649e+00, 2.71064091e+00, 2.04262256e+00,
6.85410113e-02, 1.97640287e+00, 6.93123615e-01,
3.68865470e+00, 5.36015595e+00, 1.46410993e+01,
2.46977590e+01, 2.67789491e+01, 6.17807820e+00,
9.24892508e-02, 8.25199138e-02, 1.86433879e+00,
2.06299785e+01, 3.94625655e+00, 1.67513970e+01,
2.10447611e+01, 7.72805542e+00, 6.38328645e+00,
6.07692945e+00, 7.58979454e-02, 7.45153366e-02,
2.00914887e+00, 4.97157012e+00, 6.20791344e+00,
8.78320423e+00, 5.66214682e+01, 6.33234824e-03,
7.75716297e+00, 3.54602769e-02, 4.93081954e-02,
1.77774382e+01, 2.26529534e+00, 2.97260889e+00,
4.30937328e+00, 1.91454932e+00, 3.83564785e+00,
1.32002755e+00, 2.81833885e+00, 9.36535530e-01,
9.23437131e-01, 1.58054015e-01, 7.69167098e-01,
6.24793633e-01, 7.09205538e-01, 4.30282408e-01,
1.43936851e+00, 1.24289253e+00, 8.35386782e-01,
9.41629351e-01, 3.90768904e-02, 1.88325870e-01,
9.19070998e-01, 1.34840741e+00, 9.99116769e-01,
1.66422436e+00, 1.42627012e+00, 1.15848062e+00,
4.12890645e-03, 7.20994097e-03, 2.26456765e-01,
1.45610536e+00, 2.08846695e+00, 1.95020608e+00,
2.35261800e+00, 1.78938684e+00, 6.61323613e-03,
9.81652237e-02, 1.78065458e+00, 1.17594516e+00,
1.46847607e+00, 7.85176252e-01, 7.83720874e-01,
6.55211026e-03, 1.96621523e-01, 7.55341010e-01,
1.19632044e+00, 8.63766646e-01, 1.76100698e+00,
7.16991808e-01, 7.66256342e-01, 2.25146925e-02,
5.11783557e-03, 1.81994977e-01, 6.03181275e-01,
9.37990907e-01, 4.57716277e-01, 1.21014653e+00,
1.00493828e-04, 8.31748338e-01, 3.51182631e-03,
8.32476026e-01, 4.42216505e-01, 9.52544684e-01,
8.80503489e-01, 4.42070967e-01, 7.32782656e-01,
2.98061346e-03, 2.64878736e-01, 2.58766149e-01,
2.11902989e+00, 2.47778048e+00, 4.27298883e+00,
6.56302560e+00, 5.62066856e+00, 9.64915394e+00,
4.24315359e-03, 4.57497970e-01, 1.09298863e+01,
1.08862250e+01, 8.67405090e+00, 2.36789947e+00,
4.62737330e+00, 1.28146004e-02, 1.02822432e+00,
1.45901611e+00, 2.03389029e+00, 1.14392685e+00,
2.63205051e+00, 6.04345577e-01, 1.39352412e+00,
3.35610090e+00, 1.27782159e+00, 3.82691558e+01,
3.76433434e+00]), array([ 9.21196949e+09, 1.96477017e+04, 4.33746126e+05,
8.02342293e+09, 2.18927341e+01, 4.40642631e+01,
2.60012901e+08, 1.43740766e+06, 7.29562158e+07,
7.81065842e+04, 4.39615492e+08, 1.13366799e+05,
7.18410363e+06, 2.23329365e+04, 4.64560296e+07,
1.99264966e+06, 2.33013819e+07, 1.91194588e+06,
2.25970580e+07, 3.74759001e+06, 7.15915883e+06,
1.87232766e+06, 1.26426137e+08, 5.34259014e+06,
5.11368488e+06, 2.61186774e+06, 5.03298110e+01,
5.26041902e+01, 9.63896584e+07, 9.04909459e+05,
4.26115951e+06, 4.32865721e+04, 2.79968745e+06,
2.17900202e+05, 3.92220364e+07, 4.22594331e+06,
1.48788421e+04, 1.23682209e+06, 4.16578232e+03,
7.85908069e+04, 2.96989905e+07, 5.41155519e+05,
2.03666990e+05, 2.70431025e+05, 4.46365263e+03,
4.33452658e+03, 2.17753468e+04, 2.65442064e+05,
5.26482105e+04, 7.15622415e+05, 8.91996854e+04,
5.57149540e+03, 2.07922280e+02, 4.43870782e+05,
6.52379999e+05, 5.92512469e+06, 7.51572280e+05,
1.29566248e+04, 3.46439311e+06, 1.67276923e+07,
5.12689095e+08, 1.50695965e+07, 2.40643994e+05,
9.48782968e+06, 9.77102657e+07, 2.03666990e+07,
6.43282482e+05, 9.58614155e+06, 2.24796707e+05,
3.16652280e+05, 7.17676693e+04, 4.51647692e+05,
2.03520256e+05, 3.10489446e+04, 5.33965546e+04,
4.40936099e+03, 1.76227706e+04, 4.98162415e+03,
9.01681308e+03, 7.90603562e+03, 1.63608569e+03,
3.71530850e+03, 1.47761282e+03, 1.69184467e+03,
1.05531196e+03, 1.24577287e+03, 8.16282037e+02,
1.32266156e+03, 1.45604290e+03, 4.15844561e+02,
4.15257624e+02, 3.91780161e+02, 1.73586491e+02,
4.27876761e+02, 3.11810054e+02, 2.79528542e+02,
1.26382117e+03, 5.14009702e+02, 2.86278313e+02,
1.31312384e+02, 6.66319743e+01, 1.58032672e+02,
7.68740174e+01, 9.48049297e+02, 2.62654116e+02,
3.93981174e+02, 8.42987651e+01, 1.20395364e+02,
1.10754931e+02, 4.51354223e+01, 4.14964156e+01,
1.72852820e+01, 1.11488602e+01, 3.34553846e+01,
1.97357422e+01, 1.85178488e+01, 3.50988070e+01,
1.61701025e+01, 4.72777408e+00, 1.74907098e+00,
1.41495734e+01, 1.16110727e+01, 1.74760364e+01,
3.40129743e+01, 2.22008758e+01, 2.32720351e+01,
8.82312401e-01, 8.21124263e+00, 1.56565330e+01,
2.14525317e+01, 2.13938380e+01, 1.10476136e+01,
1.30549367e+01, 1.24929449e+01, 1.53043711e+00,
1.41877243e+00, 1.19089430e+01, 1.16052034e+01,
2.31253009e+01, 1.09346283e+01, 2.70284291e+01,
9.07990876e+00, 4.24942078e-01, 9.76075518e+00,
3.78720823e+00, 2.10416761e+00, 7.06818366e-01,
3.09902510e+01, 1.66689986e+01, 5.41008785e+01,
1.87379500e+01, 7.33964183e+01, 1.35450287e+01,
1.76961376e+01, 1.07365372e+01, 8.38292158e+00,
4.03078690e-01, 7.59936126e+00, 2.57811889e+00,
1.38781152e+01, 2.04547395e+01, 5.90458191e+01,
9.58320687e+01, 1.05721950e+02, 2.40790728e+01,
7.69473845e-01, 7.29268690e-01, 7.72115060e+00,
7.74756274e+01, 1.59500013e+01, 6.50325721e+01,
8.36531348e+01, 4.32425519e+01, 3.97796261e+01,
2.55170674e+01, 1.41613121e+00, 1.27409256e+00,
8.49150485e+00, 2.10270027e+01, 2.47980701e+01,
3.48346855e+01, 2.28171592e+02, 3.09609041e-02,
3.12396990e+01, 2.61920445e-01, 3.60965991e-01,
7.07845505e+01, 1.07600147e+01, 1.23843617e+01,
1.68744264e+01, 8.05863913e+00, 1.83270944e+01,
6.68080552e+00, 1.28216294e+01, 4.53555236e+00,
4.13350080e+00, 1.54070850e+00, 4.54288906e+00,
2.31986680e+00, 2.89799932e+00, 1.57739204e+00,
6.15402995e+00, 5.70942550e+00, 3.38955870e+00,
3.65221282e+00, 1.03946467e-01, 6.42842280e-01,
3.57297638e+00, 5.92805937e+00, 3.74318798e+00,
6.44309621e+00, 5.55535465e+00, 4.39322024e+00,
1.90314183e-02, 3.36167921e-02, 8.97279284e-01,
5.65660121e+00, 9.22224088e+00, 7.56121038e+00,
9.50543777e+00, 1.15112935e+01, 5.76224979e-02,
3.51428272e-01, 7.35138056e+00, 4.52087894e+00,
5.96767759e+00, 3.44531767e+00, 3.08728637e+00,
7.41594358e-02, 9.87080579e-01, 3.05940688e+00,
5.54214858e+00, 3.43211160e+00, 7.29855627e+00,
3.19000027e+00, 3.02859271e+00, 7.32203373e-01,
5.92072266e-02, 8.03956369e-01, 2.49741511e+00,
4.69402523e+00, 1.96330283e+00, 5.38514304e+00,
2.26850985e-03, 3.36461390e+00, 7.59495923e-02,
3.89579149e+00, 1.93982537e+00, 4.69549257e+00,
3.83563049e+00, 1.88993576e+00, 2.92147678e+00,
4.77619635e-02, 1.02860634e+00, 9.46435221e-01,
8.37265019e+00, 9.21637151e+00, 1.63315101e+01,
2.46660094e+01, 2.07188610e+01, 3.55830296e+01,
1.93982537e-02, 1.50109028e+00, 4.03078690e+01,
3.99850539e+01, 3.32499568e+01, 9.95591159e+00,
1.79455857e+01, 1.97504156e-01, 4.39908960e+00,
6.65879540e+00, 1.04093201e+01, 5.61551565e+00,
1.41275633e+01, 8.82018933e+00, 7.75049743e+00,
1.45795044e+01, 5.53921389e+00, 1.58619608e+02,
1.61554291e+01]), array([ 1.73778690e+09, 3.18271870e+03, 7.96241992e+04,
1.46022684e+09, 1.77737407e-02, 8.29531205e+00,
6.27771582e+07, 2.94879450e+05, 1.53422786e+07,
9.22201177e+03, 2.81383823e+08, 1.59788223e+04,
1.25779243e+05, 1.78682101e+03, 2.15412699e+07,
1.06660439e+05, 3.43013852e+06, 5.22505691e+05,
2.93529887e+07, 1.51465920e+06, 2.48319536e+06,
8.90036599e+05, 8.73841847e+07, 4.39282658e+06,
5.13958461e+06, 1.78727087e+06, 1.08572319e+01,
1.04411167e+01, 3.92272891e+07, 6.70957588e+05,
4.95964292e+06, 1.98138297e+04, 1.43660949e+06,
7.22915752e+04, 1.03196561e+07, 1.78772072e+06,
5.11709190e+03, 5.60068520e+05, 4.32984699e+03,
1.39409827e+04, 5.17557295e+06, 5.72889365e+04,
6.69383098e+03, 1.18806503e+04, 1.56121911e+03,
9.67636455e+02, 3.07700295e+03, 1.81921052e+04,
4.07792862e+03, 3.28168663e+04, 5.75363564e+03,
1.87589215e+03, 2.25826825e+02, 1.35990935e+04,
1.53265337e+05, 3.05226097e+05, 5.50171727e+04,
3.72479305e+03, 2.31450003e+05, 4.26236885e+06,
6.07977996e+06, 1.53287830e+05, 9.53915900e+04,
7.06721000e+04, 3.13548400e+05, 3.09949566e+05,
8.28856424e+04, 2.29875513e+05, 1.02881663e+05,
7.22465898e+04, 3.63932074e+04, 1.00834826e+04,
6.70282807e+04, 1.72474113e+04, 2.50568808e+04,
5.25429744e+03, 1.06975337e+04, 4.63574787e+03,
4.62675078e+03, 7.24040388e+03, 1.85160002e+03,
4.14540675e+03, 1.44155789e+03, 1.89726023e+03,
2.98478283e+02, 7.94667502e+02, 1.08639797e+03,
1.68852786e+03, 1.23372523e+03, 5.04511522e+02,
5.50171727e+02, 4.30060646e+02, 3.32217351e+01,
2.09969463e+02, 1.88826314e+02, 3.50211520e+02,
1.04703573e+03, 1.04096269e+03, 2.77335134e+02,
1.93774711e+02, 1.16984593e+01, 4.93040239e+01,
1.84417743e+01, 5.26554379e+02, 1.03623923e+02,
1.47079841e+02, 2.49894026e+01, 4.46255399e+01,
2.90605834e+01, 2.02861766e+01, 1.11653821e+01,
4.31860063e+00, 2.29425659e+00, 8.38078436e+00,
3.99920413e+00, 4.03519247e+00, 7.37985869e+00,
4.08692570e+00, 1.27556168e+00, 5.80761814e-01,
2.90830761e+00, 2.60465601e+00, 3.69330325e+00,
7.24715169e+00, 5.55120123e+00, 5.38025662e+00,
6.06178579e-01, 1.67368267e+00, 3.13998254e+00,
4.48729597e+00, 4.44231055e+00, 2.42021577e+00,
2.57541548e+00, 2.53942714e+00, 8.63045345e-01,
5.73789074e-01, 2.86782073e+00, 2.35273764e+00,
4.77295341e+00, 2.17864405e+00, 5.34201901e+00,
2.22610367e+00, 9.01282955e-02, 1.88758836e+00,
2.34823909e+00, 1.28860745e+00, 1.35676037e-01,
6.29121144e+00, 3.22095631e+00, 1.09876896e+01,
3.54485135e+00, 1.43660949e+01, 3.39190091e+00,
4.57726682e+00, 2.13410848e+00, 1.72901474e+00,
3.22770412e-01, 1.87364288e+00, 5.16207732e-01,
2.74411082e+00, 4.01269976e+00, 1.23462494e+01,
1.89141212e+01, 2.07000425e+01, 4.83143446e+00,
9.81581936e-01, 8.59446512e-01, 3.07475368e+00,
1.52095716e+01, 2.86782073e+00, 1.27511182e+01,
1.63679463e+01, 1.01082246e+01, 1.03174068e+01,
5.15982805e+00, 1.67773136e+00, 1.18244185e+00,
2.54392569e+00, 4.38832804e+00, 5.33752047e+00,
7.38210796e+00, 5.08335283e+01, 1.31964739e-02,
6.85127996e+00, 2.46070265e-01, 3.53360500e-01,
1.57718894e+01, 2.56866767e+00, 2.54617496e+00,
4.18814291e+00, 1.69122699e+00, 4.06218372e+00,
1.42221416e+00, 2.55517204e+00, 9.00833101e-01,
8.96334559e-01, 2.25152043e+00, 9.80907154e-01,
4.30735428e-01, 6.02354818e-01, 2.89706126e-01,
1.30637669e+00, 1.30862596e+00, 7.06945927e-01,
7.19091991e-01, 1.52410614e-02, 1.16287319e-01,
7.52156277e-01, 1.23484987e+00, 7.24940096e-01,
1.28365905e+00, 1.20448471e+00, 8.65969398e-01,
1.82753282e-02, 2.70362394e-02, 3.01402336e-01,
1.14465410e+00, 1.77939842e+00, 1.53805162e+00,
1.82505862e+00, 2.01624667e+00, 4.90566041e-02,
2.06168195e-01, 1.56639244e+00, 8.88687037e-01,
1.15792480e+00, 7.16842720e-01, 6.21023768e-01,
9.35696804e-02, 3.83500733e-01, 6.17200007e-01,
1.21843019e+00, 6.86702486e-01, 1.40579448e+00,
6.58136742e-01, 6.16525226e-01, 2.58441257e-01,
7.22016044e-02, 2.97128721e-01, 5.06535866e-01,
1.01059753e+00, 4.17689655e-01, 1.17546911e+00,
1.03443981e-03, 6.80404527e-01, 1.21775541e-01,
9.59764005e-01, 4.11166769e-01, 1.00857319e+00,
7.62053070e-01, 3.89123911e-01, 5.69290531e-01,
7.86570126e-02, 2.19641329e-01, 1.84147830e-01,
1.62959696e+00, 1.78884536e+00, 3.05226097e+00,
4.68298256e+00, 4.00370267e+00, 6.84678142e+00,
1.80256591e-02, 2.59116038e-01, 7.62952779e+00,
7.66551613e+00, 6.13151319e+00, 1.89973442e+00,
3.46162832e+00, 5.53320706e-01, 1.44043325e+00,
1.37430468e+00, 2.73736301e+00, 1.26161620e+00,
3.42563998e+00, 2.55292277e+00, 1.69819973e+00,
3.10624348e+00, 1.56954142e+00, 3.42788925e+01,
3.23895048e+00]), array([ 1.64107819e+09, 2.86565909e+03, 7.95359288e+04,
1.42763229e+09, 4.66498474e-01, 8.89985180e+00,
8.10771839e+07, 2.26528648e+06, 1.64089897e+07,
2.43374924e+04, 1.46705974e+08, 1.91223155e+04,
3.51800424e+04, 2.56636886e+03, 3.15240421e+07,
7.27078106e+04, 8.28514193e+06, 1.00163656e+06,
7.47687912e+06, 1.29339973e+06, 2.30650609e+06,
5.78149857e+05, 5.26535734e+06, 1.13999108e+06,
7.19551046e+05, 2.32084335e+05, 1.03801735e+01,
1.15719579e+01, 2.27424727e+06, 6.01806330e+04,
3.09326303e+05, 4.18289450e+03, 1.20522560e+05,
1.82979232e+04, 3.30832187e+05, 8.75827139e+04,
1.61527113e+03, 2.76350614e+04, 6.14172213e+02,
2.38894532e+03, 2.33697276e+05, 3.10760028e+03,
8.31381645e+02, 6.42846726e+03, 5.19187891e+02,
5.93741623e+02, 8.13280859e+02, 1.35576679e+03,
1.41813386e+03, 9.85865580e+03, 1.21633697e+03,
1.50111073e+03, 4.26533372e+01, 1.72118761e+03,
2.60758847e+03, 5.33525147e+04, 6.25462802e+03,
3.08430224e+03, 3.89077291e+04, 2.53769435e+05,
4.21336117e+06, 1.18802089e+05, 5.41231422e+04,
5.15424361e+04, 1.79197781e+05, 1.50935465e+05,
2.78321986e+04, 4.39257687e+04, 3.07534146e+04,
1.86563546e+04, 1.37601817e+04, 4.87466711e+03,
1.63050446e+04, 7.08618888e+03, 8.14714584e+03,
2.99469439e+03, 3.93199252e+03, 1.88893351e+03,
1.01113499e+03, 2.22048256e+03, 8.75110276e+02,
1.49358367e+03, 4.37465530e+02, 4.92126320e+02,
5.09510243e+00, 1.64430407e+02, 2.08786294e+02,
6.16681233e+02, 4.56103963e+02, 1.25200090e+02,
1.46526759e+02, 1.17511736e+02, 1.16113854e+00,
2.99290223e+01, 1.35003189e+02, 1.02762284e+02,
3.82804741e+02, 2.02513744e+02, 9.61492245e+01,
4.43558864e+01, 9.47154989e-01, 4.24562000e+01,
2.57712180e+01, 2.75992182e+02, 6.65427904e+01,
7.34963597e+01, 1.69735192e+01, 2.71153358e+01,
2.39432179e+01, 8.85684003e+00, 1.05020402e+01,
3.11835323e+00, 2.38715316e+00, 6.51807511e+00,
4.51085924e+00, 3.33878854e+00, 6.95356926e+00,
2.71691005e+00, 8.77081649e-01, 3.28860814e-01,
2.27424727e+00, 2.50185121e+00, 2.98214929e+00,
6.27075744e+00, 3.46603169e+00, 3.49829052e+00,
1.24752051e-01, 1.84592174e+00, 2.91225517e+00,
4.36927883e+00, 4.15601215e+00, 1.72387584e+00,
2.45346297e+00, 2.38894532e+00, 1.87459625e-01,
3.05383557e-01, 2.42478846e+00, 2.18284726e+00,
4.60584356e+00, 2.12549823e+00, 5.52163580e+00,
1.49931857e+00, 8.05574584e-02, 1.87997272e+00,
3.65420818e-01, 2.85669830e-01, 1.39429817e-01,
6.47506334e+00, 3.27964736e+00, 1.09644167e+01,
3.76352976e+00, 1.49340445e+01, 2.61654926e+00,
3.36208658e+00, 2.04485117e+00, 1.55039504e+00,
6.34602803e-02, 1.61975152e+00, 5.44098873e-01,
2.98752576e+00, 3.94812193e+00, 1.17493815e+01,
1.83696095e+01, 1.99646293e+01, 4.57896120e+00,
6.53958099e-02, 1.18425736e-01, 1.69484290e+00,
1.44824209e+01, 4.01801605e+00, 1.21293188e+01,
1.69824800e+01, 6.23491430e+00, 5.96071427e+00,
4.65243964e+00, 1.34949424e-01, 4.41049844e-01,
2.18643157e+00, 3.82087878e+00, 5.64171032e+00,
6.90338887e+00, 5.10227106e+01, 6.83707906e-03,
6.66323983e+00, 3.41764345e-02, 8.76723218e-02,
1.51849465e+01, 2.16313353e+00, 2.45167081e+00,
3.60402778e+00, 1.47386994e+00, 3.19183166e+00,
1.01507774e+00, 2.31188256e+00, 7.43565950e-01,
7.27257321e-01, 1.78713899e-01, 5.82630249e-01,
5.04850634e-01, 5.21696911e-01, 3.51621209e-01,
1.10343108e+00, 9.49305577e-01, 6.12559272e-01,
6.69012218e-01, 4.00905527e-02, 1.60057544e-01,
6.65069473e-01, 1.05683500e+00, 7.36576538e-01,
1.33497777e+00, 1.04034715e+00, 8.42313803e-01,
1.19375580e-02, 1.39465660e-02, 3.21692186e-01,
1.04948715e+00, 1.58104093e+00, 1.41956758e+00,
1.91760802e+00, 1.69251310e+00, 2.51618847e-02,
8.02707132e-02, 1.41275738e+00, 8.59518510e-01,
1.20217893e+00, 5.75820053e-01, 5.76895347e-01,
2.40507473e-02, 2.45346297e-01, 5.79583582e-01,
1.00432480e+00, 6.57900845e-01, 1.44824209e+00,
5.17574949e-01, 5.63454169e-01, 1.22260952e-01,
1.84054527e-02, 1.78713899e-01, 4.48935335e-01,
7.95896936e-01, 3.55205523e-01, 1.02314245e+00,
3.82625525e-04, 6.16322802e-01, 1.42745307e-02,
6.98941240e-01, 3.50545914e-01, 7.73853404e-01,
6.62022806e-01, 3.29577677e-01, 5.24743577e-01,
2.68823554e-02, 2.24019628e-01, 1.86742762e-01,
1.58570054e+00, 1.73050683e+00, 3.14165127e+00,
4.69903572e+00, 3.90331800e+00, 6.73313395e+00,
5.15962008e-03, 4.18827097e-01, 8.03244779e+00,
7.71523600e+00, 7.01987907e+00, 1.75165428e+00,
3.49291405e+00, 1.21221501e-01, 1.38587503e+00,
1.30486953e+00, 2.05560411e+00, 1.01669068e+00,
2.60041985e+00, 6.60589080e-01, 1.23855972e+00,
2.84773752e+00, 1.42100131e+00, 3.24021990e+01,
2.95526694e+00])]
bin values bin_values[-1]==0
mass ranges [22.5, 30]
Out[130]:
0
<matplotlib.figure.Figure at 0xb020c48c>
Both cases below s7,s8 should deliver same number of stars
In [121]:
reload(s)
s7=s.sygma(iolevel=0,mgal=1e11,dt=7e6,tend=1e9,imf_type='salpeter',imf_bdys=[1,30],special_timesteps=-1,iniZ=-1,hardsetZ=0.0001,table='yield_tables/isotope_yield_table_h1.txt',sn1a_on=True, sn1a_table='yield_tables/sn1a_h1.txt', iniabu_table='yield_tables/iniabu/iniab1.0E-04GN93_alpha_h1.ppn',pop3_table='yield_tables/popIII_h1.txt')
s8=s.sygma(iolevel=0,mgal=1e11,dt=7e6,tend=1e9,imf_type='salpeter',imf_bdys=[1,30],special_timesteps=-1,iniZ=0.0001)
################ Start SYGMA simulation ##############
Use specified tables as input
Note that table, sn1a_table iniabu_table and pop3_table have to be specified as an input
/home/christian/NuGrid/SYGMA_PROJECT/SYGMA/yield_tables/popIII_h1.txt
10 30
Use SNIa tables from yield_tables/sn1a_t86.txt
Use initial abundance of yield_tables/iniabu/iniab1.0E-04GN93_alpha_h1.ppn
START SFRMDOT Run time: 1.08s
Use NuGrid yields
Total mass of the gas in stars:
1.000E+11
sum msssfac 36877281297.2
END SFRMDOT Run time: 1.5s
################ Simulation successful ##############
Run time: 2.18s
################ Start SYGMA simulation ##############
10 30
Use SNIa tables from yield_tables/sn1a_t86.txt
START SFRMDOT Run time: 2.54s
Use NuGrid yields
Total mass of the gas in stars:
1.000E+11
sum msssfac 36877281297.2
END SFRMDOT Run time: 3.0s
################ Simulation successful ##############
Run time: 3.99s
In [122]:
s7.plot_sn_distr(rate=True,rate_only='sn2',marker2='o',markevery=3)
s8.plot_sn_distr(rate=True,rate_only='sn2',markevery=5,marker2='x')
#plt.xlim(1e6,1e8)
In [ ]:
s7.plot_sn_distr(rate=False)
s8.plot_sn_distr(rate=False)
s7.plot_sn_distr(rate=True)
s8.plot_sn_distr(rate=True)
In [1]:
import chem_evol
reload(chem_evol)
import sygma as s
reload(s)
imf_bdys = [0.1, 100.0]
imf_yields_range = [1.0, 100.0]
table = 'yield_tables/isotope_yield_table_MESA_only.txt'
s_imf_y_0_02 = s.sygma(table=table,imf_bdys=imf_bdys,imf_yields_range=imf_yields_range,iniZ=0.02)
imf_yields_range = [1.0, 30.0]
table = 'yield_tables/isotope_yield_table_heger_GG_ertl_fullZ_CL13_56Ni.txt'
s_imf_y_0_02 = s.sygma(table=table,imf_bdys=imf_bdys,imf_yields_range=imf_yields_range,iniZ=0.02)
imf_yields_range = [1.0, 100.0]
table = 'yield_tables/isotope_yield_table_heger_GG_ertl_fullZ_CL13_56Ni.txt'
s_imf_y_0_02 = s.sygma(table=table,imf_bdys=imf_bdys,imf_yields_range=imf_yields_range,iniZ=0.02)
SYGMA run in progress..
SYGMA run completed - Run time: 0.26s
SYGMA run in progress..
SYGMA run completed - Run time: 0.27s
SYGMA run in progress..
---------------------------------------------------------------------------
IndexError Traceback (most recent call last)
<ipython-input-1-f0e06dd50fc6> in <module>()
16 imf_yields_range = [1.0, 100.0]
17 table = 'yield_tables/isotope_yield_table_heger_GG_ertl_fullZ_CL13_56Ni.txt'
---> 18 s_imf_y_0_02 = s.sygma(table=table,imf_bdys=imf_bdys,imf_yields_range=imf_yields_range,iniZ=0.02)
/Users/benoitcote/Desktop/OMEGA_SYGMA_inheritance/sygma.pyc in __init__(self, sfr, imf_type, alphaimf, imf_bdys, sn1a_rate, iniZ, dt, special_timesteps, tend, mgal, transitionmass, iolevel, ini_alpha, table, hardsetZ, sn1a_on, sn1a_table, iniabu_table, extra_source_on, extra_source_table, pop3_table, imf_bdys_pop3, imf_yields_range_pop3, starbursts, beta_pow, gauss_dtd, exp_dtd, nb_1a_per_m, direct_norm_1a, Z_trans, f_arfo, imf_yields_range, exclude_masses, netyields_on, wiersmamod, yield_interp, dt_in, ytables_in, zm_lifetime_grid_nugrid_in, isotopes_in, ytables_pop3_in, zm_lifetime_grid_pop3_in, ytables_1a_in)
164
165 # Run the simulation
--> 166 self.__run_simulation()
167
168 # Do the final update of the history class
/Users/benoitcote/Desktop/OMEGA_SYGMA_inheritance/sygma.pyc in __run_simulation(self)
191
192 # Run the timestep i
--> 193 self._evol_stars(i)
194
195 # Get the new metallicity of the gas
/Users/benoitcote/Desktop/OMEGA_SYGMA_inheritance/chem_evol.pyc in _evol_stars(self, i)
855
856 # Calculate stellar ejecta and the number of SNe
--> 857 self.__sfrmdot(i)
858
859 # If no star is forming during the current timestep ...
/Users/benoitcote/Desktop/OMEGA_SYGMA_inheritance/chem_evol.pyc in __sfrmdot(self, i)
1055 # Apply the IMF on the mass boundaries
1056 mass_bdys, mstars, yields_all, massfac, mftot = \
-> 1057 self.__apply_imf_mass_bdys( mass_bdys, mstars, yields_all )
1058
1059 # Output information
/Users/benoitcote/Desktop/OMEGA_SYGMA_inheritance/chem_evol.pyc in __apply_imf_mass_bdys(self, mass_bdys, mstars, yields_all)
1355
1356 # If the mass bin is outside of the IMF mass range ...
-> 1357 if (mass_bdys[1+w] <= imf_bd[0]) or (mass_bdys[w] >= imf_bd[1]):
1358
1359 # Skip this bin
IndexError: list index out of range
In [3]:
%matplotlib inline
s1=s.sygma(iniZ=0.02,dt=1e7,tend=2e7)
s1.plot_yield_input(fig=3,xaxis='[C/H]',yaxis='[Fe/H]',iniZ=0.0001,masses=[1,3,12,25],marker='x',color='b',shape='--')
SYGMA run in progress..
SYGMA run completed - Run time: 1.57s
---------------------------------------------------------------------------
NameError Traceback (most recent call last)
<ipython-input-3-08daa00393c3> in <module>()
1 get_ipython().magic(u'matplotlib inline')
2 s1=s.sygma(iniZ=0.02,dt=1e7,tend=2e7)
----> 3 s1.plot_yield_input(fig=3,xaxis='[C/H]',yaxis='[Fe/H]',iniZ=0.0001,masses=[1,3,12,25],marker='x',color='b',shape='--')
/astro/adampaul/nupycee/sygma.pyc in plot_yield_input(self, fig, xaxis, yaxis, iniZ, netyields, masses, label, marker, color, shape, table, fsize, fontsize, rspace, bspace, labelsize, legend_fontsize, solar_ab, netyields_iniabu)
566 ini_species=ini_elems
567 ini_species_frac=ini_elems_frac
--> 568 ini_species_frac_sol=ini_elem_frac_sol
569 isoavail=y_delay.get(M=mini,Z=Z,quantity='Isotopes')
570 yx2=0
NameError: global name 'ini_elem_frac_sol' is not defined
<matplotlib.figure.Figure at 0x3fa1cd0>
In [ ]:
Content source: NuGrid/NuPyCEE
Similar notebooks: