File name: Pandora53s1_Potchefstroom-METSI_20250220_L2Fit_fwvt1c2d20230223p1-8.txt File generation date: 20250220T223319.5Z Data description: Level 2 spectral fitting results file Data file version: fwvt1c2p1-8 Local principal investigator: Stuart Piketh Network principal investigator: Alexander Cede Instrument type: Pandora Instrument number: 53 Spectrometer number: 1 Processing software version used: BlickP v1.8.70 Instrument operation file used: Pandora53_OF_v5d20210210.txt Instrument calibration file used: Pandora53s1_CF_v2d20230223.txt Level 1 file used: Pandora53s1_Potchefstroom-METSI_20250220_L1_smca1c2d20230223p1-8.txt Full location name: Gerhard Minnebron Mooi River Ecosystem Trails for Scientific Investigations Short location name: Potchefstroom-METSI Country of location: South Africa Location latitude [deg]: -26.4907 Location longitude [deg]: 27.1322 Location altitude [m]: 1440 Local noon date: 20250220 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: Modification of session 103 "redoing cal analysis" First and last pixel inside fitting window: 1684 1833 Nominal wavelengths inside fitting window [nm]: 492.0989 492.2196 492.3403 492.4609 492.5815 492.7021 492.8227 492.9432 493.0638 493.1843 493.3048 493.4253 493.5458 493.6663 493.7867 493.9072 494.0276 494.1480 494.2684 494.3888 494.5091 494.6294 494.7498 494.8701 494.9903 495.1106 495.2309 495.3511 495.4713 495.5915 495.7117 495.8319 495.9520 496.0722 496.1923 496.3124 496.4325 496.5525 496.6726 496.7926 496.9127 497.0327 497.1526 497.2726 497.3926 497.5125 497.6324 497.7523 497.8722 497.9921 498.1119 498.2318 498.3516 498.4714 498.5912 498.7109 498.8307 498.9504 499.0702 499.1899 499.3095 499.4292 499.5489 499.6685 499.7881 499.9077 500.0273 500.1469 500.2664 500.3859 500.5055 500.6250 500.7444 500.8639 500.9834 501.1028 501.2222 501.3416 501.4610 501.5803 501.6997 501.8190 501.9383 502.0576 502.1769 502.2962 502.4154 502.5346 502.6538 502.7730 502.8922 503.0114 503.1305 503.2496 503.3687 503.4878 503.6069 503.7260 503.8450 503.9640 504.0830 504.2020 504.3210 504.4399 504.5589 504.6778 504.7967 504.9156 505.0344 505.1533 505.2721 505.3909 505.5097 505.6285 505.7473 505.8660 505.9848 506.1035 506.2222 506.3408 506.4595 506.5781 506.6968 506.8154 506.9339 507.0525 507.1711 507.2896 507.4081 507.5266 507.6451 507.7636 507.8820 508.0005 508.1189 508.2373 508.3557 508.4740 508.5924 508.7107 508.8290 508.9473 509.0656 509.1838 509.3020 509.4203 509.5385 509.6567 509.7748 509.8930 --------------------------------------------------------------------------------------- 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: Azimuth motor temperature [°C], 999=no temperature signal Column 100: Humidity in head sensor [%], -9=no humidity signal Column 101: Pressure in head sensor [hPa], -9=no pressure signal Columns 102-251: Unweighted slant column residuals for each pixel inside the fitting window multiplied by 1e5, 9e5=pixel was not used for fitting Columns 252-401: 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 20250220T153329.3Z 9182.648256 39 1 20.63 2 73.79 265.95 128.79 187.87 0.00 1 0.00 1 0 15 2.399e-03 2.387e-03 1.404e-04 1.402e-04 7.946863e+03 9.150049e+00 5.287771e+02 -6.000000e+00 1.313000e+02 284.60 -2.00 20.00 -2.00 505.27 0.00e+00 1.127643e-03 3.850740e-06 7.771034e-05 -6.000000e+00 6.706113e-05 284.60 -2.00 20.00 -2.00 500.47 0.00e+00 1.180194e-01 2.754757e-03 5.513340e-02 -6.000000e+00 4.728387e-02 227.69 -2.00 3.06 -2.00 502.06 0.00e+00 3.525 0.000 489.00 513.00 1 2.8634e+00 2.1061e-04 6.5647e-03 -6.0000e+00 3.6170e-03 -7.5385e-03 5.4752e-05 1.3100e-03 -6.0000e+00 9.5436e-04 7.0447e-02 7.1519e-04 1.2530e-05 6.7142e-05 -6.0000e+00 2.2685e-04 2.1756e-02 -2.8025e-04 12 1 1 12 1 1 99.70 18.22 1.65 0.02176 0.02672 312 142 1 4 0 65.900 21.17 15.00 18.22 26.47 53.60 19.2 883.92 -65 -71 -94 322 455 392 463 299 357 362 111 49 -20 -126 -99 1 33 -175 12 177 82 190 187 43 161 227 54 -174 42 -65 -204 3 -127 -348 -84 96 75 11 22 -9 -43 75 106 121 62 -38 -87 -105 -285 -263 -410 -438 -376 -566 -571 -290 -43 63 158 217 123 160 157 130 124 115 154 4 -94 -12 -131 25 62 -81 -110 -34 15 -45 -141 -197 -432 -445 -290 -429 -456 -449 -571 -580 -528 -288 22 81 101 220 168 25 -96 -379 -440 -316 -251 1 -98 -310 -191 -162 -84 38 39 25 11 66 144 287 484 575 488 366 311 171 195 350 284 29 172 162 121 277 243 -2 49 92 53 152 249 34 57 64 58 56 134 183 224 179 8 -147 -254 -285 -213 -163 -63 -69 -94 326 466 405 481 310 372 372 114 50 -20 -126 -98 1 33 -178 13 184 86 198 195 45 167 237 56 -181 43 -65 -204 3 -128 -358 -87 101 77 12 23 -10 -44 77 108 125 64 -39 -90 -107 -290 -264 -408 -431 -366 -551 -565 -292 -43 65 162 222 126 165 160 132 127 116 155 4 -93 -11 -129 25 61 -80 -110 -33 15 -44 -136 -190 -418 -430 -282 -422 -449 -446 -562 -570 -523 -285 22 80 102 221 166 25 -92 -359 -412 -294 -233 1 -93 -303 -189 -162 -83 37 38 24 11 66 143 287 486 579 491 372 315 171 194 345 277 29 170 160 120 273 239 -2 48 89 51 146 239 33 55 62 58 56 135 183 225 180 8 -148 -254 -282 -207 -158 SQ 20250220T153748.1Z 9182.651251 42 1 51.13 2 74.76 265.49 128.90 186.67 0.00 1 0.00 1 0 15 2.384e-03 2.365e-03 7.002e-04 6.989e-04 8.798029e+03 4.631916e+01 5.860947e+02 -6.000000e+00 1.304853e+02 284.60 -2.00 20.00 -2.00 505.27 0.00e+00 1.138135e-03 1.918364e-05 7.906730e-05 -6.000000e+00 6.664641e-05 284.60 -2.00 20.00 -2.00 500.42 0.00e+00 1.825003e-01 1.378901e-02 5.897500e-02 -6.000000e+00 4.699400e-02 227.69 -2.00 3.06 -2.00 502.01 0.00e+00 3.733 0.000 489.00 513.00 1 6.9211e+00 1.0544e-03 6.9789e-03 -6.0000e+00 3.5948e-03 1.9760e-03 2.7390e-04 1.4164e-03 -6.0000e+00 9.4859e-04 1.1800e-03 8.2179e-04 6.2065e-05 7.2729e-05 -6.0000e+00 2.2545e-04 2.2428e-02 -2.8968e-04 10 0 0 10 0 0 99.60 18.22 1.62 0.02243 0.02813 26 5 1 4 0 1966.270 21.17 15.01 18.22 26.62 54.60 19.1 884.75 28 -67 -55 410 412 325 346 216 211 340 131 158 -72 24 0 157 70 -326 -112 99 94 140 216 2 32 261 36 -90 42 -59 -59 -2 31 -248 0 -59 72 -18 -42 111 -25 72 71 138 81 -203 -283 -53 -395 -369 -469 -466 -144 -472 -657 -159 -50 10 7 231 122 20 139 140 118 113 182 -14 -70 78 -198 8 128 -339 244 1 244 -40 -82 -81 -411 -521 -187 -357 -443 -457 -586 -570 -422 -252 -34 193 79 110 264 71 -100 -486 -306 -252 -276 148 -142 -431 -271 -196 -57 64 90 -35 8 46 53 266 375 529 466 212 182 131 114 361 413 127 312 242 101 180 164 -77 -133 180 59 322 248 66 56 -29 94 26 177 132 159 212 -89 -128 -349 -326 -67 74 27 -66 -54 419 423 335 363 228 223 352 134 162 -73 24 0 157 71 -335 -118 103 98 146 225 3 34 274 38 -94 43 -59 -59 -2 32 -254 0 -63 75 -19 -43 114 -26 73 73 144 84 -212 -294 -55 -404 -375 -462 -456 -139 -462 -651 -162 -51 10 7 237 126 20 144 145 121 114 185 -14 -69 76 -196 8 125 -338 244 1 239 -40 -80 -78 -398 -500 -180 -348 -438 -448 -576 -559 -412 -250 -33 191 79 110 259 69 -96 -458 -286 -231 -252 137 -134 -421 -269 -193 -57 63 88 -34 8 45 53 266 375 531 472 214 183 132 112 356 402 125 308 235 98 175 159 -75 -130 173 56 307 237 62 55 -28 93 26 176 132 160 214 -90 -128 -346 -321 -64 71 SQ 20250220T154301.6Z 9182.654880 45 1 49.69 2 75.92 264.94 129.01 185.22 0.00 1 0.00 1 0 15 2.051e-02 2.322e-02 9.860e-04 9.838e-04 1.094186e+04 6.681836e+01 6.973066e+02 -6.000000e+00 1.123181e+03 284.60 -2.00 20.00 -2.00 505.27 0.00e+00 1.991304e-03 2.705841e-05 1.274502e-04 -6.000000e+00 5.732778e-04 284.60 -2.00 20.00 -2.00 500.44 0.00e+00 -3.211326e-01 1.943903e-02 7.750824e-02 -6.000000e+00 4.036431e-01 227.69 -2.00 3.06 -2.00 502.03 0.00e+00 4.022 0.000 489.00 513.00 1 7.5435e+00 1.4875e-03 8.6471e-03 -6.0000e+00 3.0891e-02 1.2019e-02 3.8509e-04 1.7880e-03 -6.0000e+00 8.1285e-03 6.2466e-04 -2.2453e-03 8.6597e-05 8.7639e-05 -6.0000e+00 1.9382e-03 9.5632e-03 -5.0512e-05 12 8 8 10 0 0 99.45 18.22 1.69 0.00956 -0.00601 13 2 1 4 0 3821.930 21.17 15.00 18.22 26.69 53.10 19.0 884.75 405 400 400 704 611 877 789 505 445 665 197 178 -179 -77 -69 67 190 -542 -153 -43 -216 -245 -141 -277 -91 92 -396 -351 82 -144 -273 -197 -125 -391 -145 -124 -142 -325 -117 119 -21 93 134 335 252 178 148 219 239 -9 88 -157 435 -25 -54 541 476 737 685 1127 884 1019 -1512 863 943 1014 892 917 763 992 497 778 1081 401 -23782 902 1000 807 699 697 289 309 580 423 216 430 77 61 -249 -18 338 352 458 174 273 -335 -723 -986 -1043 -906 -745 -182 -321 -543 -213 -443 139 320 250 458 464 636 422 601 645 664 810 840 844 790 682 753 672 171 417 413 135 446 429 16 15 475 -119 409 364 413 281 100 173 107 281 173 435 -97 -356 -366 -810 -662 -418 77 394 390 395 718 626 914 821 530 465 684 204 182 -181 -78 -70 68 192 -554 -160 -45 -229 -259 -149 -288 -95 95 -412 -368 83 -144 -270 -196 -127 -406 -151 -130 -148 -341 -121 122 -22 96 137 348 264 187 155 226 242 -9 88 -154 419 -25 -53 548 484 753 700 1145 901 1044 -1565 889 962 1019 900 913 751 981 490 767 1053 398 -27131 894 973 787 674 668 279 301 567 410 213 425 77 61 -247 -18 334 353 455 172 270 -325 -694 -921 -971 -828 -678 -168 -302 -527 -210 -442 137 311 241 443 455 632 424 594 639 660 813 851 855 790 669 737 651 167 408 405 133 438 420 16 15 453 -114 389 346 397 271 98 169 106 282 171 432 -97 -356 -364 -811 -647 -407 74