File name: Pandora187s1_PittsburghPA_20230114_L2Fit_fwvt1c3d20220705p1-8.txt File generation date: 20251128T091133.5Z Data description: Level 2 spectral fitting results file Data file version: fwvt1c3p1-8 Local principal investigator: Lukas Valin Network principal investigator: Alexander Cede Instrument type: Pandora Instrument number: 187 Spectrometer number: 1 Processing software version used: BlickP v1.8.75 Instrument operation file used: Pandora187_OF_v4d20210517.txt Instrument calibration file used: Pandora187s1_CF_v3d20220705.txt Level 1 file used: Pandora187s1_PittsburghPA_20230114_L1_smca1c3d20220705p1-8.txt Full location name: Pittsburgh, PA - Lawrenceville ACHD (USEPA AQS ID 42-003-0008) Short location name: PittsburghPA Country of location: United States Location latitude [deg]: 40.4655 Location longitude [deg]: -79.9608 Location altitude [m]: 265 Local noon date: 20230114 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: 1607 1755 Nominal wavelengths inside fitting window [nm]: 492.0854 492.2072 492.3291 492.4509 492.5727 492.6944 492.8162 492.9379 493.0597 493.1814 493.3031 493.4247 493.5464 493.6680 493.7896 493.9112 494.0328 494.1544 494.2760 494.3975 494.5190 494.6405 494.7620 494.8835 495.0049 495.1264 495.2478 495.3692 495.4906 495.6119 495.7333 495.8546 495.9760 496.0973 496.2185 496.3398 496.4611 496.5823 496.7035 496.8247 496.9459 497.0671 497.1882 497.3093 497.4305 497.5516 497.6726 497.7937 497.9147 498.0358 498.1568 498.2778 498.3988 498.5197 498.6407 498.7616 498.8825 499.0034 499.1243 499.2452 499.3660 499.4868 499.6076 499.7284 499.8492 499.9700 500.0907 500.2114 500.3321 500.4528 500.5735 500.6942 500.8148 500.9354 501.0560 501.1766 501.2972 501.4177 501.5383 501.6588 501.7793 501.8998 502.0202 502.1407 502.2611 502.3816 502.5020 502.6223 502.7427 502.8630 502.9834 503.1037 503.2240 503.3443 503.4645 503.5848 503.7050 503.8252 503.9454 504.0656 504.1858 504.3059 504.4260 504.5461 504.6662 504.7863 504.9064 505.0264 505.1464 505.2664 505.3864 505.5064 505.6263 505.7463 505.8662 505.9861 506.1060 506.2258 506.3457 506.4655 506.5853 506.7051 506.8249 506.9447 507.0644 507.1842 507.3039 507.4236 507.5433 507.6629 507.7826 507.9022 508.0218 508.1414 508.2610 508.3805 508.5001 508.6196 508.7391 508.8586 508.9780 509.0975 509.2169 509.3363 509.4557 509.5751 509.6945 509.8138 509.9332 --------------------------------------------------------------------------------------- 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-249: Unweighted slant column residuals for each pixel inside the fitting window multiplied by 1e5, 9e5=pixel was not used for fitting Columns 250-398: 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 20230114T220034.0Z 8414.917060 563 1 45.45 2 88.17 239.78 144.05 330.81 0.00 1 0.00 1 0 15 4.242e-03 4.252e-03 1.680e-03 1.674e-03 3.810599e+03 1.158421e+02 2.416841e+02 -6.000000e+00 2.684213e+02 274.17 -2.00 20.00 -2.00 505.40 0.00e+00 2.936357e-03 4.943150e-05 1.592822e-04 -6.000000e+00 1.262095e-04 274.17 -2.00 20.00 -2.00 501.05 0.00e+00 1.359595e+00 3.521149e-02 4.755288e-02 -6.000000e+00 9.196985e-02 223.67 -2.00 3.49 -2.00 502.56 0.00e+00 16.639 0.084 489.00 513.00 1 5.2617e+00 2.7201e-03 2.6398e-02 -6.0000e+00 7.1197e-03 -1.7895e-02 6.5170e-04 1.6733e-03 -6.0000e+00 1.7464e-03 8.8895e-04 1.5367e-03 2.7470e-04 1.0301e-04 -6.0000e+00 7.1307e-04 2.0316e-03 -3.9888e-05 11 8 0 10 0 0 99.97 18.73 3.12 0.00203 0.01269 25 8 1 4 0 1817.480 22.09 15.00 18.73 2.20 22.5 895.95 -699 -388 -123 -161 61 153 665 493 414 270 104 99 -79 -356 -283 -379 147 186 313 168 160 420 346 599 632 175 366 136 -92 -556 -799 -318 236 198 49 -105 258 -62 -216 -289 330 223 139 515 570 -105 184 -116 -165 -614 -822 -575 -724 -365 -111 -184 165 136 130 470 803 780 458 317 116 50 110 -213 -478 -306 -446 -657 -55 -228 -303 -528 -692 -505 -265 -687 -708 -714 -207 -62 -462 -345 -363 -242 -200 169 420 743 265 178 26 145 -467 -582 -602 -470 -694 -603 -406 -530 -247 -731 -618 -52 -2 -3 409 333 406 781 811 796 948 577 840 306 329 728 340 413 361 232 0 522 263 498 -19 -428 30 -246 -187 -195 -5 15 -167 38 288 -65 157 8 10 -208 -403 -635 -581 -626 -351 -114 -153 56 143 628 469 394 255 98 92 -74 -329 -262 -354 138 177 303 165 158 417 346 599 633 175 364 134 -90 -538 -774 -309 232 189 48 -102 251 -60 -211 -282 322 220 138 514 571 -106 184 -114 -160 -592 -783 -543 -686 -348 -107 -182 165 137 132 478 817 795 469 323 117 50 109 -210 -470 -298 -436 -645 -53 -225 -298 -518 -675 -489 -255 -665 -688 -698 -203 -61 -460 -346 -365 -244 -203 172 429 762 261 174 26 145 -463 -569 -583 -455 -680 -600 -392 -521 -244 -731 -618 -52 -2 -3 418 346 426 824 862 849 1011 616 889 323 345 759 356 433 381 245 0 550 277 520 -19 -447 32 -257 -197 -200 -5 15 -178 40 311 -71 170 8 11 -222 -429 -673 -616 SQ 20230114T220344.0Z 8414.919259 567 1 4.04 2 88.69 240.30 144.34 331.99 0.00 1 0.00 1 0 15 5.754e-03 5.715e-03 3.436e-03 3.423e-03 4.262251e+03 2.388009e+02 2.672111e+02 -6.000000e+00 3.641630e+02 274.17 -2.00 20.00 -2.00 505.40 0.00e+00 3.508741e-03 1.013767e-04 1.911263e-04 -6.000000e+00 1.711836e-04 274.17 -2.00 20.00 -2.00 501.10 0.00e+00 1.361741e+00 7.177248e-02 5.475713e-02 -6.000000e+00 1.247249e-01 223.67 -2.00 3.49 -2.00 502.60 0.00e+00 18.055 0.107 489.00 513.00 1 3.8149e+00 5.5417e-03 2.9805e-02 -6.0000e+00 9.6559e-03 -3.3389e-02 1.3254e-03 1.8886e-03 -6.0000e+00 2.3671e-03 3.0852e-03 8.4208e-04 5.6036e-04 1.2328e-04 -6.0000e+00 9.6734e-04 -7.0868e-03 5.3023e-05 12 9 9 12 1 1 999.00 18.73 4.30 -0.00709 -0.00125 2 11 1 4 0 1284.020 22.09 15.00 18.73 2.15 22.6 895.82 -1203 -923 -180 -364 404 727 929 442 675 626 129 282 -463 -554 -471 -567 254 581 123 -58 366 433 806 1196 1126 531 354 255 -160 -861 -685 -511 -314 -11 275 -95 -554 -161 -130 393 22 171 524 599 583 167 328 -409 -339 -828 -1185 -1049 -1236 -394 -605 -210 701 316 126 421 1332 747 620 495 267 77 -176 -30 -612 -61 -618 -304 -248 -46 -141 -563 -678 -485 -110 -1080 -883 -940 -326 -407 -389 -578 -642 -340 -194 -165 281 839 604 300 226 -29 -451 -534 -413 -871 -1223 -906 -306 -280 -660 -241 -469 -537 -86 219 443 144 771 1268 965 1477 1202 620 557 410 209 344 505 745 404 259 317 344 219 -22 208 -311 -154 -171 -327 -346 178 128 124 -58 196 -43 422 150 -121 -423 -351 -834 -623 -1074 -833 -165 -330 372 681 881 424 646 596 122 264 -434 -517 -440 -534 241 558 116 -55 352 419 781 1160 1095 518 344 245 -153 -819 -649 -488 -304 -10 264 -92 -539 -156 -126 380 21 168 518 597 585 168 328 -407 -333 -802 -1138 -999 -1184 -379 -594 -210 702 322 125 419 1317 744 618 492 274 78 -178 -30 -610 -61 -614 -302 -247 -46 -140 -561 -673 -478 -107 -1068 -875 -938 -326 -397 -381 -571 -636 -337 -194 -166 284 847 611 303 226 -28 -440 -514 -408 -865 -1188 -893 -305 -284 -679 -248 -470 -540 -86 223 456 150 811 1341 1031 1579 1249 646 578 437 223 364 535 789 417 269 330 357 227 -23 221 -319 -163 -176 -341 -366 190 138 136 -62 212 -46 458 162 -130 -454 -373 -888 -663 SQ 20230114T220852.0Z 8414.922824 575 1 45.45 2 89.55 241.15 144.78 333.95 0.00 1 0.00 1 0 15 7.262e-03 7.224e-03 9.303e-04 9.239e-04 5.402463e+03 6.500627e+01 3.238559e+02 -6.000000e+00 4.595532e+02 274.17 -2.00 20.00 -2.00 505.44 0.00e+00 4.138049e-03 2.737895e-05 2.226272e-04 -6.000000e+00 2.160511e-04 274.17 -2.00 20.00 -2.00 501.29 0.00e+00 1.234710e+00 1.916140e-02 6.940284e-02 -6.000000e+00 1.574309e-01 223.67 -2.00 3.49 -2.00 502.75 0.00e+00 20.010 0.146 489.00 513.00 1 3.7197e+00 1.4776e-03 3.5220e-02 -6.0000e+00 1.2188e-02 -4.3521e-02 3.5200e-04 2.2853e-03 -6.0000e+00 2.9888e-03 2.5941e-03 2.4159e-03 1.5118e-04 1.4112e-04 -6.0000e+00 1.2207e-03 3.3408e-03 -3.2317e-05 12 8 8 10 0 0 99.83 18.73 6.13 0.00334 0.02010 46 14 1 4 0 987.900 22.09 15.00 18.73 1.96 22.8 895.19 -1714 -1299 -641 -203 350 764 1001 1059 785 545 87 -334 -446 -528 -644 -484 -99 37 411 466 556 749 826 1021 903 706 397 -185 -642 -1023 -986 -724 -207 293 217 314 67 90 97 177 484 550 712 786 733 592 471 -7 -285 -854 -1341 -1273 -1186 -859 -455 -79 223 386 595 918 1311 1313 1052 720 514 181 -6 -205 -596 -805 -797 -542 -420 -360 -588 -790 -983 -973 -981 -1119 -1026 -1080 -706 -642 -617 -526 -521 -396 -75 207 763 818 738 486 193 43 -458 -741 -1130 -1160 -1091 -724 -525 -458 -519 -753 -748 -337 -202 217 543 762 1060 1301 1473 1556 1564 1178 1042 552 568 599 641 611 594 431 379 413 227 209 -292 -460 -425 -622 -378 -237 -23 44 162 351 309 296 360 67 -162 -504 -802 -1130 -956 -1495 -1144 -575 -178 312 692 918 976 725 502 80 -305 -404 -477 -584 -442 -91 34 390 448 515 697 773 957 850 665 373 -173 -595 -943 -909 -672 -195 279 209 305 65 88 94 173 475 545 713 761 714 578 457 -7 -272 -843 -1313 -1240 -1160 -813 -438 -77 220 385 598 924 1320 1327 1066 729 517 182 -6 -202 -586 -791 -784 -534 -415 -357 -584 -782 -968 -953 -958 -1097 -1010 -1070 -703 -643 -623 -533 -530 -404 -77 212 755 814 735 483 190 44 -457 -732 -1113 -1144 -1089 -732 -520 -461 -528 -773 -771 -347 -209 228 575 785 1102 1364 1555 1650 1662 1254 1105 585 599 631 675 646 631 460 406 442 243 222 -311 -489 -452 -667 -408 -259 -25 47 177 387 344 331 402 74 -180 -559 -886 -1249 -1057