File name: Pandora29s1_Fajardo_20220616_L2Fit_fout2c2d20220426p1-8.txt File generation date: 20260409T134842.6Z Data description: Level 2 spectral fitting results file Data file version: fout2c2p1-8 Local principal investigator: Tom Hanisco Network principal investigator: Alexander Cede Instrument type: Pandora Instrument number: 29 Spectrometer number: 1 Processing software version used: BlickP v1.8.88 Instrument operation file used: Pandora29_OF_v2d20211130.txt Instrument calibration file used: Pandora29s1_CF_v2d20220426.txt Level 1 file used: Pandora29s1_Fajardo_20220616_L1_smca1c2d20220426p1-8.txt Full location name: Cape San Juan Light, Fajardo (Puerto Rico) Short location name: Fajardo Country of location: Puerto Rico Location latitude [deg]: 18.3797 Location longitude [deg]: -65.6184 Location altitude [m]: 66 Local noon date: 20220616 Notes on s-number (L1 configuration): Corrections NOT applied although requested by the s-number are latency correction, matrix method stray light correction (replaced by simple method) Data caveats: None First and last pixel inside fitting window: 243 385 Nominal wavelengths inside fitting window [nm]: 305.0193 305.1603 305.3013 305.4423 305.5833 305.7243 305.8652 306.0062 306.1471 306.2880 306.4289 306.5699 306.7107 306.8516 306.9925 307.1334 307.2742 307.4150 307.5559 307.6967 307.8375 307.9783 308.1190 308.2598 308.4006 308.5413 308.6820 308.8228 308.9635 309.1042 309.2449 309.3855 309.5262 309.6668 309.8075 309.9481 310.0887 310.2293 310.3699 310.5105 310.6511 310.7917 310.9322 311.0727 311.2133 311.3538 311.4943 311.6348 311.7753 311.9157 312.0562 312.1966 312.3371 312.4775 312.6179 312.7583 312.8987 313.0391 313.1794 313.3198 313.4601 313.6005 313.7408 313.8811 314.0214 314.1617 314.3019 314.4422 314.5825 314.7227 314.8629 315.0031 315.1433 315.2835 315.4237 315.5639 315.7040 315.8442 315.9843 316.1244 316.2646 316.4047 316.5447 316.6848 316.8249 316.9649 317.1050 317.2450 317.3850 317.5250 317.6650 317.8050 317.9450 318.0849 318.2249 318.3648 318.5048 318.6447 318.7846 318.9245 319.0643 319.2042 319.3441 319.4839 319.6237 319.7636 319.9034 320.0432 320.1829 320.3227 320.4625 320.6022 320.7420 320.8817 321.0214 321.1611 321.3008 321.4405 321.5802 321.7198 321.8595 321.9991 322.1387 322.2783 322.4179 322.5575 322.6971 322.8367 322.9762 323.1157 323.2553 323.3948 323.5343 323.6738 323.8133 323.9527 324.0922 324.2316 324.3711 324.5105 324.6499 324.7893 324.9287 --------------------------------------------------------------------------------------- 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 Column 97: Temperature in head sensor [°C], 999=no temperature signal Column 98: Humidity in head sensor [%], -9=no humidity signal Column 99: Pressure in head sensor [hPa], -9=no pressure signal Columns 100-242: Unweighted slant column residuals for each pixel inside the fitting window multiplied by 1e5, 9e5=pixel was not used for fitting Columns 243-385: 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 20220616T141646.0Z 8202.594977 359 1 9.87 2 29.90 74.91 120.91 250.47 0.00 1 0.00 1 0 17 1.541e-02 1.461e-02 4.805e-03 3.554e-03 1.198326e-01 1.806979e-03 1.486484e-03 -6.000000e+00 8.692588e-03 228.70 -2.00 3.11 -2.00 314.45 0.00e+00 3.347758e-04 1.472964e-04 2.438588e-04 -6.000000e+00 4.023726e-04 295.00 -2.00 20.00 -2.00 313.77 0.00e+00 1.153 0.000 302.00 328.00 3 5.7635e+00 4.6241e-03 1.2613e-02 -6.0000e+00 2.1565e-02 -4.0520e-02 5.1214e-03 5.6173e-03 -6.0000e+00 2.3357e-02 4.8943e-02 2.9369e-03 2.8050e-03 -6.0000e+00 1.4009e-02 -2.5672e-02 1.2980e-03 1.1733e-03 -6.0000e+00 6.0400e-03 5.0539e-04 7.0264e-04 9.6011e-05 1.7093e-05 -6.0000e+00 3.9479e-04 1.2462e-02 5.4597e-05 12 9 9 12 41 1 99.96 22.82 0.36 0.01246 0.01774 3 3 1 5 0 2424.730 25.68 20.00 22.82 33.90 13.3 1088.31 -529 -1803 -2789 -2529 -2080 -822 901 3881 939 22 232 1090 2683 923 611 435 -1201 -672 -775 1449 557 1231 1999 2364 1612 -83 134 -974 164 503 3655 1468 86 -1440 -3877 -2354 -1473 77 -100 -521 -194 -1407 -926 -1041 81 1134 303 156 -111 -650 489 355 -464 107 877 -360 -201 -994 122 1133 2057 2406 1327 1590 -190 -1008 -207 185 1097 -550 -31 -786 -966 -199 -852 -102 -687 -1233 -389 -1681 -2500 -2335 -3581 -2312 -1848 66 -1346 -235 -386 507 -132 -322 498 973 2036 2315 3800 5143 4088 3344 3894 3355 2610 708 -942 -1904 -1991 -1689 133 -144 -267 -924 -1634 -536 -1079 -1095 -838 -950 -1040 -935 -1025 -798 778 -90 828 269 1174 1048 675 -241 -979 -361 -275 232 -199 -594 446 -134 112 -43 101 236 1053 -201 -752 -1196 -1017 -720 -283 316 1343 375 9 96 431 1121 435 296 229 -654 -367 -434 818 325 702 1134 1394 927 -50 81 -602 91 263 1884 806 51 -892 -2269 -1283 -857 47 -74 -420 -159 -1212 -814 -923 71 974 255 134 -92 -523 400 289 -399 95 794 -333 -189 -940 112 1024 1893 2352 1326 1673 -202 -974 -189 160 968 -527 -32 -825 -1023 -216 -877 -100 -600 -1061 -335 -1523 -2392 -2445 -3808 -2518 -2146 78 -1686 -303 -502 633 -153 -358 520 946 2021 2356 4182 5828 4971 4023 4372 3699 2743 788 -1110 -2281 -2489 -2191 176 -200 -386 -1284 -2285 -717 -1313 -1305 -987 -1170 -1311 -1214 -1385 -1072 1061 -119 1019 334 1360 1198 815 -294 -1251 -466 -356 308 -268 -846 647 -191 161 -63 147 331 1435 SS 20220616T142510.3Z 8202.600814 365 1 4.45 2 27.98 74.63 122.69 250.73 0.00 1 0.00 1 0 17 1.422e-02 1.434e-02 1.654e-03 1.331e-03 1.150574e-01 6.810092e-04 1.461874e-03 -6.000000e+00 8.019332e-03 228.70 -2.00 3.11 -2.00 313.95 0.00e+00 4.960220e-04 5.081249e-05 2.198053e-04 -6.000000e+00 3.712550e-04 295.00 -2.00 20.00 -2.00 313.33 0.00e+00 1.132 0.000 302.00 328.00 3 3.3915e+00 1.7331e-03 1.2350e-02 -6.0000e+00 1.9894e-02 -5.3756e-02 1.9136e-03 5.4485e-03 -6.0000e+00 2.1552e-02 4.7025e-02 1.1040e-03 2.7662e-03 -6.0000e+00 1.2926e-02 -2.4670e-02 4.8243e-04 1.2288e-03 -6.0000e+00 5.5701e-03 5.6150e-03 1.8722e-04 3.5721e-05 1.7977e-05 -6.0000e+00 3.6439e-04 1.3406e-02 -3.9333e-05 12 9 9 12 9 1 44.65 22.80 0.31 0.01341 0.01481 12 25 1 5 0 328.850 25.68 20.00 22.80 34.34 13.3 1091.74 -495 -1399 -1062 -3091 -1795 -326 2174 2149 996 222 -295 672 1002 551 184 -466 -897 372 46 848 377 644 1972 1640 697 29 -243 371 227 1830 2393 921 773 -2124 -2661 -1010 -402 107 -534 -608 -1134 -1083 -890 -585 209 -166 29 -3 -324 415 873 542 88 379 139 158 -219 -218 339 1262 1668 1857 1637 1365 -734 -944 -353 55 1130 -877 -1011 -874 -1042 -91 -558 -385 114 82 308 -1439 -1891 -2514 -3131 -2023 -1459 -1457 -1385 -1123 -363 -387 -439 -369 699 1964 2635 3177 4176 4975 4028 3386 3399 3188 2839 1533 -998 -2030 -2004 -1455 -322 -305 -726 -1107 -1035 -997 -915 -633 -366 -518 -564 -1100 -1151 -435 173 167 706 219 1165 1818 481 -155 -1000 -419 -432 -121 -545 -613 -521 -136 125 -225 235 814 1434 -221 -684 -526 -1467 -773 -136 915 942 465 110 -143 311 497 287 101 -279 -548 226 29 515 235 395 1195 1039 465 19 -161 241 141 1059 1400 569 481 -1407 -1718 -607 -259 71 -409 -500 -981 -970 -815 -518 189 -151 26 -2 -272 336 725 479 78 350 127 145 -205 -205 311 1140 1542 1760 1671 1422 -773 -932 -324 49 1035 -862 -1058 -925 -1114 -98 -565 -377 103 70 267 -1352 -1839 -2599 -3339 -2213 -1644 -1716 -1674 -1406 -456 -471 -509 -414 719 1934 2660 3347 4461 5621 4726 3927 3818 3434 3006 1652 -1138 -2478 -2473 -1859 -422 -407 -1011 -1505 -1378 -1295 -1109 -734 -433 -621 -706 -1412 -1539 -576 223 213 871 267 1369 2090 571 -189 -1278 -524 -543 -156 -726 -838 -733 -190 177 -314 327 1128 1933 SS 20220616T143009.2Z 8202.604273 368 1 52.14 2 26.84 74.41 123.74 250.88 0.00 1 0.00 1 0 17 1.421e-02 1.418e-02 1.125e-03 8.510e-04 1.151275e-01 4.337992e-04 1.454658e-03 -6.000000e+00 8.013794e-03 228.70 -2.00 3.11 -2.00 314.32 0.00e+00 4.708603e-04 3.453314e-05 2.399938e-04 -6.000000e+00 3.709833e-04 295.00 -2.00 20.00 -2.00 313.64 0.00e+00 1.120 0.000 302.00 328.00 3 4.9714e+00 1.1080e-03 1.2244e-02 -6.0000e+00 1.9880e-02 -4.4273e-02 1.2256e-03 5.4581e-03 -6.0000e+00 2.1535e-02 4.7389e-02 7.0474e-04 2.7392e-03 -6.0000e+00 1.2917e-02 -2.5475e-02 3.1107e-04 1.2274e-03 -6.0000e+00 5.5670e-03 1.1656e-03 5.1308e-04 2.2939e-05 1.7632e-05 -6.0000e+00 3.6406e-04 1.3697e-02 2.4967e-05 12 9 8 11 8 0 99.99 22.78 0.41 0.01370 0.01755 76 12 1 5 0 685.920 25.68 20.02 22.78 34.63 13.3 1093.53 -823 -1489 -1276 -2921 -2139 -129 2128 1975 1343 -416 -603 1318 1224 791 -538 113 -985 -120 -20 918 476 144 2092 1895 729 391 -153 731 245 2002 2479 1008 491 -1844 -2494 -1205 -311 277 -590 -917 -914 -1000 -1094 -728 150 -112 235 -114 46 34 976 179 -72 310 99 371 -214 -467 382 1282 1813 1908 1691 1362 -707 -935 -743 -82 901 -585 -872 -782 -946 -310 -420 -343 179 -442 -187 -1371 -1782 -2393 -3239 -1913 -1451 -1241 -1095 -900 -286 -494 -277 -331 645 1927 2403 3098 4137 5010 4092 3199 3496 3146 2910 1412 -1242 -1961 -2059 -1523 -46 -179 -498 -1060 -1123 -1178 -853 -528 -377 -561 -729 -827 -1043 -527 106 292 724 17 1026 1605 829 -145 -701 -545 -420 -197 -762 -690 -213 -125 131 -134 210 719 1420 -323 -639 -555 -1212 -825 -45 750 774 547 -182 -256 576 530 381 -271 61 -547 -66 -12 532 284 85 1208 1143 463 241 -95 451 138 1129 1312 564 292 -1160 -1454 -657 -181 184 -442 -757 -773 -877 -984 -658 133 -100 202 -94 38 28 808 149 -62 281 92 346 -193 -421 355 1176 1700 1799 1718 1388 -727 -919 -693 -72 807 -570 -890 -836 -1021 -323 -437 -321 157 -380 -162 -1253 -1718 -2541 -3489 -2109 -1621 -1495 -1328 -1132 -361 -606 -328 -374 680 1877 2406 3175 4373 5782 4818 3730 4005 3347 3094 1505 -1424 -2400 -2629 -2014 -63 -242 -704 -1507 -1519 -1534 -1052 -634 -448 -662 -931 -1093 -1371 -720 141 371 907 20 1208 1850 1014 -179 -907 -715 -553 -254 -1048 -957 -303 -182 194 -194 295 978 1967