File name: Pandora289s1_DeBilt_20251206_L2Fit_fwvt1c2d20250401p1-8.txt File generation date: 20251206T233513.5Z Data description: Level 2 spectral fitting results file Data file version: fwvt1c2p1-8 Local principal investigator: Ankie Piters Network principal investigator: Alexander Cede Instrument type: Pandora Instrument number: 289 Spectrometer number: 1 Processing software version used: BlickP v1.8.75 Instrument operation file used: Pandora289_OF_v3d20240910.txt Instrument calibration file used: Pandora289s1_CF_v2d20250401.txt Level 1 file used: Pandora289s1_DeBilt_20251206_L1_smca1c2d20250401p1-8.txt Full location name: Royal Netherlands Meteorological Institute (KNMI), De Bilt Short location name: DeBilt Country of location: Netherlands Location latitude [deg]: 52.1012 Location longitude [deg]: 5.1783 Location altitude [m]: 20 Local noon date: 20251206 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: OOB session based on clearance session 1 First and last pixel inside fitting window: 1626 1776 Nominal wavelengths inside fitting window [nm]: 492.1097 492.2301 492.3504 492.4707 492.5910 492.7113 492.8315 492.9518 493.0720 493.1922 493.3124 493.4326 493.5527 493.6729 493.7930 493.9131 494.0332 494.1533 494.2733 494.3934 494.5134 494.6334 494.7534 494.8733 494.9933 495.1132 495.2332 495.3531 495.4730 495.5928 495.7127 495.8325 495.9523 496.0721 496.1919 496.3117 496.4315 496.5512 496.6709 496.7906 496.9103 497.0300 497.1496 497.2693 497.3889 497.5085 497.6281 497.7476 497.8672 497.9867 498.1062 498.2257 498.3452 498.4647 498.5841 498.7036 498.8230 498.9424 499.0617 499.1811 499.3005 499.4198 499.5391 499.6584 499.7777 499.8969 500.0162 500.1354 500.2546 500.3738 500.4930 500.6121 500.7313 500.8504 500.9695 501.0886 501.2077 501.3267 501.4458 501.5648 501.6838 501.8028 501.9218 502.0407 502.1597 502.2786 502.3975 502.5164 502.6353 502.7541 502.8729 502.9918 503.1106 503.2293 503.3481 503.4669 503.5856 503.7043 503.8230 503.9417 504.0604 504.1790 504.2976 504.4162 504.5348 504.6534 504.7720 504.8905 505.0090 505.1276 505.2461 505.3645 505.4830 505.6014 505.7198 505.8383 505.9566 506.0750 506.1934 506.3117 506.4300 506.5483 506.6666 506.7849 506.9031 507.0214 507.1396 507.2578 507.3760 507.4941 507.6123 507.7304 507.8485 507.9666 508.0847 508.2028 508.3208 508.4388 508.5568 508.6748 508.7928 508.9108 509.0287 509.1466 509.2645 509.3824 509.5003 509.6182 509.7360 509.8538 509.9716 --------------------------------------------------------------------------------------- Column 1: Two letter code of measurement routine Column 2: UT date and time for center-time of measurement, yyyymmddThhmmssZ (ISO 8601) Column 3: Fractional days since 1-Jan-2000 midnight for center-time of measurement Column 4: Routine count (1 for the first routine of the day, 2 for the second, etc.) Column 5: Repetition count (1 for the first set in the routine, 2 for the second, etc.) Column 6: Total duration of measurement set in seconds Column 7: Data processing type index Column 8: Solar zenith angle for center-time of measurement in degree Column 9: Solar azimuth for center-time of measurement in degree, 0=north, increases clockwise Column 10: Lunar zenith angle for center-time of measurement in degree Column 11: Lunar azimuth for center-time of measurement in degree, 0=north, increases clockwise Column 12: Pointing zenith angle in degree, absolute or relative (see next column), 999=tracker not used Column 13: Zenith pointing mode: zenith angle is... 0=absolute, 1=relative to sun, 2=relative to moon Column 14: Pointing azimuth in degree, increases clockwise, absolute (0=north) or relative (see next column), 999=tracker not used Column 15: Azimuth pointing mode: like zenith angle mode but also fixed scattering angles relative to sun (3) or moon (4) Column 16: Fitting result index: 0=no error or warning, 1, 2=warning, >2=error Column 17: Number of function evaluations used, 0=linear fitting or fitting not successful or no fitting done Column 18: rms of unweighted fitting residuals, -9=fitting not successful Column 19: Normalized rms of fitting residuals weighted with independent uncertainty, -9=fitting not successful or no uncertainty used Column 20: Expected rms based on independent uncertainty, -9=fitting not successful or no uncertainty given Column 21: Expected normalized weighted rms based on independent uncertainty, -9=fitting not successful or no uncertainty given Column 22: Water vapor slant column amount [moles per square meter], -9e99=fitting not successful Column 23: Independent uncertainty of water vapor slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -3=spectral fitting was done, but no independent uncertainty could be retrieved, -5=no independent uncertainty input was given, -9=spectral fitting not successful Column 24: Structured uncertainty of water vapor slant column amount [moles per square meter], -7=not given since method "MEAS" was chosen, -9=spectral fitting not successful Column 25: Common uncertainty of water vapor slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -6=no common uncertainty input was given, -7=not given since method "MEAS" was chosen, -9=spectral fitting not successful Column 26: rms-based uncertainty of water vapor slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -3=spectral fitting was done, but no rms-based uncertainty could be retrieved, -9=spectral fitting not successful Column 27: Water vapor effective temperature [K] Column 28: Independent uncertainty of water vapor effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -2=no temperature fitting was requested and output for effective temperature and structured uncertainty of it is based on f-code, -3=spectral fitting was done, but no independent uncertainty could be retrieved, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -5=no independent uncertainty input was given, -9=spectral fitting not successful Column 29: Structured uncertainty of water vapor effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -9=spectral fitting not successful Column 30: Common uncertainty of water vapor effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -2=no temperature fitting was requested and output for effective temperature and structured uncertainty of it is based on f-code, -3=spectral fitting was done, but no common uncertainty could be retrieved, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -6=no common uncertainty input was given, -9=spectral fitting not successful Column 31: Effective water vapor fitting wavelength [nm], -9=fitting not successful Column 32: Diffuse correction applied before fitting at effective fitting wavelength for water vapor [%], 0=no diffuse correction applied or fitting not requested, >0=measured diffuse correction, <0=(negative value of) calculated diffuse correction Column 33: Nitrogen dioxide slant column amount [moles per square meter], -9e99=fitting not successful Column 34: Independent uncertainty of nitrogen dioxide slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -3=spectral fitting was done, but no independent uncertainty could be retrieved, -5=no independent uncertainty input was given, -9=spectral fitting not successful Column 35: Structured uncertainty of nitrogen dioxide slant column amount [moles per square meter], -7=not given since method "MEAS" was chosen, -9=spectral fitting not successful Column 36: Common uncertainty of nitrogen dioxide slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -6=no common uncertainty input was given, -7=not given since method "MEAS" was chosen, -9=spectral fitting not successful Column 37: rms-based uncertainty of nitrogen dioxide slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -3=spectral fitting was done, but no rms-based uncertainty could be retrieved, -9=spectral fitting not successful Column 38: Nitrogen dioxide effective temperature [K] Column 39: Independent uncertainty of nitrogen dioxide effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -2=no temperature fitting was requested and output for effective temperature and structured uncertainty of it is based on f-code, -3=spectral fitting was done, but no independent uncertainty could be retrieved, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -5=no independent uncertainty input was given, -9=spectral fitting not successful Column 40: Structured uncertainty of nitrogen dioxide effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -9=spectral fitting not successful Column 41: Common uncertainty of nitrogen dioxide effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -2=no temperature fitting was requested and output for effective temperature and structured uncertainty of it is based on f-code, -3=spectral fitting was done, but no common uncertainty could be retrieved, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -6=no common uncertainty input was given, -9=spectral fitting not successful Column 42: Effective nitrogen dioxide fitting wavelength [nm], -9=fitting not successful Column 43: Diffuse correction applied before fitting at effective fitting wavelength for nitrogen dioxide [%], 0=no diffuse correction applied or fitting not requested, >0=measured diffuse correction, <0=(negative value of) calculated diffuse correction Column 44: Ozone slant column amount [moles per square meter], -9e99=fitting not successful Column 45: Independent uncertainty of ozone slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -3=spectral fitting was done, but no independent uncertainty could be retrieved, -5=no independent uncertainty input was given, -9=spectral fitting not successful Column 46: Structured uncertainty of ozone slant column amount [moles per square meter], -7=not given since method "MEAS" was chosen, -9=spectral fitting not successful Column 47: Common uncertainty of ozone slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -6=no common uncertainty input was given, -7=not given since method "MEAS" was chosen, -9=spectral fitting not successful Column 48: rms-based uncertainty of ozone slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -3=spectral fitting was done, but no rms-based uncertainty could be retrieved, -9=spectral fitting not successful Column 49: Ozone effective temperature [K] Column 50: Independent uncertainty of ozone effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -2=no temperature fitting was requested and output for effective temperature and structured uncertainty of it is based on f-code, -3=spectral fitting was done, but no independent uncertainty could be retrieved, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -5=no independent uncertainty input was given, -9=spectral fitting not successful Column 51: Structured uncertainty of ozone effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -9=spectral fitting not successful Column 52: Common uncertainty of ozone effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -2=no temperature fitting was requested and output for effective temperature and structured uncertainty of it is based on f-code, -3=spectral fitting was done, but no common uncertainty could be retrieved, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -6=no common uncertainty input was given, -9=spectral fitting not successful Column 53: Effective ozone fitting wavelength [nm], -9=fitting not successful Column 54: Diffuse correction applied before fitting at effective fitting wavelength for ozone [%], 0=no diffuse correction applied or fitting not requested, >0=measured diffuse correction, <0=(negative value of) calculated diffuse correction Column 55: Molecular scattering air mass factor used for molecular scattering subtraction before the fitting Column 56: Estimated uncertainty of molecular scattering air mass factor, -9=molecular scattering was not subtracted before the fitting Column 57: Lower limit used for wavelength scaling [nm] Column 58: Upper limit used for wavelength scaling [nm] Column 59: Order of smoothing polynomial Column 60: Smoothing polynomial coefficient, order 0 Column 61: Independent uncertainty of smoothing polynomial coefficient, order 0, -9=fitting not successful Column 62: Structured uncertainty of smoothing polynomial coefficient, order 0, -9=fitting not successful Column 63: Common uncertainty of smoothing polynomial coefficient, order 0, -9=fitting not successful Column 64: rms-based uncertainty of smoothing polynomial coefficient, order 0, -9=fitting not successful Column 65: Smoothing polynomial coefficient, order 1 Column 66: Independent uncertainty of smoothing polynomial coefficient, order 1, -9=fitting not successful Column 67: Structured uncertainty of smoothing polynomial coefficient, order 1, -9=fitting not successful Column 68: Common uncertainty of smoothing polynomial coefficient, order 1, -9=fitting not successful Column 69: rms-based uncertainty of smoothing polynomial coefficient, order 1, -9=fitting not successful Column 70: Mean value of measured data inside fitting window [same units as measurements] Column 71: Wavelength change polynomial coefficient, order 0 Column 72: Independent uncertainty of wavelength change polynomial coefficient, order 0, -9=fitting not successful Column 73: Structured uncertainty of wavelength change polynomial coefficient, order 0, -9=fitting not successful Column 74: Common uncertainty of wavelength change polynomial coefficient, order 0, -9=fitting not successful Column 75: rms-based uncertainty of wavelength change polynomial coefficient, order 0, -9=fitting not successful Column 76: L1 based wavelength change coefficient, order 0, -9=no L1 wavelength change determination Column 77: L1 based wavelength change coefficient, order 1, -9=no L1 wavelength change determination Column 78: L2Fit data quality flag: 0=assured high quality, 1=assured medium quality, 2=assured low quality, 10=not-assured high quality, 11=not-assured medium quality, 12=not-assured low quality Column 79: Sum over 2^i using those i, for which the corresponding L2Fit data quality parameter exceeds the DQ1 limit, 0=L1 data quality above 0, 1=Spectral fitting was not successful, 2=Wavelength shift too large, 3=Normalized rms of fitting residuals weighted with independent uncertainty too large Column 80: Sum over 2^i using those i, for which the corresponding L2Fit data quality parameter exceeds the DQ2 limit (same parameters as for DQ1) Column 81: L1 data quality flag: 0=assured high quality, 1=assured medium quality, 2=assured low quality, 10=not-assured high quality, 11=not-assured medium quality, 12=not-assured low quality Column 82: Sum over 2^i using those i, for which the corresponding L1 data quality parameter exceeds the DQ1 limit, 0=Saturated data, 1=Too few dark counts measurements, 2=No temperature given or effective temperature too different from the reference temperature, 3=Dark count too high, 4=Unsuccessful dark background fitting, 5=The dark count differs significantly from the dark map for too many pixels, 6=Absolute value of estimated average residual stray light level too high, 7=Although attempted, no wavelength change could be retrieved, 8=Absolute value of retrieved wavelength shift too large, 9=Retrieved wavelength shift differs too much from the shift predicted by the effective temperature Column 83: Sum over 2^i using those i, for which the corresponding L1 data quality parameter exceeds the DQ2 limit (same parameters as for DQ1) Column 84: Atmospheric variability [%], 999=no atmospheric variability was determined Column 85: Wavelength effective temperature [°C], 999=no effective temperature given Column 86: Estimated average residual stray light level [%] (only valid for stray light correction methods 2 and higher) Column 87: Retrieved wavelength shift from L1 data [nm], -9=no wavelength change determination Column 88: Retrieved total wavelength shift [nm], -9=no wavelength change fitting Column 89: Number of bright count cycles Column 90: Number of dark count cycles Column 91: Effective position of filterwheel #1, 0=filterwheel not used, 1-9 are valid positions Column 92: Effective position of filterwheel #2, 0=filterwheel not used, 1-9 are valid positions Column 93: Sum over 2^i, 0=spectra were interpolated in time, 1=spectra are corrected for off-target signal Column 94: Integration time [ms] Column 95: Temperature at electronics board [°C], 999=no temperature signal Column 96: Spectrometer control temperature [°C], 999=no temperature signal Column 97: Auxiliary spectrometer temperature [°C], 999=no temperature signal Column 98: Temperature in head sensor [°C], 999=no temperature signal Column 99: Humidity in head sensor [%], -9=no humidity signal Column 100: Pressure in head sensor [hPa], -9=no pressure signal Columns 101-251: Unweighted slant column residuals for each pixel inside the fitting window multiplied by 1e5, 9e5=pixel was not used for fitting Columns 252-402: Normalized slant column residuals weighted with independent instrumental uncertainty for each pixel inside the fitting window multiplied by 1e5, 9e5=pixel was not used for fitting, 0=no independent instrumental uncertainty was given --------------------------------------------------------------------------------------- SO 20251206T124402.2Z 9471.530581 2279 1 2.94 2 76.37 197.45 100.85 354.75 0.00 1 0.00 1 0 15 3.055e-03 3.026e-03 1.432e-03 1.430e-03 3.837502e+03 7.887893e+01 3.154173e+02 -6.000000e+00 1.503352e+02 269.46 -2.00 20.00 -2.00 505.27 0.00e+00 6.840533e-04 3.708299e-05 3.717923e-05 -6.000000e+00 8.083873e-05 269.46 -2.00 20.00 -2.00 500.49 0.00e+00 3.672719e-01 2.719756e-02 2.217641e-02 -6.000000e+00 5.856074e-02 220.42 -2.00 3.84 -2.00 502.08 0.00e+00 4.146 0.001 489.00 513.00 1 5.3990e+00 2.0694e-03 6.2895e-03 -6.0000e+00 4.4562e-03 6.3170e-03 5.4347e-04 7.0508e-04 -6.0000e+00 1.1846e-03 4.7229e-03 6.5893e-03 1.0929e-04 4.8090e-05 -6.0000e+00 2.4166e-04 3.6460e-02 -4.9301e-04 12 1 1 12 1 1 -1612.54 24.53 1.35 0.03646 0.08217 3 5 1 4 0 732.610 33.75 19.99 24.53 14.98 0.0 947.90 -977 -171 117 178 418 164 607 416 42 346 -37 -208 -82 -286 -363 -260 171 275 151 339 360 101 82 261 345 393 88 64 -372 -341 -652 -125 84 251 185 18 45 -77 -358 -118 -288 -227 58 393 248 155 280 -8 -29 -364 -421 -617 -438 -520 -598 -234 109 155 -27 299 120 327 243 430 91 -52 170 195 -309 -424 -321 -101 65 -122 -107 110 122 -381 -219 -451 30 -142 -272 -198 20 -48 4 -175 126 -286 -128 191 435 407 374 239 223 161 -337 -371 -924 -768 -274 21 -43 65 -210 -326 -115 -473 -185 130 -11 378 129 455 541 503 204 112 170 89 -122 139 68 66 234 453 212 -21 158 118 15 -372 -245 -249 -491 3 172 98 115 133 248 215 483 282 189 4 -412 -330 -564 -949 -168 118 182 428 170 636 433 43 355 -38 -211 -84 -288 -362 -261 173 283 157 355 375 105 86 275 359 408 92 66 -380 -339 -648 -125 84 260 194 18 47 -80 -366 -120 -293 -233 60 406 259 162 290 -8 -29 -367 -417 -603 -421 -508 -591 -236 110 157 -27 304 123 332 248 438 92 -52 171 193 -305 -418 -313 -99 64 -122 -107 110 118 -369 -210 -434 29 -137 -266 -193 20 -47 4 -173 123 -281 -127 189 437 412 373 236 217 156 -319 -348 -861 -725 -264 21 -43 65 -209 -323 -112 -465 -182 130 -11 380 130 460 544 511 207 112 168 88 -120 136 66 64 230 448 210 -21 154 115 14 -354 -235 -239 -478 3 170 97 115 133 249 215 484 283 190 4 -402 -317 -548 SO 20251206T124434.1Z 9471.530951 2281 1 36.48 2 76.40 197.58 100.85 354.86 0.00 1 0.00 1 0 15 3.148e-03 3.114e-03 5.895e-04 5.886e-04 3.735758e+03 3.252827e+01 2.911788e+02 -6.000000e+00 1.549449e+02 269.46 -2.00 20.00 -2.00 505.29 0.00e+00 7.609659e-04 1.526283e-05 4.497581e-05 -6.000000e+00 8.326814e-05 269.46 -2.00 20.00 -2.00 500.54 0.00e+00 3.834569e-01 1.118628e-02 2.697325e-02 -6.000000e+00 6.009047e-02 220.42 -2.00 3.84 -2.00 502.12 0.00e+00 4.153 0.002 489.00 513.00 1 6.3931e+00 8.5113e-04 6.4398e-03 -6.0000e+00 4.5759e-03 5.4573e-03 2.2336e-04 7.5721e-04 -6.0000e+00 1.2101e-03 1.7426e-03 8.4602e-04 4.5023e-05 4.3794e-05 -6.0000e+00 2.4820e-04 -7.2259e-04 2.6684e-05 10 0 0 10 0 0 99.99 24.53 1.42 -0.00072 0.00515 96 10 1 4 0 379.870 33.75 19.99 24.53 15.00 0.0 947.90 -730 -404 -200 60 304 335 550 591 282 271 173 -104 -184 -421 -261 -210 -21 98 221 331 274 154 202 270 409 429 535 62 -243 -461 -524 -335 -10 79 278 218 -42 -118 -211 -222 -307 -139 183 141 257 162 284 84 9 -236 -624 -856 -711 -613 -267 -217 69 171 101 143 189 304 346 215 287 244 122 -32 -349 -471 -357 -310 66 131 50 5 -154 -287 -290 -234 -245 -252 -227 -140 -128 -89 -124 -205 -254 -255 -87 211 343 250 394 361 241 21 -556 -623 -859 -558 -372 -83 41 107 -127 -274 -291 -357 -144 87 99 150 417 527 475 483 295 340 16 101 201 111 94 307 279 274 416 352 107 20 -30 -220 -259 -332 -161 -14 337 127 124 273 108 305 253 225 -55 -191 -550 -504 -422 -702 -400 -201 61 313 345 573 619 295 278 178 -106 -186 -419 -258 -209 -22 100 232 349 288 162 213 280 425 445 552 65 -247 -461 -515 -332 -10 81 291 228 -43 -123 -218 -228 -316 -144 187 147 269 169 295 86 9 -237 -622 -830 -680 -594 -262 -222 71 174 102 144 190 309 354 220 290 248 122 -32 -343 -457 -350 -305 64 130 50 5 -149 -279 -281 -226 -238 -247 -218 -136 -126 -89 -123 -204 -252 -253 -87 209 346 253 393 351 234 20 -531 -581 -793 -522 -353 -81 41 108 -127 -267 -280 -346 -142 87 98 150 418 529 481 494 301 341 16 100 196 109 94 306 277 271 411 349 106 19 -29 -211 -248 -315 -157 -13 332 127 125 275 107 304 253 225 -54 -187 -538 -485 -410 SO 20251206T124849.3Z 9471.533904 2308 1 36.77 2 76.60 198.56 100.91 355.77 0.00 1 0.00 1 0 15 3.448e-03 3.415e-03 9.285e-04 9.270e-04 3.798830e+03 5.122201e+01 2.984379e+02 -6.000000e+00 1.697050e+02 269.46 -2.00 20.00 -2.00 505.29 0.00e+00 7.336135e-04 2.400784e-05 4.397400e-05 -6.000000e+00 9.119130e-05 269.46 -2.00 20.00 -2.00 500.52 0.00e+00 4.099552e-01 1.760758e-02 2.673419e-02 -6.000000e+00 6.578041e-02 220.42 -2.00 3.84 -2.00 502.11 0.00e+00 4.212 0.002 489.00 513.00 1 7.2802e+00 1.3394e-03 6.5091e-03 -6.0000e+00 5.0096e-03 8.8565e-03 3.5176e-04 7.6270e-04 -6.0000e+00 1.3240e-03 7.1058e-04 3.5111e-04 7.0821e-05 4.5232e-05 -6.0000e+00 2.7170e-04 -4.5289e-03 8.0790e-05 10 0 0 10 0 0 99.94 24.53 1.45 -0.00453 -0.00209 11 1 1 4 0 3341.590 33.75 20.00 24.53 15.06 0.0 948.30 -815 -67 -378 30 177 564 557 641 -60 223 105 -95 -72 -530 -74 -251 22 148 442 133 324 103 135 261 601 356 372 178 -80 -317 -600 -300 -52 146 186 -32 -60 -210 -423 -104 -264 -20 146 166 476 152 34 281 -78 -435 -538 -906 -840 -827 -485 -190 -54 173 155 201 564 460 479 348 186 225 22 -208 -446 -668 -116 -450 69 233 63 224 -150 -117 -114 -33 -243 -385 -357 -58 -156 -2 -78 -615 -16 -290 -13 306 255 361 382 444 428 -108 -624 -298 -904 -500 -322 -35 -324 218 -199 -18 -313 -443 -169 371 219 -66 281 558 354 612 88 331 -84 65 300 -42 123 -20 144 351 198 466 -134 -5 71 -219 -105 -302 -481 39 296 232 225 316 229 461 303 187 64 -300 -486 -536 -426 -791 -65 -383 30 183 579 582 672 -63 229 109 -97 -73 -531 -73 -250 22 153 465 140 341 109 142 277 623 369 393 186 -81 -317 -591 -298 -52 153 195 -34 -61 -220 -437 -106 -272 -21 153 172 496 159 35 286 -81 -437 -535 -879 -805 -804 -476 -193 -54 176 157 207 580 465 487 354 188 228 21 -205 -437 -648 -113 -442 69 230 63 222 -145 -113 -110 -32 -235 -376 -342 -56 -153 -2 -78 -613 -15 -287 -13 301 255 363 379 440 413 -105 -594 -277 -833 -466 -304 -35 -326 218 -197 -18 -307 -438 -166 368 221 -66 279 555 357 619 89 329 -84 64 290 -41 121 -19 143 345 194 457 -132 -5 68 -210 -100 -291 -466 38 289 229 224 316 230 464 307 189 64 -291 -471 -513 -411