File name: Pandora240s1_Thessaloniki_20230302_L2Fit_fwvt1c10d20221220p1-8.txt File generation date: 20251121T091548.3Z Data description: Level 2 spectral fitting results file Data file version: fwvt1c10p1-8 Local principal investigator: Alkiviadis Bais Network principal investigator: Alexander Cede Instrument type: Pandora Instrument number: 240 Spectrometer number: 1 Processing software version used: BlickP v1.8.75 Instrument operation file used: Pandora240_OF_v2d20220707.txt Instrument calibration file used: Pandora240s1_CF_v10d20221220.txt Level 1 file used: Pandora240s1_Thessaloniki_20230302_L1_smca1c10d20221220p1-8.txt Full location name: Thessaloniki Short location name: Thessaloniki Country of location: Greece Location latitude [deg]: 40.6336 Location longitude [deg]: 22.9561 Location altitude [m]: 60 Local noon date: 20230302 Notes on s-number (L1 configuration): Corrections NOT applied although requested by the s-number are blind pixel subtraction, 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 1 for testing the CF, absolute calibration done with an uncalibrated lamp, certificate of FEL1441 (GSFC) used First and last pixel inside fitting window: 1620 1769 Nominal wavelengths inside fitting window [nm]: 492.0731 492.1939 492.3148 492.4356 492.5563 492.6771 492.7979 492.9186 493.0393 493.1600 493.2807 493.4013 493.5220 493.6426 493.7632 493.8838 494.0044 494.1249 494.2455 494.3660 494.4865 494.6070 494.7274 494.8479 494.9683 495.0887 495.2091 495.3295 495.4499 495.5702 495.6906 495.8109 495.9312 496.0514 496.1717 496.2919 496.4122 496.5324 496.6526 496.7727 496.8929 497.0130 497.1331 497.2532 497.3733 497.4934 497.6134 497.7335 497.8535 497.9735 498.0934 498.2134 498.3333 498.4533 498.5732 498.6931 498.8129 498.9328 499.0526 499.1724 499.2922 499.4120 499.5318 499.6515 499.7713 499.8910 500.0107 500.1303 500.2500 500.3696 500.4893 500.6089 500.7285 500.8480 500.9676 501.0871 501.2066 501.3261 501.4456 501.5651 501.6845 501.8039 501.9233 502.0427 502.1621 502.2815 502.4008 502.5201 502.6394 502.7587 502.8780 502.9972 503.1165 503.2357 503.3549 503.4740 503.5932 503.7123 503.8315 503.9506 504.0697 504.1887 504.3078 504.4268 504.5458 504.6648 504.7838 504.9028 505.0217 505.1406 505.2596 505.3784 505.4973 505.6162 505.7350 505.8538 505.9726 506.0914 506.2102 506.3289 506.4476 506.5664 506.6851 506.8037 506.9224 507.0410 507.1596 507.2782 507.3968 507.5154 507.6339 507.7525 507.8710 507.9895 508.1079 508.2264 508.3448 508.4633 508.5817 508.7000 508.8184 508.9368 509.0551 509.1734 509.2917 509.4100 509.5282 509.6465 509.7647 509.8829 --------------------------------------------------------------------------------------- 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-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 20230302T100317.5Z 8461.418953 342 1 54.64 2 48.62 167.70 100.56 40.63 0.00 1 0.00 1 0 15 4.382e-03 4.343e-03 1.041e-03 1.036e-03 1.844710e+03 5.618751e+01 1.526759e+02 -6.000000e+00 2.198206e+02 276.14 -2.00 20.00 -2.00 505.20 0.00e+00 3.646067e-04 2.721791e-05 2.208669e-05 -6.000000e+00 1.177853e-04 276.14 -2.00 20.00 -2.00 500.15 0.00e+00 -6.674593e-02 2.011874e-02 1.093223e-02 -6.000000e+00 8.426894e-02 225.58 -2.00 3.50 -2.00 501.77 0.00e+00 1.510 0.000 489.00 513.00 1 7.2997e+00 1.5308e-03 2.4121e-03 -6.0000e+00 6.4130e-03 1.4719e-02 4.1072e-04 3.1678e-04 -6.0000e+00 1.7182e-03 1.0678e-03 -4.3787e-04 6.9979e-05 1.7016e-05 -6.0000e+00 3.0372e-04 4.9036e-04 -1.0348e-05 11 8 0 10 0 0 99.42 23.67 2.24 0.00049 -0.00255 29 3 1 2 0 1881.253 20.00 23.67 18.85 28.00 28.00 18.3 1020.91 -1090 -814 -204 -297 111 453 742 926 718 358 414 117 -308 -182 -464 -178 -386 18 99 -22 188 61 -46 445 351 285 531 508 -425 -484 -673 -554 -867 170 91 96 191 252 -34 -266 -105 -100 335 189 512 360 546 604 409 -147 -349 -966 -967 -702 -561 -281 -47 160 67 75 -59 305 274 311 49 360 79 -48 -120 -688 -115 -196 -296 50 70 356 -41 -306 -440 -631 -487 -170 -463 -159 -77 9 -186 -32 -192 -82 -347 251 604 576 655 478 641 -373 -412 -663 -1004 -580 -859 176 168 186 297 -555 -516 -282 -220 -129 319 97 281 558 586 896 255 471 202 -202 124 346 2 70 277 -34 465 748 78 199 -30 -362 -1249 -605 -193 -120 -74 250 156 143 163 312 747 348 90 -370 -493 -795 -1094 -827 -214 -314 119 491 821 1034 801 392 446 125 -324 -187 -479 -185 -405 19 106 -24 201 65 -49 475 380 307 567 542 -442 -497 -688 -572 -903 177 99 104 203 266 -36 -274 -108 -104 350 201 548 387 580 632 423 -149 -354 -963 -949 -690 -547 -284 -49 166 70 77 -61 316 286 324 50 365 79 -47 -118 -663 -109 -188 -287 49 70 353 -39 -291 -415 -592 -457 -163 -447 -152 -76 9 -185 -31 -190 -81 -333 242 598 574 638 460 607 -360 -385 -607 -920 -536 -821 168 167 183 288 -532 -484 -266 -212 -126 315 97 281 537 568 873 247 464 197 -194 117 331 2 67 267 -32 441 713 75 185 -28 -331 -1129 -548 -181 -115 -72 234 148 135 155 298 713 333 85 -337 -454 -732 SQ 20230302T134313.1Z 8461.571679 477 1 54.76 2 63.61 232.46 66.85 74.02 0.00 1 0.00 1 0 15 4.762e-03 4.760e-03 1.602e-03 1.594e-03 2.491292e+03 8.825919e+01 1.995196e+02 -6.000000e+00 2.388522e+02 276.15 -2.00 20.00 -2.00 505.21 0.00e+00 4.871936e-04 4.183218e-05 2.961983e-05 -6.000000e+00 1.279898e-04 276.15 -2.00 20.00 -2.00 500.15 0.00e+00 5.561317e-02 3.100498e-02 1.551003e-02 -6.000000e+00 9.157038e-02 225.58 -2.00 3.50 -2.00 501.78 0.00e+00 2.237 0.000 489.00 513.00 1 8.0292e+00 2.3590e-03 3.5276e-03 -6.0000e+00 6.9685e-03 1.4176e-02 6.3207e-04 4.4027e-04 -6.0000e+00 1.8673e-03 4.5885e-04 -3.1610e-06 1.0741e-04 2.2510e-05 -6.0000e+00 3.3004e-04 1.8873e-03 -2.4670e-05 12 8 8 10 0 0 98.57 23.72 2.12 0.00189 0.00187 14 2 1 2 0 3910.608 20.00 23.72 19.18 27.75 28.00 18.1 1021.25 -1034 -1288 -220 -64 80 861 461 1548 266 523 187 142 -444 -215 -573 -124 -347 -34 361 -17 353 177 76 354 329 341 572 391 -382 -179 -776 -681 -831 158 444 -12 -68 -42 -263 -487 146 -281 253 314 463 402 651 511 786 -16 -611 -574 -952 -1166 -728 -306 -283 262 27 167 -123 241 78 237 293 31 148 -51 -57 -850 -524 -170 -296 317 493 840 119 -116 -559 -722 -573 -355 -231 -121 -30 144 -66 -262 -31 -193 -394 255 362 674 799 516 292 2 -371 -997 -697 -740 -555 161 195 128 109 -611 -600 -521 306 -517 51 521 225 158 711 778 7 441 393 53 143 508 67 -144 185 155 334 892 -19 106 -141 -448 -1244 -854 -57 -44 268 148 136 8 400 458 313 108 -188 -334 -236 -242 -1039 -1320 -221 -67 86 931 513 1720 285 572 202 153 -468 -221 -593 -129 -364 -37 385 -18 377 190 81 378 356 368 612 418 -397 -184 -797 -707 -869 164 483 -13 -73 -45 -275 -502 150 -294 265 334 497 433 693 536 812 -17 -596 -572 -939 -1159 -714 -311 -294 274 28 172 -127 251 82 248 304 32 148 -50 -56 -823 -499 -164 -289 311 492 832 115 -111 -529 -680 -541 -342 -223 -117 -30 144 -63 -259 -31 -184 -378 246 359 670 776 495 276 2 -346 -913 -636 -685 -529 153 193 126 105 -586 -562 -493 294 -508 51 519 225 152 688 759 7 435 384 51 135 485 65 -138 179 149 317 852 -18 98 -129 -411 -1129 -780 -54 -41 248 140 129 7 381 439 301 104 -177 -305 -219 -223 SQ 20230302T143736.5Z 8461.609450 517 1 34.57 2 72.35 243.57 57.07 81.48 0.00 1 0.00 1 0 15 4.231e-03 4.193e-03 2.267e-04 2.260e-04 3.111102e+03 1.257451e+01 2.452667e+02 -6.000000e+00 2.122221e+02 276.15 -2.00 20.00 -2.00 505.23 0.00e+00 6.056594e-04 5.936025e-06 3.585805e-05 -6.000000e+00 1.137274e-04 276.15 -2.00 20.00 -2.00 500.28 0.00e+00 1.644233e-01 4.364300e-03 1.986296e-02 -6.000000e+00 8.136799e-02 225.59 -2.00 3.50 -2.00 501.89 0.00e+00 3.254 0.001 489.00 513.00 1 3.6877e+00 3.3219e-04 5.0121e-03 -6.0000e+00 6.1920e-03 5.8554e-03 8.8683e-05 5.7348e-04 -6.0000e+00 1.6594e-03 3.0195e-02 -1.0494e-04 1.5173e-05 2.8422e-05 -6.0000e+00 2.9329e-04 2.6785e-03 -3.5024e-05 12 9 1 12 1 1 99.81 23.72 1.93 0.00268 0.00195 399 70 1 2 0 83.610 20.00 23.72 18.86 27.75 28.00 17.9 1017.39 -998 -842 -115 -210 210 515 688 606 652 465 498 250 -270 -137 -404 -277 -296 45 136 77 256 109 80 298 301 463 582 530 -455 -378 -635 -713 -852 253 163 166 203 -27 -220 -357 -108 -117 164 284 485 326 665 423 459 -193 -405 -780 -984 -890 -540 -218 -84 44 10 7 -43 339 224 501 271 135 199 47 -195 -569 -353 -228 -260 115 173 16 -109 -259 -443 -687 -311 -239 -365 -259 -28 44 -136 -103 -225 -161 -244 242 442 473 650 498 432 47 -511 -849 -962 -869 -701 209 47 124 110 -378 -463 -373 -144 -166 74 285 349 482 607 757 378 516 239 18 47 443 70 152 410 291 420 733 82 -22 -76 -372 -950 -511 -214 -122 -46 106 40 116 238 428 510 456 172 -422 -568 -751 -983 -842 -116 -217 220 546 728 648 694 486 525 263 -279 -142 -409 -282 -304 46 146 83 274 117 86 319 319 488 624 568 -477 -382 -641 -729 -877 264 174 177 212 -28 -226 -371 -112 -120 170 297 513 338 697 439 468 -198 -399 -773 -962 -874 -535 -218 -87 45 10 7 -44 349 232 517 279 136 197 47 -195 -558 -342 -223 -257 112 172 16 -107 -251 -426 -658 -298 -233 -349 -247 -27 43 -135 -101 -222 -159 -237 236 444 478 640 483 422 44 -485 -788 -893 -818 -664 202 47 121 108 -368 -442 -359 -141 -161 73 284 347 480 610 750 372 508 233 17 45 432 68 146 395 280 399 700 79 -21 -71 -349 -881 -476 -203 -118 -44 102 39 112 232 421 501 448 166 -405 -531 -700