File name: Pandora162s1_Brussels-Uccle_20241030_L2Fit_fwvt1c13d20231120p1-8.txt File generation date: 20241031T000441.7Z Data description: Level 2 spectral fitting results file Data file version: fwvt1c13p1-8 Local principal investigator: Michel Van Roozendael Network principal investigator: Alexander Cede Instrument type: Pandora Instrument number: 162 Spectrometer number: 1 Processing software version used: BlickP v1.8.60 Instrument operation file used: Pandora162_OF_v8d20191031.txt Instrument calibration file used: Pandora162s1_CF_v13d20231120.txt Level 1 file used: Pandora162s1_Brussels-Uccle_20241030_L1_smca1c13d20231120p1-8.txt Full location name: Uccle (BIRA-IASB, KMI-IRM) Short location name: Brussels-Uccle Country of location: Belgium Location latitude [deg]: 50.7980 Location longitude [deg]: 4.3580 Location altitude [m]: 107 Local noon date: 20241030 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 12 for fieldcal First and last pixel inside fitting window: 1648 1798 Nominal wavelengths inside fitting window [nm]: 492.0019 492.1226 492.2433 492.3640 492.4847 492.6053 492.7260 492.8466 492.9672 493.0878 493.2084 493.3289 493.4495 493.5700 493.6905 493.8110 493.9314 494.0519 494.1723 494.2927 494.4131 494.5335 494.6539 494.7742 494.8946 495.0149 495.1352 495.2555 495.3758 495.4960 495.6162 495.7365 495.8567 495.9768 496.0970 496.2172 496.3373 496.4574 496.5775 496.6976 496.8177 496.9377 497.0577 497.1778 497.2978 497.4177 497.5377 497.6577 497.7776 497.8975 498.0174 498.1373 498.2572 498.3770 498.4968 498.6166 498.7364 498.8562 498.9760 499.0957 499.2155 499.3352 499.4549 499.5746 499.6942 499.8139 499.9335 500.0531 500.1727 500.2923 500.4118 500.5314 500.6509 500.7704 500.8899 501.0094 501.1289 501.2483 501.3677 501.4871 501.6065 501.7259 501.8453 501.9646 502.0839 502.2032 502.3225 502.4418 502.5611 502.6803 502.7995 502.9187 503.0379 503.1571 503.2762 503.3954 503.5145 503.6336 503.7527 503.8718 503.9908 504.1099 504.2289 504.3479 504.4669 504.5858 504.7048 504.8237 504.9426 505.0615 505.1804 505.2993 505.4181 505.5370 505.6558 505.7746 505.8934 506.0121 506.1309 506.2496 506.3683 506.4870 506.6057 506.7244 506.8430 506.9616 507.0803 507.1988 507.3174 507.4360 507.5545 507.6731 507.7916 507.9101 508.0285 508.1470 508.2654 508.3839 508.5023 508.6207 508.7390 508.8574 508.9757 509.0941 509.2124 509.3306 509.4489 509.5672 509.6854 509.8036 509.9218 --------------------------------------------------------------------------------------- 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: Spectrometer control temperature [°C], 999=no temperature signal Column 96: Auxiliary spectrometer temperature [°C], 999=no temperature signal Column 97: Temperature in head sensor [°C], 999=no temperature signal Column 98: Temperature at electronics board 1 [°C], 999=no temperature signal Column 99: Temperature at electronics board 2 [°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-252: Unweighted slant column residuals for each pixel inside the fitting window multiplied by 1e5, 9e5=pixel was not used for fitting Columns 253-403: 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 20241030T123347.4Z 9069.523465 303 1 49.74 2 66.50 197.90 67.20 221.76 0.00 1 0.00 1 0 15 1.933e-03 1.929e-03 6.928e-04 6.905e-04 2.691866e+03 4.189335e+01 2.064441e+02 -6.000000e+00 1.064155e+02 274.35 -2.00 20.00 -2.00 505.25 0.00e+00 7.842713e-04 1.880124e-05 4.547832e-05 -6.000000e+00 5.358407e-05 274.35 -2.00 20.00 -2.00 500.33 0.00e+00 2.634470e-02 1.392482e-02 2.026170e-02 -6.000000e+00 3.876510e-02 222.76 -2.00 3.79 -2.00 501.95 0.00e+00 2.490 0.000 489.00 513.00 1 6.8065e+00 1.0688e-03 4.0569e-03 -6.0000e+00 2.9683e-03 1.8574e-02 2.7330e-04 5.2802e-04 -6.0000e+00 7.6768e-04 1.5051e-03 1.2525e-03 7.3695e-05 3.0708e-05 -6.0000e+00 2.1345e-04 5.6488e-04 -5.4975e-05 10 0 0 10 0 0 99.57 18.78 2.01 0.00056 0.00925 38 8 1 2 0 1308.370 15.00 18.78 23.50 21.17 21.79 1.3 975.96 900000 -223 17 -8 -98 78 305 251 57 187 176 281 82 -117 107 -178 -215 -150 236 -58 -146 86 60 216 401 151 -31 217 29 -397 -175 -196 -157 -78 332 -166 -308 93 103 -172 -17 49 -114 -94 195 277 88 232 155 121 -199 -170 -224 -318 -56 -111 -126 -271 -107 -66 -30 137 247 31 62 -55 -13 -95 170 -322 -372 -129 -149 -256 285 -48 -248 -331 62 52 -78 -7 -24 -331 -3 78 -78 -73 -150 20 186 108 227 209 3 121 111 43 -209 -63 -302 -165 -338 -235 -18 -380 -279 -250 -23 -14 148 -23 -127 -130 -22 180 238 328 184 121 365 175 -29 302 271 71 13 369 44 302 281 -182 -228 158 -152 -72 422 -180 -49 -110 32 88 61 274 113 120 -14 -172 -397 -301 49 900000 -222 17 -8 -100 82 325 271 62 201 187 293 84 -119 108 -180 -219 -156 250 -62 -158 94 65 235 437 165 -34 234 30 -411 -177 -197 -160 -82 354 -172 -320 95 110 -182 -18 51 -116 -96 201 286 91 237 156 126 -202 -167 -217 -308 -55 -112 -130 -271 -107 -66 -30 138 248 32 62 -55 -12 -97 170 -320 -366 -126 -147 -255 285 -48 -247 -325 60 50 -75 -7 -24 -325 -3 74 -74 -70 -144 19 179 105 223 207 3 118 106 40 -202 -60 -282 -155 -324 -221 -17 -373 -274 -243 -22 -14 142 -22 -126 -131 -22 182 241 333 187 122 361 171 -28 291 261 69 12 357 43 291 269 -173 -215 146 -140 -67 393 -172 -48 -108 32 87 61 275 114 120 -14 -168 -380 -284 46 SQ 20241030T130310.7Z 9069.543874 320 1 7.61 2 68.22 205.40 70.53 228.54 0.00 1 0.00 1 0 15 2.402e-03 2.397e-03 1.508e-03 1.505e-03 2.719622e+03 9.046543e+01 2.083695e+02 -6.000000e+00 1.322664e+02 274.34 -2.00 20.00 -2.00 505.24 0.00e+00 7.670044e-04 4.087955e-05 4.458855e-05 -6.000000e+00 6.659906e-05 274.34 -2.00 20.00 -2.00 500.42 0.00e+00 1.141589e-01 3.015047e-02 2.033058e-02 -6.000000e+00 4.816507e-02 222.76 -2.00 3.79 -2.00 502.01 0.00e+00 2.672 0.000 489.00 513.00 1 6.5914e+00 2.3090e-03 4.2798e-03 -6.0000e+00 3.6883e-03 1.4601e-02 5.9386e-04 5.3550e-04 -6.0000e+00 9.5342e-04 1.8076e-03 7.5562e-04 1.6088e-04 3.2595e-05 -6.0000e+00 2.6529e-04 -8.4066e-04 -1.7975e-05 12 1 1 12 1 1 -1054.86 18.79 2.07 -0.00084 0.00440 4 7 1 2 0 1518.780 15.00 18.79 24.93 21.17 21.79 1.2 982.57 900000 -303 -41 190 -94 86 313 273 512 157 194 525 -125 -235 -63 -378 -402 -46 441 176 -142 -26 150 81 281 158 75 68 94 -272 -107 -123 -16 -178 378 -663 -30 94 186 83 -81 -245 -28 -205 207 260 43 213 -120 279 -236 -72 -311 -324 -204 -14 -206 -179 31 -227 215 266 268 -243 243 -19 -150 -243 82 -347 -271 -7 91 150 112 -40 -292 -216 -171 -184 -305 -464 -45 -294 103 -165 48 159 -11 12 31 -34 -84 463 21 0 198 163 -101 -230 -385 -75 -516 -374 73 -462 -114 -165 -123 254 227 8 -191 -138 3 153 176 495 -162 -11 125 307 62 322 452 115 154 510 54 423 431 -236 -377 148 -101 -183 284 175 -236 115 -42 235 324 136 163 108 -173 -272 -369 -94 -278 900000 -296 -41 190 -97 90 324 286 536 164 200 547 -130 -240 -64 -383 -412 -48 453 184 -150 -27 158 86 298 168 79 72 97 -274 -109 -124 -16 -181 393 -707 -31 98 193 85 -84 -253 -29 -215 219 268 44 224 -124 282 -233 -71 -304 -316 -201 -14 -208 -183 32 -235 221 274 277 -245 251 -19 -151 -242 80 -337 -268 -7 90 146 109 -40 -285 -213 -167 -178 -296 -452 -44 -282 99 -161 47 156 -10 12 31 -34 -83 457 20 0 195 157 -98 -219 -362 -71 -484 -362 72 -454 -112 -160 -122 249 224 8 -189 -138 3 155 179 502 -165 -11 124 300 60 310 435 112 150 493 52 407 413 -226 -366 141 -95 -174 273 167 -229 113 -42 235 324 133 160 106 -172 -266 -354 -89 -262 SQ 20241030T152350.7Z 9069.641559 402 1 45.22 2 82.91 237.29 90.13 256.93 0.00 1 0.00 1 0 22 2.544e-03 2.518e-03 5.086e-04 5.078e-04 6.986908e+03 3.326863e+01 5.066430e+02 -6.000000e+00 1.400681e+02 274.33 -2.00 20.00 -2.00 505.31 0.00e+00 1.887159e-03 1.387502e-05 1.074865e-04 -6.000000e+00 7.052982e-05 274.33 -2.00 20.00 -2.00 500.56 0.00e+00 5.000979e-01 1.015491e-02 5.131777e-02 -6.000000e+00 5.101987e-02 222.75 -2.00 3.79 -2.00 502.14 0.00e+00 7.453 0.009 489.00 513.00 1 5.4679e+00 7.7745e-04 1.1727e-02 -6.0000e+00 3.9068e-03 1.0232e-03 1.9739e-04 1.3755e-03 -6.0000e+00 1.0102e-03 2.6980e-03 1.4496e-03 5.3323e-05 7.4765e-05 -6.0000e+00 2.8095e-04 -1.0148e-03 -3.4595e-05 10 0 0 10 0 0 96.92 18.80 1.62 -0.00101 0.00904 42 8 1 2 0 1076.260 15.00 18.80 21.60 21.17 21.79 1.3 959.39 900000 -224 -141 147 -85 117 589 485 478 252 287 256 13 -228 -117 -222 -311 -52 165 48 -99 231 62 164 425 271 94 280 1 -479 -321 -236 -318 -107 257 -211 -190 75 92 -161 -38 -10 -55 48 112 270 107 279 79 147 -205 -387 -395 -526 -125 -86 5 -203 65 -57 25 254 326 32 134 99 -11 -163 -75 -182 -400 -258 -272 -79 156 21 -327 -362 -263 -281 -432 -275 -235 -496 -92 -23 -67 -90 -26 6 37 85 384 437 259 161 148 98 -99 -398 -637 -543 -550 -234 -136 -374 -304 -220 -313 66 220 109 -40 134 140 333 411 552 331 139 132 106 170 242 428 195 197 221 39 230 373 -97 -107 191 -208 -122 292 -100 -72 -8 107 172 192 162 216 173 -32 -89 -441 -369 -275 900000 -209 -135 145 -87 117 594 496 489 256 289 254 13 -231 -117 -224 -317 -52 167 50 -103 240 65 171 445 285 99 291 1 -480 -314 -229 -313 -108 266 -223 -201 78 95 -166 -39 -10 -57 50 114 277 109 283 82 150 -205 -379 -381 -508 -122 -87 5 -205 65 -58 25 258 331 33 136 100 -11 -167 -75 -182 -396 -254 -270 -79 158 22 -329 -358 -256 -272 -418 -267 -229 -490 -92 -22 -65 -87 -25 6 36 83 379 434 256 156 140 94 -93 -369 -586 -502 -521 -219 -132 -371 -303 -217 -306 65 216 108 -41 137 144 344 409 552 332 139 131 108 170 241 427 196 198 223 39 231 373 -97 -105 184 -201 -118 287 -100 -71 -8 108 174 196 165 221 177 -33 -89 -433 -357 -267