File name: Pandora163s1_Tsukuba-NIES-West_20241006_L2Fit_fwvt1c12d20240101p1-8.txt File generation date: 20241007T171531.6Z Data description: Level 2 spectral fitting results file Data file version: fwvt1c12p1-8 Local principal investigator: Jinsoo Park Network principal investigator: Alexander Cede Instrument type: Pandora Instrument number: 163 Spectrometer number: 1 Processing software version used: BlickP v1.8.60 Instrument operation file used: Pandora163_OF_v9d20190605.txt Instrument calibration file used: Pandora163s1_CF_v12d20240101.txt Level 1 file used: Pandora163s1_Tsukuba-NIES-West_20241006_L1_smca1c12d20240101p1-8.txt Full location name: National Institute for Environmental Studies (NIES), lidar building Short location name: Tsukuba-NIES-West Country of location: Japan Location latitude [deg]: 36.0506 Location longitude [deg]: 140.1202 Location altitude [m]: 30 Local noon date: 20241006 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 11 for fieldcal for VP from 20240101 First and last pixel inside fitting window: 1645 1795 Nominal wavelengths inside fitting window [nm]: 492.0139 492.1349 492.2559 492.3768 492.4978 492.6187 492.7397 492.8606 492.9815 493.1023 493.2232 493.3440 493.4649 493.5857 493.7065 493.8272 493.9480 494.0687 494.1895 494.3102 494.4309 494.5515 494.6722 494.7928 494.9135 495.0341 495.1547 495.2752 495.3958 495.5163 495.6368 495.7574 495.8778 495.9983 496.1188 496.2392 496.3596 496.4800 496.6004 496.7208 496.8412 496.9615 497.0818 497.2021 497.3224 497.4427 497.5629 497.6832 497.8034 497.9236 498.0438 498.1640 498.2841 498.4042 498.5244 498.6445 498.7646 498.8846 499.0047 499.1247 499.2447 499.3647 499.4847 499.6047 499.7246 499.8446 499.9645 500.0844 500.2043 500.3241 500.4440 500.5638 500.6836 500.8034 500.9232 501.0430 501.1627 501.2824 501.4022 501.5219 501.6415 501.7612 501.8808 502.0005 502.1201 502.2397 502.3592 502.4788 502.5983 502.7179 502.8374 502.9569 503.0763 503.1958 503.3152 503.4347 503.5541 503.6735 503.7928 503.9122 504.0315 504.1508 504.2701 504.3894 504.5087 504.6280 504.7472 504.8664 504.9856 505.1048 505.2240 505.3431 505.4622 505.5814 505.7005 505.8195 505.9386 506.0576 506.1767 506.2957 506.4147 506.5337 506.6526 506.7716 506.8905 507.0094 507.1283 507.2472 507.3660 507.4849 507.6037 507.7225 507.8413 507.9600 508.0788 508.1975 508.3163 508.4350 508.5536 508.6723 508.7910 508.9096 509.0282 509.1468 509.2654 509.3840 509.5025 509.6210 509.7395 509.8580 509.9765 --------------------------------------------------------------------------------------- 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 --------------------------------------------------------------------------------------- SQ 20241006T013041.8Z 9045.062984 197 1 50.50 2 43.41 159.12 73.31 133.81 0.00 1 0.00 1 0 15 2.501e-03 2.477e-03 1.438e-03 1.428e-03 3.909323e+03 9.327259e+01 2.668824e+02 -6.000000e+00 1.413088e+02 288.52 -2.00 20.00 -2.00 505.29 0.00e+00 3.603656e-04 3.974058e-05 2.574669e-05 -6.000000e+00 7.087626e-05 288.52 -2.00 20.00 -2.00 500.23 0.00e+00 -3.437368e-02 2.829867e-02 2.684038e-02 -6.000000e+00 4.887673e-02 227.45 -2.00 3.41 -2.00 501.87 0.00e+00 1.375 0.000 489.00 513.00 1 7.6634e+00 2.1724e-03 3.0541e-03 -6.0000e+00 3.7607e-03 6.0385e-03 5.5519e-04 6.4105e-04 -6.0000e+00 9.6517e-04 7.5068e-04 2.5906e-04 1.3585e-04 2.4557e-05 -6.0000e+00 2.4681e-04 9.7755e-03 3.4601e-05 10 0 0 10 0 0 99.07 21.43 1.54 0.00978 0.01157 32 6 1 4 0 1578.030 25.07 20.00 21.43 32.56 10.6 994.82 70 -583 -319 13 83 137 283 253 333 242 27 308 -189 -154 -573 68 -53 -23 131 -115 302 101 402 46 186 262 102 260 73 -169 -136 238 -190 257 -10 -303 -361 219 -239 -112 111 -242 -47 221 332 -182 28 -175 -126 50 -364 353 -593 -396 15 -245 -356 -288 -445 -2 97 153 274 86 326 403 108 426 442 73 -30 -27 -286 2 -17 -265 -275 -372 -292 -470 -24 -262 -223 -173 95 -535 -266 -55 -51 162 320 185 154 73 155 149 -294 14 -306 -50 -554 -378 154 -107 -197 -434 -38 -141 -52 111 -190 181 -36 166 175 273 172 239 224 -63 253 -252 334 400 300 -15 154 142 282 98 -162 -224 131 378 -216 200 323 -84 26 400 -66 59 303 -177 105 -97 -109 -380 51 -239 -378 66 -558 -314 13 87 145 306 277 364 262 29 324 -196 -158 -585 69 -55 -25 142 -126 332 111 442 51 206 290 113 284 78 -175 -138 240 -196 273 -11 -337 -400 238 -258 -120 118 -259 -50 241 365 -202 31 -190 -134 52 -369 346 -580 -388 15 -253 -354 -290 -448 -2 97 154 276 86 326 396 105 406 444 73 -30 -25 -268 1 -16 -254 -259 -344 -268 -429 -22 -241 -207 -162 90 -513 -255 -53 -49 155 309 181 153 72 153 144 -278 13 -276 -44 -488 -337 141 -102 -194 -431 -37 -136 -50 106 -184 178 -36 167 176 277 176 244 228 -63 249 -246 321 385 291 -15 150 138 273 94 -155 -212 121 346 -199 185 305 -81 26 395 -66 59 287 -169 100 -92 -101 -368 48 -225 -360 SQ 20241006T013508.2Z 9045.066068 200 1 51.49 2 43.10 160.67 72.69 134.65 0.00 1 0.00 1 0 15 3.037e-03 3.011e-03 1.743e-03 1.732e-03 4.095239e+03 1.107875e+02 2.802777e+02 -6.000000e+00 1.715897e+02 288.52 -2.00 20.00 -2.00 505.19 0.00e+00 3.877167e-04 4.792950e-05 2.697680e-05 -6.000000e+00 8.606476e-05 288.52 -2.00 20.00 -2.00 500.20 0.00e+00 -8.363187e-03 3.405084e-02 2.713491e-02 -6.000000e+00 5.935370e-02 227.45 -2.00 3.41 -2.00 501.82 0.00e+00 1.368 0.000 489.00 513.00 1 7.8832e+00 2.6130e-03 3.0862e-03 -6.0000e+00 4.5668e-03 4.8151e-03 6.7139e-04 6.4943e-04 -6.0000e+00 1.1721e-03 6.0150e-04 3.7873e-04 1.6452e-04 2.9515e-05 -6.0000e+00 2.9970e-04 1.0085e-02 2.8687e-05 10 0 0 10 0 0 99.52 21.44 1.30 0.01009 0.01271 19 3 1 4 0 2709.190 25.07 20.01 21.44 32.73 10.6 996.03 -493 -435 -223 -394 240 377 262 169 75 -193 -19 104 58 -405 -291 -8 -191 28 50 -77 271 256 367 1073 144 222 -73 276 -36 201 -174 -112 -101 9 -91 -29 -8 274 -416 197 140 -32 -312 -342 438 301 105 240 -139 -425 -424 -242 -510 -143 -123 -422 -479 -145 -295 302 -50 287 239 -79 346 453 124 469 327 -7 429 345 -261 -1 -120 -222 -515 -410 114 16 -194 -644 -18 -272 57 -481 -78 157 -132 -111 394 -246 183 -65 563 200 64 -191 -249 -375 -547 -163 253 -104 -150 -654 -343 -42 413 78 -90 44 316 572 222 20 123 -12 10 54 29 118 128 278 14 -28 355 103 283 652 182 -574 -160 462 -10 739 -39 -50 -137 174 99 -148 491 -376 -4 -181 -180 -675 -82 134 -505 -482 -430 -227 -415 259 415 295 192 85 -217 -21 113 62 -431 -308 -8 -204 30 53 -82 289 272 390 1137 155 239 -78 291 -37 201 -171 -110 -101 9 -97 -31 -9 288 -435 203 144 -33 -328 -363 466 322 112 252 -144 -431 -417 -232 -482 -136 -119 -423 -493 -151 -307 312 -52 297 248 -82 357 460 124 462 319 -7 409 330 -254 -1 -119 -219 -503 -393 107 15 -183 -615 -17 -264 55 -477 -78 155 -130 -110 392 -249 187 -66 569 198 62 -182 -231 -342 -497 -150 239 -102 -152 -672 -348 -42 406 77 -91 42 304 554 217 20 122 -12 10 52 27 111 120 259 13 -26 333 97 264 606 169 -529 -145 409 -9 659 -36 -46 -130 167 96 -144 478 -370 -4 -177 -173 -636 -76 122 -466 SQ 20241006T015011.3Z 9045.076519 212 1 50.52 2 42.23 166.05 70.67 137.56 0.00 1 0.00 1 0 15 3.193e-03 3.267e-03 1.610e-03 1.602e-03 3.935908e+03 1.027769e+02 2.664526e+02 -6.000000e+00 1.803845e+02 288.52 -2.00 20.00 -2.00 505.24 0.00e+00 4.508288e-04 4.442753e-05 2.885477e-05 -6.000000e+00 9.047488e-05 288.52 -2.00 20.00 -2.00 500.27 0.00e+00 -5.766385e-02 3.162805e-02 2.741984e-02 -6.000000e+00 6.238861e-02 227.45 -2.00 3.41 -2.00 501.87 0.00e+00 1.349 0.000 489.00 513.00 1 7.8005e+00 2.4364e-03 3.0788e-03 -6.0000e+00 4.8004e-03 5.3983e-03 6.2039e-04 6.4973e-04 -6.0000e+00 1.2319e-03 6.5714e-04 5.1740e-04 1.5280e-04 2.3992e-05 -6.0000e+00 3.1505e-04 6.8852e-03 5.6984e-05 10 0 0 10 0 0 95.72 21.46 1.49 0.00689 0.01047 13 2 1 4 0 3885.080 25.07 20.01 21.46 33.29 10.6 1000.62 -403 -582 -68 -368 83 75 165 -33 8 205 -53 240 -235 -462 -419 194 -51 -150 -20 -157 83 263 776 2178 34 386 115 -51 -86 -150 -46 14 -157 -207 -516 11 -363 279 -359 66 -183 311 -297 34 396 38 126 232 -337 -101 -413 -129 -96 -432 -243 -98 -596 -185 288 387 -124 188 327 137 213 61 345 384 180 16 242 72 -348 478 -26 -386 -298 -121 38 28 -158 -475 8 -372 79 -343 -531 7 -205 34 108 19 241 70 343 140 17 -319 -295 -232 -188 -198 252 13 -117 -525 -85 -106 229 -55 -181 544 213 301 -140 126 -53 426 -128 -74 72 -39 311 460 352 75 109 215 -56 278 125 -333 188 333 96 138 287 100 -33 183 -127 -278 147 -306 -363 -106 -340 -520 77 168 -215 -399 -585 -70 -394 85 78 176 -35 8 217 -55 247 -239 -466 -419 193 -51 -155 -21 -169 90 284 832 2311 37 418 124 -55 -90 -152 -46 13 -159 -217 -558 12 -394 297 -380 69 -192 325 -316 36 426 42 136 247 -354 -104 -411 -125 -92 -416 -239 -99 -622 -195 302 404 -130 198 344 144 223 63 350 384 178 16 234 70 -342 474 -26 -388 -294 -117 36 26 -151 -459 8 -365 78 -344 -535 7 -206 34 102 18 234 68 330 133 16 -307 -278 -214 -173 -185 227 12 -112 -512 -82 -100 214 -52 -173 524 208 296 -140 125 -53 428 -128 -73 70 -38 294 435 335 72 104 205 -53 263 118 -310 172 300 87 125 266 95 -32 178 -125 -275 146 -306 -362 -105 -331 -496 72 154 -201