File name: Pandora123s1_StonyPlain_20211103_L2Fit_fout2c1d20210901p1-8.txt File generation date: 20250129T030912.8Z Data description: Level 2 spectral fitting results file Data file version: fout2c1p1-8 Local principal investigator: Vitali Fioletov Network principal investigator: Alexander Cede Instrument type: Pandora Instrument number: 123 Spectrometer number: 1 Processing software version used: BlickP v1.8.70 Instrument operation file used: Pandora123_OF_v6d20161109.txt Instrument calibration file used: Pandora123s1_CF_v1d20210901.txt Level 1 file used: Pandora123s1_StonyPlain_20211103_L1_smca1c1d20210901p1-8.txt Full location name: Stony Plain Upper Air Station (CWSE) Short location name: StonyPlain Country of location: Canada Location latitude [deg]: 53.5500 Location longitude [deg]: -114.1100 Location altitude [m]: 766 Local noon date: 20211103 Notes on s-number (L1 configuration): Corrections NOT applied although requested by the s-number are latency correction, radiometric temperature correction, matrix method stray light correction (replaced by simple method) Data caveats: FieldCal for VP from 20210920 based on Session 3 First and last pixel inside fitting window: 218 360 Nominal wavelengths inside fitting window [nm]: 305.0517 305.1930 305.3343 305.4756 305.6168 305.7581 305.8993 306.0405 306.1817 306.3229 306.4641 306.6052 306.7464 306.8875 307.0287 307.1698 307.3109 307.4520 307.5931 307.7341 307.8752 308.0162 308.1573 308.2983 308.4393 308.5803 308.7213 308.8622 309.0032 309.1441 309.2851 309.4260 309.5669 309.7078 309.8487 309.9895 310.1304 310.2713 310.4121 310.5529 310.6937 310.8345 310.9753 311.1161 311.2568 311.3976 311.5383 311.6791 311.8198 311.9605 312.1012 312.2418 312.3825 312.5232 312.6638 312.8044 312.9450 313.0856 313.2262 313.3668 313.5074 313.6479 313.7885 313.9290 314.0695 314.2100 314.3505 314.4910 314.6315 314.7719 314.9124 315.0528 315.1932 315.3336 315.4740 315.6144 315.7548 315.8951 316.0355 316.1758 316.3162 316.4565 316.5968 316.7371 316.8773 317.0176 317.1578 317.2981 317.4383 317.5785 317.7187 317.8589 317.9991 318.1393 318.2794 318.4196 318.5597 318.6998 318.8399 318.9800 319.1201 319.2601 319.4002 319.5402 319.6803 319.8203 319.9603 320.1003 320.2403 320.3803 320.5202 320.6602 320.8001 320.9400 321.0799 321.2198 321.3597 321.4996 321.6395 321.7793 321.9192 322.0590 322.1988 322.3386 322.4784 322.6182 322.7579 322.8977 323.0374 323.1772 323.3169 323.4566 323.5963 323.7360 323.8756 324.0153 324.1549 324.2946 324.4342 324.5738 324.7134 324.8530 324.9926 --------------------------------------------------------------------------------------- Column 1: Two letter code of measurement routine Column 2: UT date and time for center-time of measurement, yyyymmddThhmmssZ (ISO 8601) Column 3: Fractional days since 1-Jan-2000 midnight for center-time of measurement Column 4: Routine count (1 for the first routine of the day, 2 for the second, etc.) Column 5: Repetition count (1 for the first set in the routine, 2 for the second, etc.) Column 6: Total duration of measurement set in seconds Column 7: Data processing type index Column 8: Solar zenith angle for center-time of measurement in degree Column 9: Solar azimuth for center-time of measurement in degree, 0=north, increases clockwise Column 10: Lunar zenith angle for center-time of measurement in degree Column 11: Lunar azimuth for center-time of measurement in degree, 0=north, increases clockwise Column 12: Pointing zenith angle in degree, absolute or relative (see next column), 999=tracker not used Column 13: Zenith pointing mode: zenith angle is... 0=absolute, 1=relative to sun, 2=relative to moon Column 14: Pointing azimuth in degree, increases clockwise, absolute (0=north) or relative (see next column), 999=tracker not used Column 15: Azimuth pointing mode: like zenith angle mode but also fixed scattering angles relative to sun (3) or moon (4) Column 16: Fitting result index: 0=no error or warning, 1, 2=warning, >2=error Column 17: Number of function evaluations used, 0=linear fitting or fitting not successful or no fitting done Column 18: rms of unweighted fitting residuals, -9=fitting not successful Column 19: Normalized rms of fitting residuals weighted with independent uncertainty, -9=fitting not successful or no uncertainty used Column 20: Expected rms based on independent uncertainty, -9=fitting not successful or no uncertainty given Column 21: Expected normalized weighted rms based on independent uncertainty, -9=fitting not successful or no uncertainty given Column 22: Ozone slant column amount [moles per square meter], -9e99=fitting not successful Column 23: Independent uncertainty of ozone slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -3=spectral fitting was done, but no independent uncertainty could be retrieved, -5=no independent uncertainty input was given, -9=spectral fitting not successful Column 24: Structured uncertainty of ozone slant column amount [moles per square meter], -7=not given since method "MEAS" was chosen, -9=spectral fitting not successful Column 25: Common uncertainty of ozone slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -6=no common uncertainty input was given, -7=not given since method "MEAS" was chosen, -9=spectral fitting not successful Column 26: rms-based uncertainty of ozone slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -3=spectral fitting was done, but no rms-based uncertainty could be retrieved, -9=spectral fitting not successful Column 27: Ozone effective temperature [K] Column 28: Independent uncertainty of ozone effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -2=no temperature fitting was requested and output for effective temperature and structured uncertainty of it is based on f-code, -3=spectral fitting was done, but no independent uncertainty could be retrieved, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -5=no independent uncertainty input was given, -9=spectral fitting not successful Column 29: Structured uncertainty of ozone effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -9=spectral fitting not successful Column 30: Common uncertainty of ozone effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -2=no temperature fitting was requested and output for effective temperature and structured uncertainty of it is based on f-code, -3=spectral fitting was done, but no common uncertainty could be retrieved, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -6=no common uncertainty input was given, -9=spectral fitting not successful Column 31: Effective ozone fitting wavelength [nm], -9=fitting not successful Column 32: Diffuse correction applied before fitting at effective fitting wavelength for ozone [%], 0=no diffuse correction applied or fitting not requested, >0=measured diffuse correction, <0=(negative value of) calculated diffuse correction Column 33: Sulfur dioxide slant column amount [moles per square meter], -9e99=fitting not successful Column 34: Independent uncertainty of sulfur dioxide slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -3=spectral fitting was done, but no independent uncertainty could be retrieved, -5=no independent uncertainty input was given, -9=spectral fitting not successful Column 35: Structured uncertainty of sulfur dioxide slant column amount [moles per square meter], -7=not given since method "MEAS" was chosen, -9=spectral fitting not successful Column 36: Common uncertainty of sulfur dioxide slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -6=no common uncertainty input was given, -7=not given since method "MEAS" was chosen, -9=spectral fitting not successful Column 37: rms-based uncertainty of sulfur dioxide slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -3=spectral fitting was done, but no rms-based uncertainty could be retrieved, -9=spectral fitting not successful Column 38: Sulfur dioxide effective temperature [K] Column 39: Independent uncertainty of sulfur dioxide effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -2=no temperature fitting was requested and output for effective temperature and structured uncertainty of it is based on f-code, -3=spectral fitting was done, but no independent uncertainty could be retrieved, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -5=no independent uncertainty input was given, -9=spectral fitting not successful Column 40: Structured uncertainty of sulfur dioxide effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -9=spectral fitting not successful Column 41: Common uncertainty of sulfur dioxide effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -2=no temperature fitting was requested and output for effective temperature and structured uncertainty of it is based on f-code, -3=spectral fitting was done, but no common uncertainty could be retrieved, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -6=no common uncertainty input was given, -9=spectral fitting not successful Column 42: Effective sulfur dioxide fitting wavelength [nm], -9=fitting not successful Column 43: Diffuse correction applied before fitting at effective fitting wavelength for sulfur dioxide [%], 0=no diffuse correction applied or fitting not requested, >0=measured diffuse correction, <0=(negative value of) calculated diffuse correction Column 44: Molecular scattering air mass factor used for molecular scattering subtraction before the fitting Column 45: Estimated uncertainty of molecular scattering air mass factor, -9=molecular scattering was not subtracted before the fitting Column 46: Lower limit used for wavelength scaling [nm] Column 47: Upper limit used for wavelength scaling [nm] Column 48: Order of smoothing polynomial Column 49: Smoothing polynomial coefficient, order 0 Column 50: Independent uncertainty of smoothing polynomial coefficient, order 0, -9=fitting not successful Column 51: Structured uncertainty of smoothing polynomial coefficient, order 0, -9=fitting not successful Column 52: Common uncertainty of smoothing polynomial coefficient, order 0, -9=fitting not successful Column 53: rms-based uncertainty of smoothing polynomial coefficient, order 0, -9=fitting not successful Column 54: Smoothing polynomial coefficient, order 1 Column 55: Independent uncertainty of smoothing polynomial coefficient, order 1, -9=fitting not successful Column 56: Structured uncertainty of smoothing polynomial coefficient, order 1, -9=fitting not successful Column 57: Common uncertainty of smoothing polynomial coefficient, order 1, -9=fitting not successful Column 58: rms-based uncertainty of smoothing polynomial coefficient, order 1, -9=fitting not successful Column 59: Smoothing polynomial coefficient, order 2 Column 60: Independent uncertainty of smoothing polynomial coefficient, order 2, -9=fitting not successful Column 61: Structured uncertainty of smoothing polynomial coefficient, order 2, -9=fitting not successful Column 62: Common uncertainty of smoothing polynomial coefficient, order 2, -9=fitting not successful Column 63: rms-based uncertainty of smoothing polynomial coefficient, order 2, -9=fitting not successful Column 64: Smoothing polynomial coefficient, order 3 Column 65: Independent uncertainty of smoothing polynomial coefficient, order 3, -9=fitting not successful Column 66: Structured uncertainty of smoothing polynomial coefficient, order 3, -9=fitting not successful Column 67: Common uncertainty of smoothing polynomial coefficient, order 3, -9=fitting not successful Column 68: rms-based uncertainty of smoothing polynomial coefficient, order 3, -9=fitting not successful Column 69: Mean value of measured data inside fitting window [same units as measurements] Column 70: Wavelength change polynomial coefficient, order 0 Column 71: Independent uncertainty of wavelength change polynomial coefficient, order 0, -9=fitting not successful Column 72: Structured uncertainty of wavelength change polynomial coefficient, order 0, -9=fitting not successful Column 73: Common uncertainty of wavelength change polynomial coefficient, order 0, -9=fitting not successful Column 74: rms-based uncertainty of wavelength change polynomial coefficient, order 0, -9=fitting not successful Column 75: L1 based wavelength change coefficient, order 0, -9=no L1 wavelength change determination Column 76: L1 based wavelength change coefficient, order 1, -9=no L1 wavelength change determination Column 77: L2Fit data quality flag: 0=assured high quality, 1=assured medium quality, 2=assured low quality, 10=not-assured high quality, 11=not-assured medium quality, 12=not-assured low quality Column 78: Sum over 2^i using those i, for which the corresponding L2Fit data quality parameter exceeds the DQ1 limit, 0=L1 data quality above 0, 1=Spectral fitting was not successful, 2=Wavelength shift too large, 3=Normalized rms of fitting residuals weighted with independent uncertainty too large Column 79: Sum over 2^i using those i, for which the corresponding L2Fit data quality parameter exceeds the DQ2 limit (same parameters as for DQ1) Column 80: L1 data quality flag: 0=assured high quality, 1=assured medium quality, 2=assured low quality, 10=not-assured high quality, 11=not-assured medium quality, 12=not-assured low quality Column 81: Sum over 2^i using those i, for which the corresponding L1 data quality parameter exceeds the DQ1 limit, 0=Saturated data, 1=Too few dark counts measurements, 2=No temperature given or effective temperature too different from the reference temperature, 3=Dark count too high, 4=Unsuccessful dark background fitting, 5=The dark count differs significantly from the dark map for too many pixels, 6=Absolute value of estimated average residual stray light level too high, 7=Although attempted, no wavelength change could be retrieved, 8=Absolute value of retrieved wavelength shift too large, 9=Retrieved wavelength shift differs too much from the shift predicted by the effective temperature Column 82: Sum over 2^i using those i, for which the corresponding L1 data quality parameter exceeds the DQ2 limit (same parameters as for DQ1) Column 83: Atmospheric variability [%], 999=no atmospheric variability was determined Column 84: Wavelength effective temperature [°C], 999=no effective temperature given Column 85: Estimated average residual stray light level [%] (only valid for stray light correction methods 2 and higher) Column 86: Retrieved wavelength shift from L1 data [nm], -9=no wavelength change determination Column 87: Retrieved total wavelength shift [nm], -9=no wavelength change fitting Column 88: Number of bright count cycles Column 89: Number of dark count cycles Column 90: Effective position of filterwheel #1, 0=filterwheel not used, 1-9 are valid positions Column 91: Effective position of filterwheel #2, 0=filterwheel not used, 1-9 are valid positions Column 92: Sum over 2^i, 0=spectra were interpolated in time, 1=spectra are corrected for off-target signal Column 93: Integration time [ms] Column 94: Temperature at electronics board [°C], 999=no temperature signal Column 95: Spectrometer control temperature [°C], 999=no temperature signal Column 96: Auxiliary spectrometer temperature [°C], 999=no temperature signal Columns 97-239: Unweighted slant column residuals for each pixel inside the fitting window multiplied by 1e5, 9e5=pixel was not used for fitting Columns 240-382: Normalized slant column residuals weighted with independent instrumental uncertainty for each pixel inside the fitting window multiplied by 1e5, 9e5=pixel was not used for fitting, 0=no independent instrumental uncertainty was given --------------------------------------------------------------------------------------- SS 20211103T220148.8Z 7977.917926 144 1 49.56 2 77.11 219.93 76.82 234.07 0.00 1 0.00 1 0 17 3.018e-01 2.380e-02 5.098e-02 3.516e-03 5.247669e-01 1.754239e-03 6.623227e-03 -6.000000e+00 1.682522e-01 221.61 -2.00 3.88 -2.00 319.97 0.00e+00 5.596426e-03 5.252920e-04 4.188534e-04 -6.000000e+00 7.269746e-03 268.08 -2.00 20.00 -2.00 319.53 0.00e+00 4.368 0.002 302.00 328.00 3 4.2701e+00 4.9096e-03 4.8902e-02 -6.0000e+00 4.1411e-01 1.8947e-01 5.8605e-03 2.9470e-02 -6.0000e+00 4.5415e-01 -1.7925e-01 5.1074e-03 1.5639e-02 -6.0000e+00 2.6823e-01 7.0889e-02 2.6762e-03 7.7382e-03 -6.0000e+00 1.1242e-01 7.5133e-05 -1.7062e-04 9.8206e-05 7.1443e-05 -6.0000e+00 8.1639e-03 2.6619e-03 0.0000e+00 12 9 8 11 32 0 71.64 21.22 0.09 0.00266 0.00138 23 5 1 5 0 2154.181 26.51 20.00 21.22 -81459 -24403 -103740 -112223 900000 -94614 900000 -55482 -88003 5637 -159196 12400 -42354 -39913 31157 80687 -22904 -8302 36700 -3741 -13965 -14655 -71897 91037 29488 3341 24547 -25687 21235 28537 101658 40033 -2285 31621 -1272 2231 30539 4652 10646 12314 -4297 4193 8225 20152 4537 11390 4115 14031 1684 2998 -2584 11662 1568 9526 11834 5584 4925 9826 1459 2693 10903 7365 -1108 2623 -1447 73 -51 2865 1054 2642 -1223 737 -556 3003 13 -1081 -5913 3310 -5697 1599 -2909 -464 -4975 -5229 256 -255 -561 159 1074 657 -1314 -2783 -2283 -4846 -977 -2659 670 -91 -1767 1084 549 -1881 672 1815 1343 -2890 321 -3233 1028 135 1106 436 532 129 -1830 386 173 1061 801 -1541 -38 1662 1357 2042 2381 316 2029 963 542 -753 877 -207 854 459 -408 314 455 -173 -736 -689 -659 -1105 -699 -2284 -449 -4460 -5035 900000 -3007 900000 -1383 -3345 88 -13094 201 -1465 -1716 815 1659 -1334 -453 1361 -217 -909 -945 -8333 2402 1526 232 1380 -2097 1172 1545 3066 2634 -237 2688 -143 236 2638 617 1463 1941 -884 843 1655 3679 981 2342 922 2898 404 745 -741 3201 515 3078 3886 2008 1808 3389 532 960 3730 2882 -532 1336 -788 38 -25 1399 581 1499 -784 478 -357 1753 8 -586 -3126 1615 -3403 955 -2096 -366 -4414 -4527 224 -239 -552 159 1051 617 -1168 -2512 -1951 -4391 -950 -2785 776 -111 -2327 1362 687 -2161 710 1960 1569 -3717 420 -4332 1442 211 1726 704 840 207 -2797 554 253 1621 1305 -2479 -60 2647 2090 3013 3357 434 2669 1353 855 -1243 1526 -380 1613 864 -824 655 965 -373 -1578 -1415 -1273 -2167 -1356 SS 20211103T221238.5Z 7977.925445 147 1 50.74 2 78.17 222.37 78.14 236.36 0.00 1 0.00 1 0 17 3.689e-01 3.738e-02 9.659e-02 4.462e-03 5.447836e-01 2.269895e-03 7.422907e-03 -6.000000e+00 2.057357e-01 221.61 -2.00 3.88 -2.00 320.20 0.00e+00 4.553053e-03 7.205336e-04 3.334276e-03 -6.000000e+00 8.926422e-03 268.08 -2.00 20.00 -2.00 319.77 0.00e+00 4.731 0.002 302.00 328.00 3 4.3350e+00 6.4090e-03 5.3107e-02 -6.0000e+00 5.0628e-01 2.5973e-01 7.7910e-03 3.2560e-02 -6.0000e+00 5.5532e-01 -3.1354e-01 6.8616e-03 2.0114e-02 -6.0000e+00 3.2796e-01 1.4082e-01 3.5477e-03 8.4512e-03 -6.0000e+00 1.3780e-01 5.0374e-05 -4.4242e-04 1.2635e-04 1.0099e-04 -6.0000e+00 9.9865e-03 2.1497e-03 0.0000e+00 12 9 8 11 32 0 3.40 21.24 0.44 0.00215 -0.00117 19 4 1 5 0 2669.786 26.51 20.00 21.24 -61299 14888 88704 -134535 900000 -75462 900000 -41733 -83593 56185 -186000 -112737 -22291 -4121 55759 168977 4888 8779 82921 9250 21814 -9476 -84805 900000 68772 18378 41743 -10651 70384 48862 900000 70334 9856 49856 -22 8535 73095 29750 22721 27730 3703 16858 24900 39101 10264 31796 13985 23709 10266 9141 -1664 24496 501 18084 27598 12399 7579 24856 8027 7595 18379 14183 -1517 2983 -2127 1404 1217 1833 1921 4850 -1339 471 12 5983 820 5232 -6762 3470 -8094 1460 -4697 815 -5664 -479 -282 -281 -251 1062 614 -15046 -1294 -4938 -2853 -7014 -3436 -4637 -42 -817 -2892 1865 156 -2893 275 3036 1945 -2677 1551 -3650 1261 -72 1503 96 993 -47 -2178 943 253 1594 1499 -1713 1371 3382 2759 3489 3327 1289 2724 1102 -6 -875 1030 75 639 -136 -377 594 472 -830 -815 -842 -1154 -775 -1712 -1577 205 666 -8189 900000 -2297 900000 -1039 -3469 600 -20689 -7130 -696 -135 1234 1403 230 427 2032 490 1026 -596 -11430 900000 2229 1017 1815 -749 2356 1931 900000 3333 880 3114 -2 739 3941 2658 2631 3562 667 2825 3996 5533 1951 4939 2620 4095 2077 1966 -435 5432 153 4925 7112 3828 2490 6769 2503 2345 5281 4681 -660 1368 -1056 645 538 822 955 2671 -782 279 7 3363 442 2413 -3271 1535 -4498 866 -3142 579 -4633 -393 -246 -261 -244 1044 599 -15167 -1144 -4192 -2264 -6005 -3173 -4983 -46 -1005 -3877 2346 184 -3149 296 3278 2111 -3214 1879 -4958 1787 -106 2369 159 1593 -73 -3198 1391 381 2317 2319 -2827 2171 5053 4000 5219 4783 1806 3664 1594 -8 -1490 1847 142 1168 -265 -783 1271 1030 -1852 -1798 -1777 -2294 -1449 -3202 SS 20211103T222123.1Z 7977.931517 152 1 49.72 2 79.06 224.31 79.23 238.19 0.00 1 0.00 1 0 174 1.559e+00 6.712e-01 2.007e-01 5.516e-03 -9.912734e-01 2.750868e-03 2.045617e-01 -6.000000e+00 8.705157e-01 221.61 -2.00 3.88 -2.00 319.26 0.00e+00 -8.788736e-02 8.407240e-04 1.915728e-01 -6.000000e+00 3.772089e-02 268.08 -2.00 20.00 -2.00 318.55 0.00e+00 5.086 0.003 302.00 328.00 3 7.6944e+00 7.6519e-03 7.2666e-01 -6.0000e+00 2.1422e+00 -3.5573e+00 9.7640e-03 1.3183e+00 -6.0000e+00 2.3488e+00 2.9041e+00 7.2875e-03 1.0397e+00 -6.0000e+00 1.3883e+00 -1.0898e+00 3.8225e-03 3.5651e-01 -6.0000e+00 5.8426e-01 3.5333e-05 -2.8885e-03 1.0757e-04 4.4531e-03 -6.0000e+00 4.2279e-02 -6.8399e-03 0.0000e+00 12 8 8 10 0 0 65.27 21.25 -0.39 -0.00684 -0.02855 13 3 1 5 0 3823.434 26.51 20.00 21.25 -61365 197961 -141315 -146255 900000 17399 900000 218319 177038 422497 -83597 168165 243680 116836 900000 301177 45598 43456 275538 67719 107852 123185 32248 900000 381394 327308 379054 264697 408193 319737 900000 294514 97335 218903 85482 143406 225095 79513 215228 268301 220294 249677 258393 291326 223138 232060 187575 196632 146544 146643 106650 156726 120143 150941 174468 168088 186435 227662 214894 211624 227875 208815 169007 159319 138106 128214 112043 103672 92264 101071 104760 118424 135358 164777 175268 188843 161755 156003 98481 98265 74723 76098 65812 67549 89947 97216 105559 114046 121008 -118772 125998 119360 106500 81878 62656 40967 31215 18065 9007 15808 21101 32613 58735 75017 84478 82476 87687 74573 67384 48692 35621 22213 13038 3012 -8859 -11039 -14134 -10728 -1956 9567 32121 50020 56066 58119 55309 51493 47373 29249 8435 -10237 -20066 -31578 -38507 -42497 -41863 -33123 -25523 -19903 -12147 1600 16282 30674 35402 -2061 1064 -8296 -8888 900000 304 900000 2347 5709 2271 -68231 12050 5715 5547 900000 1619 2226 1946 1481 3044 4110 5357 4723 900000 5777 12953 10371 19528 5484 6092 900000 3957 7752 5241 7977 8408 7698 16096 19349 21105 36272 34279 32319 24624 34451 29819 32626 26131 27422 25005 27252 26112 31698 32891 35034 43855 49826 47273 55016 56868 51660 56725 60933 59295 56472 48005 40495 36402 38099 45357 54097 62426 69716 72322 75043 74550 65418 56936 47943 45506 45581 48941 50314 51579 69128 75197 85097 91939 96001 -526841 100672 91679 78272 64503 53233 37582 30894 20545 11286 18122 22733 33059 52473 67981 81723 90459 96790 94308 87333 66170 51457 33725 19309 4274 -12260 -14200 -18497 -14525 -2807 14844 47424 68819 74813 79243 73422 67163 59146 39589 11962 -16411 -33609 -56077 -66087 -77929 -82976 -66685 -52143 -41682 -25378 3154 30289 57298 66294