diff --git a/test_problems/ChemEquil_ionizedGas/.cvsignore b/test_problems/ChemEquil_ionizedGas/.cvsignore index 59f637038..4b2391680 100644 --- a/test_problems/ChemEquil_ionizedGas/.cvsignore +++ b/test_problems/ChemEquil_ionizedGas/.cvsignore @@ -7,3 +7,6 @@ ct2ctml.log diff_test.out .depends *.d +diff_csv.txt +diff_out.txt +table.csv diff --git a/test_problems/ChemEquil_ionizedGas/ionizedGasEquil.cpp b/test_problems/ChemEquil_ionizedGas/ionizedGasEquil.cpp index d99b2536c..796d2947a 100644 --- a/test_problems/ChemEquil_ionizedGas/ionizedGasEquil.cpp +++ b/test_problems/ChemEquil_ionizedGas/ionizedGasEquil.cpp @@ -46,6 +46,18 @@ int main(int argc, char **argv) { IndVar1[6] = 1.0E-2; int nj = 6; int ni = 7; + FILE *FF = fopen("table.csv","w"); + int kk = gas->nSpecies(); + std::vector Xmol(kk, 0.0); + const std::vector &snames = gas->speciesNames(); + fprintf(FF,"Temperature, Pressure,"); + for (int k = 0; k < kk; k++) { + fprintf(FF, "%11s", snames[k].c_str()); + if (k < kk-1) { + fprintf(FF,","); + } + } + fprintf(FF,"\n"); for (int j=0; jtemperature(); pres = gas->pressure(); printf("Final T = %g, pres = %g atm\n", tkelvin, pres/OneAtm); - + gas->getMoleFractions(DATA_PTR(Xmol)); + fprintf(FF,"%10.4g, %10.4g,", tkelvin, pres); + for (int k = 0; k < kk; k++) { + if (fabs(Xmol[k]) < 1.0E-130) { + fprintf(FF," %10.4g", 0.0); + } else { + fprintf(FF," %10.4g", Xmol[k]); + } + if (k < kk-1) { + fprintf(FF,","); + } + } + fprintf(FF,"\n"); } } delete gas; + fclose(FF); } catch (CanteraError) { diff --git a/test_problems/ChemEquil_ionizedGas/output_blessed.txt b/test_problems/ChemEquil_ionizedGas/output_blessed.txt index d0daf3cbd..1fe49c8b6 100644 --- a/test_problems/ChemEquil_ionizedGas/output_blessed.txt +++ b/test_problems/ChemEquil_ionizedGas/output_blessed.txt @@ -89,7 +89,7 @@ Initial T = 207.284, pres = 5.8944e-07 atm N 9.55809e-114 4.63943e-114 -18.9446 O 1.18094e-57 6.54771e-58 -20.4172 NO 1.82556e-23 1.89829e-23 -39.3618 - NO+ 2.2952e-228 2.38659e-228 12.2863 + NO+ 2.2952e-228 2.3866e-228 12.2863 Electron 2.2952e-228 4.33484e-233 -541.209 N+ 0 0 O+ 0 0 @@ -122,8 +122,8 @@ Initial T = 207.284, pres = 5.8944e-06 atm N 3.02253e-114 1.46712e-114 -17.7933 O 3.73447e-58 2.07057e-58 -19.266 NO 1.82556e-23 1.89829e-23 -37.0592 - NO+ 1.2588e-229 1.30893e-229 11.6857 - Electron 1.2588e-229 2.37745e-234 -541.81 + NO+ 1.25881e-229 1.30893e-229 11.6857 + Electron 1.25881e-229 2.37745e-234 -541.81 N+ 0 0 O+ 0 0 N2+ 0 0 @@ -188,8 +188,8 @@ Initial T = 207.284, pres = 0.00058944 atm N 3.02253e-115 1.46712e-115 -15.4907 O 3.73447e-59 2.07057e-59 -16.9634 NO 1.82556e-23 1.89829e-23 -32.454 - NO+ 1.26648e-231 1.31691e-231 11.6917 - Electron 1.26648e-231 2.39194e-236 -541.804 + NO+ 1.2656e-231 1.316e-231 11.691 + Electron 1.2656e-231 2.39028e-236 -541.804 N+ 0 0 O+ 0 0 N2+ 0 0 @@ -221,8 +221,8 @@ Initial T = 207.284, pres = 0.0058944 atm N 9.55809e-116 4.63943e-116 -14.3394 O 1.18094e-59 6.54771e-60 -15.8121 NO 1.82556e-23 1.89829e-23 -30.1515 - NO+ 2.40463e-232 2.50039e-232 12.3329 - Electron 2.40463e-232 4.54152e-237 -541.163 + NO+ 2.40459e-232 2.50035e-232 12.3329 + Electron 2.40459e-232 4.54145e-237 -541.163 N+ 0 0 O+ 0 0 N2+ 0 0 @@ -254,12 +254,12 @@ Initial T = 414.143, pres = 1.17767e-08 atm N 2.06855e-53 1.00406e-53 -20.8551 O 1.8737e-25 1.03887e-25 -22.3285 NO 5.58575e-12 5.80828e-12 -43.1837 - NO+ 4.00599e-102 4.16552e-102 11.996 - Electron 4.00599e-102 7.56595e-107 -254.374 - N+ 7.16723e-207 3.47879e-207 34.3246 - O+ 3.0296e-168 1.6797e-168 32.8512 - N2+ 6.19807e-167 6.01689e-167 13.4694 - O2+ 3.22979e-125 3.58144e-125 10.5226 + NO+ 4.00605e-102 4.16558e-102 11.996 + Electron 4.00605e-102 7.56606e-107 -254.374 + N+ 7.16734e-207 3.47884e-207 34.3246 + O+ 3.02965e-168 1.67972e-168 32.8512 + N2+ 6.19817e-167 6.01698e-167 13.4694 + O2+ 3.22984e-125 3.58149e-125 10.5226 Final T = 414.143, pres = 1.17767e-08 atm Initial T = 414.143, pres = 1.17767e-07 atm @@ -287,12 +287,12 @@ Initial T = 414.143, pres = 1.17767e-07 atm N 6.54134e-54 3.17512e-54 -19.7038 O 5.92517e-26 3.28519e-26 -21.1772 NO 5.58575e-12 5.80828e-12 -40.8811 - NO+ 5.13663e-103 5.34118e-103 12.2446 - Electron 5.13663e-103 9.70133e-108 -254.125 - N+ 2.90616e-208 1.41057e-208 33.4219 - O+ 1.22844e-169 6.81081e-170 31.9485 - N2+ 7.9474e-168 7.71507e-168 13.718 - O2+ 4.14136e-126 4.59225e-126 10.7712 + NO+ 5.13662e-103 5.34116e-103 12.2446 + Electron 5.13662e-103 9.7013e-108 -254.125 + N+ 2.90615e-208 1.41057e-208 33.4219 + O+ 1.22844e-169 6.81079e-170 31.9485 + N2+ 7.94737e-168 7.71505e-168 13.718 + O2+ 4.14134e-126 4.59223e-126 10.7712 Final T = 414.143, pres = 1.17767e-07 atm Initial T = 414.143, pres = 1.17767e-06 atm @@ -320,12 +320,12 @@ Initial T = 414.143, pres = 1.17767e-06 atm N 2.06855e-54 1.00406e-54 -18.5525 O 1.8737e-26 1.03887e-26 -20.026 NO 5.58575e-12 5.80828e-12 -38.5785 - NO+ 2.82554e-104 2.93805e-104 11.6469 - Electron 2.82554e-104 5.33647e-109 -254.723 - N+ 5.05525e-210 2.45369e-210 31.6729 - O+ 2.13686e-171 1.18474e-171 30.1995 - N2+ 4.37167e-169 4.24387e-169 13.1203 - O2+ 2.27806e-127 2.52609e-127 10.1735 + NO+ 2.82555e-104 2.93806e-104 11.6469 + Electron 2.82555e-104 5.33649e-109 -254.723 + N+ 5.05526e-210 2.45369e-210 31.6729 + O+ 2.13687e-171 1.18474e-171 30.1995 + N2+ 4.37169e-169 4.24389e-169 13.1203 + O2+ 2.27807e-127 2.52609e-127 10.1735 Final T = 414.143, pres = 1.17767e-06 atm Initial T = 414.143, pres = 1.17767e-05 atm @@ -353,12 +353,12 @@ Initial T = 414.143, pres = 1.17767e-05 atm N 6.54134e-55 3.17512e-55 -17.4013 O 5.92517e-27 3.28519e-27 -18.8747 NO 5.58575e-12 5.80828e-12 -36.2759 - NO+ 3.94044e-105 4.09735e-105 11.9795 - Electron 3.94044e-105 7.44214e-110 -254.39 + NO+ 3.94045e-105 4.09736e-105 11.9795 + Electron 3.94045e-105 7.44215e-110 -254.39 N+ 2.22939e-211 1.08209e-211 30.8542 - O+ 9.42367e-173 5.22475e-173 29.3808 - N2+ 6.09665e-170 5.91843e-170 13.4529 - O2+ 3.17694e-128 3.52283e-128 10.5061 + O+ 9.42368e-173 5.22475e-173 29.3808 + N2+ 6.09666e-170 5.91843e-170 13.4529 + O2+ 3.17694e-128 3.52284e-128 10.5061 Final T = 414.143, pres = 1.17767e-05 atm Initial T = 414.143, pres = 0.000117767 atm @@ -390,7 +390,7 @@ Initial T = 414.143, pres = 0.000117767 atm Electron 2.91676e-106 5.50875e-111 -254.691 N+ 5.21845e-213 2.5329e-213 29.4021 O+ 2.20585e-174 1.22298e-174 27.9287 - N2+ 4.51281e-171 4.38089e-171 13.1521 + N2+ 4.51281e-171 4.38088e-171 13.1521 O2+ 2.35161e-129 2.60764e-129 10.2053 Final T = 414.143, pres = 0.000117767 atm @@ -419,12 +419,12 @@ Initial T = 414.143, pres = 0.00117767 atm N 6.54134e-56 3.17512e-56 -15.0987 O 5.92517e-28 3.28519e-28 -16.5721 NO 5.58575e-12 5.80828e-12 -31.6707 - NO+ 3.85738e-107 4.01099e-107 11.9582 - Electron 3.85738e-107 7.28527e-112 -254.411 - N+ 2.1824e-214 1.05928e-214 28.5303 - O+ 9.22503e-176 5.11462e-176 27.0569 - N2+ 5.96814e-172 5.79367e-172 13.4316 - O2+ 3.10997e-130 3.44857e-130 10.4848 + NO+ 3.85557e-107 4.00911e-107 11.9577 + Electron 3.85557e-107 7.28185e-112 -254.412 + N+ 2.18137e-214 1.05878e-214 28.5298 + O+ 9.22071e-176 5.11222e-176 27.0564 + N2+ 5.96534e-172 5.79096e-172 13.4311 + O2+ 3.10852e-130 3.44696e-130 10.4843 Final T = 414.143, pres = 0.00117767 atm Initial T = 414.143, pres = 0.0117767 atm @@ -452,12 +452,12 @@ Initial T = 414.143, pres = 0.0117767 atm N 2.06855e-56 1.00406e-56 -13.9474 O 1.8737e-28 1.03887e-28 -15.4208 NO 5.58575e-12 5.80828e-12 -29.3682 - NO+ 6.14689e-108 6.39166e-108 12.4242 - Electron 6.14689e-108 1.16094e-112 -253.946 - N+ 1.09976e-215 5.33793e-216 27.845 - O+ 4.64869e-177 2.57736e-177 26.3716 - N2+ 9.51046e-173 9.23244e-173 13.8976 - O2+ 4.95586e-131 5.49544e-131 10.9508 + NO+ 2.85055e-102 2.96406e-102 25.4712 + Electron 2.85055e-102 5.38371e-107 -240.898 + N+ 5.1e-210 2.47541e-210 40.892 + O+ 2.15578e-171 1.19522e-171 39.4186 + N2+ 4.41037e-167 4.28145e-167 26.9447 + O2+ 2.29823e-125 2.54845e-125 23.9978 Final T = 414.143, pres = 0.0117767 atm Initial T = 1153.16, pres = 3.27918e-08 atm @@ -485,12 +485,12 @@ Initial T = 1153.16, pres = 3.27918e-08 atm N 2.86312e-15 1.38975e-15 -21.3563 O 2.23425e-05 1.23879e-05 -22.8775 NO 0.000134494 0.000139854 -44.2338 - NO+ 9.03856e-22 9.39859e-22 11.6755 - Electron 9.03856e-22 1.70709e-26 -69.7353 - N+ 1.59896e-54 7.76101e-55 34.553 - O+ 3.05445e-41 1.69349e-41 33.0319 - N2+ 1.21258e-44 1.17715e-44 13.1967 - O2+ 3.22543e-30 3.57664e-30 10.1544 + NO+ 7.02819e-19 7.30815e-19 18.3317 + Electron 7.02819e-19 1.3274e-23 -63.0791 + N+ 1.24332e-51 6.0348e-52 41.2092 + O+ 2.37508e-38 1.31683e-38 39.688 + N2+ 9.42879e-42 9.15326e-42 19.8529 + O2+ 2.50803e-27 2.78112e-27 16.8106 Final T = 1152.46, pres = 3.27723e-08 atm Initial T = 1152.46, pres = 3.27723e-07 atm @@ -518,12 +518,12 @@ Initial T = 1152.46, pres = 3.27723e-07 atm N 9.11146e-16 4.42265e-16 -20.2051 O 7.08903e-06 3.93051e-06 -21.7263 NO 0.000134663 0.000140028 -41.9314 - NO+ 9.05883e-23 9.4196e-23 11.6643 - Electron 9.05883e-23 1.71091e-27 -69.7332 - N+ 5.12845e-56 2.48922e-56 33.3906 - O+ 9.75583e-43 5.40893e-43 31.8694 - N2+ 1.22376e-45 1.18799e-45 13.1855 - O2+ 3.24058e-31 3.59342e-31 10.1431 + NO+ 9.67397e-23 1.00592e-22 11.73 + Electron 9.67398e-23 1.82709e-27 -69.6675 + N+ 5.4767e-56 2.65826e-56 33.4563 + O+ 1.04183e-42 5.77622e-43 31.9351 + N2+ 1.30686e-45 1.26866e-45 13.2512 + O2+ 3.46064e-31 3.83743e-31 10.2088 Final T = 1152.61, pres = 3.27763e-07 atm Initial T = 1152.61, pres = 3.27763e-06 atm @@ -551,12 +551,12 @@ Initial T = 1152.61, pres = 3.27763e-06 atm N 2.88709e-16 1.40138e-16 -19.0539 O 2.24413e-06 1.24425e-06 -20.575 NO 0.000134716 0.000140083 -39.6289 - NO+ 9.06539e-24 9.4264e-24 11.6607 - Electron 9.06539e-24 1.71214e-28 -69.7325 - N+ 1.62909e-57 7.90717e-58 32.2358 - O+ 3.09489e-44 1.7159e-44 30.7146 - N2+ 1.22735e-46 1.19147e-46 13.1819 - O2+ 3.24545e-32 3.5988e-32 10.1396 + NO+ 9.16628e-20 9.5313e-20 20.8821 + Electron 9.16628e-20 1.7312e-24 -60.5111 + N+ 1.64721e-53 7.99516e-54 41.4572 + O+ 3.12933e-40 1.73499e-40 39.936 + N2+ 1.24101e-42 1.20473e-42 22.4033 + O2+ 3.28157e-28 3.63885e-28 19.361 Final T = 1152.66, pres = 3.27775e-06 atm Initial T = 1152.66, pres = 3.27775e-05 atm @@ -584,12 +584,12 @@ Initial T = 1152.66, pres = 3.27775e-05 atm N 9.1356e-17 4.43436e-17 -17.9026 O 7.09896e-07 3.936e-07 -19.4238 NO 0.000134733 0.000140101 -37.3264 - NO+ 9.06748e-25 9.42856e-25 11.6596 - Electron 9.06748e-25 1.71254e-29 -69.7322 - N+ 5.15899e-59 2.50404e-59 31.0834 - O+ 9.79676e-46 5.4316e-46 29.5622 - N2+ 1.22849e-47 1.19258e-47 13.1808 - O2+ 3.247e-33 3.60052e-33 10.1384 + NO+ 1.27093e-23 1.32154e-23 14.2998 + Electron 1.27093e-23 2.40034e-28 -67.092 + N+ 7.231e-58 3.50974e-58 33.7236 + O+ 1.37315e-44 7.61311e-45 32.2024 + N2+ 1.72189e-46 1.67155e-46 15.821 + O2+ 4.5511e-32 5.0466e-32 12.7787 Final T = 1152.67, pres = 3.27779e-05 atm Initial T = 1152.67, pres = 0.000327779 atm @@ -598,7 +598,7 @@ Initial T = 1152.67, pres = 0.000327779 atm temperature 1152.68 K pressure 33.2124 Pa - density 1e-04 kg/m^3 + density 0.0001 kg/m^3 mean mol. weight 28.8565 amu 1 kg 1 kmol @@ -617,12 +617,12 @@ Initial T = 1152.67, pres = 0.000327779 atm N 2.88951e-17 1.40255e-17 -16.7513 O 2.24513e-07 1.24481e-07 -18.2725 NO 0.000134738 0.000140106 -35.0238 - NO+ 9.06814e-26 9.42925e-26 11.6592 - Electron 9.06814e-26 1.71266e-30 -69.7322 - N+ 1.63215e-60 7.92205e-61 29.9317 - O+ 3.099e-47 1.71817e-47 28.4106 - N2+ 1.22885e-48 1.19293e-48 13.1804 - O2+ 3.24749e-34 3.60106e-34 10.1381 + NO+ 5.59478e-26 5.81757e-26 11.1763 + Electron 5.59478e-26 1.05666e-30 -70.2151 + N+ 1.00699e-60 4.88768e-61 29.4488 + O+ 1.91199e-47 1.06006e-47 27.9276 + N2+ 7.58165e-49 7.36002e-49 12.6975 + O2+ 2.00361e-34 2.22175e-34 9.65516 Final T = 1152.68, pres = 0.00032778 atm Initial T = 1152.68, pres = 0.0032778 atm @@ -650,12 +650,12 @@ Initial T = 1152.68, pres = 0.0032778 atm N 9.13803e-18 4.43553e-18 -15.6 O 7.09996e-08 3.93655e-08 -17.1212 NO 0.00013474 0.000140108 -32.7212 - NO+ 9.06835e-27 9.42946e-27 11.6591 - Electron 9.06835e-27 1.7127e-31 -69.7322 - N+ 5.16206e-62 2.50553e-62 28.7803 - O+ 9.80088e-49 5.43388e-49 27.2592 - N2+ 1.22896e-49 1.19304e-49 13.1803 - O2+ 3.24764e-35 3.60123e-35 10.138 + NO+ 2.9014e-21 3.01693e-21 24.335 + Electron 2.9014e-21 5.47974e-26 -57.0562 + N+ 1.65159e-56 8.01638e-57 41.4562 + O+ 3.13577e-43 1.73856e-43 39.9351 + N2+ 3.93204e-44 3.81709e-44 25.8562 + O2+ 1.03908e-29 1.1522e-29 22.8139 Final T = 1152.68, pres = 0.00327781 atm Initial T = 1152.68, pres = 0.0327781 atm @@ -683,12 +683,12 @@ Initial T = 1152.68, pres = 0.0327781 atm N 2.88976e-18 1.40267e-18 -14.4487 O 2.24523e-08 1.24486e-08 -15.9699 NO 0.00013474 0.000140109 -30.4186 - NO+ 9.8979e-28 1.0292e-27 11.7466 - Electron 9.8979e-28 1.86937e-32 -69.6446 - N+ 1.78178e-63 8.6483e-64 27.7165 - O+ 3.38291e-50 1.87558e-50 26.1954 - N2+ 1.34142e-50 1.3022e-50 13.2678 - O2+ 3.54476e-36 3.93069e-36 10.2255 + NO+ 9.1752e-22 9.54057e-22 25.4863 + Electron 9.1752e-22 1.73288e-26 -55.9049 + N+ 1.65168e-57 8.01684e-58 41.4562 + O+ 3.13591e-44 1.73863e-44 39.935 + N2+ 1.24347e-44 1.20712e-44 27.0075 + O2+ 3.28593e-30 3.64369e-30 23.9652 Final T = 1152.68, pres = 0.0327781 atm Initial T = 3008.94, pres = 8.55633e-08 atm @@ -1060,7 +1060,7 @@ Initial T = 3318.59, pres = 0.00114605 atm temperature 3417.09 K pressure 119.127 Pa - density 1e-04 kg/m^3 + density 0.0001 kg/m^3 mean mol. weight 23.8496 amu 1 kg 1 kmol diff --git a/test_problems/ChemEquil_ionizedGas/runtest b/test_problems/ChemEquil_ionizedGas/runtest index be204380e..fbc3d54b4 100755 --- a/test_problems/ChemEquil_ionizedGas/runtest +++ b/test_problems/ChemEquil_ionizedGas/runtest @@ -2,7 +2,7 @@ # # temp_success="1" -/bin/rm -f output.txt outputa.txt +/bin/rm -f output.txt outputa.txt table.csv diff_out.txt diff_csv.txt testName="ChemEquil_ionizedGas" ################################################################# @@ -19,14 +19,34 @@ then echo "ERROR: $testName: ionizedGasEquil returned with bad status, $retnStat, check output" fi -diff -w output.txt output_blessed.txt > diff_test.out -retnStat=$? -if [ $retnStat = "0" ] -then - echo "successful diff comparison on $testName test" - temp_success="0" -else - echo "FAILED: unsuccessful diff comparison on $testName test" - echo "FAILED" > csvCode.txt -fi +diff -w output.txt output_blessed.txt > diff_out.txt +retnStat_txt=$? +csvdiff -a 1.0E-50 table.csv table_blessed.csv > diff_csv.txt +retnStat_csv=$? + +eCode=1 +if test $retnStat_csv = "1" +then + eCode=0 + echo "Successful test comparison on " $testName + if [ $retnStat_txt != "0" ] + then + echo " But, text output files have differences. See diff_out.txt" + fi +else + echo "Unsuccessful test comparison on " $testName " test" + if test $retnStat_csv != "1" + then + echo " csv files are different - see diff_csv.txt" + fi + if test $retnStat_txt != "0" + then + echo " And, text output files have differences. See diff_out.txt" + fi +fi +# +# +exit 0 + + diff --git a/test_problems/ChemEquil_ionizedGas/table_blessed.csv b/test_problems/ChemEquil_ionizedGas/table_blessed.csv new file mode 100644 index 000000000..158cc1895 --- /dev/null +++ b/test_problems/ChemEquil_ionizedGas/table_blessed.csv @@ -0,0 +1,43 @@ +Temperature, Pressure, N2, O2, N, O, NO, NO+, Electron, N+, O+, N2+, O2+ + 207.3, 0.0005973, 0.7885, 0.2115, 9.558e-113, 1.181e-56, 1.826e-23, 0, 0, 0, 0, 0, 0 + 207.3, 0.005973, 0.7885, 0.2115, 3.023e-113, 3.734e-57, 1.826e-23, 0, 0, 0, 0, 0, 0 + 207.3, 0.05973, 0.7885, 0.2115, 9.558e-114, 1.181e-57, 1.826e-23, 0, 0, 0, 0, 0, 0 + 207.3, 0.5973, 0.7885, 0.2115, 3.023e-114, 3.734e-58, 1.826e-23, 0, 0, 0, 0, 0, 0 + 207.3, 5.973, 0.7885, 0.2115, 9.558e-115, 1.181e-58, 1.826e-23, 0, 0, 0, 0, 0, 0 + 207.3, 59.73, 0.7885, 0.2115, 3.023e-115, 3.734e-59, 1.826e-23, 0, 0, 0, 0, 0, 0 + 207.3, 597.3, 0.7885, 0.2115, 9.558e-116, 1.181e-59, 1.826e-23, 0, 0, 0, 0, 0, 0 + 414.1, 0.001193, 0.7885, 0.2115, 2.069e-53, 1.874e-25, 5.586e-12, 4.006e-102, 4.006e-102, 0, 0, 0, 3.23e-125 + 414.1, 0.01193, 0.7885, 0.2115, 6.541e-54, 5.925e-26, 5.586e-12, 5.137e-103, 5.137e-103, 0, 0, 0, 4.141e-126 + 414.1, 0.1193, 0.7885, 0.2115, 2.069e-54, 1.874e-26, 5.586e-12, 2.826e-104, 2.826e-104, 0, 0, 0, 2.278e-127 + 414.1, 1.193, 0.7885, 0.2115, 6.541e-55, 5.925e-27, 5.586e-12, 3.94e-105, 3.94e-105, 0, 0, 0, 3.177e-128 + 414.1, 11.93, 0.7885, 0.2115, 2.069e-55, 1.874e-27, 5.586e-12, 2.917e-106, 2.917e-106, 0, 0, 0, 2.352e-129 + 414.1, 119.3, 0.7885, 0.2115, 6.541e-56, 5.925e-28, 5.586e-12, 3.856e-107, 3.856e-107, 0, 0, 0, 3.109e-130 + 414.1, 1193, 0.7885, 0.2115, 2.069e-56, 1.874e-28, 5.586e-12, 2.851e-102, 2.851e-102, 0, 0, 0, 2.298e-125 + 1152, 0.003321, 0.7884, 0.2114, 2.863e-15, 2.234e-05, 0.0001345, 7.028e-19, 7.028e-19, 1.243e-51, 2.375e-38, 9.429e-42, 2.508e-27 + 1153, 0.03321, 0.7884, 0.2115, 9.111e-16, 7.089e-06, 0.0001347, 9.674e-23, 9.674e-23, 5.477e-56, 1.042e-42, 1.307e-45, 3.461e-31 + 1153, 0.3321, 0.7884, 0.2115, 2.887e-16, 2.244e-06, 0.0001347, 9.166e-20, 9.166e-20, 1.647e-53, 3.129e-40, 1.241e-42, 3.282e-28 + 1153, 3.321, 0.7884, 0.2115, 9.136e-17, 7.099e-07, 0.0001347, 1.271e-23, 1.271e-23, 7.231e-58, 1.373e-44, 1.722e-46, 4.551e-32 + 1153, 33.21, 0.7884, 0.2115, 2.89e-17, 2.245e-07, 0.0001347, 5.595e-26, 5.595e-26, 1.007e-60, 1.912e-47, 7.582e-49, 2.004e-34 + 1153, 332.1, 0.7884, 0.2115, 9.138e-18, 7.1e-08, 0.0001347, 2.901e-21, 2.901e-21, 1.652e-56, 3.136e-43, 3.932e-44, 1.039e-29 + 1153, 3321, 0.7884, 0.2115, 2.89e-18, 2.245e-08, 0.0001347, 9.175e-22, 9.175e-22, 1.652e-57, 3.136e-44, 1.243e-44, 3.286e-30 + 1764, 0.005446, 0.7345, 0.1293, 7.047e-08, 0.1335, 0.00278, 5.61e-11, 5.61e-11, 1.279e-29, 2.263e-21, 1.383e-25, 1.154e-16 + 1879, 0.05766, 0.7384, 0.1357, 1.608e-07, 0.1218, 0.004184, 1.472e-10, 1.472e-10, 4.308e-28, 2.092e-20, 3.125e-24, 6.605e-16 + 2007, 0.6114, 0.7426, 0.1427, 3.489e-07, 0.1084, 0.006245, 3.68e-10, 3.68e-10, 1.256e-26, 1.732e-19, 6.362e-23, 3.535e-15 + 2147, 6.489, 0.7472, 0.1506, 7.059e-07, 0.09302, 0.009208, 8.567e-10, 8.567e-10, 2.983e-25, 1.225e-18, 1.116e-21, 1.716e-14 + 2299, 68.85, 0.752, 0.1591, 1.288e-06, 0.07561, 0.01332, 1.793e-09, 1.793e-09, 5.239e-24, 6.863e-18, 1.573e-20, 7.206e-14 + 2457, 728.7, 0.7568, 0.1679, 2.009e-06, 0.05658, 0.01867, 3.186e-09, 3.186e-09, 5.818e-23, 2.697e-17, 1.587e-19, 2.417e-13 + 2611, 7665, 0.7613, 0.1765, 2.472e-06, 0.0373, 0.02495, 4.397e-09, 4.398e-09, 3.215e-22, 6.159e-17, 9.61e-19, 5.738e-13 + 2946, 0.01041, 0.6308, 1.393e-06, 0.02415, 0.3449, 0.0001029, 2.912e-06, 2.913e-06, 7.396e-12, 8.35e-10, 3.216e-12, 3.032e-12 + 3080, 0.1085, 0.6361, 5.853e-06, 0.01774, 0.3459, 0.0002489, 3.28e-06, 3.281e-06, 6.298e-12, 8.265e-10, 4.524e-12, 9.624e-12 + 3207, 1.126, 0.6411, 2.749e-05, 0.0117, 0.3466, 0.0006234, 3.384e-06, 3.384e-06, 3.787e-12, 6.566e-10, 4.893e-12, 2.838e-11 + 3319, 11.61, 0.6451, 0.0001486, 0.006702, 0.3464, 0.001629, 3.11e-06, 3.11e-06, 1.453e-12, 3.921e-10, 3.762e-12, 7.586e-11 + 3417, 119.1, 0.6477, 0.0008763, 0.003474, 0.3435, 0.004358, 2.623e-06, 2.623e-06, 4.057e-13, 1.906e-10, 2.273e-12, 1.873e-10 + 3555, 1230, 0.6499, 0.00421, 0.002098, 0.3329, 0.01082, 2.485e-06, 2.485e-06, 1.882e-13, 1.254e-10, 2.028e-12, 4.968e-10 + 3800, 1.297e+04, 0.6534, 0.01249, 0.001863, 0.3095, 0.02275, 3.185e-06, 3.186e-06, 3.119e-13, 1.785e-10, 4.815e-12, 1.621e-09 + 3294, 0.0125, 0.5182, 1.58e-07, 0.1606, 0.3211, 4.646e-05, 1.403e-05, 1.408e-05, 4.794e-09, 5.129e-08, 4.044e-10, 1.185e-11 + 3516, 0.1328, 0.5258, 5.179e-07, 0.1514, 0.3227, 0.0001044, 1.954e-05, 1.963e-05, 9.133e-09, 8.456e-08, 1.06e-09, 4.501e-11 + 3767, 1.414, 0.5345, 1.729e-06, 0.1408, 0.3244, 0.0002361, 2.703e-05, 2.719e-05, 1.669e-08, 1.359e-07, 2.705e-09, 1.701e-10 + 4050, 15.11, 0.5445, 5.91e-06, 0.1286, 0.3262, 0.0005385, 3.701e-05, 3.726e-05, 2.88e-08, 2.11e-07, 6.656e-09, 6.382e-10 + 4372, 161.8, 0.556, 2.081e-05, 0.1146, 0.3281, 0.001241, 4.987e-05, 5.025e-05, 4.594e-08, 3.117e-07, 1.555e-08, 2.363e-09 + 4734, 1736, 0.5688, 7.603e-05, 0.09874, 0.3294, 0.002891, 6.556e-05, 6.609e-05, 6.571e-08, 4.295e-07, 3.37e-08, 8.558e-09 + 5140, 1.863e+04, 0.5827, 0.0002894, 0.08126, 0.3288, 0.006811, 8.306e-05, 8.378e-05, 8.126e-08, 5.379e-07, 6.601e-08, 2.99e-08