File name: Pandora108s1_Toronto-West_20221206_L2Fit_fwvt1c6d20211104p1-8.txt File generation date: 20240820T140611.5Z Data description: Level 2 spectral fitting results file Data file version: fwvt1c6p1-8 Local principal investigator: Vitali Fioletov Network principal investigator: Alexander Cede Instrument type: Pandora Instrument number: 108 Spectrometer number: 1 Processing software version used: BlickP v1.8.59 Instrument operation file used: Pandora108_OF_v3d20210713.txt Instrument calibration file used: Pandora108s1_CF_v6d20211104.txt Level 1 file used: Pandora108s1_Toronto-West_20221206_L1_smca1c6d20211104p1-8.txt Full location name: Ontario Ministry of the Environment, Toronto West (Etobicoke) Short location name: Toronto-West Country of location: Canada Location latitude [deg]: 43.7094 Location longitude [deg]: -79.5435 Location altitude [m]: 141 Local noon date: 20221206 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: None First and last pixel inside fitting window: 1653 1802 Nominal wavelengths inside fitting window [nm]: 492.1141 492.2350 492.3559 492.4768 492.5976 492.7185 492.8393 492.9601 493.0809 493.2016 493.3224 493.4431 493.5638 493.6845 493.8052 493.9259 494.0466 494.1672 494.2878 494.4084 494.5290 494.6496 494.7701 494.8907 495.0112 495.1317 495.2522 495.3726 495.4931 495.6135 495.7339 495.8543 495.9747 496.0951 496.2154 496.3358 496.4561 496.5764 496.6967 496.8169 496.9372 497.0574 497.1776 497.2978 497.4180 497.5382 497.6583 497.7785 497.8986 498.0187 498.1388 498.2588 498.3789 498.4989 498.6189 498.7389 498.8589 498.9789 499.0988 499.2187 499.3387 499.4585 499.5784 499.6983 499.8181 499.9380 500.0578 500.1776 500.2973 500.4171 500.5368 500.6566 500.7763 500.8960 501.0156 501.1353 501.2549 501.3746 501.4942 501.6138 501.7333 501.8529 501.9724 502.0920 502.2115 502.3310 502.4504 502.5699 502.6893 502.8087 502.9281 503.0475 503.1669 503.2862 503.4056 503.5249 503.6442 503.7635 503.8827 504.0020 504.1212 504.2404 504.3596 504.4788 504.5980 504.7171 504.8363 504.9554 505.0745 505.1935 505.3126 505.4316 505.5507 505.6697 505.7887 505.9076 506.0266 506.1455 506.2645 506.3834 506.5023 506.6211 506.7400 506.8588 506.9776 507.0964 507.2152 507.3340 507.4527 507.5715 507.6902 507.8089 507.9276 508.0462 508.1649 508.2835 508.4021 508.5207 508.6393 508.7578 508.8764 508.9949 509.1134 509.2319 509.3504 509.4688 509.5873 509.7057 509.8241 509.9425 --------------------------------------------------------------------------------------- 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 20221206T181924.0Z 8375.763472 427 1 50.52 2 68.19 197.43 108.56 33.81 0.00 1 0.00 1 0 15 2.362e-03 2.362e-03 1.075e-03 1.069e-03 3.143597e+03 6.430270e+01 2.287628e+02 -6.000000e+00 1.274721e+02 275.08 -2.00 20.00 -2.00 505.23 0.00e+00 7.164851e-04 2.886913e-05 4.395372e-05 -6.000000e+00 6.507382e-05 275.08 -2.00 20.00 -2.00 500.26 0.00e+00 1.916729e-01 2.142625e-02 2.484775e-02 -6.000000e+00 4.674865e-02 222.33 -2.00 3.57 -2.00 501.86 0.00e+00 2.668 0.000 489.00 513.00 1 7.1458e+00 1.6432e-03 4.4183e-03 -6.0000e+00 3.5760e-03 8.9664e-03 4.2335e-04 6.2173e-04 -6.0000e+00 9.2942e-04 1.0335e-03 7.4715e-04 1.0556e-04 3.2105e-05 -6.0000e+00 2.4061e-04 3.9646e-03 -2.7151e-05 10 0 0 10 0 0 99.76 18.18 0.87 0.00396 0.00915 21 4 1 4 0 2405.330 23.64 15.00 18.18 11.83 26.7 1016.23 -501 -324 114 149 37 408 450 333 536 97 158 -36 -133 -263 -231 -336 -69 44 139 -152 68 163 95 288 298 182 119 47 -278 -339 -407 -152 -237 176 125 -73 -1 -98 -108 -226 -70 135 255 191 361 498 78 367 -28 -179 -458 -302 -553 -305 -85 -96 -39 86 213 -63 245 99 158 130 82 -97 -81 18 -255 -316 -250 -228 -93 -32 -26 -174 -318 -245 -97 -373 -256 -338 -215 -319 61 -52 -259 124 -132 -128 -75 242 344 214 202 22 -72 -124 -291 -451 -324 -559 -38 59 -310 -164 7 -69 -126 67 -75 52 65 357 337 534 440 288 280 285 207 188 226 42 140 318 91 179 80 107 45 29 -27 -137 -59 140 -101 50 66 110 53 155 295 -95 166 -93 -85 -409 -355 -289 -470 -311 112 151 38 432 483 360 573 103 165 -37 -135 -265 -232 -341 -71 47 149 -165 74 177 104 315 326 199 129 50 -289 -343 -408 -154 -247 189 136 -80 -1 -105 -114 -239 -74 143 274 207 394 542 85 390 -30 -182 -453 -293 -536 -300 -86 -100 -41 91 226 -67 259 97 156 128 87 -100 -83 18 -254 -312 -245 -225 -93 -33 -26 -174 -314 -238 -93 -360 -248 -331 -212 -318 61 -53 -243 116 -124 -121 -72 232 332 206 192 20 -71 -120 -274 -421 -303 -537 -35 56 -299 -157 7 -64 -118 63 -72 50 64 351 333 530 438 287 276 277 197 178 213 40 133 302 87 170 76 101 42 27 -24 -123 -53 128 -95 48 64 108 52 153 292 -94 165 -91 -81 -384 -329 -268 SQ 20221206T182404.6Z 8375.766720 430 1 49.43 2 68.45 198.55 108.09 34.80 0.00 1 0.00 1 0 15 2.695e-03 2.680e-03 1.411e-03 1.406e-03 3.722634e+03 8.612989e+01 2.670390e+02 -6.000000e+00 1.454399e+02 275.07 -2.00 20.00 -2.00 505.27 0.00e+00 7.105034e-04 3.793644e-05 4.473855e-05 -6.000000e+00 7.424654e-05 275.07 -2.00 20.00 -2.00 500.38 0.00e+00 2.038115e-01 2.792267e-02 2.942246e-02 -6.000000e+00 5.333912e-02 222.33 -2.00 3.57 -2.00 501.98 0.00e+00 2.698 0.000 489.00 513.00 1 7.5570e+00 2.1342e-03 4.6529e-03 -6.0000e+00 4.0801e-03 9.6458e-03 5.5033e-04 7.2455e-04 -6.0000e+00 1.0605e-03 6.8088e-04 6.8227e-04 1.3836e-04 3.7053e-05 -6.0000e+00 2.7453e-04 3.9002e-03 -2.8011e-05 10 0 0 10 0 0 98.94 18.19 0.90 0.00390 0.00863 18 4 1 4 0 2745.170 23.64 15.01 18.19 11.78 26.7 1015.77 -263 -561 145 -92 338 400 488 106 468 56 -10 -16 -369 -70 -391 -119 -121 271 -56 -69 175 37 242 193 466 -57 428 -203 -55 -470 -357 -233 -333 268 51 1 -28 70 -282 -12 -34 176 223 100 441 334 353 386 73 -166 -335 -174 -602 -137 -278 9 -213 222 -37 270 104 290 83 13 -79 -243 43 -201 20 -444 42 -305 -143 -90 -114 -240 -532 14 -319 -285 -502 -76 -372 -169 -22 -138 35 -139 5 -280 376 -55 391 200 226 32 -38 -143 -413 -303 -794 -248 -339 239 -410 67 -16 -227 -88 -271 102 40 331 337 577 154 588 465 166 558 25 133 -110 231 -42 384 47 219 -68 155 129 -97 224 -233 106 -33 96 173 157 44 63 258 -7 187 -83 -34 -210 -320 -568 -354 -253 -551 145 -94 352 423 522 113 499 59 -10 -16 -377 -71 -396 -121 -125 284 -60 -74 187 40 259 207 500 -61 454 -213 -56 -474 -357 -235 -343 282 54 1 -30 74 -293 -12 -35 183 234 106 467 353 371 401 74 -166 -328 -168 -580 -133 -277 9 -219 228 -39 277 107 298 85 13 -80 -245 43 -198 19 -431 40 -296 -139 -89 -113 -235 -514 14 -302 -269 -478 -72 -358 -163 -22 -135 34 -136 4 -275 370 -55 392 200 223 31 -36 -133 -379 -275 -726 -230 -324 233 -407 66 -16 -220 -85 -264 101 40 330 338 581 156 595 471 167 552 25 129 -107 224 -41 373 46 212 -66 150 123 -91 208 -216 99 -32 92 168 154 44 63 256 -7 187 -83 -33 -204 -305 -536 -335 SQ 20221206T183631.7Z 8375.775367 439 1 50.41 2 69.22 201.52 106.79 37.39 0.00 1 0.00 1 0 15 2.630e-03 2.616e-03 1.304e-03 1.300e-03 3.322630e+03 7.886068e+01 2.434413e+02 -6.000000e+00 1.418951e+02 275.07 -2.00 20.00 -2.00 505.27 0.00e+00 6.609964e-04 3.506536e-05 4.104228e-05 -6.000000e+00 7.243664e-05 275.07 -2.00 20.00 -2.00 500.38 0.00e+00 2.273708e-01 2.580481e-02 2.560559e-02 -6.000000e+00 5.203768e-02 222.33 -2.00 3.57 -2.00 501.98 0.00e+00 2.791 0.000 489.00 513.00 1 7.4283e+00 1.9724e-03 4.5878e-03 -6.0000e+00 3.9806e-03 1.1096e-02 5.0905e-04 6.4906e-04 -6.0000e+00 1.0346e-03 7.6410e-04 7.5768e-04 1.2819e-04 3.4261e-05 -6.0000e+00 2.6783e-04 3.1520e-03 -2.1167e-05 10 0 0 10 0 0 98.46 18.16 0.81 0.00315 0.00841 19 4 1 4 0 2652.430 23.64 14.98 18.16 11.95 26.8 1017.02 -725 -332 -141 120 163 434 587 318 503 -107 137 -212 -126 -216 -174 -311 2 83 69 154 -32 236 29 411 147 309 148 -38 -179 -531 -209 -480 -73 -17 242 -221 18 102 53 -130 -66 300 98 310 274 584 -107 408 52 -78 -316 -515 -440 -515 -63 -278 45 22 84 99 107 246 -36 226 -83 -35 -112 14 -235 -217 -53 -340 -106 -98 226 -174 -231 -390 -47 -482 -135 -207 -180 53 -133 -81 -122 114 -102 -19 176 156 406 49 223 27 147 -438 30 -587 -405 -592 -189 150 -314 -25 -219 -206 -403 80 -52 182 178 268 507 231 627 217 211 88 355 -219 290 124 139 421 137 289 -148 225 -150 -95 -78 41 -229 109 46 -17 255 -75 269 -120 281 37 124 -99 -199 -301 -570 -73 -700 -326 -142 124 170 460 627 341 536 -113 143 -219 -128 -219 -176 -317 2 87 73 165 -34 253 31 440 158 331 157 -40 -184 -535 -209 -484 -75 -18 258 -236 19 107 55 -135 -69 312 103 328 290 617 -112 422 53 -78 -310 -497 -423 -501 -63 -283 46 23 86 101 110 253 -37 233 -85 -35 -112 14 -229 -210 -51 -329 -104 -96 223 -171 -223 -372 -44 -457 -128 -198 -173 51 -129 -79 -120 112 -100 -19 174 156 407 49 220 26 140 -410 28 -535 -369 -552 -181 147 -312 -25 -215 -200 -390 77 -52 181 178 269 511 234 634 220 212 87 347 -213 281 120 136 409 133 280 -143 217 -144 -90 -73 38 -213 103 44 -17 249 -74 266 -119 280 37 123 -98 -193 -287 -538 -69