File name: Pandora134s1_BristolPA_20240418_L2Fit_fwvt1c5d20230501p1-8.txt File generation date: 20240825T110256.9Z Data description: Level 2 spectral fitting results file Data file version: fwvt1c5p1-8 Local principal investigator: Lukas Valin Network principal investigator: Alexander Cede Instrument type: Pandora Instrument number: 134 Spectrometer number: 1 Processing software version used: BlickP v1.8.59 Instrument operation file used: Pandora134_OF_v6d20200325.txt Instrument calibration file used: Pandora134s1_CF_v5d20230501.txt Level 1 file used: Pandora134s1_BristolPA_20240418_L1_smca1c5d20230501p1-8.txt Full location name: Bristol, Pennsylvania; Philadelphia-Camden-Wilmington, PA-NJ-DE-MD; A420170012 Short location name: BristolPA Country of location: United States Location latitude [deg]: 40.1074 Location longitude [deg]: -74.8824 Location altitude [m]: 10 Local noon date: 20240418 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 3 for VP from 20230501 First and last pixel inside fitting window: 1648 1797 Nominal wavelengths inside fitting window [nm]: 492.0138 492.1355 492.2572 492.3789 492.5006 492.6223 492.7439 492.8656 492.9872 493.1088 493.2304 493.3519 493.4735 493.5950 493.7165 493.8380 493.9595 494.0810 494.2024 494.3239 494.4453 494.5667 494.6881 494.8094 494.9308 495.0521 495.1734 495.2947 495.4160 495.5373 495.6585 495.7798 495.9010 496.0222 496.1434 496.2646 496.3857 496.5068 496.6280 496.7491 496.8701 496.9912 497.1123 497.2333 497.3543 497.4753 497.5963 497.7173 497.8382 497.9592 498.0801 498.2010 498.3219 498.4427 498.5636 498.6844 498.8052 498.9260 499.0468 499.1676 499.2883 499.4091 499.5298 499.6505 499.7712 499.8919 500.0125 500.1331 500.2538 500.3744 500.4949 500.6155 500.7361 500.8566 500.9771 501.0976 501.2181 501.3386 501.4590 501.5794 501.6999 501.8203 501.9406 502.0610 502.1813 502.3017 502.4220 502.5423 502.6626 502.7828 502.9031 503.0233 503.1435 503.2637 503.3839 503.5041 503.6242 503.7443 503.8644 503.9845 504.1046 504.2247 504.3447 504.4648 504.5848 504.7048 504.8247 504.9447 505.0646 505.1846 505.3045 505.4244 505.5442 505.6641 505.7839 505.9038 506.0236 506.1433 506.2631 506.3829 506.5026 506.6223 506.7420 506.8617 506.9814 507.1010 507.2207 507.3403 507.4599 507.5795 507.6991 507.8186 507.9381 508.0577 508.1772 508.2966 508.4161 508.5355 508.6550 508.7744 508.8938 509.0132 509.1325 509.2519 509.3712 509.4905 509.6098 509.7291 509.8483 509.9676 --------------------------------------------------------------------------------------- Column 1: Two letter code of measurement routine Column 2: UT date and time for center-time of measurement, yyyymmddThhmmssZ (ISO 8601) Column 3: Fractional days since 1-Jan-2000 midnight for center-time of measurement Column 4: Routine count (1 for the first routine of the day, 2 for the second, etc.) Column 5: Repetition count (1 for the first set in the routine, 2 for the second, etc.) Column 6: Total duration of measurement set in seconds Column 7: Data processing type index Column 8: Solar zenith angle for center-time of measurement in degree Column 9: Solar azimuth for center-time of measurement in degree, 0=north, increases clockwise Column 10: Lunar zenith angle for center-time of measurement in degree Column 11: Lunar azimuth for center-time of measurement in degree, 0=north, increases clockwise Column 12: Pointing zenith angle in degree, absolute or relative (see next column), 999=tracker not used Column 13: Zenith pointing mode: zenith angle is... 0=absolute, 1=relative to sun, 2=relative to moon Column 14: Pointing azimuth in degree, increases clockwise, absolute (0=north) or relative (see next column), 999=tracker not used Column 15: Azimuth pointing mode: like zenith angle mode but also fixed scattering angles relative to sun (3) or moon (4) Column 16: Fitting result index: 0=no error or warning, 1, 2=warning, >2=error Column 17: Number of function evaluations used, 0=linear fitting or fitting not successful or no fitting done Column 18: rms of unweighted fitting residuals, -9=fitting not successful Column 19: Normalized rms of fitting residuals weighted with independent uncertainty, -9=fitting not successful or no uncertainty used Column 20: Expected rms based on independent uncertainty, -9=fitting not successful or no uncertainty given Column 21: Expected normalized weighted rms based on independent uncertainty, -9=fitting not successful or no uncertainty given Column 22: Water vapor slant column amount [moles per square meter], -9e99=fitting not successful Column 23: Independent uncertainty of water vapor slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -3=spectral fitting was done, but no independent uncertainty could be retrieved, -5=no independent uncertainty input was given, -9=spectral fitting not successful Column 24: Structured uncertainty of water vapor slant column amount [moles per square meter], -7=not given since method "MEAS" was chosen, -9=spectral fitting not successful Column 25: Common uncertainty of water vapor slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -6=no common uncertainty input was given, -7=not given since method "MEAS" was chosen, -9=spectral fitting not successful Column 26: rms-based uncertainty of water vapor slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -3=spectral fitting was done, but no rms-based uncertainty could be retrieved, -9=spectral fitting not successful Column 27: Water vapor effective temperature [K] Column 28: Independent uncertainty of water vapor effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -2=no temperature fitting was requested and output for effective temperature and structured uncertainty of it is based on f-code, -3=spectral fitting was done, but no independent uncertainty could be retrieved, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -5=no independent uncertainty input was given, -9=spectral fitting not successful Column 29: Structured uncertainty of water vapor effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -9=spectral fitting not successful Column 30: Common uncertainty of water vapor effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -2=no temperature fitting was requested and output for effective temperature and structured uncertainty of it is based on f-code, -3=spectral fitting was done, but no common uncertainty could be retrieved, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -6=no common uncertainty input was given, -9=spectral fitting not successful Column 31: Effective water vapor fitting wavelength [nm], -9=fitting not successful Column 32: Diffuse correction applied before fitting at effective fitting wavelength for water vapor [%], 0=no diffuse correction applied or fitting not requested, >0=measured diffuse correction, <0=(negative value of) calculated diffuse correction Column 33: Nitrogen dioxide slant column amount [moles per square meter], -9e99=fitting not successful Column 34: Independent uncertainty of nitrogen dioxide slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -3=spectral fitting was done, but no independent uncertainty could be retrieved, -5=no independent uncertainty input was given, -9=spectral fitting not successful Column 35: Structured uncertainty of nitrogen dioxide slant column amount [moles per square meter], -7=not given since method "MEAS" was chosen, -9=spectral fitting not successful Column 36: Common uncertainty of nitrogen dioxide slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -6=no common uncertainty input was given, -7=not given since method "MEAS" was chosen, -9=spectral fitting not successful Column 37: rms-based uncertainty of nitrogen dioxide slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -3=spectral fitting was done, but no rms-based uncertainty could be retrieved, -9=spectral fitting not successful Column 38: Nitrogen dioxide effective temperature [K] Column 39: Independent uncertainty of nitrogen dioxide effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -2=no temperature fitting was requested and output for effective temperature and structured uncertainty of it is based on f-code, -3=spectral fitting was done, but no independent uncertainty could be retrieved, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -5=no independent uncertainty input was given, -9=spectral fitting not successful Column 40: Structured uncertainty of nitrogen dioxide effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -9=spectral fitting not successful Column 41: Common uncertainty of nitrogen dioxide effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -2=no temperature fitting was requested and output for effective temperature and structured uncertainty of it is based on f-code, -3=spectral fitting was done, but no common uncertainty could be retrieved, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -6=no common uncertainty input was given, -9=spectral fitting not successful Column 42: Effective nitrogen dioxide fitting wavelength [nm], -9=fitting not successful Column 43: Diffuse correction applied before fitting at effective fitting wavelength for nitrogen dioxide [%], 0=no diffuse correction applied or fitting not requested, >0=measured diffuse correction, <0=(negative value of) calculated diffuse correction Column 44: Ozone slant column amount [moles per square meter], -9e99=fitting not successful Column 45: Independent uncertainty of ozone slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -3=spectral fitting was done, but no independent uncertainty could be retrieved, -5=no independent uncertainty input was given, -9=spectral fitting not successful Column 46: Structured uncertainty of ozone slant column amount [moles per square meter], -7=not given since method "MEAS" was chosen, -9=spectral fitting not successful Column 47: Common uncertainty of ozone slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -6=no common uncertainty input was given, -7=not given since method "MEAS" was chosen, -9=spectral fitting not successful Column 48: rms-based uncertainty of ozone slant column amount [moles per square meter], -1=cross section is zero in this wavelength range, -3=spectral fitting was done, but no rms-based uncertainty could be retrieved, -9=spectral fitting not successful Column 49: Ozone effective temperature [K] Column 50: Independent uncertainty of ozone effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -2=no temperature fitting was requested and output for effective temperature and structured uncertainty of it is based on f-code, -3=spectral fitting was done, but no independent uncertainty could be retrieved, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -5=no independent uncertainty input was given, -9=spectral fitting not successful Column 51: Structured uncertainty of ozone effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -9=spectral fitting not successful Column 52: Common uncertainty of ozone effective temperature [K], -1=temperature fitting was requested, but cross section is zero in this wavelength range, -2=no temperature fitting was requested and output for effective temperature and structured uncertainty of it is based on f-code, -3=spectral fitting was done, but no common uncertainty could be retrieved, -4=temperature fitting was requested, but differential optical depth is too small to retrieve the temperature, -6=no common uncertainty input was given, -9=spectral fitting not successful Column 53: Effective ozone fitting wavelength [nm], -9=fitting not successful Column 54: Diffuse correction applied before fitting at effective fitting wavelength for ozone [%], 0=no diffuse correction applied or fitting not requested, >0=measured diffuse correction, <0=(negative value of) calculated diffuse correction Column 55: Molecular scattering air mass factor used for molecular scattering subtraction before the fitting Column 56: Estimated uncertainty of molecular scattering air mass factor, -9=molecular scattering was not subtracted before the fitting Column 57: Lower limit used for wavelength scaling [nm] Column 58: Upper limit used for wavelength scaling [nm] Column 59: Order of smoothing polynomial Column 60: Smoothing polynomial coefficient, order 0 Column 61: Independent uncertainty of smoothing polynomial coefficient, order 0, -9=fitting not successful Column 62: Structured uncertainty of smoothing polynomial coefficient, order 0, -9=fitting not successful Column 63: Common uncertainty of smoothing polynomial coefficient, order 0, -9=fitting not successful Column 64: rms-based uncertainty of smoothing polynomial coefficient, order 0, -9=fitting not successful Column 65: Smoothing polynomial coefficient, order 1 Column 66: Independent uncertainty of smoothing polynomial coefficient, order 1, -9=fitting not successful Column 67: Structured uncertainty of smoothing polynomial coefficient, order 1, -9=fitting not successful Column 68: Common uncertainty of smoothing polynomial coefficient, order 1, -9=fitting not successful Column 69: rms-based uncertainty of smoothing polynomial coefficient, order 1, -9=fitting not successful Column 70: Mean value of measured data inside fitting window [same units as measurements] Column 71: Wavelength change polynomial coefficient, order 0 Column 72: Independent uncertainty of wavelength change polynomial coefficient, order 0, -9=fitting not successful Column 73: Structured uncertainty of wavelength change polynomial coefficient, order 0, -9=fitting not successful Column 74: Common uncertainty of wavelength change polynomial coefficient, order 0, -9=fitting not successful Column 75: rms-based uncertainty of wavelength change polynomial coefficient, order 0, -9=fitting not successful Column 76: L1 based wavelength change coefficient, order 0, -9=no L1 wavelength change determination Column 77: L1 based wavelength change coefficient, order 1, -9=no L1 wavelength change determination Column 78: L2Fit data quality flag: 0=assured high quality, 1=assured medium quality, 2=assured low quality, 10=not-assured high quality, 11=not-assured medium quality, 12=not-assured low quality Column 79: Sum over 2^i using those i, for which the corresponding L2Fit data quality parameter exceeds the DQ1 limit, 0=L1 data quality above 0, 1=Spectral fitting was not successful, 2=Wavelength shift too large, 3=Normalized rms of fitting residuals weighted with independent uncertainty too large Column 80: Sum over 2^i using those i, for which the corresponding L2Fit data quality parameter exceeds the DQ2 limit (same parameters as for DQ1) Column 81: L1 data quality flag: 0=assured high quality, 1=assured medium quality, 2=assured low quality, 10=not-assured high quality, 11=not-assured medium quality, 12=not-assured low quality Column 82: Sum over 2^i using those i, for which the corresponding L1 data quality parameter exceeds the DQ1 limit, 0=Saturated data, 1=Too few dark counts measurements, 2=No temperature given or effective temperature too different from the reference temperature, 3=Dark count too high, 4=Unsuccessful dark background fitting, 5=The dark count differs significantly from the dark map for too many pixels, 6=Absolute value of estimated average residual stray light level too high, 7=Although attempted, no wavelength change could be retrieved, 8=Absolute value of retrieved wavelength shift too large, 9=Retrieved wavelength shift differs too much from the shift predicted by the effective temperature Column 83: Sum over 2^i using those i, for which the corresponding L1 data quality parameter exceeds the DQ2 limit (same parameters as for DQ1) Column 84: Atmospheric variability [%], 999=no atmospheric variability was determined Column 85: Wavelength effective temperature [°C], 999=no effective temperature given Column 86: Estimated average residual stray light level [%] (only valid for stray light correction methods 2 and higher) Column 87: Retrieved wavelength shift from L1 data [nm], -9=no wavelength change determination Column 88: Retrieved total wavelength shift [nm], -9=no wavelength change fitting Column 89: Number of bright count cycles Column 90: Number of dark count cycles Column 91: Effective position of filterwheel #1, 0=filterwheel not used, 1-9 are valid positions Column 92: Effective position of filterwheel #2, 0=filterwheel not used, 1-9 are valid positions Column 93: Sum over 2^i, 0=spectra were interpolated in time, 1=spectra are corrected for off-target signal Column 94: Integration time [ms] Column 95: Temperature at electronics board [°C], 999=no temperature signal Column 96: Spectrometer control temperature [°C], 999=no temperature signal Column 97: Auxiliary spectrometer temperature [°C], 999=no temperature signal Column 98: Temperature in head sensor [°C], 999=no temperature signal Column 99: Humidity in head sensor [%], -9=no humidity signal Column 100: Pressure in head sensor [hPa], -9=no pressure signal Columns 101-250: Unweighted slant column residuals for each pixel inside the fitting window multiplied by 1e5, 9e5=pixel was not used for fitting Columns 251-400: Normalized slant column residuals weighted with independent instrumental uncertainty for each pixel inside the fitting window multiplied by 1e5, 9e5=pixel was not used for fitting, 0=no independent instrumental uncertainty was given --------------------------------------------------------------------------------------- SQ 20240418T203512.2Z 8874.857780 555 1 8.99 2 55.60 254.43 70.77 88.82 0.00 1 0.00 1 0 15 2.469e-03 2.466e-03 9.998e-04 9.980e-04 1.545366e+03 5.836632e+01 1.202793e+02 -6.000000e+00 1.376363e+02 281.34 -2.00 20.00 -2.00 505.26 0.00e+00 1.009691e-03 2.740084e-05 5.561124e-05 -6.000000e+00 6.933121e-05 281.34 -2.00 20.00 -2.00 500.42 0.00e+00 9.454730e-02 1.977834e-02 1.266389e-02 -6.000000e+00 4.911526e-02 229.04 -2.00 3.49 -2.00 502.02 0.00e+00 1.764 0.000 489.00 513.00 1 6.0909e+00 1.5170e-03 3.0460e-03 -6.0000e+00 3.7679e-03 1.5523e-02 3.8999e-04 3.0111e-04 -6.0000e+00 9.6972e-04 3.3695e-03 1.9160e-04 1.0096e-04 2.4803e-05 -6.0000e+00 2.5532e-04 -1.7859e-04 1.5257e-05 12 1 1 12 1 1 93.38 21.55 0.28 -0.00018 0.00115 9 10 1 4 0 896.576 27.15 20.00 21.55 16.72 30.6 1023.75 -275 -185 -182 -94 -43 278 463 490 421 387 212 256 17 5 -334 -297 -176 -314 -178 -154 -269 -93 225 214 153 188 -39 -121 -44 -99 -124 44 -192 -114 -281 -123 -180 212 139 -89 -34 -24 265 227 94 383 276 320 42 149 -204 -86 -457 -158 -23 -163 -42 -68 140 201 251 258 89 93 55 83 148 -249 -187 -113 -176 -182 -144 -565 -193 -388 -304 -180 -151 -293 -498 -517 -403 -27 -128 -36 221 190 287 475 466 598 407 209 144 26 -173 -74 -119 -199 -214 -168 -262 -384 -489 -301 56 24 -45 42 30 -10 72 140 488 493 438 9 55 -62 55 -135 174 212 -80 -178 123 -190 -39 -47 -64 -282 -67 -61 171 115 73 56 356 306 311 259 295 432 -51 -311 -174 -373 -333 -299 -270 -183 -178 -95 -44 283 480 513 442 403 218 258 17 5 -339 -302 -181 -317 -184 -161 -284 -98 237 226 162 200 -41 -126 -45 -98 -125 44 -197 -116 -294 -130 -190 220 143 -91 -35 -24 273 237 99 404 289 332 43 148 -205 -84 -445 -154 -23 -162 -43 -69 143 205 256 264 91 96 56 83 146 -243 -188 -111 -173 -179 -143 -551 -189 -377 -301 -176 -146 -283 -484 -507 -398 -27 -124 -35 215 185 279 463 459 596 411 211 144 26 -167 -72 -113 -188 -202 -161 -251 -380 -492 -301 55 24 -44 41 30 -10 71 139 485 492 439 9 54 -61 55 -133 171 209 -79 -177 121 -188 -38 -46 -62 -269 -63 -57 160 109 71 56 355 297 303 254 290 425 -50 -302 -171 -360 -317 -285 SQ 20240418T211411.3Z 8874.884853 590 1 14.08 2 62.88 261.65 63.57 95.16 0.00 1 0.00 1 0 15 2.560e-03 2.557e-03 1.270e-03 1.268e-03 1.833255e+03 7.499604e+01 1.364014e+02 -6.000000e+00 1.426716e+02 281.34 -2.00 20.00 -2.00 505.27 0.00e+00 1.226523e-03 3.487153e-05 6.731081e-05 -6.000000e+00 7.186651e-05 281.34 -2.00 20.00 -2.00 500.45 0.00e+00 1.679903e-01 2.511290e-02 1.620781e-02 -6.000000e+00 5.090347e-02 229.04 -2.00 3.49 -2.00 502.04 0.00e+00 2.182 0.000 489.00 513.00 1 6.8316e+00 1.9259e-03 3.7654e-03 -6.0000e+00 3.9051e-03 1.3166e-02 4.9413e-04 3.7496e-04 -6.0000e+00 1.0049e-03 1.4998e-03 8.8088e-05 1.2804e-04 2.9068e-05 -6.0000e+00 2.6467e-04 4.3904e-03 6.4292e-07 12 1 1 12 1 1 98.80 21.56 0.34 0.00439 0.00500 6 4 1 4 0 2011.258 27.15 20.00 21.56 16.52 30.8 1023.04 -309 -241 -176 235 -108 198 308 315 390 226 83 -87 -176 71 -387 -158 -271 -113 -329 -130 -275 122 85 192 275 105 -12 -102 -92 88 -9 -51 -145 -22 -137 -115 -144 112 165 -112 -89 24 413 219 328 363 500 115 475 333 -288 -77 -411 -252 -194 -363 -160 44 -5 95 225 538 192 74 35 110 -73 -162 -420 92 -285 -166 -243 -277 -174 -282 -206 -231 -442 -455 -387 -157 -359 -66 -53 -77 111 340 376 404 578 503 207 242 -31 7 -22 -120 -10 -170 -343 -170 -285 -387 -597 -508 34 99 26 185 33 -107 300 310 528 575 389 -3 -26 -3 -112 99 123 293 -177 -88 -47 -201 -141 111 -223 -136 4 8 188 211 -5 -132 451 65 289 204 365 174 -87 -290 -258 -175 -379 -407 -300 -236 -174 237 -111 206 322 331 409 236 86 -90 -179 72 -390 -159 -275 -115 -342 -136 -288 128 89 201 288 111 -12 -106 -94 88 -9 -51 -146 -22 -143 -120 -150 116 170 -114 -91 24 425 228 342 378 518 119 484 335 -286 -75 -400 -247 -192 -366 -163 45 -6 97 230 548 197 76 36 111 -73 -161 -415 90 -277 -162 -239 -274 -172 -279 -201 -224 -425 -437 -374 -152 -348 -65 -52 -76 110 335 370 398 571 501 208 243 -31 7 -21 -115 -9 -160 -323 -162 -276 -383 -596 -506 33 97 26 180 33 -106 300 312 532 581 393 -3 -27 -3 -110 97 121 287 -174 -86 -46 -197 -138 109 -216 -131 3 7 179 203 -4 -129 444 64 286 203 363 174 -87 -287 -251 -169 -363 -390 SQ 20240418T211812.0Z 8874.887639 593 1 5.54 2 63.64 262.35 62.83 95.84 0.00 1 0.00 1 0 15 3.533e-03 3.532e-03 2.804e-03 2.800e-03 2.051875e+03 1.663862e+02 1.707870e+02 -6.000000e+00 1.969139e+02 281.34 -2.00 20.00 -2.00 505.27 0.00e+00 1.347179e-03 7.702363e-05 7.121647e-05 -6.000000e+00 9.918950e-05 281.34 -2.00 20.00 -2.00 500.46 0.00e+00 2.197939e-01 5.540996e-02 1.461100e-02 -6.000000e+00 7.026103e-02 229.04 -2.00 3.49 -2.00 502.05 0.00e+00 2.239 0.000 489.00 513.00 1 6.9991e+00 4.2499e-03 3.8428e-03 -6.0000e+00 5.3901e-03 1.2452e-02 1.0894e-03 3.6835e-04 -6.0000e+00 1.3871e-03 1.2505e-03 1.9111e-04 2.8223e-04 3.3609e-05 -6.0000e+00 3.6527e-04 -7.8741e-04 2.6372e-05 12 1 1 12 1 1 999.00 21.56 0.28 -0.00079 0.00054 1 3 1 4 0 2766.858 27.15 20.00 21.56 16.50 30.7 1022.98 -327 -365 -237 -131 -53 210 472 87 482 734 -133 -307 -229 227 -387 -331 -501 -68 -428 188 -252 -141 -336 693 227 266 244 473 363 76 -232 -319 -190 -288 -213 -34 -330 487 -85 -315 -247 170 343 -71 1015 226 222 105 125 258 -329 -124 -184 -258 198 -202 329 -180 106 -271 551 -22 119 98 -300 -10 -277 545 221 -326 -352 162 -192 -149 -329 -188 -271 -445 -297 -388 -588 -391 -186 496 -311 -60 4 47 -84 544 582 816 527 471 452 360 -28 -37 -319 -494 29 -137 -572 -189 -697 -824 169 134 363 360 266 -415 162 41 517 -84 211 -173 150 -30 -144 131 373 218 404 -145 -194 -128 -809 360 -285 -224 -378 -167 -19 542 -140 209 206 722 147 536 283 466 -30 -326 -333 -233 -560 -604 -319 -357 -234 -132 -54 217 492 91 504 764 -138 -315 -233 229 -390 -333 -509 -70 -444 196 -263 -147 -352 722 238 279 255 491 372 77 -232 -319 -192 -295 -221 -36 -344 503 -87 -323 -253 174 353 -73 1051 235 230 108 127 259 -327 -122 -179 -252 196 -203 334 -184 108 -277 561 -23 122 100 -306 -10 -278 540 218 -320 -345 159 -189 -148 -327 -186 -266 -433 -287 -375 -569 -379 -181 483 -305 -59 4 47 -83 535 575 810 527 472 450 355 -28 -36 -305 -468 28 -130 -556 -186 -694 -819 167 131 354 352 262 -414 162 41 520 -85 213 -175 150 -30 -142 129 365 214 398 -143 -191 -126 -797 351 -278 -216 -362 -159 -19 521 -136 204 204 713 146 532 282 464 -30 -322 -325 -224 -538 -581