File name: Pandora145s1_Toronto-Scarborough_20241104_L2Fit_fwvt1c8d20221111p1-8.txt File generation date: 20241105T123344.1Z Data description: Level 2 spectral fitting results file Data file version: fwvt1c8p1-8 Local principal investigator: Vitali Fioletov Network principal investigator: Alexander Cede Instrument type: Pandora Instrument number: 145 Spectrometer number: 1 Processing software version used: BlickP v1.8.60 Instrument operation file used: Pandora145_OF_v5d20190718.txt Instrument calibration file used: Pandora145s1_CF_v8d20221111.txt Level 1 file used: Pandora145s1_Toronto-Scarborough_20241104_L1_smca1c8d20221111p1-8.txt Full location name: University of Toronto - Scarborough Short location name: Toronto-Scarborough Country of location: Canada Location latitude [deg]: 43.7843 Location longitude [deg]: -79.1874 Location altitude [m]: 137 Local noon date: 20241104 Notes on s-number (L1 configuration): Corrections NOT applied although requested by the s-number are dark map for dark correction (replaced by immediate dark method), latency correction, matrix method stray light correction (replaced by simple method) Data caveats: based on session 4 for field cal VP 20221111 First and last pixel inside fitting window: 1647 1796 Nominal wavelengths inside fitting window [nm]: 492.0874 492.2085 492.3296 492.4507 492.5718 492.6928 492.8138 492.9348 493.0558 493.1768 493.2978 493.4187 493.5397 493.6606 493.7815 493.9024 494.0232 494.1441 494.2649 494.3857 494.5065 494.6273 494.7481 494.8688 494.9895 495.1102 495.2309 495.3516 495.4723 495.5929 495.7136 495.8342 495.9548 496.0754 496.1959 496.3165 496.4370 496.5575 496.6780 496.7985 496.9190 497.0394 497.1598 497.2803 497.4007 497.5210 497.6414 497.7618 497.8821 498.0024 498.1227 498.2430 498.3632 498.4835 498.6037 498.7239 498.8441 498.9643 499.0845 499.2046 499.3248 499.4449 499.5650 499.6850 499.8051 499.9252 500.0452 500.1652 500.2852 500.4052 500.5251 500.6451 500.7650 500.8849 501.0048 501.1247 501.2446 501.3644 501.4842 501.6040 501.7238 501.8436 501.9634 502.0831 502.2029 502.3226 502.4423 502.5619 502.6816 502.8012 502.9209 503.0405 503.1601 503.2796 503.3992 503.5187 503.6383 503.7578 503.8773 503.9967 504.1162 504.2356 504.3551 504.4745 504.5939 504.7132 504.8326 504.9519 505.0712 505.1905 505.3098 505.4291 505.5484 505.6676 505.7868 505.9060 506.0252 506.1444 506.2635 506.3827 506.5018 506.6209 506.7400 506.8590 506.9781 507.0971 507.2161 507.3351 507.4541 507.5731 507.6920 507.8110 507.9299 508.0488 508.1677 508.2865 508.4054 508.5242 508.6430 508.7618 508.8806 508.9993 509.1181 509.2368 509.3555 509.4742 509.5929 509.7116 509.8302 509.9488 --------------------------------------------------------------------------------------- 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-250: Unweighted slant column residuals for each pixel inside the fitting window multiplied by 1e5, 9e5=pixel was not used for fitting Columns 251-400: 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 20241104T174445.6Z 9074.739417 270 1 31.19 2 60.33 192.32 76.38 156.56 0.00 1 0.00 1 0 22 3.855e-03 3.828e-03 1.518e-03 1.513e-03 4.655098e+03 9.048861e+01 3.485031e+02 -6.000000e+00 2.009617e+02 278.69 -2.00 20.00 -2.00 505.24 0.00e+00 5.932739e-04 4.040222e-05 3.972424e-05 -6.000000e+00 1.052477e-04 278.69 -2.00 20.00 -2.00 500.32 0.00e+00 1.060014e-01 2.943622e-02 2.968725e-02 -6.000000e+00 7.458125e-02 223.48 -2.00 3.57 -2.00 501.93 0.00e+00 2.011 0.000 489.00 513.00 1 7.9061e+00 2.2452e-03 3.9386e-03 -6.0000e+00 5.6952e-03 8.5145e-03 5.8930e-04 7.7105e-04 -6.0000e+00 1.5025e-03 5.3204e-04 4.2546e-04 1.2419e-04 4.6853e-05 -6.0000e+00 3.2838e-04 -1.5377e-02 -3.5623e-06 11 8 0 10 0 0 28.15 18.21 0.57 -0.01538 -0.01243 8 8 1 4 0 3898.530 22.90 15.00 18.21 19.71 24.7 1081.00 -1255 -453 -455 299 543 605 530 704 557 279 255 144 -151 -367 -514 -239 -252 10 58 338 268 285 126 151 287 389 96 -88 -100 -609 -1011 -264 83 -127 368 370 51 -126 -101 -120 -365 65 300 193 502 357 282 -114 -289 -142 -738 -834 -798 -499 -573 353 262 220 390 252 56 413 457 272 313 181 94 -100 -339 -244 -577 -541 -64 139 296 52 -358 -672 -624 -402 -599 -287 34 -371 10 -82 110 -187 -143 -323 398 85 87 446 617 144 222 13 -267 -307 -991 -719 -472 -117 499 120 4 -275 -342 -264 4 262 98 146 714 286 340 290 226 392 -53 -9 -126 192 240 205 15 187 -42 240 158 34 -503 -246 -253 -31 -29 363 103 299 434 181 424 263 569 -140 -87 -672 -408 -813 -1244 -453 -454 308 572 624 558 747 589 291 262 152 -157 -376 -526 -246 -262 10 61 360 285 303 135 162 309 417 103 -93 -103 -627 -1032 -269 83 -133 396 398 54 -133 -105 -123 -377 67 314 205 536 382 299 -120 -299 -143 -720 -826 -781 -493 -563 356 270 229 403 260 58 429 476 284 323 184 94 -99 -332 -234 -552 -522 -62 137 294 51 -346 -640 -589 -377 -567 -273 33 -359 10 -81 108 -184 -141 -320 394 86 86 438 616 141 212 12 -252 -285 -927 -659 -450 -117 488 118 4 -271 -335 -251 4 256 97 146 709 286 345 294 228 387 -51 -9 -120 184 230 197 14 179 -40 230 149 32 -480 -231 -239 -29 -27 350 101 295 429 180 424 255 550 -135 -85 -639 -380 -762 SQ 20241104T182611.5Z 9074.768188 298 1 31.16 2 62.63 203.37 74.06 165.40 0.00 1 0.00 1 0 22 3.733e-03 3.720e-03 1.145e-03 1.142e-03 3.993550e+03 6.717777e+01 3.032767e+02 -6.000000e+00 1.945858e+02 278.68 -2.00 20.00 -2.00 505.26 0.00e+00 5.895198e-04 3.055152e-05 3.944930e-05 -6.000000e+00 1.019095e-04 278.68 -2.00 20.00 -2.00 500.38 0.00e+00 1.200411e-01 2.211361e-02 2.582446e-02 -6.000000e+00 7.220264e-02 223.48 -2.00 3.57 -2.00 501.98 0.00e+00 2.164 0.000 489.00 513.00 1 7.3825e+00 1.6881e-03 3.8940e-03 -6.0000e+00 5.5138e-03 8.7897e-03 4.4279e-04 6.7782e-04 -6.0000e+00 1.4543e-03 8.7938e-04 8.2780e-05 9.3866e-05 4.1292e-05 -6.0000e+00 3.1794e-04 -1.4591e-02 6.0078e-06 10 0 0 10 0 0 96.15 18.25 0.56 -0.01459 -0.01402 10 10 1 4 0 3115.110 22.90 14.99 18.25 20.77 24.7 1086.52 -917 -375 -288 192 524 720 562 850 605 443 90 39 -397 -133 -615 -299 -71 71 178 254 -5 -36 193 123 328 530 324 -160 -23 -560 -723 -572 -23 -148 383 627 63 -134 -204 -370 -407 29 293 271 376 302 334 254 -268 -396 -600 -933 -628 -825 -375 42 223 107 409 72 381 247 388 424 276 26 -3 -222 -154 -109 -387 -314 -219 -44 -138 155 -340 -659 -332 -319 -361 -176 -144 -527 -289 13 -150 -252 -138 -114 109 -37 472 418 524 206 340 -46 -507 -652 -952 -564 -316 -117 177 131 10 -280 -397 -166 233 166 264 370 430 491 580 375 340 -89 185 117 -84 -149 372 -145 217 114 363 -82 234 35 -239 -331 -248 -306 -35 411 206 318 113 362 259 666 114 195 -134 -460 -431 -685 -906 -366 -292 196 545 744 593 903 642 463 93 41 -409 -135 -624 -304 -73 73 189 272 -6 -39 206 132 345 555 339 -171 -24 -571 -728 -579 -23 -155 403 659 65 -138 -212 -383 -422 30 308 281 394 317 347 261 -278 -400 -602 -915 -607 -809 -369 43 225 109 414 73 385 252 396 433 278 26 -3 -220 -151 -105 -371 -304 -214 -44 -138 154 -331 -631 -321 -308 -342 -168 -139 -514 -286 12 -149 -250 -137 -114 106 -36 474 420 515 199 328 -43 -478 -605 -887 -535 -303 -115 178 131 10 -271 -381 -161 229 166 260 367 427 489 584 379 340 -88 183 114 -82 -146 365 -143 213 111 354 -81 226 33 -223 -307 -230 -290 -34 394 201 313 112 358 258 664 112 193 -130 -435 -409 -653 SQ 20241104T191751.4Z 9074.804067 331 1 29.64 2 67.27 216.11 72.71 176.84 0.00 1 0.00 1 0 15 3.823e-03 3.793e-03 6.226e-04 6.214e-04 4.130114e+03 3.649453e+01 3.090023e+02 -6.000000e+00 1.992956e+02 278.68 -2.00 20.00 -2.00 505.27 0.00e+00 7.177658e-04 1.662953e-05 4.584158e-05 -6.000000e+00 1.043746e-04 278.68 -2.00 20.00 -2.00 500.45 0.00e+00 1.104928e-01 1.200178e-02 2.781340e-02 -6.000000e+00 7.397620e-02 223.48 -2.00 3.57 -2.00 502.04 0.00e+00 2.568 0.000 489.00 513.00 1 6.1098e+00 9.1630e-04 4.4703e-03 -6.0000e+00 5.6489e-03 -2.7850e-03 2.3982e-04 7.2125e-04 -6.0000e+00 1.4905e-03 2.9629e-03 4.5906e-04 5.1086e-05 4.3271e-05 -6.0000e+00 3.2569e-04 -1.2563e-02 -3.3286e-05 11 8 0 10 0 0 99.91 18.27 0.58 -0.01256 -0.00938 33 33 1 4 0 898.060 22.90 15.00 18.27 21.87 24.6 1091.91 -847 -429 -341 17 625 581 690 742 563 502 97 -220 -145 -437 -542 -294 -109 -86 165 363 223 27 99 115 309 466 361 -18 25 -483 -885 -528 -166 -24 375 428 168 -63 -173 -277 -183 28 217 326 322 500 330 134 -109 -305 -734 -993 -802 -676 -250 -19 162 156 138 124 194 303 440 428 295 159 56 -96 -299 -235 -572 -421 23 17 117 138 -193 -488 -571 -449 -564 -241 -146 -293 -180 -94 -140 -217 -291 -162 68 44 362 648 589 212 185 -164 -486 -777 -754 -541 -375 -99 311 264 -263 -257 -437 -309 -32 315 258 444 499 560 474 443 406 221 -4 26 -20 -89 208 165 139 284 307 22 246 22 -491 -256 -251 -267 -19 275 241 322 273 408 331 465 274 231 -180 -495 -793 -761 -825 -425 -340 17 638 605 712 773 585 528 101 -227 -147 -436 -540 -294 -111 -90 172 381 234 28 104 121 328 493 380 -18 25 -486 -882 -528 -169 -25 398 443 177 -66 -178 -283 -187 29 225 343 333 517 339 140 -112 -305 -733 -969 -772 -659 -251 -20 167 157 142 128 199 306 445 432 304 162 56 -95 -292 -232 -564 -406 22 17 117 136 -186 -479 -555 -434 -550 -230 -140 -285 -178 -94 -140 -215 -289 -161 68 44 364 648 579 211 179 -155 -460 -726 -704 -516 -361 -97 311 263 -260 -249 -420 -300 -31 316 254 441 496 559 478 448 407 218 -4 25 -20 -87 206 163 137 280 302 21 240 21 -464 -239 -235 -254 -19 266 236 319 271 407 332 468 277 231 -176 -472 -765 -735