File name: Pandora32s1_GreenbeltMD_20231207_L2Fit_fwvt1c12d20221109p1-8.txt File generation date: 20240223T220257.0Z Data description: Level 2 spectral fitting results file Data file version: fwvt1c12p1-8 Local principal investigator: Thomas Hanisco Network principal investigator: Alexander Cede Instrument type: Pandora Instrument number: 32 Spectrometer number: 1 Processing software version used: BlickP v1.8.59 Instrument operation file used: Pandora32_OF_v1d20220415.txt Instrument calibration file used: Pandora32s1_CF_v12d20221109.txt Level 1 file used: Pandora32s1_GreenbeltMD_20231207_L1_smca1c12d20221109p1-8.txt Full location name: Goddard Space Flight Center Short location name: GreenbeltMD Country of location: USA Location latitude [deg]: 38.9926 Location longitude [deg]: -76.8396 Location altitude [m]: 90 Local noon date: 20231207 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: Based on session 11 with new SO2 and HCHO field cal First and last pixel inside fitting window: 1680 1830 Nominal wavelengths inside fitting window [nm]: 492.0337 492.1546 492.2756 492.3965 492.5174 492.6382 492.7591 492.8799 493.0008 493.1216 493.2424 493.3631 493.4839 493.6046 493.7254 493.8461 493.9668 494.0874 494.2081 494.3287 494.4493 494.5700 494.6905 494.8111 494.9317 495.0522 495.1727 495.2933 495.4137 495.5342 495.6547 495.7751 495.8955 496.0159 496.1363 496.2567 496.3771 496.4974 496.6177 496.7380 496.8583 496.9786 497.0989 497.2191 497.3393 497.4595 497.5797 497.6999 497.8201 497.9402 498.0603 498.1804 498.3005 498.4206 498.5406 498.6607 498.7807 498.9007 499.0207 499.1407 499.2606 499.3806 499.5005 499.6204 499.7403 499.8602 499.9800 500.0999 500.2197 500.3395 500.4593 500.5790 500.6988 500.8185 500.9383 501.0580 501.1776 501.2973 501.4170 501.5366 501.6562 501.7758 501.8954 502.0150 502.1345 502.2541 502.3736 502.4931 502.6126 502.7321 502.8515 502.9709 503.0904 503.2098 503.3291 503.4485 503.5679 503.6872 503.8065 503.9258 504.0451 504.1644 504.2836 504.4028 504.5221 504.6413 504.7604 504.8796 504.9988 505.1179 505.2370 505.3561 505.4752 505.5942 505.7133 505.8323 505.9513 506.0703 506.1893 506.3083 506.4272 506.5462 506.6651 506.7840 506.9028 507.0217 507.1405 507.2594 507.3782 507.4970 507.6157 507.7345 507.8533 507.9720 508.0907 508.2094 508.3281 508.4467 508.5653 508.6840 508.8026 508.9212 509.0397 509.1583 509.2768 509.3954 509.5139 509.6324 509.7508 509.8693 509.9877 --------------------------------------------------------------------------------------- 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: Water vapor slant column amount [moles per square meter], -9e99=fitting not successful Column 23: Independent uncertainty of water vapor 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 water vapor 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 water vapor 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 water vapor 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: Water vapor effective temperature [K] Column 28: Independent uncertainty of water vapor 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 water vapor 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 water vapor 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 water vapor fitting wavelength [nm], -9=fitting not successful Column 32: Diffuse correction applied before fitting at effective fitting wavelength for water vapor [%], 0=no diffuse correction applied or fitting not requested, >0=measured diffuse correction, <0=(negative value of) calculated diffuse correction Column 33: Nitrogen dioxide slant column amount [moles per square meter], -9e99=fitting not successful Column 34: Independent uncertainty of nitrogen 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 nitrogen 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 nitrogen 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 nitrogen 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: Nitrogen dioxide effective temperature [K] Column 39: Independent uncertainty of nitrogen 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 nitrogen 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 nitrogen 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 nitrogen dioxide fitting wavelength [nm], -9=fitting not successful Column 43: Diffuse correction applied before fitting at effective fitting wavelength for nitrogen dioxide [%], 0=no diffuse correction applied or fitting not requested, >0=measured diffuse correction, <0=(negative value of) calculated diffuse correction Column 44: Ozone slant column amount [moles per square meter], -9e99=fitting not successful Column 45: 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 46: 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 47: 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 48: 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 49: Ozone effective temperature [K] Column 50: 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 51: 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 52: 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 53: Effective ozone fitting wavelength [nm], -9=fitting not successful Column 54: 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 55: Molecular scattering air mass factor used for molecular scattering subtraction before the fitting Column 56: Estimated uncertainty of molecular scattering air mass factor, -9=molecular scattering was not subtracted before the fitting Column 57: Lower limit used for wavelength scaling [nm] Column 58: Upper limit used for wavelength scaling [nm] Column 59: Order of smoothing polynomial Column 60: Smoothing polynomial coefficient, order 0 Column 61: Independent uncertainty of smoothing polynomial coefficient, order 0, -9=fitting not successful Column 62: Structured uncertainty of smoothing polynomial coefficient, order 0, -9=fitting not successful Column 63: Common uncertainty of smoothing polynomial coefficient, order 0, -9=fitting not successful Column 64: rms-based uncertainty of smoothing polynomial coefficient, order 0, -9=fitting not successful Column 65: Smoothing polynomial coefficient, order 1 Column 66: Independent uncertainty of smoothing polynomial coefficient, order 1, -9=fitting not successful Column 67: Structured uncertainty of smoothing polynomial coefficient, order 1, -9=fitting not successful Column 68: Common uncertainty of smoothing polynomial coefficient, order 1, -9=fitting not successful Column 69: rms-based uncertainty of smoothing polynomial coefficient, order 1, -9=fitting not successful Column 70: Mean value of measured data inside fitting window [same units as measurements] Column 71: Wavelength change polynomial coefficient, order 0 Column 72: Independent uncertainty of wavelength change polynomial coefficient, order 0, -9=fitting not successful Column 73: Structured uncertainty of wavelength change polynomial coefficient, order 0, -9=fitting not successful Column 74: Common uncertainty of wavelength change polynomial coefficient, order 0, -9=fitting not successful Column 75: rms-based uncertainty of wavelength change polynomial coefficient, order 0, -9=fitting not successful Column 76: L1 based wavelength change coefficient, order 0, -9=no L1 wavelength change determination Column 77: L1 based wavelength change coefficient, order 1, -9=no L1 wavelength change determination Column 78: 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 79: 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 80: 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 81: 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 82: 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 83: 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 84: Atmospheric variability [%], 999=no atmospheric variability was determined Column 85: Wavelength effective temperature [°C], 999=no effective temperature given Column 86: Estimated average residual stray light level [%] (only valid for stray light correction methods 2 and higher) Column 87: Retrieved wavelength shift from L1 data [nm], -9=no wavelength change determination Column 88: Retrieved total wavelength shift [nm], -9=no wavelength change fitting Column 89: Number of bright count cycles Column 90: Number of dark count cycles Column 91: Effective position of filterwheel #1, 0=filterwheel not used, 1-9 are valid positions Column 92: Effective position of filterwheel #2, 0=filterwheel not used, 1-9 are valid positions Column 93: Sum over 2^i, 0=spectra were interpolated in time, 1=spectra are corrected for off-target signal Column 94: Integration time [ms] Column 95: Temperature at electronics board [°C], 999=no temperature signal Column 96: Spectrometer control temperature [°C], 999=no temperature signal Column 97: Auxiliary spectrometer temperature [°C], 999=no temperature signal Column 98: Temperature in head sensor [°C], 999=no temperature signal Column 99: Humidity in head sensor [%], -9=no humidity signal Column 100: Pressure in head sensor [hPa], -9=no pressure signal Columns 101-251: Unweighted slant column residuals for each pixel inside the fitting window multiplied by 1e5, 9e5=pixel was not used for fitting Columns 252-402: 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 --------------------------------------------------------------------------------------- SQ 20231207T212846.9Z 8741.894988 534 1 7.27 2 88.14 238.53 122.86 292.13 0.00 1 0.00 1 0 22 6.240e-03 6.188e-03 8.441e-04 8.419e-04 1.484871e+04 6.065481e+01 9.576492e+02 -6.000000e+00 3.547417e+02 279.14 -2.00 20.00 -2.00 505.41 0.00e+00 4.946953e-03 2.362390e-05 2.848873e-04 -6.000000e+00 1.761157e-04 279.14 -2.00 20.00 -2.00 500.98 0.00e+00 8.947662e-01 1.652762e-02 1.177168e-01 -6.000000e+00 1.248590e-01 223.25 -2.00 3.47 -2.00 502.49 0.00e+00 16.487 0.079 489.00 513.00 1 3.2994e+00 1.2672e-03 2.8551e-02 -6.0000e+00 9.5998e-03 -3.2597e-02 3.1187e-04 3.0165e-03 -6.0000e+00 2.4382e-03 6.0188e-03 -1.5242e-04 8.7755e-05 1.5156e-04 -6.0000e+00 7.1721e-04 5.8430e-03 -5.5185e-05 12 9 9 12 1 1 75.06 21.00 5.92 0.00584 0.00479 15 25 1 4 0 454.400 25.68 19.98 21.00 10.85 3.8 918.01 -1038 -829 -302 107 577 902 1314 1349 1194 975 482 294 -152 -267 -411 -521 -381 -86 305 413 325 469 434 469 538 576 382 346 -188 -700 -1041 -877 -914 -154 189 52 10 -52 -122 -322 -79 -31 128 373 578 752 528 426 210 -227 -668 -1107 -1374 -1253 -731 -353 10 397 263 582 563 604 717 663 509 380 73 -99 -305 -488 -688 -694 -524 -553 -269 -421 -644 -827 -843 -829 -870 -955 -794 -863 -764 -469 -626 -432 -565 -263 -314 132 426 712 390 383 52 -228 -530 -998 -1186 -1275 -887 -480 -199 -245 -160 -413 -411 -203 54 309 441 734 966 1222 1339 1433 1019 858 464 533 444 717 515 505 572 335 281 313 241 39 -194 -253 -512 -313 -173 45 102 401 375 283 305 357 424 288 9 -558 -743 -665 -906 -934 -752 -274 100 543 868 1258 1305 1160 943 463 279 -145 -253 -389 -495 -360 -82 298 401 318 462 429 457 527 568 377 338 -185 -678 -994 -836 -886 -149 188 51 10 -51 -120 -316 -78 -30 129 373 584 764 539 432 210 -228 -660 -1076 -1352 -1215 -721 -350 11 400 267 592 575 619 739 685 524 388 74 -99 -308 -489 -685 -692 -525 -550 -269 -421 -638 -824 -830 -811 -851 -939 -784 -861 -756 -466 -625 -429 -560 -260 -312 132 429 714 387 377 50 -216 -503 -936 -1111 -1190 -851 -469 -198 -245 -161 -417 -414 -205 55 314 453 758 986 1255 1390 1502 1059 888 484 548 451 727 526 519 592 348 293 325 250 40 -198 -257 -524 -323 -179 47 107 426 397 303 329 387 458 310 9 -591 -791 -705 -955 SQ 20231207T213136.0Z 8741.896944 538 1 48.60 2 88.61 238.98 123.36 292.68 0.00 1 0.00 1 0 22 6.482e-03 6.433e-03 1.530e-04 1.525e-04 1.577474e+04 1.107130e+01 1.015356e+03 -6.000000e+00 3.685328e+02 279.14 -2.00 20.00 -2.00 505.42 0.00e+00 5.344037e-03 4.286848e-06 3.053761e-04 -6.000000e+00 1.829781e-04 279.14 -2.00 20.00 -2.00 501.01 0.00e+00 9.404954e-01 2.993703e-03 1.246991e-01 -6.000000e+00 1.297462e-01 223.25 -2.00 3.47 -2.00 502.51 0.00e+00 17.749 0.098 489.00 513.00 1 1.6593e+00 2.2962e-04 3.1478e-02 -6.0000e+00 9.9752e-03 -3.2657e-02 5.6372e-05 3.2320e-03 -6.0000e+00 2.5344e-03 2.5704e-02 3.5307e-04 1.5769e-05 1.6094e-04 -6.0000e+00 7.4512e-04 2.0144e-03 -5.1987e-05 12 8 8 10 0 0 89.29 21.00 6.56 0.00201 0.00446 608 143 1 4 0 79.920 25.68 20.00 21.00 10.83 3.8 918.19 -1098 -863 -314 206 637 932 1243 1379 1221 919 519 247 -73 -351 -549 -554 -431 -113 291 322 311 360 466 571 682 649 468 316 -43 -707 -1049 -967 -754 -283 155 140 -29 -22 -136 -275 -184 -12 222 418 630 718 638 466 199 -213 -678 -1035 -1242 -1143 -735 -241 112 273 356 426 494 688 671 673 562 359 99 -169 -328 -581 -759 -765 -626 -482 -300 -383 -688 -881 -864 -953 -969 -898 -859 -788 -693 -588 -584 -581 -523 -371 -249 81 462 623 543 403 119 -214 -617 -1067 -1333 -1324 -977 -457 -246 -222 -324 -410 -403 -271 58 242 557 878 1162 1347 1457 1483 1184 854 595 432 445 621 678 557 541 420 280 389 257 -21 -241 -380 -510 -396 -126 -95 65 279 370 344 344 452 411 256 14 -452 -775 -895 -832 -979 -776 -289 191 592 887 1178 1320 1172 879 493 237 -70 -330 -516 -522 -411 -107 282 309 302 351 456 561 674 632 456 313 -42 -678 -993 -915 -724 -278 153 140 -29 -22 -135 -273 -183 -12 222 423 632 725 646 469 202 -213 -667 -1021 -1217 -1125 -722 -242 113 280 360 432 503 704 690 693 577 366 100 -172 -330 -581 -755 -762 -627 -487 -299 -382 -694 -877 -850 -932 -946 -880 -847 -784 -696 -583 -580 -577 -517 -365 -247 80 462 622 535 394 115 -202 -583 -996 -1244 -1254 -933 -445 -243 -221 -325 -412 -404 -273 59 245 569 903 1201 1403 1507 1551 1248 899 619 444 461 641 691 572 560 437 292 405 267 -22 -247 -387 -523 -410 -131 -101 68 298 399 369 371 489 445 277 14 -487 -827 -955 -894 SQ 20231207T213610.2Z 8741.900118 543 1 48.32 2 89.37 239.71 124.16 293.58 0.00 1 0.00 1 0 22 9.400e-03 9.298e-03 4.784e-04 4.758e-04 1.660766e+04 3.473718e+01 1.045142e+03 -6.000000e+00 5.343631e+02 279.14 -2.00 20.00 -2.00 505.46 0.00e+00 7.090587e-03 1.339050e-05 3.999056e-04 -6.000000e+00 2.653175e-04 279.14 -2.00 20.00 -2.00 501.17 0.00e+00 8.165899e-01 9.276208e-03 1.412108e-01 -6.000000e+00 1.881394e-01 223.25 -2.00 3.47 -2.00 502.64 0.00e+00 19.576 0.132 489.00 513.00 1 3.6077e+00 7.1081e-04 3.6776e-02 -6.0000e+00 1.4464e-02 -4.1687e-02 1.7391e-04 3.5774e-03 -6.0000e+00 3.6752e-03 2.7880e-03 9.5254e-05 4.9459e-05 1.9080e-04 -6.0000e+00 1.0804e-03 5.4961e-03 -7.7391e-05 12 8 8 10 0 0 98.78 21.00 10.71 0.00550 0.00616 84 20 1 4 0 575.170 25.68 19.98 21.00 11.08 3.8 918.76 -1782 -1469 -696 -42 857 1267 1794 1956 1647 1361 810 279 -147 -522 -718 -865 -654 -260 437 482 404 577 623 798 983 942 772 511 -115 -987 -1478 -1468 -1197 -447 257 259 163 89 -175 -310 -212 -11 371 637 942 1135 1035 812 374 -211 -805 -1560 -1794 -1608 -971 -352 201 421 633 704 744 994 1023 1034 776 615 173 -214 -497 -828 -1091 -1138 -914 -656 -531 -467 -1003 -1291 -1324 -1500 -1406 -1368 -1162 -1074 -963 -879 -835 -752 -763 -548 -312 117 676 862 831 568 212 -293 -937 -1488 -1876 -1843 -1485 -650 -215 -147 -408 -453 -607 -367 -9 384 931 1307 1449 1963 2018 1951 1668 1276 928 567 648 873 959 815 727 645 361 578 347 8 -360 -627 -759 -647 -334 -125 56 404 416 416 468 657 639 312 53 -623 -1182 -1374 -1291 -1550 -1287 -623 -38 785 1185 1657 1825 1544 1270 750 256 -134 -485 -666 -784 -598 -241 413 462 380 545 591 761 944 911 748 493 -110 -928 -1374 -1365 -1129 -431 247 252 160 87 -170 -301 -207 -11 370 630 942 1143 1046 818 374 -208 -783 -1539 -1759 -1585 -948 -352 200 425 645 720 763 1025 1061 1075 806 634 177 -217 -500 -827 -1086 -1133 -914 -661 -538 -473 -1011 -1286 -1306 -1472 -1376 -1342 -1145 -1067 -965 -887 -843 -756 -764 -547 -312 117 684 872 829 554 205 -277 -871 -1405 -1727 -1719 -1424 -629 -215 -147 -414 -459 -616 -375 -9 393 959 1356 1515 2015 2098 2056 1773 1356 977 590 667 897 990 850 764 682 384 613 367 8 -377 -654 -796 -686 -352 -134 60 435 454 458 519 730 709 346 59 -680 -1281 -1489 -1411