├── .github └── workflows │ └── testing_in_conda.yml ├── .gitignore ├── .travis.yml ├── 1-magnetic-inversion-raglan-reproduce.ipynb ├── 2-linear-inversion-app.ipynb ├── 3-magnetic-inversion-raglan-lp.ipynb ├── LICENSE ├── README.md ├── data └── Raglan_1997 │ ├── inv.inp │ ├── maginv3d.sus │ ├── mesh.msh │ └── obs.mag ├── environment.yml ├── environment_for_test.yml ├── figures ├── mag3disco.gif ├── magobs.gif ├── map.gif └── osmand-section.jpg ├── pdfs ├── inversion-for-geologists-part1.pdf └── inversion-for-geologists-part2.pdf ├── requirements_dev.txt └── tests └── test_notebooks.py /.github/workflows/testing_in_conda.yml: -------------------------------------------------------------------------------- 1 | name: TestingInConda 2 | 3 | on: 4 | push: 5 | branches: 6 | - '*' 7 | pull_request: 8 | branches: 9 | - '*' 10 | jobs: 11 | setup-build: 12 | name: Ex1 (${{ matrix.python-version }}, ${{ matrix.os }}) 13 | runs-on: ${{ matrix.os }} 14 | strategy: 15 | fail-fast: false 16 | matrix: 17 | os: ["ubuntu-latest"] 18 | python-version: [3.7, 3.8] 19 | 20 | steps: 21 | - uses: actions/checkout@v2 22 | - name: Setup Conda 23 | uses: s-weigand/setup-conda@v1 24 | with: 25 | update-conda: true 26 | conda-channels: conda-forge 27 | python-version: ${{ matrix.python-version }} 28 | 29 | - name: Install Env 30 | shell: bash 31 | run: | 32 | python --version 33 | conda env create -f environment_for_test.yml 34 | source activate t21-tue-simpeg 35 | pip install -r requirements_dev.txt 36 | 37 | - name: Run Tests 38 | shell: bash 39 | run: | 40 | source activate t21-tue-simpeg 41 | conda list 42 | pytest . 43 | -------------------------------------------------------------------------------- /.gitignore: -------------------------------------------------------------------------------- 1 | # Byte-compiled / optimized / DLL files 2 | __pycache__/ 3 | *.py[cod] 4 | *$py.class 5 | 6 | # C extensions 7 | *.so 8 | 9 | # Distribution / packaging 10 | .Python 11 | build/ 12 | develop-eggs/ 13 | dist/ 14 | downloads/ 15 | eggs/ 16 | .eggs/ 17 | lib/ 18 | lib64/ 19 | parts/ 20 | sdist/ 21 | var/ 22 | wheels/ 23 | *.egg-info/ 24 | .installed.cfg 25 | *.egg 26 | MANIFEST 27 | 28 | # PyInstaller 29 | # Usually these files are written by a python script from a template 30 | # before PyInstaller builds the exe, so as to inject date/other infos into it. 31 | *.manifest 32 | *.spec 33 | 34 | # Installer logs 35 | pip-log.txt 36 | pip-delete-this-directory.txt 37 | 38 | # Unit test / coverage reports 39 | htmlcov/ 40 | .tox/ 41 | .coverage 42 | .coverage.* 43 | .cache 44 | nosetests.xml 45 | coverage.xml 46 | *.cover 47 | .hypothesis/ 48 | .pytest_cache/ 49 | 50 | # Translations 51 | *.mo 52 | *.pot 53 | 54 | # Django stuff: 55 | *.log 56 | local_settings.py 57 | db.sqlite3 58 | 59 | # Flask stuff: 60 | instance/ 61 | .webassets-cache 62 | 63 | # Scrapy stuff: 64 | .scrapy 65 | 66 | # Sphinx documentation 67 | docs/_build/ 68 | 69 | # PyBuilder 70 | target/ 71 | 72 | # Jupyter Notebook 73 | .ipynb_checkpoints 74 | 75 | # pyenv 76 | .python-version 77 | 78 | # celery beat schedule file 79 | celerybeat-schedule 80 | 81 | # SageMath parsed files 82 | *.sage.py 83 | 84 | # Environments 85 | .env 86 | .venv 87 | env/ 88 | venv/ 89 | ENV/ 90 | env.bak/ 91 | venv.bak/ 92 | 93 | # Spyder project settings 94 | .spyderproject 95 | .spyproject 96 | 97 | # Rope project settings 98 | .ropeproject 99 | 100 | # mkdocs documentation 101 | /site 102 | 103 | # mypy 104 | .mypy_cache/ 105 | 106 | # repo specific 107 | *.plot 108 | *.gmt 109 | *.ps 110 | -------------------------------------------------------------------------------- /.travis.yml: -------------------------------------------------------------------------------- 1 | language: python 2 | python: 3 | - 3.6 4 | 5 | env: 6 | matrix: 7 | - TEST=tests/test_notebooks.py 8 | 9 | before_install: 10 | - wget http://repo.continuum.io/miniconda/Miniconda3-latest-Linux-x86_64.sh -O miniconda.sh; 11 | - chmod +x miniconda.sh 12 | - ./miniconda.sh -b -p $HOME/miniconda 13 | - export PATH=$HOME/anaconda/bin:$HOME/miniconda/bin:$PATH 14 | - conda update --yes conda 15 | 16 | install: 17 | - conda env create -f environment.yml 18 | - source activate t21-tue-simpeg 19 | - pip install -r requirements_dev.txt 20 | 21 | script: 22 | - travis_wait 45 pytest $TEST -v -s 23 | -------------------------------------------------------------------------------- /2-linear-inversion-app.ipynb: -------------------------------------------------------------------------------- 1 | { 2 | "cells": [ 3 | { 4 | "cell_type": "code", 5 | "execution_count": 35, 6 | "metadata": {}, 7 | "outputs": [], 8 | "source": [ 9 | "from geoscilabs.inversion.LinearInversionDirect import LinearInversionDirectApp\n", 10 | "from ipywidgets import interact, FloatSlider, ToggleButtons, IntSlider, FloatText, IntText" 11 | ] 12 | }, 13 | { 14 | "cell_type": "code", 15 | "execution_count": 36, 16 | "metadata": { 17 | "iooxa": { 18 | "id": { 19 | "block": "XoBdBeGBCevcAMVqqZtM", 20 | "project": "VNMrkxzChhdveZyf6lmb", 21 | "version": 1 22 | }, 23 | "outputId": { 24 | "block": "nzHSDMZI4foQGe88I1D5", 25 | "project": "VNMrkxzChhdveZyf6lmb", 26 | "version": 1 27 | } 28 | } 29 | }, 30 | "outputs": [], 31 | "source": [ 32 | "app = LinearInversionDirectApp()" 33 | ] 34 | }, 35 | { 36 | "cell_type": "markdown", 37 | "metadata": { 38 | "iooxa": { 39 | "id": { 40 | "block": "GAKaomYwcnkDpl2fpbAI", 41 | "project": "VNMrkxzChhdveZyf6lmb", 42 | "version": 1 43 | } 44 | } 45 | }, 46 | "source": [ 47 | "# Linear Inversion App\n", 48 | "\n", 49 | "This app is based upon the inversion tutorial: \"INVERSION FOR APPLIED GEOPHYSICS\" by Oldenburg and Li (2005).\n", 50 | "\n", 51 | "Douglas W. Oldenburg and Yaoguo Li (2005) 5. Inversion for Applied Geophysics: A Tutorial. Near-Surface Geophysics: pp. 89-150.\n", 52 | "eISBN: 978-1-56080-171-9 \n", 53 | "print ISBN: 978-1-56080-130-6 \n", 54 | "https://doi.org/10.1190/1.9781560801719.ch5 \n" 55 | ] 56 | }, 57 | { 58 | "cell_type": "markdown", 59 | "metadata": { 60 | "iooxa": { 61 | "id": { 62 | "block": "wzGxtmJjFucP0mYtZzVD", 63 | "project": "VNMrkxzChhdveZyf6lmb", 64 | "version": 2 65 | } 66 | } 67 | }, 68 | "source": [ 69 | "## Purpose\n", 70 | "\n", 71 | "By using a simple decaying and oscillating kernel function, which emulates the physics of electromagnetic (EM) survey, we understand basic concepts of inverting data. Three items that we are going to explore are:\n", 72 | "\n", 73 | "- Step1: Create a model ($\\mathbf{m}$)\n", 74 | "- Step2: Generate a sensitivity kernel (or matrix), $\\mathbf{G}$\n", 75 | "- Step3: Simulate data ($\\mathbf{d} = \\mathbf{G}\\mathbf{m}$)\n", 76 | "- Step4: All three steps together\n", 77 | "- Step5: Invert the data, and explore inversion results" 78 | ] 79 | }, 80 | { 81 | "cell_type": "markdown", 82 | "metadata": { 83 | "iooxa": { 84 | "id": { 85 | "block": "TVYiH721oymtWnC2zQ1w", 86 | "project": "VNMrkxzChhdveZyf6lmb", 87 | "version": 2 88 | } 89 | } 90 | }, 91 | "source": [ 92 | "## Forward problem\n", 93 | "\n", 94 | "\n", 95 | "Let $g_j(x)$ denote the kernel function for $j$th datum. With a given model $m(x)$, $j$th datum can be computed by solving following integral equation:\n", 96 | "\n", 97 | " $$ d_j = \\int_a^{b} g_j(x) m(x) dx $$\n", 98 | "\n", 99 | "where \n", 100 | "\n", 101 | "$$ g_j(x) = e^{p_jx} cos (2 \\pi q_jx) $$ \n", 102 | "\n", 103 | "By discretizing $g_j(x)$ we obtain" 104 | ] 105 | }, 106 | { 107 | "cell_type": "markdown", 108 | "metadata": { 109 | "iooxa": { 110 | "id": { 111 | "block": "s9H8M9Mtzk7wkXBPF8iZ", 112 | "project": "VNMrkxzChhdveZyf6lmb", 113 | "version": 2 114 | } 115 | } 116 | }, 117 | "source": [ 118 | "$$ \\mathbf{g}_j(\\mathbf{x}) = e^{p_j\\mathbf{x}} cos (2 \\pi q_j \\mathbf{x})$$\n", 119 | "\n", 120 | "where\n", 121 | "\n", 122 | "- $\\mathbf{g}_j$: $j$th row vector for the sensitivty matrix ($1 \\times M$)\n", 123 | "- $\\mathbf{x}$: model location ($1 \\times M$)\n", 124 | "- $p_j$: decaying constant (<0)\n", 125 | "- $q_j$: oscillating constant (>0)\n", 126 | "\n", 127 | "By stacking multiple rows of $\\mathbf{g}_j$, we obtain sensitivity matrix, $\\mathbf{G}$: \n", 128 | "\n", 129 | "\\begin{align}\n", 130 | " \\mathbf{G} = \n", 131 | " \\begin{bmatrix}\n", 132 | " \\mathbf{g}_1\\\\\n", 133 | " \\vdots\\\\\n", 134 | " \\mathbf{g}_{N}\n", 135 | " \\end{bmatrix}\n", 136 | "\\end{align}\n", 137 | "\n", 138 | "Here, the size of the matrix $\\mathbf{G}$ is $(N \\times M)$. \n", 139 | "Finally data, $\\mathbf{d}$, can be written as a linear equation:\n", 140 | "\n", 141 | "$$ \\mathbf{d} = \\mathbf{G}\\mathbf{m}$$\n", 142 | "\n", 143 | "where $\\mathbf{m}$ is an inversion model; this is a column vector ($M \\times 1$). \n", 144 | "\n", 145 | "In real measurments, there will be various noises source, and hence observation, $\\mathbf{d}^{obs}$, can be written as \n", 146 | "\n", 147 | "$$ \\mathbf{d}^{obs} = \\mathbf{G}\\mathbf{m} + \\mathbf{noise}$$" 148 | ] 149 | }, 150 | { 151 | "cell_type": "markdown", 152 | "metadata": { 153 | "iooxa": { 154 | "id": { 155 | "block": "UY71pbGgRaPr7PBxH3OT", 156 | "project": "VNMrkxzChhdveZyf6lmb", 157 | "version": 2 158 | } 159 | } 160 | }, 161 | "source": [ 162 | "## Step1: Create a model, $\\mathbf{m}$\n", 163 | "\n", 164 | "The model $m$ is a function defined on the interval (-2,2). Here we generate a model that is the sum of a: (a) background $m_{ref}$, (b) box car $m_1$ and (c) Gaussian $m_2$. The box car is defined by\n", 165 | "- `m$_{background}$` : amplitude of the background\n", 166 | "- `m1` : amplitude\n", 167 | "- `$m1_{center}$` : center\n", 168 | "- `$m1_{width}$` : width\n", 169 | "the Gaussian is defined by \n", 170 | "- `m2` : amplitude\n", 171 | "- `$m2_{center}$` : center\n", 172 | "- `$m2_{sigma}$` : width of Gaussian (as defined by a standard deviation)\n", 173 | "- `M`: # of model parameters" 174 | ] 175 | }, 176 | { 177 | "cell_type": "code", 178 | "execution_count": 38, 179 | "metadata": { 180 | "iooxa": { 181 | "id": { 182 | "block": "BoaOYiDvY7yFMo9QvOpF", 183 | "project": "FesuZy9nERkRH9r67vBe", 184 | "version": 1 185 | }, 186 | "outputId": { 187 | "block": "vuewWXXAWRcE2KNdwnEk", 188 | "project": "FesuZy9nERkRH9r67vBe", 189 | "version": 1 190 | } 191 | } 192 | }, 193 | "outputs": [ 194 | { 195 | "data": { 196 | "application/vnd.jupyter.widget-view+json": { 197 | "model_id": "b8ec921437494d54ab190cabc88af603", 198 | "version_major": 2, 199 | "version_minor": 0 200 | }, 201 | "text/plain": [ 202 | "interactive(children=(FloatSlider(value=0.0, continuous_update=False, description='m$_{background}$', max=2.0,…" 203 | ] 204 | }, 205 | "metadata": {}, 206 | "output_type": "display_data" 207 | } 208 | ], 209 | "source": [ 210 | "app.reset_to_defaults()\n", 211 | "Q_model = app.interact_plot_model()" 212 | ] 213 | }, 214 | { 215 | "cell_type": "markdown", 216 | "metadata": { 217 | "iooxa": { 218 | "id": { 219 | "block": "0zUv3gNqulptgg5RQABq", 220 | "project": "VNMrkxzChhdveZyf6lmb", 221 | "version": 4 222 | } 223 | } 224 | }, 225 | "source": [ 226 | "## Step2: Generate a sensitivity kernel (or matrix), $\\mathbf{G}$\n", 227 | "\n", 228 | "By using the following app, we explore each row vector of the sensitivity matrix, $\\mathbf{g}_j$. Parameters of the apps are:\n", 229 | "\n", 230 | "- `M`: # of model parameters\n", 231 | "- `N`: # of data\n", 232 | "- `p`: decaying constant (<0)\n", 233 | "- `q`: oscillating constant (>0)\n", 234 | "- `ymin`: maximum limit for y-axis\n", 235 | "- `ymax`: minimum limit for y-axis\n", 236 | "- `show_singular`: show singualr values" 237 | ] 238 | }, 239 | { 240 | "cell_type": "code", 241 | "execution_count": 39, 242 | "metadata": { 243 | "iooxa": { 244 | "id": { 245 | "block": "x9PxRUeZbvYukM3fIzi0", 246 | "project": "VNMrkxzChhdveZyf6lmb", 247 | "version": 27 248 | }, 249 | "outputId": { 250 | "block": "xjP1IRGiznblm7UECCWU", 251 | "project": "VNMrkxzChhdveZyf6lmb", 252 | "version": 27 253 | } 254 | } 255 | }, 256 | "outputs": [ 257 | { 258 | "data": { 259 | "application/vnd.jupyter.widget-view+json": { 260 | "model_id": "7912e3d288204167ae717f5fa298c456", 261 | "version_major": 2, 262 | "version_minor": 0 263 | }, 264 | "text/plain": [ 265 | "interactive(children=(IntSlider(value=20, continuous_update=False, description='N', min=1), IntSlider(value=10…" 266 | ] 267 | }, 268 | "metadata": {}, 269 | "output_type": "display_data" 270 | } 271 | ], 272 | "source": [ 273 | "app.reset_to_defaults()\n", 274 | "Q_kernel = app.interact_plot_G()" 275 | ] 276 | }, 277 | { 278 | "cell_type": "markdown", 279 | "metadata": { 280 | "iooxa": { 281 | "id": { 282 | "block": "B3aFbgUWKTzmqttCPDu3", 283 | "project": "VNMrkxzChhdveZyf6lmb", 284 | "version": 1 285 | } 286 | } 287 | }, 288 | "source": [ 289 | "## Step3: Simulate data" 290 | ] 291 | }, 292 | { 293 | "cell_type": "markdown", 294 | "metadata": { 295 | "iooxa": { 296 | "id": { 297 | "block": "fOWoINnfZAyspfxAKRnr", 298 | "project": "VNMrkxzChhdveZyf6lmb", 299 | "version": 3 300 | } 301 | } 302 | }, 303 | "source": [ 304 | "The $j$-th datum is the inner product of the $j$-th kernel $g_j(x)$ and the model $m(x)$. In discrete form it can be written as the dot product of the vector $g_j$ and the model vector $m$." 305 | ] 306 | }, 307 | { 308 | "cell_type": "markdown", 309 | "metadata": { 310 | "iooxa": { 311 | "id": { 312 | "block": "xeKrzyj4YLEYLfraRG1y", 313 | "project": "VNMrkxzChhdveZyf6lmb", 314 | "version": 1 315 | } 316 | } 317 | }, 318 | "source": [ 319 | "### $$ d_j = \\mathbf{g}_j \\mathbf{m} $$\n", 320 | "\n", 321 | "If there are $N$ data, these data can be written as a column vector, $\\mathbf{d}$:\n", 322 | "\n", 323 | "\\begin{align}\n", 324 | " \\mathbf{d} = \\mathbf{G}\\mathbf{m} = \n", 325 | " \\begin{bmatrix}\n", 326 | " d_1\\\\\n", 327 | " \\vdots\\\\\n", 328 | " d_{N}\n", 329 | " \\end{bmatrix}\n", 330 | "\\end{align}" 331 | ] 332 | }, 333 | { 334 | "cell_type": "markdown", 335 | "metadata": { 336 | "iooxa": { 337 | "id": { 338 | "block": "k5SdPFYbjSTR0Unuprld", 339 | "project": "VNMrkxzChhdveZyf6lmb", 340 | "version": 1 341 | } 342 | } 343 | }, 344 | "source": [ 345 | "### Adding Noise\n", 346 | "\n", 347 | "Observational data are always contaminated with noise. Here we add Gaussian noise $N(0,\\epsilon)$ (zero mean and standard deviation $\\sigma$). Here we choose \n", 348 | "\n", 349 | "$$ \\epsilon = \\% |d| + \\text{floor} $$\n", 350 | " " 351 | ] 352 | }, 353 | { 354 | "cell_type": "code", 355 | "execution_count": 40, 356 | "metadata": { 357 | "iooxa": { 358 | "id": { 359 | "block": "lSiXnfA4hZvPpbkCW4st", 360 | "project": "FesuZy9nERkRH9r67vBe", 361 | "version": 1 362 | }, 363 | "outputId": { 364 | "block": "NsHNyjWDNLM2ZFx0quKY", 365 | "project": "FesuZy9nERkRH9r67vBe", 366 | "version": 1 367 | } 368 | } 369 | }, 370 | "outputs": [ 371 | { 372 | "data": { 373 | "application/vnd.jupyter.widget-view+json": { 374 | "model_id": "b666c7a20012473d9ee8a647df9ffc13", 375 | "version_major": 2, 376 | "version_minor": 0 377 | }, 378 | "text/plain": [ 379 | "interactive(children=(Checkbox(value=False, description='add_noise'), FloatText(value=0.0, description='percen…" 380 | ] 381 | }, 382 | "metadata": {}, 383 | "output_type": "display_data" 384 | } 385 | ], 386 | "source": [ 387 | "app.reset_to_defaults()\n", 388 | "Q_data = app.interact_plot_data()" 389 | ] 390 | }, 391 | { 392 | "cell_type": "markdown", 393 | "metadata": { 394 | "iooxa": { 395 | "id": { 396 | "block": "Kfas1vohMX9B1w8haaeF", 397 | "project": "VNMrkxzChhdveZyf6lmb", 398 | "version": 1 399 | } 400 | } 401 | }, 402 | "source": [ 403 | "## Step4: All three steps together" 404 | ] 405 | }, 406 | { 407 | "cell_type": "code", 408 | "execution_count": 41, 409 | "metadata": { 410 | "iooxa": { 411 | "id": { 412 | "block": "gpeg40wOdccGxZMMBIxU", 413 | "project": "VNMrkxzChhdveZyf6lmb", 414 | "version": 26 415 | }, 416 | "outputId": { 417 | "block": "yFInpwVxV9MsRBuL5pNr", 418 | "project": "VNMrkxzChhdveZyf6lmb", 419 | "version": 26 420 | } 421 | } 422 | }, 423 | "outputs": [ 424 | { 425 | "data": { 426 | "application/vnd.jupyter.widget-view+json": { 427 | "model_id": "e7467249c0c14dd7af088033b97d9a7c", 428 | "version_major": 2, 429 | "version_minor": 0 430 | }, 431 | "text/plain": [ 432 | "VBox(children=(HBox(children=(VBox(children=(Button(description='Model', layout=Layout(height='30px', width='1…" 433 | ] 434 | }, 435 | "metadata": {}, 436 | "output_type": "display_data" 437 | } 438 | ], 439 | "source": [ 440 | "app.interact_plot_all_three_together()" 441 | ] 442 | }, 443 | { 444 | "cell_type": "markdown", 445 | "metadata": { 446 | "iooxa": { 447 | "id": { 448 | "block": "I5cRNN1g9FcrGvbH7UKL", 449 | "project": "VNMrkxzChhdveZyf6lmb", 450 | "version": 1 451 | } 452 | } 453 | }, 454 | "source": [ 455 | "## Inverse Problem\n", 456 | "\n", 457 | "In the inverse problem we attempt to find the model $\\mathbf{m}$ that gave rise to the observational data $\\mathbf{d}^{obs}$. The inverse problem is formulated as an optimization problem: \n", 458 | "\n", 459 | "\n", 460 | "$$\\text{minimize} \\ \\ \\ \\phi(\\mathbf{m}) = \\phi_d(\\mathbf{m}) + \\beta \\phi_m(\\mathbf{m}) $$\n", 461 | "\n", 462 | "where \n", 463 | "\n", 464 | "- $\\phi_d$: data misfit\n", 465 | "- $\\phi_m$: model regularization\n", 466 | "- $\\beta$: trade-off (or Tikhonov) parameter $0<\\beta<\\infty$\n", 467 | "\n", 468 | "Data misfit is defined as \n", 469 | "\n", 470 | "$$ \\phi_d = \\sum_{j=1}^{N}\\Big(\\frac{\\mathbf{g}_j\\mathbf{m}-d^{obs}_j}{\\epsilon_j}\\Big)^2$$\n", 471 | "\n", 472 | "where $\\epsilon_j$ is an estimate of the standard deviation of the $j$th datum.\n", 473 | "\n", 474 | "\n", 475 | "\n" 476 | ] 477 | }, 478 | { 479 | "cell_type": "markdown", 480 | "metadata": { 481 | "iooxa": { 482 | "id": { 483 | "block": "DPEqdTQucMMV6G5ONR7c", 484 | "project": "VNMrkxzChhdveZyf6lmb", 485 | "version": 1 486 | } 487 | } 488 | }, 489 | "source": [ 490 | "The model regularization term, $\\phi_m$, can be written as \n", 491 | "\n", 492 | "$$ \\phi_m(\\mathbf{m}) = \\alpha_s \\int (\\mathbf{m}-\\mathbf{m}_{ref}) dx + \\alpha_x \\int (\\frac{d \\mathbf{m}}{dx}) dx$$\n", 493 | "\n", 494 | "The first term is referred to as the \"smallness\" term. Minimizing this generates a model that is close to a reference model $m_{ref}$. The second term penalizes roughness of the model. It is generically referred to as a \"flattest\" or \"smoothness\" term. " 495 | ] 496 | }, 497 | { 498 | "cell_type": "markdown", 499 | "metadata": { 500 | "iooxa": { 501 | "id": { 502 | "block": "j7hVPYIA73jbxv0bAinV", 503 | "project": "VNMrkxzChhdveZyf6lmb", 504 | "version": 6 505 | } 506 | } 507 | }, 508 | "source": [ 509 | "## Step5: Invert the data, and explore inversion results\n", 510 | "\n" 511 | ] 512 | }, 513 | { 514 | "cell_type": "markdown", 515 | "metadata": { 516 | "iooxa": { 517 | "id": { 518 | "block": "8k53qHfi1hcwT29TMxNC", 519 | "project": "VNMrkxzChhdveZyf6lmb", 520 | "version": 3 521 | } 522 | } 523 | }, 524 | "source": [ 525 | "In the inverse problem we define parameters needed to evaluate the data misfit and the model regularization terms. We then deal with parameters associated with the inversion.\n", 526 | "\n", 527 | "### Parameters\n", 528 | "\n", 529 | "- `mode`: `Run` or `Explore`\n", 530 | " - `Run`: Each click of the app, will run `n_beta` times of inversion\n", 531 | " - `Explore`: Not running inversions, but explore result of the inversions\n", 532 | "\n", 533 | "\n", 534 | "- `noise option`: `error contaminated` or `clean data`\n", 535 | "\n", 536 | "#### Misfit\n", 537 | "- `percent`: percentage of the uncertainty (%)\n", 538 | "\n", 539 | "- `floor`: floor of the uncertainty (%)\n", 540 | "\n", 541 | "- `chifact`: chi factor for stopping criteria (when $\\phi_d^{\\ast}=N \\rightarrow$ `chifact=1`)\n", 542 | "\n", 543 | "#### Model norm\n", 544 | "- `mref`: reference model\n", 545 | "\n", 546 | "- `alpha_s`: $\\alpha_s$ for smallness\n", 547 | "\n", 548 | "- `alpha_x`: $\\alpha_x$ for smoothness\n", 549 | "\n", 550 | "#### Beta\n", 551 | "- `beta_min`: minimum $\\beta$\n", 552 | "\n", 553 | "- `beta_max`: maximum $\\beta$\n", 554 | "\n", 555 | "- `n_beta`: the number of $\\beta$\n", 556 | "\n", 557 | "#### Plotting options\n", 558 | "\n", 559 | "- `data`: `obs & pred` or `normalized misfit`\n", 560 | " - `obs & pred`: show observed and predicted data\n", 561 | " - `normalized misfit`: show normalized misfit\n", 562 | "\n", 563 | "\n", 564 | "- `tikhonov`: `phi_d & phi_m` or `phi_d vs phi_m`\n", 565 | " - `phi_d & phi_m`: show $\\phi_d$ and $\\phi_m$ as a function of $\\beta$\n", 566 | " - `phi_d vs phi_m`: show tikhonov curve\n", 567 | " \n", 568 | "- `i_beta`: i-th $\\beta$ value\n", 569 | "\n", 570 | "- `scale`: `linear` or `log`\n", 571 | " - `linear`: linear scale for plotting the third panel\n", 572 | " - `log`: log scale for plotting the third panel " 573 | ] 574 | }, 575 | { 576 | "cell_type": "markdown", 577 | "metadata": {}, 578 | "source": [ 579 | "### $$\\text{minimize} \\ \\ \\ \\phi(\\mathbf{m}) = \\phi_d(\\mathbf{m}) + \\beta \\phi_m(\\mathbf{m}) $$" 580 | ] 581 | }, 582 | { 583 | "cell_type": "code", 584 | "execution_count": 43, 585 | "metadata": { 586 | "iooxa": { 587 | "id": { 588 | "block": "zYsKlEza7On1cZcJUO1x", 589 | "project": "VNMrkxzChhdveZyf6lmb", 590 | "version": 21 591 | }, 592 | "outputId": { 593 | "block": "XHsuSTVmGPtaRtV3WO9m", 594 | "project": "VNMrkxzChhdveZyf6lmb", 595 | "version": 21 596 | } 597 | }, 598 | "scrolled": false 599 | }, 600 | "outputs": [ 601 | { 602 | "data": { 603 | "application/vnd.jupyter.widget-view+json": { 604 | "model_id": "87337f8a284040a2ad713455eae8302e", 605 | "version_major": 2, 606 | "version_minor": 0 607 | }, 608 | "text/plain": [ 609 | "VBox(children=(HBox(children=(RadioButtons(description='mode', options=('Run', 'Explore'), value='Run'), Radio…" 610 | ] 611 | }, 612 | "metadata": {}, 613 | "output_type": "display_data" 614 | } 615 | ], 616 | "source": [ 617 | "app.interact_plot_inversion()" 618 | ] 619 | }, 620 | { 621 | "cell_type": "code", 622 | "execution_count": null, 623 | "metadata": {}, 624 | "outputs": [], 625 | "source": [] 626 | } 627 | ], 628 | "metadata": { 629 | "iooxa": { 630 | "id": { 631 | "block": "pwZt1IzpzYFhHsoFuTHA", 632 | "project": "FesuZy9nERkRH9r67vBe", 633 | "version": 1 634 | } 635 | }, 636 | "kernelspec": { 637 | "display_name": "Python 3", 638 | "language": "python", 639 | "name": "python3" 640 | }, 641 | "language_info": { 642 | "codemirror_mode": { 643 | "name": "ipython", 644 | "version": 3 645 | }, 646 | "file_extension": ".py", 647 | "mimetype": "text/x-python", 648 | "name": "python", 649 | "nbconvert_exporter": "python", 650 | "pygments_lexer": "ipython3", 651 | "version": "3.7.3" 652 | } 653 | }, 654 | "nbformat": 4, 655 | "nbformat_minor": 2 656 | } 657 | -------------------------------------------------------------------------------- /LICENSE: -------------------------------------------------------------------------------- 1 | MIT License 2 | 3 | Copyright (c) 2020 SimPEG 4 | 5 | Permission is hereby granted, free of charge, to any person obtaining a copy 6 | of this software and associated documentation files (the "Software"), to deal 7 | in the Software without restriction, including without limitation the rights 8 | to use, copy, modify, merge, publish, distribute, sublicense, and/or sell 9 | copies of the Software, and to permit persons to whom the Software is 10 | furnished to do so, subject to the following conditions: 11 | 12 | The above copyright notice and this permission notice shall be included in all 13 | copies or substantial portions of the Software. 14 | 15 | THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR 16 | IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, 17 | FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE 18 | AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER 19 | LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, 20 | OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE 21 | SOFTWARE. 22 | -------------------------------------------------------------------------------- /README.md: -------------------------------------------------------------------------------- 1 | **[summary](#summary) | [prerequisites](#prerequisites) | [setup](#setup) | [resources](#resources) | [license](#license)** 2 | 3 | # Transform 2021: Inversion for Geologists 4 | 5 | [![Binder](https://mybinder.org/badge_logo.svg)](https://mybinder.org/v2/gh/simpeg/transform-2021-simpeg/HEAD) 6 | [![License](https://img.shields.io/github/license/simpeg/transform-2021-simpeg.svg)](https://github.com/simpeg/transform-2021-simpeg/blob/master/LICENSE) 7 | 8 | | | Info | 9 | |--------:|:-----| 10 | | When | Tuesday, April 20 • 20:00 - 22:00 UTC | 11 | | Slack | [Software Underground](https://softwareunderground.org/) channel `t21-tue-inversion-for-geologists` | 12 | | YouTube | https://youtu.be/5MiaebDwWUQ | 13 | | conda environment | `t21-tue-simpeg` | 14 | | Slides | [SimPEG Transform presentation](http://bit.ly/transform-2021-slides) | 15 | 16 | **Team** 17 | - [Seogi Kang](https://github.com/sgkang) (Instructor) 18 | - [Lindsey Heagy](http://github.com/lheagy) 19 | - [Doug Oldenburg](http://github.com/dougoldenburg) 20 | - [Dom Fournier](https://github.com/fourndo) 21 | - [Joe Capriotti](https://github.com/jcapriot) 22 | - [Dieter Werthmüller](http://github.com/prisae) 23 | - and the [SimPEG contributors](https://github.com/simpeg/simpeg/graphs/contributors) 24 | 25 | ## Before the tutorial 26 | 27 | Make sure you've done these things before the tutorial on Tuesday: 28 | 29 | 1. Sign-up for the [Software Underground Slack](https://softwareunderground.org/slack) 30 | 1. Join the channel `t21-tue-inversion-for-geologists`. This is where **all communication will happen**. 31 | 1. Set up your computer ([instructions below](#usage)). We will not have time to 32 | solve many computer issues during the tutorial so make sure you do this 33 | ahead of time. If you need any help, ask at the `t21-tue-inversion-for-geologists` channel on 34 | Slack. 35 | 1. If you have some data you'd like to process, please have it ready and make 36 | sure you can load it with pandas or numpy. You'll have some time at the end 37 | of the tutorial to work on your own data. 38 | 39 | ## Summary 40 | 41 | This repo contains the notebooks and tutorial resources for the Transform 2021 tutorial on Inversion for Geologists. 42 | In this tutorial, we will provide fundamental concepts about the inversion framework. As a motivating example, 43 | we will use the magnetic data acquired at [Raglan deposits](https://en.wikipedia.org/wiki/Raglan_Mine) in Northern Quebec, Candada. 44 | About 20 years ago, this magnetic data set was inverted in 3D (sort of the first time), and resulting susceptiblity model 45 | played an imporant role for spotting a drill hole, which found mineralized zones. 46 | 47 | 48 | Starting from field data in a text file we will learn how to 49 | - load those data into SimPEG 50 | - construct a survey object that contains the geometry of the sources and receivers 51 | - set up and run a forward simulation 52 | - define the inverse problem consisting of a data misfit and regularization 53 | - run an inversion and discuss inversion strategies 54 | 55 | Then, we will work with a synthetic example to 56 | - demonstrate how to explore aspects of the physics with SimPEG 57 | - explore the role and influence of parameters used in an inversion 58 | 59 | ## Prerequisites 60 | 61 | **Software** 62 | 63 | * Some knowledge of Python is assumed. 64 | * All coding will be done in Jupyter notebooks. I'll explain how they work 65 | briefly but it will help if you've used them before. 66 | * We'll use [numpy](https://numpy.org/), [matplotlib](https://matplotlib.org/), and 67 | [ipywidgets](https://ipywidgets.readthedocs.io/) 68 | You don't need to be an expert in these tools but some familiarity will help. 69 | 70 | **Geophysical Inversions** 71 | 72 | * This tutorial will focus on providing fundamental concepts of the geophysical inversion. 73 | I do not assume an extensive background in inversions, but it will help if you have been 74 | introduced to some concepts. The [Inversion Module](https://curvenote.com/@geosci/inversion-module/inverse-theory-overview) 75 | provides a nice overview. 76 | * As a motivating example, I will use magnetic mehod. 77 | I'll explain the basic principles, but if these are new to you, then I would recommend 78 | taking a read through the [Magnetics](https://gpg.geosci.xyz/content/magnetics/index.html). 79 | 80 | ## Usage 81 | 82 | There are a few things you'll need to follow the tutorial: 83 | 84 | 1. A working Python installation ([Anaconda](https://www.anaconda.com/products/individual) or [Miniconda](https://docs.conda.io/en/latest/miniconda.html)) 85 | 2. The SimPEG *conda environment* installed 86 | 3. A web browser that works with Jupyter 87 | (basically anything except Internet Explorer) 88 | 89 | To get things setup, please do the following. 90 | 91 | **If you have any trouble**, please ask for help in the 92 | `t21-tue-inversion-for-geologists` channel on the Software Underground slack. 93 | 94 | **Windows users:** When you see "*terminal*" in the instructions, 95 | this means the "*Anaconda Prompt*" program for you. 96 | 97 | ### Step 1: Python 98 | 99 | **Follow the general instructions for Transform2021:** https://www.youtube.com/playlist?list=PLgLft9vxdduAW-jmhYqXvtfGYJS6v2FjM 100 | 101 | This will get you a working Python 3 installation with the `conda` package 102 | manager. If you already have one, you can skip this step. 103 | 104 | ### Step 2: Download the SimPEG tutorials 105 | 106 | To access the notebooks, there are 3 options (in order of preference): 107 | 1. Use git to clone this repository 108 | 2. From GitHub, you can use the `download` option to download this repository as a zip file (follow all instructions below, replacing the `git clone` step with download and unzip the zip file with the repository contents. 109 | 3. You can run the notebooks online with binder through: https://mybinder.org/v2/gh/simpeg/transform-2021-simpeg/master 110 | 111 | To clone this repository, open up a terminal and navigate to where you want this repository stored on your computer. 112 | 113 | Then run 114 | ``` 115 | git clone https://github.com/simpeg/transform-2021-simpeg.git 116 | ``` 117 | to clone the repository, and `cd` into the `transform-2021-simpeg` directory 118 | ``` 119 | cd transform-2021-simpeg 120 | ``` 121 | 122 | ### Step 3: Create the SimPEG tutorial conda environment 123 | 124 | From inside of the `transform-2021-simpeg` repository, create the `t21-tue-simpeg` conda environment 125 | ``` 126 | conda env create -f environment.yml 127 | ``` 128 | and activate the environment 129 | ``` 130 | conda activate t21-tue-simpeg 131 | ``` 132 | 133 | ### Step 4: Launching the notebooks 134 | 135 | Once you have activated the conda environment, you can launch the notebooks 136 | ``` 137 | jupyter notebook 138 | ``` 139 | Jupyter will then launch in your web-browser. 140 | 141 | If you are able to open any one of the notebooks and run the first cell, then you should be good to go! 142 | If you run into issues, please post in the #t21-tue-simpeg slack channel. 143 | 144 | ## Resources 145 | 146 | **Resources on Geophysics and Inversions** 147 | - [Geophysics for Practicing Geoscientists](https://gpg.geosci.xyz/) 148 | - [Inversion Module](https://curvenote.com/@geosci/inversion-module/inverse-theory-overview) 149 | - [SimPEG](https://www.simepg.xyz) 150 | 151 | **Resources on SimPEG** 152 | - [Docs](http://docs.simpeg.xyz/) 153 | - [Discourse](http://simpeg.discourse.group/) 154 | - [Slack](http://slack.simpeg.xyz/) 155 | - [This Tutorial](https://github.com/simpeg/transform-2021-simpeg/) 156 | 157 | 158 | ## License 159 | 160 | All code and text in this repository is free software: you can redistribute it and/or 161 | modify it under the terms of the MIT License. 162 | A copy of this license is provided in [LICENSE](LICENSE). 163 | -------------------------------------------------------------------------------- /data/Raglan_1997/inv.inp: -------------------------------------------------------------------------------- 1 | 1 60 !! irest, maxit 2 | 50 !! ndpt, number of subspace vectors 3 | 27.0 !! chifact 4 | mesh !! mesh 5 | null !! topo.idx; default=null 6 | maginv3d.mtx !! sensitivity 7 | obs.mag !! obs.mag 8 | maginv3d.wtw !! weighting matrix 9 | null !! reference model 10 | null !! initial model 11 | 0.0005 0.02 12 | 0 13 | 1 14 | -------------------------------------------------------------------------------- /data/Raglan_1997/mesh.msh: -------------------------------------------------------------------------------- 1 | 40 40 10 2 | 500 39000 0 3 | 40*100.0 4 | 40*100.0 5 | 10*100.0 6 | -------------------------------------------------------------------------------- /data/Raglan_1997/obs.mag: -------------------------------------------------------------------------------- 1 | 83.00 -32.00 60000.00 !! incl, decl, geomag 2 | 83.00 -32.00 !! aincl, adecl: anomaly 3 | 1638 !! # of data 4 | 509.00 42958.00 40.00 0.429000E+02 0.585800E+01 5 | 504.00 42858.00 40.00 0.883200E+02 0.676640E+01 6 | 499.30 42758.00 40.00 0.861100E+02 0.672220E+01 7 | 498.10 42658.00 40.00 0.563100E+02 0.612620E+01 8 | 498.60 42558.00 40.00 0.217100E+02 0.543420E+01 9 | 500.80 42458.00 40.00 0.500000E+01 0.510000E+01 10 | 503.00 42358.00 40.00 0.556000E+01 0.511120E+01 11 | 505.20 42258.00 40.00 0.207000E+01 0.504140E+01 12 | 508.00 42158.00 40.00 0.200000E+00 0.500400E+01 13 | 510.90 42058.00 40.00 -0.900000E-01 0.500180E+01 14 | 514.80 41958.00 40.00 -0.430000E+00 0.500860E+01 15 | 517.00 41858.00 40.00 -0.399000E+01 0.507980E+01 16 | 517.00 41758.00 40.00 -0.740000E+01 0.514800E+01 17 | 516.00 41658.00 40.00 -0.100300E+02 0.520060E+01 18 | 514.80 41558.00 40.00 -0.117900E+02 0.523580E+01 19 | 511.10 41458.00 40.00 -0.120800E+02 0.524160E+01 20 | 507.30 41358.00 40.00 -0.104700E+02 0.520940E+01 21 | 505.10 41258.00 40.00 -0.139000E+01 0.502780E+01 22 | 503.00 41159.00 40.00 0.858000E+01 0.517160E+01 23 | 503.00 41059.00 40.00 0.446000E+01 0.508920E+01 24 | 500.80 40959.00 40.00 0.906400E+02 0.681280E+01 25 | 498.50 40859.00 40.00 0.196330E+03 0.892660E+01 26 | 497.20 40759.00 40.00 0.269790E+03 0.103958E+02 27 | 496.60 40659.00 40.00 0.324000E+03 0.114800E+02 28 | 496.70 40559.00 40.00 0.372460E+03 0.124492E+02 29 | 497.30 40459.00 40.00 0.428640E+03 0.135728E+02 30 | 498.30 40359.00 40.00 0.506040E+03 0.151208E+02 31 | 499.60 40259.00 40.00 0.618150E+03 0.173630E+02 32 | 501.20 40159.00 40.00 0.782430E+03 0.206486E+02 33 | 504.60 40059.00 40.00 0.997450E+03 0.249490E+02 34 | 564.00 40041.00 40.00 0.105820E+04 0.261640E+02 35 | 564.80 40141.00 40.00 0.916060E+03 0.233212E+02 36 | 568.50 40241.00 40.00 0.708360E+03 0.191672E+02 37 | 573.50 40341.00 40.00 0.514910E+03 0.152982E+02 38 | 578.70 40440.00 40.00 0.353900E+03 0.120780E+02 39 | 579.00 40540.00 40.00 0.224980E+03 0.949960E+01 40 | 579.00 40640.00 40.00 0.146290E+03 0.792580E+01 41 | 579.40 40740.00 40.00 0.922500E+02 0.684500E+01 42 | 581.50 40840.00 40.00 0.567300E+02 0.613460E+01 43 | 583.50 40940.00 40.00 0.400400E+02 0.580080E+01 44 | 589.30 41040.00 40.00 0.378300E+02 0.575660E+01 45 | 597.00 41140.00 40.00 0.168000E+02 0.533600E+01 46 | 607.60 41239.00 40.00 0.141100E+02 0.528220E+01 47 | 619.00 41339.00 40.00 0.179000E+01 0.503580E+01 48 | 628.00 41438.00 40.00 -0.374000E+01 0.507480E+01 49 | 632.60 41538.00 40.00 -0.705000E+01 0.514100E+01 50 | 634.00 41638.00 40.00 -0.629000E+01 0.512580E+01 51 | 631.90 41738.00 40.00 -0.625000E+01 0.512500E+01 52 | 630.50 41838.00 40.00 -0.481000E+01 0.509620E+01 53 | 626.10 41938.00 40.00 -0.262000E+01 0.505240E+01 54 | 620.30 42038.00 40.00 0.390000E+00 0.500780E+01 55 | 615.40 42138.00 40.00 -0.154000E+01 0.503080E+01 56 | 611.80 42238.00 40.00 -0.400000E-01 0.500080E+01 57 | 608.90 42337.00 40.00 0.370000E+01 0.507400E+01 58 | 606.10 42437.00 40.00 0.886000E+01 0.517720E+01 59 | 603.60 42537.00 40.00 0.186400E+02 0.537280E+01 60 | 598.10 42637.00 40.00 0.319400E+02 0.563880E+01 61 | 590.90 42737.00 40.00 0.772400E+02 0.654480E+01 62 | 585.70 42837.00 40.00 0.146950E+03 0.793900E+01 63 | 581.70 42937.00 40.00 0.150930E+03 0.801860E+01 64 | 651.00 42958.00 40.00 0.809000E+02 0.661800E+01 65 | 654.90 42858.00 40.00 0.114990E+03 0.729980E+01 66 | 662.50 42758.00 40.00 0.871400E+02 0.674280E+01 67 | 671.60 42659.00 40.00 0.457400E+02 0.591480E+01 68 | 679.10 42559.00 40.00 0.285700E+02 0.557140E+01 69 | 686.50 42459.00 40.00 0.129700E+02 0.525940E+01 70 | 693.00 42359.00 40.00 0.672000E+01 0.513440E+01 71 | 699.30 42260.00 40.00 0.212000E+01 0.504240E+01 72 | 702.80 42160.00 40.00 -0.690000E+00 0.501380E+01 73 | 704.30 42060.00 40.00 -0.164000E+01 0.503280E+01 74 | 706.00 41960.00 40.00 -0.370000E+01 0.507400E+01 75 | 708.50 41860.00 40.00 -0.274000E+01 0.505480E+01 76 | 711.00 41760.00 40.00 -0.301000E+01 0.506020E+01 77 | 715.20 41660.00 40.00 -0.273000E+01 0.505460E+01 78 | 720.00 41560.00 40.00 -0.178000E+01 0.503560E+01 79 | 723.10 41460.00 40.00 0.960000E+00 0.501920E+01 80 | 719.90 41360.00 40.00 0.632000E+01 0.512640E+01 81 | 714.20 41260.00 40.00 0.150800E+02 0.530160E+01 82 | 709.50 41160.00 40.00 0.283900E+02 0.556780E+01 83 | 704.00 41061.00 40.00 0.404400E+02 0.580880E+01 84 | 703.00 40961.00 40.00 0.564900E+02 0.612980E+01 85 | 704.60 40861.00 40.00 0.799500E+02 0.659900E+01 86 | 703.10 40761.00 40.00 0.112000E+03 0.724000E+01 87 | 698.50 40661.00 40.00 0.163120E+03 0.826240E+01 88 | 696.00 40561.00 40.00 0.238730E+03 0.977460E+01 89 | 696.00 40461.00 40.00 0.339210E+03 0.117842E+02 90 | 697.00 40361.00 40.00 0.478750E+03 0.145750E+02 91 | 697.90 40261.00 40.00 0.655240E+03 0.181048E+02 92 | 699.90 40161.00 40.00 0.834650E+03 0.216930E+02 93 | 702.30 40061.00 40.00 0.976050E+03 0.245210E+02 94 | 759.00 40029.00 40.00 0.961600E+03 0.242320E+02 95 | 760.90 40129.00 40.00 0.890860E+03 0.228172E+02 96 | 767.00 40228.00 40.00 0.749950E+03 0.199990E+02 97 | 777.50 40328.00 40.00 0.590480E+03 0.168096E+02 98 | 791.30 40427.00 40.00 0.444150E+03 0.138830E+02 99 | 803.30 40526.00 40.00 0.335660E+03 0.117132E+02 100 | 811.10 40626.00 40.00 0.245360E+03 0.990720E+01 101 | 813.00 40725.00 40.00 0.176890E+03 0.853780E+01 102 | 809.80 40825.00 40.00 0.125540E+03 0.751080E+01 103 | 803.20 40925.00 40.00 0.903100E+02 0.680620E+01 104 | 796.10 41025.00 40.00 0.612600E+02 0.622520E+01 105 | 791.30 41125.00 40.00 0.424500E+02 0.584900E+01 106 | 791.00 41225.00 40.00 0.259600E+02 0.551920E+01 107 | 794.90 41325.00 40.00 0.137100E+02 0.527420E+01 108 | 802.50 41424.00 40.00 0.627000E+01 0.512540E+01 109 | 810.40 41524.00 40.00 0.327000E+01 0.506540E+01 110 | 815.50 41624.00 40.00 0.224000E+01 0.504480E+01 111 | 817.20 41724.00 40.00 0.700000E-01 0.500140E+01 112 | 819.10 41824.00 40.00 -0.480000E+00 0.500960E+01 113 | 822.00 41924.00 40.00 -0.137000E+01 0.502740E+01 114 | 824.00 42024.00 40.00 -0.220000E+01 0.504400E+01 115 | 824.10 42124.00 40.00 -0.900000E-01 0.500180E+01 116 | 820.60 42224.00 40.00 0.344000E+01 0.506880E+01 117 | 813.30 42323.00 40.00 0.645000E+01 0.512900E+01 118 | 804.30 42423.00 40.00 0.853000E+01 0.517060E+01 119 | 798.60 42523.00 40.00 0.816000E+01 0.516320E+01 120 | 797.00 42623.00 40.00 0.217000E+01 0.504340E+01 121 | 796.10 42723.00 40.00 0.579000E+01 0.511580E+01 122 | 792.00 42823.00 40.00 0.470200E+02 0.594040E+01 123 | 789.90 42923.00 40.00 0.708800E+02 0.641760E+01 124 | 890.00 42969.00 40.00 0.783000E+02 0.656600E+01 125 | 889.00 42869.00 40.00 0.612000E+02 0.622400E+01 126 | 889.00 42769.00 40.00 0.246600E+02 0.549320E+01 127 | 891.00 42669.00 40.00 0.207700E+02 0.541540E+01 128 | 891.50 42569.00 40.00 -0.879000E+01 0.517580E+01 129 | 891.00 42469.00 40.00 0.419000E+01 0.508380E+01 130 | 891.00 42369.00 40.00 0.547000E+01 0.510940E+01 131 | 892.60 42269.00 40.00 0.607000E+01 0.512140E+01 132 | 898.30 42169.00 40.00 0.439000E+01 0.508780E+01 133 | 904.10 42069.00 40.00 0.163000E+01 0.503260E+01 134 | 910.10 41969.00 40.00 -0.390000E+00 0.500780E+01 135 | 916.80 41870.00 40.00 0.298000E+01 0.505960E+01 136 | 921.60 41770.00 40.00 0.232000E+01 0.504640E+01 137 | 923.00 41670.00 40.00 0.249000E+01 0.504980E+01 138 | 923.00 41570.00 40.00 0.427000E+01 0.508540E+01 139 | 922.60 41470.00 40.00 0.902000E+01 0.518040E+01 140 | 920.30 41370.00 40.00 0.184900E+02 0.536980E+01 141 | 916.00 41270.00 40.00 0.275900E+02 0.555180E+01 142 | 908.70 41170.00 40.00 0.440400E+02 0.588080E+01 143 | 901.70 41070.00 40.00 0.651300E+02 0.630260E+01 144 | 898.10 40970.00 40.00 0.925400E+02 0.685080E+01 145 | 896.00 40871.00 40.00 0.126710E+03 0.753420E+01 146 | 895.30 40771.00 40.00 0.169150E+03 0.838300E+01 147 | 897.00 40671.00 40.00 0.228120E+03 0.956240E+01 148 | 892.80 40571.00 40.00 0.298760E+03 0.109752E+02 149 | 887.50 40471.00 40.00 0.386330E+03 0.127266E+02 150 | 882.30 40371.00 40.00 0.497550E+03 0.149510E+02 151 | 880.60 40271.00 40.00 0.617160E+03 0.173432E+02 152 | 884.00 40171.00 40.00 0.707310E+03 0.191462E+02 153 | 891.20 40071.00 40.00 0.748550E+03 0.199710E+02 154 | 992.00 40028.00 40.00 0.675900E+03 0.185180E+02 155 | 994.00 40128.00 40.00 0.641420E+03 0.178284E+02 156 | 996.20 40228.00 40.00 0.617160E+03 0.173432E+02 157 | 999.60 40328.00 40.00 0.556560E+03 0.161312E+02 158 | 1004.50 40427.00 40.00 0.477710E+03 0.145542E+02 159 | 1009.50 40527.00 40.00 0.397190E+03 0.129438E+02 160 | 1016.10 40627.00 40.00 0.318220E+03 0.113644E+02 161 | 1023.00 40727.00 40.00 0.251620E+03 0.100324E+02 162 | 1025.00 40827.00 40.00 0.195020E+03 0.890040E+01 163 | 1023.70 40927.00 40.00 0.148760E+03 0.797520E+01 164 | 1022.00 41027.00 40.00 0.106040E+03 0.712080E+01 165 | 1018.60 41127.00 40.00 0.718000E+02 0.643600E+01 166 | 1014.60 41227.00 40.00 0.505800E+02 0.601160E+01 167 | 1010.60 41327.00 40.00 0.289600E+02 0.557920E+01 168 | 1004.00 41426.00 40.00 0.177400E+02 0.535480E+01 169 | 999.10 41526.00 40.00 0.106500E+02 0.521300E+01 170 | 998.00 41626.00 40.00 0.626000E+01 0.512520E+01 171 | 996.60 41726.00 40.00 0.271000E+01 0.505420E+01 172 | 996.00 41826.00 40.00 0.289000E+01 0.505780E+01 173 | 996.00 41926.00 40.00 0.141000E+01 0.502820E+01 174 | 994.30 42026.00 40.00 0.181000E+01 0.503620E+01 175 | 991.80 42126.00 40.00 0.334000E+01 0.506680E+01 176 | 991.00 42226.00 40.00 0.605000E+01 0.512100E+01 177 | 991.90 42326.00 40.00 0.778000E+01 0.515560E+01 178 | 994.60 42426.00 40.00 0.843000E+01 0.516860E+01 179 | 997.30 42526.00 40.00 0.125100E+02 0.525020E+01 180 | 1000.00 42626.00 40.00 0.313800E+02 0.562760E+01 181 | 1004.50 42726.00 40.00 0.214500E+02 0.542900E+01 182 | 1008.50 42826.00 40.00 0.311600E+02 0.562320E+01 183 | 1011.50 42926.00 40.00 0.739500E+02 0.647900E+01 184 | 1109.00 42963.00 40.00 0.563000E+02 0.612600E+01 185 | 1108.70 42863.00 40.00 0.375400E+02 0.575080E+01 186 | 1108.00 42763.00 40.00 0.276200E+02 0.555240E+01 187 | 1108.00 42663.00 40.00 0.138300E+02 0.527660E+01 188 | 1108.90 42563.00 40.00 0.117900E+02 0.523580E+01 189 | 1110.00 42463.00 40.00 0.960000E+01 0.519200E+01 190 | 1111.00 42363.00 40.00 0.886000E+01 0.517720E+01 191 | 1111.00 42263.00 40.00 0.669000E+01 0.513380E+01 192 | 1111.00 42163.00 40.00 0.607000E+01 0.512140E+01 193 | 1109.40 42063.00 40.00 0.519000E+01 0.510380E+01 194 | 1107.10 41963.00 40.00 0.505000E+01 0.510100E+01 195 | 1104.00 41863.00 40.00 0.411000E+01 0.508220E+01 196 | 1100.50 41763.00 40.00 0.642000E+01 0.512840E+01 197 | 1097.10 41663.00 40.00 0.924000E+01 0.518480E+01 198 | 1092.70 41563.00 40.00 0.135200E+02 0.527040E+01 199 | 1089.10 41463.00 40.00 0.209300E+02 0.541860E+01 200 | 1086.90 41363.00 40.00 0.436500E+02 0.587300E+01 201 | 1086.10 41263.00 40.00 0.903500E+02 0.680700E+01 202 | 1089.80 41163.00 40.00 0.728500E+02 0.645700E+01 203 | 1093.60 41063.00 40.00 0.964100E+02 0.692820E+01 204 | 1095.00 40963.00 40.00 0.138250E+03 0.776500E+01 205 | 1097.60 40863.00 40.00 0.184910E+03 0.869820E+01 206 | 1102.50 40764.00 40.00 0.237750E+03 0.975500E+01 207 | 1104.80 40664.00 40.00 0.300730E+03 0.110146E+02 208 | 1108.10 40564.00 40.00 0.372420E+03 0.124484E+02 209 | 1110.50 40464.00 40.00 0.441660E+03 0.138332E+02 210 | 1112.60 40364.00 40.00 0.507090E+03 0.151418E+02 211 | 1114.00 40264.00 40.00 0.544610E+03 0.158922E+02 212 | 1114.00 40164.00 40.00 0.565300E+03 0.163060E+02 213 | 1115.30 40064.00 40.00 0.608110E+03 0.171622E+02 214 | 1181.00 40026.00 40.00 0.605300E+03 0.171060E+02 215 | 1184.20 40126.00 40.00 0.583350E+03 0.166670E+02 216 | 1190.10 40226.00 40.00 0.528540E+03 0.155708E+02 217 | 1197.40 40325.00 40.00 0.506890E+03 0.151378E+02 218 | 1204.30 40425.00 40.00 0.477410E+03 0.145482E+02 219 | 1212.70 40525.00 40.00 0.425500E+03 0.135100E+02 220 | 1219.00 40624.00 40.00 0.366670E+03 0.123334E+02 221 | 1217.60 40724.00 40.00 0.306280E+03 0.111256E+02 222 | 1214.30 40824.00 40.00 0.244610E+03 0.989220E+01 223 | 1209.00 40924.00 40.00 0.190760E+03 0.881520E+01 224 | 1202.80 41024.00 40.00 0.144880E+03 0.789760E+01 225 | 1198.50 41124.00 40.00 0.105100E+03 0.710200E+01 226 | 1194.40 41224.00 40.00 0.846100E+02 0.669220E+01 227 | 1191.90 41324.00 40.00 0.320030E+03 0.114006E+02 228 | 1193.50 41424.00 40.00 0.100950E+03 0.701900E+01 229 | 1199.00 41524.00 40.00 0.298500E+02 0.559700E+01 230 | 1203.30 41623.00 40.00 0.138800E+02 0.527760E+01 231 | 1207.00 41723.00 40.00 0.118200E+02 0.523640E+01 232 | 1211.80 41823.00 40.00 0.775000E+01 0.515500E+01 233 | 1214.60 41923.00 40.00 0.704000E+01 0.514080E+01 234 | 1217.40 42023.00 40.00 0.104300E+02 0.520860E+01 235 | 1221.10 42123.00 40.00 0.706000E+01 0.514120E+01 236 | 1225.80 42223.00 40.00 0.566000E+01 0.511320E+01 237 | 1229.50 42323.00 40.00 0.750000E+01 0.515000E+01 238 | 1231.30 42423.00 40.00 0.115200E+02 0.523040E+01 239 | 1232.00 42523.00 40.00 0.120500E+02 0.524100E+01 240 | 1232.00 42623.00 40.00 0.998000E+01 0.519960E+01 241 | 1230.40 42723.00 40.00 0.123000E+02 0.524600E+01 242 | 1225.50 42823.00 40.00 0.172500E+02 0.534500E+01 243 | 1217.50 42922.00 40.00 0.231200E+02 0.546240E+01 244 | 1308.00 42970.00 40.00 0.238000E+02 0.547600E+01 245 | 1311.50 42870.00 40.00 0.246000E+01 0.504920E+01 246 | 1314.40 42770.00 40.00 0.424000E+01 0.508480E+01 247 | 1315.00 42670.00 40.00 0.793000E+01 0.515860E+01 248 | 1314.10 42570.00 40.00 0.788000E+01 0.515760E+01 249 | 1311.40 42470.00 40.00 0.102000E+02 0.520400E+01 250 | 1307.60 42370.00 40.00 0.995000E+01 0.519900E+01 251 | 1302.30 42270.00 40.00 0.728000E+01 0.514560E+01 252 | 1297.30 42170.00 40.00 0.770000E+01 0.515400E+01 253 | 1292.50 42070.00 40.00 0.115100E+02 0.523020E+01 254 | 1287.60 41970.00 40.00 0.114600E+02 0.522920E+01 255 | 1283.00 41871.00 40.00 0.987000E+01 0.519740E+01 256 | 1279.10 41771.00 40.00 0.114000E+02 0.522800E+01 257 | 1276.90 41671.00 40.00 0.167100E+02 0.533420E+01 258 | 1276.70 41571.00 40.00 0.250000E+02 0.550000E+01 259 | 1276.00 41471.00 40.00 0.672300E+02 0.634460E+01 260 | 1277.00 41371.00 40.00 0.212210E+03 0.924420E+01 261 | 1280.40 41271.00 40.00 0.232040E+03 0.964080E+01 262 | 1285.00 41171.00 40.00 0.102900E+03 0.705800E+01 263 | 1290.80 41071.00 40.00 0.131610E+03 0.763220E+01 264 | 1295.80 40971.00 40.00 0.172620E+03 0.845240E+01 265 | 1301.10 40871.00 40.00 0.223000E+03 0.946000E+01 266 | 1305.50 40771.00 40.00 0.279960E+03 0.105992E+02 267 | 1310.90 40672.00 40.00 0.336860E+03 0.117372E+02 268 | 1315.20 40572.00 40.00 0.394290E+03 0.128858E+02 269 | 1316.80 40472.00 40.00 0.434220E+03 0.136844E+02 270 | 1315.50 40372.00 40.00 0.465280E+03 0.143056E+02 271 | 1312.20 40272.00 40.00 0.477710E+03 0.145542E+02 272 | 1307.80 40172.00 40.00 0.507050E+03 0.151410E+02 273 | 1305.60 40072.00 40.00 0.527750E+03 0.155550E+02 274 | 1378.00 40047.00 40.00 0.465300E+03 0.143060E+02 275 | 1383.50 40147.00 40.00 0.461480E+03 0.142296E+02 276 | 1391.10 40246.00 40.00 0.440090E+03 0.138018E+02 277 | 1399.20 40346.00 40.00 0.437600E+03 0.137520E+02 278 | 1403.30 40446.00 40.00 0.433880E+03 0.136776E+02 279 | 1405.00 40546.00 40.00 0.413530E+03 0.132706E+02 280 | 1405.00 40646.00 40.00 0.375650E+03 0.125130E+02 281 | 1404.00 40746.00 40.00 0.323690E+03 0.114738E+02 282 | 1402.70 40846.00 40.00 0.264230E+03 0.102846E+02 283 | 1399.70 40946.00 40.00 0.213350E+03 0.926700E+01 284 | 1396.90 41046.00 40.00 0.169860E+03 0.839720E+01 285 | 1396.00 41146.00 40.00 0.132230E+03 0.764460E+01 286 | 1395.00 41246.00 40.00 0.101790E+03 0.703580E+01 287 | 1396.80 41346.00 40.00 0.225450E+03 0.950900E+01 288 | 1398.00 41446.00 40.00 0.781600E+02 0.656320E+01 289 | 1397.00 41546.00 40.00 0.399000E+02 0.579800E+01 290 | 1394.60 41646.00 40.00 0.272200E+02 0.554440E+01 291 | 1392.80 41746.00 40.00 0.210900E+02 0.542180E+01 292 | 1390.10 41846.00 40.00 0.153600E+02 0.530720E+01 293 | 1387.40 41945.00 40.00 0.125400E+02 0.525080E+01 294 | 1384.60 42045.00 40.00 0.103400E+02 0.520680E+01 295 | 1384.00 42145.00 40.00 0.108800E+02 0.521760E+01 296 | 1385.00 42245.00 40.00 0.918000E+01 0.518360E+01 297 | 1387.70 42345.00 40.00 0.931000E+01 0.518620E+01 298 | 1391.60 42445.00 40.00 0.977000E+01 0.519540E+01 299 | 1395.50 42545.00 40.00 0.921000E+01 0.518420E+01 300 | 1399.00 42645.00 40.00 0.101400E+02 0.520280E+01 301 | 1401.00 42745.00 40.00 0.749000E+01 0.514980E+01 302 | 1403.90 42845.00 40.00 0.129100E+02 0.525820E+01 303 | 1404.00 42945.00 40.00 0.183600E+02 0.536720E+01 304 | 1502.00 42967.00 40.00 0.187000E+02 0.537400E+01 305 | 1504.60 42867.00 40.00 0.151500E+02 0.530300E+01 306 | 1506.20 42767.00 40.00 0.112300E+02 0.522460E+01 307 | 1507.80 42667.00 40.00 0.109600E+02 0.521920E+01 308 | 1510.50 42567.00 40.00 0.115100E+02 0.523020E+01 309 | 1513.00 42467.00 40.00 0.956000E+01 0.519120E+01 310 | 1515.00 42367.00 40.00 0.115000E+02 0.523000E+01 311 | 1515.90 42267.00 40.00 0.125200E+02 0.525040E+01 312 | 1513.50 42167.00 40.00 0.136800E+02 0.527360E+01 313 | 1510.50 42067.00 40.00 0.164800E+02 0.532960E+01 314 | 1505.80 41967.00 40.00 0.206700E+02 0.541340E+01 315 | 1498.90 41867.00 40.00 0.248700E+02 0.549740E+01 316 | 1492.50 41768.00 40.00 0.277600E+02 0.555520E+01 317 | 1490.00 41668.00 40.00 0.318900E+02 0.563780E+01 318 | 1487.80 41568.00 40.00 0.403400E+02 0.580680E+01 319 | 1485.50 41468.00 40.00 0.526500E+02 0.605300E+01 320 | 1484.00 41368.00 40.00 0.134540E+03 0.769080E+01 321 | 1485.20 41268.00 40.00 0.105920E+03 0.711840E+01 322 | 1488.80 41168.00 40.00 0.131780E+03 0.763560E+01 323 | 1493.80 41068.00 40.00 0.165020E+03 0.830040E+01 324 | 1497.10 40968.00 40.00 0.207190E+03 0.914380E+01 325 | 1501.30 40868.00 40.00 0.250560E+03 0.100112E+02 326 | 1504.70 40768.00 40.00 0.298990E+03 0.109798E+02 327 | 1508.00 40668.00 40.00 0.342240E+03 0.118448E+02 328 | 1511.40 40568.00 40.00 0.375320E+03 0.125064E+02 329 | 1515.00 40468.00 40.00 0.395180E+03 0.129036E+02 330 | 1516.90 40368.00 40.00 0.397330E+03 0.129466E+02 331 | 1517.00 40268.00 40.00 0.377270E+03 0.125454E+02 332 | 1517.30 40168.00 40.00 0.359410E+03 0.121882E+02 333 | 1519.70 40068.00 40.00 0.394370E+03 0.128874E+02 334 | 1575.00 40053.00 40.00 0.333600E+03 0.116720E+02 335 | 1581.30 40153.00 40.00 0.349660E+03 0.119932E+02 336 | 1588.50 40252.00 40.00 0.329970E+03 0.115994E+02 337 | 1595.90 40352.00 40.00 0.350980E+03 0.120196E+02 338 | 1600.60 40452.00 40.00 0.365310E+03 0.123062E+02 339 | 1602.40 40552.00 40.00 0.366930E+03 0.123386E+02 340 | 1605.50 40652.00 40.00 0.348830E+03 0.119766E+02 341 | 1609.60 40752.00 40.00 0.316970E+03 0.113394E+02 342 | 1612.00 40852.00 40.00 0.279200E+03 0.105840E+02 343 | 1611.60 40952.00 40.00 0.239140E+03 0.978280E+01 344 | 1609.60 41052.00 40.00 0.200210E+03 0.900420E+01 345 | 1606.80 41152.00 40.00 0.163460E+03 0.826920E+01 346 | 1602.80 41252.00 40.00 0.127020E+03 0.754040E+01 347 | 1598.20 41351.00 40.00 0.106260E+03 0.712520E+01 348 | 1594.60 41451.00 40.00 0.185770E+03 0.871540E+01 349 | 1591.60 41551.00 40.00 0.814100E+02 0.662820E+01 350 | 1590.00 41651.00 40.00 0.482000E+02 0.596400E+01 351 | 1589.90 41751.00 40.00 0.376400E+02 0.575280E+01 352 | 1592.70 41851.00 40.00 0.343800E+02 0.568760E+01 353 | 1597.30 41951.00 40.00 0.270600E+02 0.554120E+01 354 | 1601.60 42051.00 40.00 0.219000E+02 0.543800E+01 355 | 1606.50 42151.00 40.00 0.169600E+02 0.533920E+01 356 | 1610.00 42251.00 40.00 0.144300E+02 0.528860E+01 357 | 1612.10 42351.00 40.00 0.111800E+02 0.522360E+01 358 | 1610.60 42451.00 40.00 0.884000E+01 0.517680E+01 359 | 1606.90 42551.00 40.00 0.106800E+02 0.521360E+01 360 | 1604.10 42651.00 40.00 0.114300E+02 0.522860E+01 361 | 1601.50 42751.00 40.00 0.106800E+02 0.521360E+01 362 | 1600.00 42851.00 40.00 0.118400E+02 0.523680E+01 363 | 1599.70 42951.00 40.00 0.377000E+01 0.507540E+01 364 | 1712.00 42967.00 40.00 0.142000E+02 0.528400E+01 365 | 1709.30 42867.00 40.00 0.131600E+02 0.526320E+01 366 | 1706.90 42767.00 40.00 0.123500E+02 0.524700E+01 367 | 1705.00 42667.00 40.00 0.120900E+02 0.524180E+01 368 | 1704.50 42567.00 40.00 0.110500E+02 0.522100E+01 369 | 1704.00 42467.00 40.00 0.100800E+02 0.520160E+01 370 | 1704.00 42367.00 40.00 0.114900E+02 0.522980E+01 371 | 1703.00 42267.00 40.00 0.143700E+02 0.528740E+01 372 | 1703.00 42167.00 40.00 0.201000E+02 0.540200E+01 373 | 1704.00 42067.00 40.00 0.258800E+02 0.551760E+01 374 | 1705.00 41967.00 40.00 0.382800E+02 0.576560E+01 375 | 1703.60 41867.00 40.00 0.417500E+02 0.583500E+01 376 | 1701.00 41767.00 40.00 0.483800E+02 0.596760E+01 377 | 1700.00 41667.00 40.00 0.569300E+02 0.613860E+01 378 | 1701.00 41567.00 40.00 0.843600E+02 0.668720E+01 379 | 1701.00 41467.00 40.00 0.168000E+03 0.836000E+01 380 | 1701.00 41367.00 40.00 0.132210E+03 0.764420E+01 381 | 1701.00 41267.00 40.00 0.133000E+03 0.766000E+01 382 | 1703.00 41167.00 40.00 0.162200E+03 0.824400E+01 383 | 1704.60 41067.00 40.00 0.196800E+03 0.893600E+01 384 | 1706.00 40967.00 40.00 0.233900E+03 0.967800E+01 385 | 1708.30 40867.00 40.00 0.270380E+03 0.104076E+02 386 | 1709.60 40767.00 40.00 0.302960E+03 0.110592E+02 387 | 1710.00 40667.00 40.00 0.327660E+03 0.115532E+02 388 | 1707.70 40567.00 40.00 0.341630E+03 0.118326E+02 389 | 1705.00 40467.00 40.00 0.343140E+03 0.118628E+02 390 | 1705.60 40367.00 40.00 0.331770E+03 0.116354E+02 391 | 1705.10 40267.00 40.00 0.307660E+03 0.111532E+02 392 | 1701.50 40167.00 40.00 0.296430E+03 0.109286E+02 393 | 1700.00 40067.00 40.00 0.291820E+03 0.108364E+02 394 | 1789.00 40040.00 40.00 0.220300E+03 0.940600E+01 395 | 1789.00 40140.00 40.00 0.269210E+03 0.103842E+02 396 | 1791.00 40240.00 40.00 0.267100E+03 0.103420E+02 397 | 1792.30 40340.00 40.00 0.289290E+03 0.107858E+02 398 | 1793.00 40440.00 40.00 0.317620E+03 0.113524E+02 399 | 1794.80 40540.00 40.00 0.328130E+03 0.115626E+02 400 | 1798.00 40640.00 40.00 0.327770E+03 0.115554E+02 401 | 1803.20 40739.00 40.00 0.315190E+03 0.113038E+02 402 | 1806.10 40839.00 40.00 0.287550E+03 0.107510E+02 403 | 1808.00 40939.00 40.00 0.259000E+03 0.101800E+02 404 | 1809.00 41039.00 40.00 0.225310E+03 0.950620E+01 405 | 1809.80 41139.00 40.00 0.194760E+03 0.889520E+01 406 | 1810.90 41239.00 40.00 0.163960E+03 0.827920E+01 407 | 1812.00 41339.00 40.00 0.137370E+03 0.774740E+01 408 | 1811.80 41439.00 40.00 0.112330E+03 0.724660E+01 409 | 1809.10 41539.00 40.00 0.914800E+02 0.682960E+01 410 | 1805.50 41639.00 40.00 0.699100E+02 0.639820E+01 411 | 1802.90 41739.00 40.00 0.596200E+02 0.619240E+01 412 | 1801.00 41839.00 40.00 0.489900E+02 0.597980E+01 413 | 1798.00 41939.00 40.00 0.464000E+02 0.592800E+01 414 | 1797.00 42039.00 40.00 0.341500E+02 0.568300E+01 415 | 1798.00 42139.00 40.00 0.247000E+02 0.549400E+01 416 | 1798.80 42239.00 40.00 0.183700E+02 0.536740E+01 417 | 1801.60 42339.00 40.00 0.138800E+02 0.527760E+01 418 | 1802.90 42439.00 40.00 0.100000E+02 0.520000E+01 419 | 1804.50 42539.00 40.00 0.950000E+01 0.519000E+01 420 | 1806.30 42639.00 40.00 0.103000E+02 0.520600E+01 421 | 1808.00 42739.00 40.00 0.133100E+02 0.526620E+01 422 | 1806.90 42839.00 40.00 0.122000E+02 0.524400E+01 423 | 1804.10 42939.00 40.00 0.111400E+02 0.522280E+01 424 | 1909.00 42970.00 40.00 0.930000E+01 0.518600E+01 425 | 1910.80 42870.00 40.00 0.940000E+01 0.518800E+01 426 | 1911.00 42770.00 40.00 0.118800E+02 0.523760E+01 427 | 1910.00 42670.00 40.00 0.107800E+02 0.521560E+01 428 | 1908.40 42570.00 40.00 0.104000E+02 0.520800E+01 429 | 1905.90 42470.00 40.00 0.114200E+02 0.522840E+01 430 | 1903.50 42370.00 40.00 0.130800E+02 0.526160E+01 431 | 1900.10 42270.00 40.00 0.169700E+02 0.533940E+01 432 | 1897.60 42170.00 40.00 0.225000E+02 0.545000E+01 433 | 1895.20 42070.00 40.00 0.298300E+02 0.559660E+01 434 | 1894.00 41970.00 40.00 0.412500E+02 0.582500E+01 435 | 1894.60 41870.00 40.00 0.493400E+02 0.598680E+01 436 | 1896.00 41770.00 40.00 0.581700E+02 0.616340E+01 437 | 1896.30 41670.00 40.00 0.706300E+02 0.641260E+01 438 | 1896.90 41570.00 40.00 0.899600E+02 0.679920E+01 439 | 1899.00 41470.00 40.00 0.115020E+03 0.730040E+01 440 | 1903.50 41370.00 40.00 0.134820E+03 0.769640E+01 441 | 1907.50 41270.00 40.00 0.159960E+03 0.819920E+01 442 | 1911.80 41170.00 40.00 0.188810E+03 0.877620E+01 443 | 1916.40 41070.00 40.00 0.217410E+03 0.934820E+01 444 | 1922.50 40971.00 40.00 0.248430E+03 0.996860E+01 445 | 1927.50 40871.00 40.00 0.275960E+03 0.105192E+02 446 | 1931.00 40771.00 40.00 0.295640E+03 0.109128E+02 447 | 1934.20 40671.00 40.00 0.307750E+03 0.111550E+02 448 | 1936.00 40571.00 40.00 0.307390E+03 0.111478E+02 449 | 1935.00 40471.00 40.00 0.304000E+03 0.110800E+02 450 | 1928.10 40371.00 40.00 0.290990E+03 0.108198E+02 451 | 1920.40 40271.00 40.00 0.266560E+03 0.103312E+02 452 | 1913.40 40172.00 40.00 0.255590E+03 0.101118E+02 453 | 1906.00 40072.00 40.00 0.205530E+03 0.911060E+01 454 | 1972.00 40029.00 40.00 0.140400E+03 0.780800E+01 455 | 1977.00 40129.00 40.00 0.197620E+03 0.895240E+01 456 | 1985.10 40228.00 40.00 0.265880E+03 0.103176E+02 457 | 1992.60 40328.00 40.00 0.279880E+03 0.105976E+02 458 | 1998.50 40428.00 40.00 0.282540E+03 0.106508E+02 459 | 2002.00 40528.00 40.00 0.297840E+03 0.109568E+02 460 | 2003.00 40628.00 40.00 0.302700E+03 0.110540E+02 461 | 2004.00 40728.00 40.00 0.302710E+03 0.110542E+02 462 | 2003.00 40828.00 40.00 0.292840E+03 0.108568E+02 463 | 2002.00 40928.00 40.00 0.272890E+03 0.104578E+02 464 | 2001.00 41028.00 40.00 0.246540E+03 0.993080E+01 465 | 2000.00 41128.00 40.00 0.216890E+03 0.933780E+01 466 | 1997.60 41228.00 40.00 0.187360E+03 0.874720E+01 467 | 1995.90 41328.00 40.00 0.161240E+03 0.822480E+01 468 | 1995.00 41428.00 40.00 0.137120E+03 0.774240E+01 469 | 1994.90 41528.00 40.00 0.115700E+03 0.731400E+01 470 | 1997.20 41628.00 40.00 0.913500E+02 0.682700E+01 471 | 2001.10 41727.00 40.00 0.745400E+02 0.649080E+01 472 | 2002.00 41827.00 40.00 0.600000E+02 0.620000E+01 473 | 2001.00 41927.00 40.00 0.468400E+02 0.593680E+01 474 | 2002.00 42027.00 40.00 0.360400E+02 0.572080E+01 475 | 2004.40 42127.00 40.00 0.276500E+02 0.555300E+01 476 | 2009.10 42227.00 40.00 0.212900E+02 0.542580E+01 477 | 2011.00 42327.00 40.00 0.158400E+02 0.531680E+01 478 | 2009.50 42427.00 40.00 0.144800E+02 0.528960E+01 479 | 2005.70 42527.00 40.00 0.115500E+02 0.523100E+01 480 | 2002.00 42627.00 40.00 0.100000E+02 0.520000E+01 481 | 1998.60 42727.00 40.00 0.991000E+01 0.519820E+01 482 | 1994.40 42827.00 40.00 0.999000E+01 0.519980E+01 483 | 1990.00 42927.00 40.00 0.150200E+02 0.530040E+01 484 | 2092.00 42963.00 40.00 0.900000E+01 0.518000E+01 485 | 2093.60 42863.00 40.00 0.132300E+02 0.526460E+01 486 | 2095.00 42763.00 40.00 0.108200E+02 0.521640E+01 487 | 2095.00 42663.00 40.00 0.714000E+01 0.514280E+01 488 | 2092.50 42563.00 40.00 0.752000E+01 0.515040E+01 489 | 2089.00 42463.00 40.00 0.950000E+01 0.519000E+01 490 | 2085.40 42363.00 40.00 0.130300E+02 0.526060E+01 491 | 2083.90 42263.00 40.00 0.146600E+02 0.529320E+01 492 | 2082.80 42163.00 40.00 0.202800E+02 0.540560E+01 493 | 2085.10 42063.00 40.00 0.297500E+02 0.559500E+01 494 | 2088.30 41963.00 40.00 0.396900E+02 0.579380E+01 495 | 2094.60 41863.00 40.00 0.539600E+02 0.607920E+01 496 | 2103.20 41764.00 40.00 0.691600E+02 0.638320E+01 497 | 2108.30 41664.00 40.00 0.860700E+02 0.672140E+01 498 | 2113.30 41564.00 40.00 0.107000E+03 0.714000E+01 499 | 2118.30 41464.00 40.00 0.127240E+03 0.754480E+01 500 | 2123.20 41364.00 40.00 0.156960E+03 0.813920E+01 501 | 2128.00 41264.00 40.00 0.187010E+03 0.874020E+01 502 | 2125.20 41164.00 40.00 0.215130E+03 0.930260E+01 503 | 2122.30 41064.00 40.00 0.241540E+03 0.983080E+01 504 | 2120.00 40964.00 40.00 0.265680E+03 0.103136E+02 505 | 2117.60 40864.00 40.00 0.285250E+03 0.107050E+02 506 | 2115.30 40764.00 40.00 0.298050E+03 0.109610E+02 507 | 2113.10 40664.00 40.00 0.297140E+03 0.109428E+02 508 | 2111.00 40564.00 40.00 0.288070E+03 0.107614E+02 509 | 2109.70 40465.00 40.00 0.274900E+03 0.104980E+02 510 | 2107.50 40365.00 40.00 0.264330E+03 0.102866E+02 511 | 2104.00 40265.00 40.00 0.248640E+03 0.997280E+01 512 | 2105.00 40165.00 40.00 0.193460E+03 0.886920E+01 513 | 2106.00 40065.00 40.00 0.136000E+03 0.772000E+01 514 | 2171.00 40043.00 40.00 0.114900E+03 0.729800E+01 515 | 2177.00 40142.00 40.00 0.152000E+03 0.804000E+01 516 | 2182.70 40242.00 40.00 0.221450E+03 0.942900E+01 517 | 2185.20 40342.00 40.00 0.238250E+03 0.976500E+01 518 | 2187.90 40442.00 40.00 0.253140E+03 0.100628E+02 519 | 2191.10 40542.00 40.00 0.273710E+03 0.104742E+02 520 | 2191.30 40642.00 40.00 0.290130E+03 0.108026E+02 521 | 2188.00 40742.00 40.00 0.299860E+03 0.109972E+02 522 | 2191.10 40842.00 40.00 0.299300E+03 0.109860E+02 523 | 2195.20 40942.00 40.00 0.283430E+03 0.106686E+02 524 | 2198.40 41042.00 40.00 0.264220E+03 0.102844E+02 525 | 2198.70 41142.00 40.00 0.240290E+03 0.980580E+01 526 | 2195.80 41242.00 40.00 0.213050E+03 0.926100E+01 527 | 2193.90 41342.00 40.00 0.183990E+03 0.867980E+01 528 | 2191.00 41442.00 40.00 0.152730E+03 0.805460E+01 529 | 2188.10 41542.00 40.00 0.127310E+03 0.754620E+01 530 | 2185.80 41642.00 40.00 0.105570E+03 0.711140E+01 531 | 2184.00 41742.00 40.00 0.847800E+02 0.669560E+01 532 | 2185.00 41842.00 40.00 0.694800E+02 0.638960E+01 533 | 2188.50 41942.00 40.00 0.546600E+02 0.609320E+01 534 | 2193.00 42041.00 40.00 0.401900E+02 0.580380E+01 535 | 2198.50 42141.00 40.00 0.289600E+02 0.557920E+01 536 | 2203.00 42241.00 40.00 0.191400E+02 0.538280E+01 537 | 2206.70 42341.00 40.00 0.155900E+02 0.531180E+01 538 | 2208.00 42441.00 40.00 0.105300E+02 0.521060E+01 539 | 2208.00 42541.00 40.00 0.691000E+01 0.513820E+01 540 | 2207.00 42641.00 40.00 0.610000E+01 0.512200E+01 541 | 2204.30 42741.00 40.00 0.872000E+01 0.517440E+01 542 | 2201.90 42841.00 40.00 0.920000E+01 0.518400E+01 543 | 2202.00 42941.00 40.00 0.548000E+01 0.510960E+01 544 | 2286.00 42966.00 40.00 0.770000E+01 0.515400E+01 545 | 2284.00 42866.00 40.00 0.973000E+01 0.519460E+01 546 | 2285.00 42766.00 40.00 0.701000E+01 0.514020E+01 547 | 2282.80 42666.00 40.00 0.575000E+01 0.511500E+01 548 | 2280.30 42566.00 40.00 0.607000E+01 0.512140E+01 549 | 2279.00 42466.00 40.00 0.104600E+02 0.520920E+01 550 | 2280.60 42366.00 40.00 0.151300E+02 0.530260E+01 551 | 2284.00 42266.00 40.00 0.177500E+02 0.535500E+01 552 | 2287.40 42166.00 40.00 0.237100E+02 0.547420E+01 553 | 2289.80 42066.00 40.00 0.370400E+02 0.574080E+01 554 | 2292.20 41966.00 40.00 0.551600E+02 0.610320E+01 555 | 2295.30 41866.00 40.00 0.733500E+02 0.646700E+01 556 | 2300.20 41766.00 40.00 0.886300E+02 0.677260E+01 557 | 2304.60 41666.00 40.00 0.107520E+03 0.715040E+01 558 | 2308.00 41566.00 40.00 0.125880E+03 0.751760E+01 559 | 2311.50 41466.00 40.00 0.148340E+03 0.796680E+01 560 | 2312.90 41366.00 40.00 0.178940E+03 0.857880E+01 561 | 2313.10 41267.00 40.00 0.210540E+03 0.921080E+01 562 | 2314.00 41167.00 40.00 0.242700E+03 0.985400E+01 563 | 2315.00 41067.00 40.00 0.269560E+03 0.103912E+02 564 | 2315.00 40967.00 40.00 0.294920E+03 0.108984E+02 565 | 2315.00 40867.00 40.00 0.308040E+03 0.111608E+02 566 | 2312.90 40767.00 40.00 0.310000E+03 0.112000E+02 567 | 2312.60 40667.00 40.00 0.296590E+03 0.109318E+02 568 | 2310.20 40567.00 40.00 0.274790E+03 0.104958E+02 569 | 2308.00 40467.00 40.00 0.247660E+03 0.995320E+01 570 | 2305.70 40367.00 40.00 0.220710E+03 0.941420E+01 571 | 2304.30 40267.00 40.00 0.212930E+03 0.925860E+01 572 | 2304.00 40167.00 40.00 0.162280E+03 0.824560E+01 573 | 2305.50 40067.00 40.00 0.113480E+03 0.726960E+01 574 | 2401.00 40039.00 40.00 0.802000E+02 0.660400E+01 575 | 2395.30 40139.00 40.00 0.114180E+03 0.728360E+01 576 | 2390.40 40238.00 40.00 0.175790E+03 0.851580E+01 577 | 2385.90 40338.00 40.00 0.194480E+03 0.888960E+01 578 | 2386.00 40438.00 40.00 0.215720E+03 0.931440E+01 579 | 2387.00 40538.00 40.00 0.255120E+03 0.101024E+02 580 | 2387.00 40638.00 40.00 0.288340E+03 0.107668E+02 581 | 2387.00 40738.00 40.00 0.309790E+03 0.111958E+02 582 | 2389.80 40838.00 40.00 0.322480E+03 0.114496E+02 583 | 2392.90 40938.00 40.00 0.317830E+03 0.113566E+02 584 | 2394.90 41038.00 40.00 0.299430E+03 0.109886E+02 585 | 2396.00 41138.00 40.00 0.272800E+03 0.104560E+02 586 | 2396.00 41238.00 40.00 0.241530E+03 0.983060E+01 587 | 2396.90 41338.00 40.00 0.208810E+03 0.917620E+01 588 | 2399.00 41438.00 40.00 0.176750E+03 0.853500E+01 589 | 2401.10 41538.00 40.00 0.148840E+03 0.797680E+01 590 | 2402.00 41638.00 40.00 0.125610E+03 0.751220E+01 591 | 2403.00 41738.00 40.00 0.108220E+03 0.716440E+01 592 | 2404.00 41838.00 40.00 0.981300E+02 0.696260E+01 593 | 2405.90 41938.00 40.00 0.840200E+02 0.668040E+01 594 | 2407.00 42038.00 40.00 0.572700E+02 0.614540E+01 595 | 2407.00 42138.00 40.00 0.312600E+02 0.562520E+01 596 | 2407.00 42238.00 40.00 0.184500E+02 0.536900E+01 597 | 2408.00 42338.00 40.00 0.122700E+02 0.524540E+01 598 | 2408.00 42438.00 40.00 0.871000E+01 0.517420E+01 599 | 2405.70 42538.00 40.00 0.821000E+01 0.516420E+01 600 | 2402.80 42638.00 40.00 0.513000E+01 0.510260E+01 601 | 2399.80 42738.00 40.00 0.361000E+01 0.507220E+01 602 | 2398.00 42838.00 40.00 0.464000E+01 0.509280E+01 603 | 2398.00 42938.00 40.00 0.695000E+01 0.513900E+01 604 | 2496.00 42965.00 40.00 -0.400000E+00 0.500800E+01 605 | 2502.50 42865.00 40.00 -0.183000E+01 0.503660E+01 606 | 2506.50 42765.00 40.00 0.171000E+01 0.503420E+01 607 | 2508.40 42665.00 40.00 0.213000E+01 0.504260E+01 608 | 2511.20 42565.00 40.00 0.539000E+01 0.510780E+01 609 | 2512.90 42465.00 40.00 0.755000E+01 0.515100E+01 610 | 2514.60 42365.00 40.00 0.112100E+02 0.522420E+01 611 | 2515.00 42265.00 40.00 0.165000E+02 0.533000E+01 612 | 2514.50 42165.00 40.00 0.267400E+02 0.553480E+01 613 | 2513.10 42065.00 40.00 0.463700E+02 0.592740E+01 614 | 2510.70 41965.00 40.00 0.812300E+02 0.662460E+01 615 | 2508.30 41865.00 40.00 0.114820E+03 0.729640E+01 616 | 2504.90 41765.00 40.00 0.122500E+03 0.745000E+01 617 | 2501.30 41665.00 40.00 0.125690E+03 0.751380E+01 618 | 2501.00 41565.00 40.00 0.146950E+03 0.793900E+01 619 | 2502.30 41465.00 40.00 0.176500E+03 0.853000E+01 620 | 2506.20 41365.00 40.00 0.207240E+03 0.914480E+01 621 | 2510.00 41266.00 40.00 0.243340E+03 0.986680E+01 622 | 2513.50 41166.00 40.00 0.281900E+03 0.106380E+02 623 | 2516.80 41066.00 40.00 0.316100E+03 0.113220E+02 624 | 2519.20 40966.00 40.00 0.339350E+03 0.117870E+02 625 | 2519.30 40866.00 40.00 0.351170E+03 0.120234E+02 626 | 2518.00 40766.00 40.00 0.338680E+03 0.117736E+02 627 | 2515.50 40666.00 40.00 0.315360E+03 0.113072E+02 628 | 2511.90 40566.00 40.00 0.280310E+03 0.106062E+02 629 | 2507.30 40466.00 40.00 0.236480E+03 0.972960E+01 630 | 2502.70 40366.00 40.00 0.201790E+03 0.903580E+01 631 | 2497.10 40266.00 40.00 0.171880E+03 0.843760E+01 632 | 2490.20 40166.00 40.00 0.123000E+03 0.746000E+01 633 | 2485.00 40067.00 40.00 0.740100E+02 0.648020E+01 634 | 2589.00 40032.00 40.00 0.734000E+02 0.646800E+01 635 | 2590.50 40132.00 40.00 0.631300E+02 0.626260E+01 636 | 2593.00 40232.00 40.00 0.112360E+03 0.724720E+01 637 | 2594.90 40332.00 40.00 0.164290E+03 0.828580E+01 638 | 2593.00 40432.00 40.00 0.205780E+03 0.911560E+01 639 | 2591.30 40532.00 40.00 0.255520E+03 0.101104E+02 640 | 2589.00 40632.00 40.00 0.302440E+03 0.110488E+02 641 | 2589.00 40732.00 40.00 0.340660E+03 0.118132E+02 642 | 2589.00 40832.00 40.00 0.364500E+03 0.122900E+02 643 | 2590.00 40932.00 40.00 0.372380E+03 0.124476E+02 644 | 2591.30 41032.00 40.00 0.359990E+03 0.121998E+02 645 | 2592.10 41131.00 40.00 0.330080E+03 0.116016E+02 646 | 2593.10 41231.00 40.00 0.288570E+03 0.107714E+02 647 | 2594.00 41331.00 40.00 0.246420E+03 0.992840E+01 648 | 2595.00 41431.00 40.00 0.208920E+03 0.917840E+01 649 | 2596.00 41531.00 40.00 0.171500E+03 0.843000E+01 650 | 2597.00 41631.00 40.00 0.142970E+03 0.785940E+01 651 | 2597.00 41731.00 40.00 0.129700E+03 0.759400E+01 652 | 2597.00 41831.00 40.00 0.129540E+03 0.759080E+01 653 | 2598.00 41931.00 40.00 0.107060E+03 0.714120E+01 654 | 2598.00 42031.00 40.00 0.646700E+02 0.629340E+01 655 | 2599.30 42131.00 40.00 0.367600E+02 0.573520E+01 656 | 2601.00 42231.00 40.00 0.204500E+02 0.540900E+01 657 | 2601.60 42331.00 40.00 0.114100E+02 0.522820E+01 658 | 2602.00 42431.00 40.00 0.564000E+01 0.511280E+01 659 | 2601.00 42531.00 40.00 0.131000E+01 0.502620E+01 660 | 2599.60 42631.00 40.00 -0.167000E+01 0.503340E+01 661 | 2598.00 42731.00 40.00 -0.361000E+01 0.507220E+01 662 | 2596.20 42831.00 40.00 -0.300000E+01 0.506000E+01 663 | 2594.50 42931.00 40.00 -0.347000E+01 0.506940E+01 664 | 2708.00 42967.00 40.00 -0.820000E+01 0.516400E+01 665 | 2708.00 42867.00 40.00 -0.800000E+01 0.516000E+01 666 | 2708.00 42767.00 40.00 -0.109300E+02 0.521860E+01 667 | 2708.90 42667.00 40.00 -0.958000E+01 0.519160E+01 668 | 2710.30 42567.00 40.00 -0.625000E+01 0.512500E+01 669 | 2711.00 42467.00 40.00 -0.200000E+01 0.504000E+01 670 | 2709.80 42367.00 40.00 0.287000E+01 0.505740E+01 671 | 2708.30 42267.00 40.00 0.932000E+01 0.518640E+01 672 | 2705.80 42167.00 40.00 0.202200E+02 0.540440E+01 673 | 2703.30 42067.00 40.00 0.350000E+02 0.570000E+01 674 | 2700.70 41967.00 40.00 0.678400E+02 0.635680E+01 675 | 2699.10 41867.00 40.00 0.109130E+03 0.718260E+01 676 | 2697.50 41767.00 40.00 0.118720E+03 0.737440E+01 677 | 2697.00 41667.00 40.00 0.131890E+03 0.763780E+01 678 | 2697.50 41567.00 40.00 0.162580E+03 0.825160E+01 679 | 2700.00 41467.00 40.00 0.200320E+03 0.900640E+01 680 | 2701.00 41367.00 40.00 0.239360E+03 0.978720E+01 681 | 2701.90 41267.00 40.00 0.278050E+03 0.105610E+02 682 | 2703.00 41167.00 40.00 0.325870E+03 0.115174E+02 683 | 2703.00 41067.00 40.00 0.374210E+03 0.124842E+02 684 | 2703.00 40967.00 40.00 0.402270E+03 0.130454E+02 685 | 2703.00 40867.00 40.00 0.405000E+03 0.131000E+02 686 | 2703.40 40767.00 40.00 0.391450E+03 0.128290E+02 687 | 2705.80 40667.00 40.00 0.354480E+03 0.120896E+02 688 | 2704.00 40567.00 40.00 0.300710E+03 0.110142E+02 689 | 2699.10 40467.00 40.00 0.240120E+03 0.980240E+01 690 | 2694.50 40367.00 40.00 0.184780E+03 0.869560E+01 691 | 2693.00 40267.00 40.00 0.124100E+03 0.748200E+01 692 | 2690.50 40167.00 40.00 0.721400E+02 0.644280E+01 693 | 2689.00 40067.00 40.00 0.469400E+02 0.593880E+01 694 | 2802.00 40043.00 40.00 0.246000E+02 0.549200E+01 695 | 2797.80 40143.00 40.00 0.193300E+02 0.538660E+01 696 | 2792.40 40242.00 40.00 0.561600E+02 0.612320E+01 697 | 2788.20 40342.00 40.00 0.136380E+03 0.772760E+01 698 | 2783.90 40442.00 40.00 0.207860E+03 0.915720E+01 699 | 2780.60 40542.00 40.00 0.275840E+03 0.105168E+02 700 | 2780.40 40642.00 40.00 0.350780E+03 0.120156E+02 701 | 2782.60 40742.00 40.00 0.407820E+03 0.131564E+02 702 | 2784.70 40842.00 40.00 0.439200E+03 0.137840E+02 703 | 2788.00 40942.00 40.00 0.442170E+03 0.138434E+02 704 | 2787.10 41042.00 40.00 0.425230E+03 0.135046E+02 705 | 2787.80 41142.00 40.00 0.381900E+03 0.126380E+02 706 | 2789.00 41242.00 40.00 0.324820E+03 0.114964E+02 707 | 2790.00 41342.00 40.00 0.278420E+03 0.105684E+02 708 | 2790.00 41442.00 40.00 0.231430E+03 0.962860E+01 709 | 2788.90 41542.00 40.00 0.186790E+03 0.873580E+01 710 | 2790.00 41642.00 40.00 0.145910E+03 0.791820E+01 711 | 2790.00 41742.00 40.00 0.135000E+03 0.770000E+01 712 | 2789.00 41842.00 40.00 0.158840E+03 0.817680E+01 713 | 2789.00 41942.00 40.00 0.855700E+02 0.671140E+01 714 | 2790.00 42042.00 40.00 0.409900E+02 0.581980E+01 715 | 2790.50 42142.00 40.00 0.225800E+02 0.545160E+01 716 | 2794.40 42242.00 40.00 0.106900E+02 0.521380E+01 717 | 2799.30 42342.00 40.00 -0.360000E+00 0.500720E+01 718 | 2804.30 42442.00 40.00 -0.643000E+01 0.512860E+01 719 | 2808.80 42542.00 40.00 -0.899000E+01 0.517980E+01 720 | 2808.60 42642.00 40.00 -0.133500E+02 0.526700E+01 721 | 2804.20 42741.00 40.00 -0.197800E+02 0.539560E+01 722 | 2801.10 42841.00 40.00 -0.190700E+02 0.538140E+01 723 | 2799.60 42941.00 40.00 -0.147000E+02 0.529400E+01 724 | 2914.00 42963.00 40.00 -0.156000E+02 0.531200E+01 725 | 2913.00 42863.00 40.00 -0.201500E+02 0.540300E+01 726 | 2913.00 42763.00 40.00 -0.203000E+02 0.540600E+01 727 | 2914.00 42663.00 40.00 -0.201000E+02 0.540200E+01 728 | 2910.30 42563.00 40.00 -0.182400E+02 0.536480E+01 729 | 2902.90 42463.00 40.00 -0.142700E+02 0.528540E+01 730 | 2897.90 42363.00 40.00 -0.932000E+01 0.518640E+01 731 | 2895.90 42263.00 40.00 -0.208000E+01 0.504160E+01 732 | 2895.90 42163.00 40.00 0.970000E+01 0.519400E+01 733 | 2897.20 42063.00 40.00 0.201600E+02 0.540320E+01 734 | 2899.80 41963.00 40.00 0.362800E+02 0.572560E+01 735 | 2902.00 41863.00 40.00 0.204480E+03 0.908960E+01 736 | 2903.00 41763.00 40.00 0.337680E+03 0.117536E+02 737 | 2904.80 41663.00 40.00 0.132830E+03 0.765660E+01 738 | 2908.80 41564.00 40.00 0.171160E+03 0.842320E+01 739 | 2913.00 41464.00 40.00 0.218420E+03 0.936840E+01 740 | 2917.50 41364.00 40.00 0.271850E+03 0.104370E+02 741 | 2921.30 41264.00 40.00 0.331340E+03 0.116268E+02 742 | 2924.00 41164.00 40.00 0.400320E+03 0.130064E+02 743 | 2925.00 41064.00 40.00 0.457770E+03 0.141554E+02 744 | 2922.70 40964.00 40.00 0.507470E+03 0.151494E+02 745 | 2921.00 40864.00 40.00 0.536950E+03 0.157390E+02 746 | 2917.30 40764.00 40.00 0.534280E+03 0.156856E+02 747 | 2912.40 40664.00 40.00 0.494420E+03 0.148884E+02 748 | 2908.60 40564.00 40.00 0.419160E+03 0.133832E+02 749 | 2904.40 40464.00 40.00 0.331680E+03 0.116336E+02 750 | 2901.10 40364.00 40.00 0.237160E+03 0.974320E+01 751 | 2898.80 40264.00 40.00 0.121460E+03 0.742920E+01 752 | 2898.00 40164.00 40.00 0.353400E+02 0.570680E+01 753 | 2898.00 40064.00 40.00 0.240400E+02 0.548080E+01 754 | 3003.00 40036.00 40.00 -0.370000E+02 0.574000E+01 755 | 3001.00 40136.00 40.00 0.404100E+02 0.580820E+01 756 | 3001.00 40236.00 40.00 0.172890E+03 0.845780E+01 757 | 3000.00 40336.00 40.00 0.355720E+03 0.121144E+02 758 | 3000.00 40436.00 40.00 0.477340E+03 0.145468E+02 759 | 3001.00 40536.00 40.00 0.536630E+03 0.157326E+02 760 | 3001.00 40636.00 40.00 0.611250E+03 0.172250E+02 761 | 3001.00 40736.00 40.00 0.650780E+03 0.180156E+02 762 | 3000.00 40836.00 40.00 0.641160E+03 0.178232E+02 763 | 2999.00 40936.00 40.00 0.591940E+03 0.168388E+02 764 | 2999.90 41036.00 40.00 0.522990E+03 0.154598E+02 765 | 3002.40 41136.00 40.00 0.452060E+03 0.140412E+02 766 | 3005.50 41236.00 40.00 0.381000E+03 0.126200E+02 767 | 3008.20 41336.00 40.00 0.310850E+03 0.112170E+02 768 | 3010.00 41435.00 40.00 0.247010E+03 0.994020E+01 769 | 3012.60 41535.00 40.00 0.196710E+03 0.893420E+01 770 | 3014.00 41635.00 40.00 0.150170E+03 0.800340E+01 771 | 3016.00 41735.00 40.00 0.115000E+03 0.730000E+01 772 | 3017.00 41835.00 40.00 0.177900E+03 0.855800E+01 773 | 3017.90 41935.00 40.00 0.390800E+02 0.578160E+01 774 | 3016.50 42035.00 40.00 0.113300E+02 0.522660E+01 775 | 3013.80 42135.00 40.00 -0.213000E+01 0.504260E+01 776 | 3012.00 42235.00 40.00 -0.155700E+02 0.531140E+01 777 | 3011.00 42335.00 40.00 -0.240900E+02 0.548180E+01 778 | 3011.00 42435.00 40.00 -0.281500E+02 0.556300E+01 779 | 3011.00 42535.00 40.00 -0.278200E+02 0.555640E+01 780 | 3011.00 42635.00 40.00 -0.271600E+02 0.554320E+01 781 | 3012.00 42735.00 40.00 -0.229500E+02 0.545900E+01 782 | 3010.70 42835.00 40.00 -0.179600E+02 0.535920E+01 783 | 3010.00 42935.00 40.00 -0.192000E+02 0.538400E+01 784 | 3126.00 42963.00 40.00 -0.810000E+01 0.516200E+01 785 | 3124.30 42863.00 40.00 -0.140500E+02 0.528100E+01 786 | 3123.00 42763.00 40.00 -0.247200E+02 0.549440E+01 787 | 3121.10 42663.00 40.00 -0.340900E+02 0.568180E+01 788 | 3120.00 42563.00 40.00 -0.387100E+02 0.577420E+01 789 | 3117.90 42463.00 40.00 -0.468400E+02 0.593680E+01 790 | 3117.10 42363.00 40.00 -0.493900E+02 0.598780E+01 791 | 3114.80 42263.00 40.00 -0.506200E+02 0.601240E+01 792 | 3114.00 42163.00 40.00 -0.407200E+02 0.581440E+01 793 | 3111.90 42063.00 40.00 -0.209200E+02 0.541840E+01 794 | 3112.00 41963.00 40.00 0.722000E+01 0.514440E+01 795 | 3111.00 41863.00 40.00 0.334500E+02 0.566900E+01 796 | 3110.00 41763.00 40.00 0.684500E+02 0.636900E+01 797 | 3109.50 41663.00 40.00 0.107170E+03 0.714340E+01 798 | 3109.00 41563.00 40.00 0.155780E+03 0.811560E+01 799 | 3108.00 41463.00 40.00 0.209460E+03 0.918920E+01 800 | 3107.00 41363.00 40.00 0.271960E+03 0.104392E+02 801 | 3106.00 41263.00 40.00 0.345860E+03 0.119172E+02 802 | 3105.00 41163.00 40.00 0.429570E+03 0.135914E+02 803 | 3104.00 41063.00 40.00 0.534740E+03 0.156948E+02 804 | 3102.10 40963.00 40.00 0.647450E+03 0.179490E+02 805 | 3097.30 40863.00 40.00 0.749810E+03 0.199962E+02 806 | 3092.70 40763.00 40.00 0.820750E+03 0.214150E+02 807 | 3089.10 40663.00 40.00 0.837330E+03 0.217466E+02 808 | 3084.30 40563.00 40.00 0.784610E+03 0.206922E+02 809 | 3080.80 40464.00 40.00 0.741240E+03 0.198248E+02 810 | 3077.40 40364.00 40.00 0.792010E+03 0.208402E+02 811 | 3075.00 40264.00 40.00 0.921510E+03 0.234302E+02 812 | 3075.20 40164.00 40.00 0.765690E+03 0.203138E+02 813 | 3077.80 40064.00 40.00 0.199770E+03 0.899540E+01 814 | 3211.00 42962.00 40.00 -0.320000E+01 0.506400E+01 815 | 3214.40 42862.00 40.00 -0.105200E+02 0.521040E+01 816 | 3215.00 42762.00 40.00 -0.235000E+02 0.547000E+01 817 | 3213.00 42662.00 40.00 -0.406400E+02 0.581280E+01 818 | 3211.70 42562.00 40.00 -0.609600E+02 0.621920E+01 819 | 3210.30 42462.00 40.00 -0.791200E+02 0.658240E+01 820 | 3210.00 42362.00 40.00 -0.978300E+02 0.695660E+01 821 | 3210.00 42262.00 40.00 -0.107340E+03 0.714680E+01 822 | 3211.00 42162.00 40.00 -0.802400E+02 0.660480E+01 823 | 3211.50 42062.00 40.00 -0.494300E+02 0.598860E+01 824 | 3212.00 41962.00 40.00 -0.160700E+02 0.532140E+01 825 | 3212.00 41862.00 40.00 0.129700E+02 0.525940E+01 826 | 3212.00 41762.00 40.00 0.470900E+02 0.594180E+01 827 | 3212.10 41662.00 40.00 0.816800E+02 0.663360E+01 828 | 3213.70 41562.00 40.00 0.133150E+03 0.766300E+01 829 | 3216.10 41462.00 40.00 0.189950E+03 0.879900E+01 830 | 3217.00 41362.00 40.00 0.257300E+03 0.101460E+02 831 | 3217.00 41262.00 40.00 0.337200E+03 0.117440E+02 832 | 3217.00 41162.00 40.00 0.430980E+03 0.136196E+02 833 | 3216.00 41062.00 40.00 0.549270E+03 0.159854E+02 834 | 3215.60 40962.00 40.00 0.695980E+03 0.189196E+02 835 | 3213.30 40862.00 40.00 0.863190E+03 0.222638E+02 836 | 3211.00 40762.00 40.00 0.102468E+04 0.254936E+02 837 | 3208.50 40662.00 40.00 0.110991E+04 0.271982E+02 838 | 3204.90 40562.00 40.00 0.110791E+04 0.271582E+02 839 | 3201.50 40462.00 40.00 0.109120E+04 0.268240E+02 840 | 3196.00 40362.00 40.00 0.130824E+04 0.311648E+02 841 | 3190.10 40263.00 40.00 0.324699E+04 0.699398E+02 842 | 3183.10 40163.00 40.00 0.438011E+04 0.926022E+02 843 | 3177.20 40063.00 40.00 0.843220E+03 0.218644E+02 844 | 3298.00 40045.00 40.00 0.443970E+04 0.937940E+02 845 | 3300.00 40145.00 40.00 0.358490E+04 0.766980E+02 846 | 3301.80 40245.00 40.00 0.438304E+04 0.926608E+02 847 | 3303.00 40345.00 40.00 0.318398E+04 0.686796E+02 848 | 3303.00 40445.00 40.00 0.179320E+04 0.408640E+02 849 | 3303.00 40545.00 40.00 0.138114E+04 0.326228E+02 850 | 3300.40 40645.00 40.00 0.126703E+04 0.303406E+02 851 | 3295.60 40745.00 40.00 0.116376E+04 0.282752E+02 852 | 3290.50 40844.00 40.00 0.994670E+03 0.248934E+02 853 | 3286.50 40944.00 40.00 0.798210E+03 0.209642E+02 854 | 3282.20 41044.00 40.00 0.613890E+03 0.172778E+02 855 | 3278.70 41144.00 40.00 0.472270E+03 0.144454E+02 856 | 3277.10 41244.00 40.00 0.364660E+03 0.122932E+02 857 | 3278.10 41344.00 40.00 0.284660E+03 0.106932E+02 858 | 3282.70 41444.00 40.00 0.207690E+03 0.915380E+01 859 | 3290.30 41544.00 40.00 0.144590E+03 0.789180E+01 860 | 3299.20 41643.00 40.00 0.905000E+02 0.681000E+01 861 | 3313.80 41742.00 40.00 0.441000E+02 0.588200E+01 862 | 3320.20 41842.00 40.00 0.446000E+01 0.508920E+01 863 | 3329.10 41941.00 40.00 -0.234800E+02 0.546960E+01 864 | 3336.30 42041.00 40.00 -0.573000E+02 0.614600E+01 865 | 3341.10 42141.00 40.00 0.486000E+01 0.509720E+01 866 | 3344.00 42241.00 40.00 0.833110E+03 0.216622E+02 867 | 3344.50 42341.00 40.00 0.340300E+02 0.568060E+01 868 | 3342.90 42441.00 40.00 -0.125040E+03 0.750080E+01 869 | 3340.30 42541.00 40.00 -0.960500E+02 0.692100E+01 870 | 3336.60 42641.00 40.00 -0.534100E+02 0.606820E+01 871 | 3332.30 42741.00 40.00 -0.209100E+02 0.541820E+01 872 | 3331.00 42841.00 40.00 -0.709000E+01 0.514180E+01 873 | 3332.50 42941.00 40.00 -0.225000E+01 0.504500E+01 874 | 3407.00 42968.00 40.00 0.118000E+02 0.523600E+01 875 | 3406.00 42868.00 40.00 0.149500E+02 0.529900E+01 876 | 3405.00 42768.00 40.00 0.139900E+02 0.527980E+01 877 | 3405.00 42668.00 40.00 -0.110100E+02 0.522020E+01 878 | 3405.00 42568.00 40.00 -0.593800E+02 0.618760E+01 879 | 3403.90 42468.00 40.00 -0.965900E+02 0.693180E+01 880 | 3403.10 42368.00 40.00 0.394200E+02 0.578840E+01 881 | 3400.80 42268.00 40.00 0.423240E+03 0.134648E+02 882 | 3400.00 42168.00 40.00 0.243460E+03 0.986920E+01 883 | 3398.80 42068.00 40.00 -0.783900E+02 0.656780E+01 884 | 3396.80 41968.00 40.00 -0.483500E+02 0.596700E+01 885 | 3393.60 41868.00 40.00 -0.203000E+02 0.540600E+01 886 | 3391.50 41768.00 40.00 0.812000E+01 0.516240E+01 887 | 3389.40 41668.00 40.00 0.442700E+02 0.588540E+01 888 | 3388.00 41568.00 40.00 0.936300E+02 0.687260E+01 889 | 3388.00 41468.00 40.00 0.148090E+03 0.796180E+01 890 | 3388.00 41368.00 40.00 0.207700E+03 0.915400E+01 891 | 3389.00 41268.00 40.00 0.307640E+03 0.111528E+02 892 | 3389.00 41168.00 40.00 0.381660E+03 0.126332E+02 893 | 3389.40 41068.00 40.00 0.493840E+03 0.148768E+02 894 | 3390.00 40968.00 40.00 0.642500E+03 0.178500E+02 895 | 3391.00 40868.00 40.00 0.836680E+03 0.217336E+02 896 | 3393.10 40768.00 40.00 0.104567E+04 0.259134E+02 897 | 3393.90 40668.00 40.00 0.120945E+04 0.291890E+02 898 | 3395.80 40568.00 40.00 0.133593E+04 0.317186E+02 899 | 3396.00 40468.00 40.00 0.152270E+04 0.354540E+02 900 | 3397.00 40368.00 40.00 0.176137E+04 0.402274E+02 901 | 3397.00 40268.00 40.00 0.251949E+04 0.553898E+02 902 | 3397.80 40168.00 40.00 0.303898E+04 0.657796E+02 903 | 3399.00 40068.00 40.00 0.234637E+04 0.519274E+02 904 | 3506.00 40032.00 40.00 0.807900E+03 0.211580E+02 905 | 3508.80 40132.00 40.00 0.199725E+04 0.449450E+02 906 | 3510.10 40232.00 40.00 0.387288E+04 0.824576E+02 907 | 3507.60 40332.00 40.00 0.157968E+04 0.365936E+02 908 | 3504.90 40432.00 40.00 0.768600E+03 0.203720E+02 909 | 3499.60 40531.00 40.00 0.976000E+03 0.245200E+02 910 | 3489.90 40631.00 40.00 0.106260E+04 0.262520E+02 911 | 3474.70 40730.00 40.00 0.101689E+04 0.253378E+02 912 | 3467.90 40829.00 40.00 0.894130E+03 0.228826E+02 913 | 3458.60 40929.00 40.00 0.717740E+03 0.193548E+02 914 | 3454.70 41029.00 40.00 0.553410E+03 0.160682E+02 915 | 3455.40 41129.00 40.00 0.428650E+03 0.135730E+02 916 | 3462.80 41229.00 40.00 0.329410E+03 0.115882E+02 917 | 3478.10 41327.00 40.00 0.255780E+03 0.101156E+02 918 | 3492.20 41426.00 40.00 0.183980E+03 0.867960E+01 919 | 3503.60 41526.00 40.00 0.109680E+03 0.719360E+01 920 | 3509.00 41625.00 40.00 0.537600E+02 0.607520E+01 921 | 3505.80 41725.00 40.00 0.112700E+02 0.522540E+01 922 | 3488.40 41824.00 40.00 -0.151500E+02 0.530300E+01 923 | 3471.80 41922.00 40.00 -0.382200E+02 0.576440E+01 924 | 3454.10 42021.00 40.00 -0.615900E+02 0.623180E+01 925 | 3438.70 42120.00 40.00 -0.871300E+02 0.674260E+01 926 | 3431.90 42219.00 40.00 0.402070E+03 0.130414E+02 927 | 3432.10 42319.00 40.00 0.204990E+03 0.909980E+01 928 | 3446.60 42418.00 40.00 0.100100E+02 0.520020E+01 929 | 3459.00 42517.00 40.00 -0.759600E+02 0.651920E+01 930 | 3474.50 42616.00 40.00 -0.282500E+02 0.556500E+01 931 | 3493.90 42714.00 40.00 0.384200E+02 0.576840E+01 932 | 3514.80 42812.00 40.00 0.741200E+02 0.648240E+01 933 | 3530.50 42911.00 40.00 0.622300E+02 0.624460E+01 934 | 3578.00 42960.00 40.00 0.357000E+02 0.571400E+01 935 | 3580.40 42860.00 40.00 0.591700E+02 0.618340E+01 936 | 3581.00 42760.00 40.00 0.719800E+02 0.643960E+01 937 | 3581.00 42660.00 40.00 0.452200E+02 0.590440E+01 938 | 3581.00 42560.00 40.00 -0.105400E+02 0.521080E+01 939 | 3579.50 42460.00 40.00 0.677500E+02 0.635500E+01 940 | 3578.00 42360.00 40.00 0.548390E+03 0.159678E+02 941 | 3577.00 42260.00 40.00 0.179500E+03 0.859000E+01 942 | 3578.00 42160.00 40.00 -0.121340E+03 0.742680E+01 943 | 3578.00 42060.00 40.00 -0.935000E+02 0.687000E+01 944 | 3578.80 41960.00 40.00 -0.580000E+02 0.616000E+01 945 | 3580.90 41860.00 40.00 -0.420600E+02 0.584120E+01 946 | 3582.10 41760.00 40.00 -0.186200E+02 0.537240E+01 947 | 3586.10 41660.00 40.00 0.140600E+02 0.528120E+01 948 | 3587.00 41560.00 40.00 0.556400E+02 0.611280E+01 949 | 3586.60 41460.00 40.00 0.110980E+03 0.721960E+01 950 | 3589.00 41360.00 40.00 0.189820E+03 0.879640E+01 951 | 3592.40 41260.00 40.00 0.246880E+03 0.993760E+01 952 | 3594.80 41160.00 40.00 0.321100E+03 0.114220E+02 953 | 3598.10 41060.00 40.00 0.402680E+03 0.130536E+02 954 | 3600.50 40960.00 40.00 0.491420E+03 0.148284E+02 955 | 3603.80 40860.00 40.00 0.591270E+03 0.168254E+02 956 | 3607.30 40760.00 40.00 0.713760E+03 0.192752E+02 957 | 3611.10 40660.00 40.00 0.866790E+03 0.223358E+02 958 | 3615.80 40561.00 40.00 0.107296E+04 0.264592E+02 959 | 3619.20 40461.00 40.00 0.186552E+04 0.423104E+02 960 | 3621.00 40361.00 40.00 0.178934E+04 0.407868E+02 961 | 3621.00 40261.00 40.00 0.145591E+04 0.341182E+02 962 | 3618.80 40161.00 40.00 0.870330E+03 0.224066E+02 963 | 3616.10 40061.00 40.00 0.669050E+03 0.183810E+02 964 | 3678.00 40050.00 40.00 0.632000E+03 0.176400E+02 965 | 3678.00 40150.00 40.00 0.427400E+03 0.135480E+02 966 | 3680.80 40250.00 40.00 0.461660E+03 0.142332E+02 967 | 3682.60 40350.00 40.00 0.800900E+03 0.210180E+02 968 | 3685.50 40450.00 40.00 0.241123E+04 0.532246E+02 969 | 3688.30 40550.00 40.00 0.126451E+04 0.302902E+02 970 | 3691.80 40650.00 40.00 0.970020E+03 0.244004E+02 971 | 3695.00 40749.00 40.00 0.699110E+03 0.189822E+02 972 | 3697.90 40849.00 40.00 0.537130E+03 0.157426E+02 973 | 3700.30 40949.00 40.00 0.455730E+03 0.141146E+02 974 | 3697.60 41049.00 40.00 0.371740E+03 0.124348E+02 975 | 3695.00 41149.00 40.00 0.293060E+03 0.108612E+02 976 | 3693.00 41249.00 40.00 0.226810E+03 0.953620E+01 977 | 3690.60 41349.00 40.00 0.169840E+03 0.839680E+01 978 | 3688.00 41449.00 40.00 0.109360E+03 0.718720E+01 979 | 3683.30 41549.00 40.00 0.548200E+02 0.609640E+01 980 | 3680.90 41649.00 40.00 0.104800E+02 0.520960E+01 981 | 3680.90 41749.00 40.00 -0.202900E+02 0.540580E+01 982 | 3683.10 41849.00 40.00 -0.472300E+02 0.594460E+01 983 | 3685.80 41949.00 40.00 -0.695100E+02 0.639020E+01 984 | 3689.60 42049.00 40.00 -0.815300E+02 0.663060E+01 985 | 3693.80 42149.00 40.00 -0.114110E+03 0.728220E+01 986 | 3697.50 42249.00 40.00 -0.125050E+03 0.750100E+01 987 | 3700.00 42349.00 40.00 0.622300E+02 0.624460E+01 988 | 3702.50 42449.00 40.00 0.581090E+03 0.166218E+02 989 | 3705.00 42549.00 40.00 -0.727100E+02 0.645420E+01 990 | 3706.00 42649.00 40.00 -0.409700E+02 0.581940E+01 991 | 3707.00 42749.00 40.00 0.446200E+02 0.589240E+01 992 | 3706.00 42849.00 40.00 0.412000E+02 0.582400E+01 993 | 3702.70 42949.00 40.00 0.220100E+02 0.544020E+01 994 | 3829.00 42969.00 40.00 -0.102000E+02 0.520400E+01 995 | 3829.00 42869.00 40.00 -0.149400E+02 0.529880E+01 996 | 3830.10 42769.00 40.00 0.120000E+02 0.524000E+01 997 | 3831.70 42669.00 40.00 0.293000E+01 0.505860E+01 998 | 3830.80 42569.00 40.00 -0.115060E+03 0.730120E+01 999 | 3827.50 42469.00 40.00 0.792600E+02 0.658520E+01 1000 | 3821.50 42369.00 40.00 0.245300E+03 0.990600E+01 1001 | 3814.10 42269.00 40.00 -0.118780E+03 0.737560E+01 1002 | 3807.00 42170.00 40.00 -0.121040E+03 0.742080E+01 1003 | 3799.80 42070.00 40.00 -0.104220E+03 0.708440E+01 1004 | 3794.40 41970.00 40.00 -0.868400E+02 0.673680E+01 1005 | 3791.80 41870.00 40.00 -0.666200E+02 0.633240E+01 1006 | 3792.90 41770.00 40.00 -0.442500E+02 0.588500E+01 1007 | 3795.00 41670.00 40.00 -0.187900E+02 0.537580E+01 1008 | 3797.00 41570.00 40.00 0.168000E+02 0.533600E+01 1009 | 3796.10 41470.00 40.00 0.587100E+02 0.617420E+01 1010 | 3796.00 41370.00 40.00 0.100860E+03 0.701720E+01 1011 | 3797.00 41270.00 40.00 0.153560E+03 0.807120E+01 1012 | 3798.00 41170.00 40.00 0.209100E+03 0.918200E+01 1013 | 3799.30 41070.00 40.00 0.262520E+03 0.102504E+02 1014 | 3801.70 40970.00 40.00 0.327590E+03 0.115518E+02 1015 | 3803.10 40870.00 40.00 0.403430E+03 0.130686E+02 1016 | 3805.00 40770.00 40.00 0.526960E+03 0.155392E+02 1017 | 3805.00 40670.00 40.00 0.647460E+03 0.179492E+02 1018 | 3805.00 40570.00 40.00 0.625150E+03 0.175030E+02 1019 | 3805.80 40470.00 40.00 0.495000E+03 0.149000E+02 1020 | 3806.00 40370.00 40.00 0.295020E+03 0.109004E+02 1021 | 3807.00 40270.00 40.00 0.216860E+03 0.933720E+01 1022 | 3807.00 40170.00 40.00 0.303250E+03 0.110650E+02 1023 | 3807.90 40070.00 40.00 0.447730E+03 0.139546E+02 1024 | 3877.00 40043.00 40.00 0.555400E+03 0.161080E+02 1025 | 3877.00 40143.00 40.00 0.370580E+03 0.124116E+02 1026 | 3877.90 40243.00 40.00 0.202470E+03 0.904940E+01 1027 | 3882.50 40343.00 40.00 0.120300E+03 0.740600E+01 1028 | 3887.20 40443.00 40.00 0.128640E+03 0.757280E+01 1029 | 3889.00 40542.00 40.00 0.203470E+03 0.906940E+01 1030 | 3889.00 40642.00 40.00 0.333250E+03 0.116650E+02 1031 | 3889.00 40742.00 40.00 0.403720E+03 0.130744E+02 1032 | 3889.90 40842.00 40.00 0.353460E+03 0.120692E+02 1033 | 3891.30 40942.00 40.00 0.282800E+03 0.106560E+02 1034 | 3893.10 41042.00 40.00 0.231990E+03 0.963980E+01 1035 | 3894.00 41142.00 40.00 0.187980E+03 0.875960E+01 1036 | 3895.00 41242.00 40.00 0.142490E+03 0.784980E+01 1037 | 3895.00 41342.00 40.00 0.986100E+02 0.697220E+01 1038 | 3894.00 41442.00 40.00 0.785200E+02 0.657040E+01 1039 | 3894.00 41542.00 40.00 0.257600E+02 0.551520E+01 1040 | 3894.00 41642.00 40.00 -0.342000E+01 0.506840E+01 1041 | 3892.30 41742.00 40.00 -0.331600E+02 0.566320E+01 1042 | 3888.10 41842.00 40.00 -0.530600E+02 0.606120E+01 1043 | 3883.70 41942.00 40.00 -0.686600E+02 0.637320E+01 1044 | 3878.20 42042.00 40.00 -0.822700E+02 0.664540E+01 1045 | 3875.90 42142.00 40.00 -0.970600E+02 0.694120E+01 1046 | 3876.00 42242.00 40.00 -0.111210E+03 0.722420E+01 1047 | 3879.30 42342.00 40.00 -0.103140E+03 0.706280E+01 1048 | 3883.30 42442.00 40.00 0.150200E+03 0.800400E+01 1049 | 3888.40 42542.00 40.00 -0.121790E+03 0.743580E+01 1050 | 3893.70 42642.00 40.00 -0.990000E+02 0.698000E+01 1051 | 3900.60 42741.00 40.00 -0.404400E+02 0.580880E+01 1052 | 3904.40 42841.00 40.00 -0.200700E+02 0.540140E+01 1053 | 3905.00 42941.00 40.00 -0.627000E+01 0.512540E+01 1054 | 4001.00 42965.00 40.00 0.366000E+02 0.573200E+01 1055 | 3997.60 42865.00 40.00 0.355600E+02 0.571120E+01 1056 | 3994.50 42765.00 40.00 -0.329000E+01 0.506580E+01 1057 | 3988.40 42665.00 40.00 -0.577900E+02 0.615580E+01 1058 | 3982.50 42565.00 40.00 -0.915700E+02 0.683140E+01 1059 | 3978.90 42465.00 40.00 -0.100460E+03 0.700920E+01 1060 | 3976.90 42365.00 40.00 -0.957100E+02 0.691420E+01 1061 | 3977.00 42265.00 40.00 -0.962200E+02 0.692440E+01 1062 | 3978.00 42165.00 40.00 -0.854200E+02 0.670840E+01 1063 | 3978.00 42065.00 40.00 -0.718600E+02 0.643720E+01 1064 | 3978.00 41965.00 40.00 -0.569400E+02 0.613880E+01 1065 | 3978.00 41865.00 40.00 -0.423200E+02 0.584640E+01 1066 | 3978.30 41765.00 40.00 -0.279000E+02 0.555800E+01 1067 | 3980.60 41665.00 40.00 -0.723000E+01 0.514460E+01 1068 | 3983.90 41565.00 40.00 0.184600E+02 0.536920E+01 1069 | 3986.00 41465.00 40.00 0.582700E+02 0.616540E+01 1070 | 3986.00 41365.00 40.00 0.771000E+02 0.654200E+01 1071 | 3986.90 41265.00 40.00 0.101640E+03 0.703280E+01 1072 | 3988.90 41165.00 40.00 0.134750E+03 0.769500E+01 1073 | 3991.40 41066.00 40.00 0.163840E+03 0.827680E+01 1074 | 3996.10 40966.00 40.00 0.180010E+03 0.860020E+01 1075 | 4000.80 40866.00 40.00 0.182080E+03 0.864160E+01 1076 | 4005.50 40766.00 40.00 0.173530E+03 0.847060E+01 1077 | 4009.30 40666.00 40.00 0.132590E+03 0.765180E+01 1078 | 4013.00 40566.00 40.00 0.546100E+02 0.609220E+01 1079 | 4015.30 40466.00 40.00 -0.125700E+02 0.525140E+01 1080 | 4017.00 40366.00 40.00 -0.469100E+02 0.593820E+01 1081 | 4018.00 40266.00 40.00 -0.217900E+02 0.543580E+01 1082 | 4018.00 40166.00 40.00 0.139110E+03 0.778220E+01 1083 | 4016.90 40066.00 40.00 0.450360E+03 0.140072E+02 1084 | 4094.00 40042.00 40.00 0.794700E+03 0.208940E+02 1085 | 4091.10 40142.00 40.00 0.281710E+03 0.106342E+02 1086 | 4091.00 40242.00 40.00 -0.813700E+02 0.662740E+01 1087 | 4091.00 40342.00 40.00 -0.162200E+03 0.824400E+01 1088 | 4091.00 40442.00 40.00 -0.948800E+02 0.689760E+01 1089 | 4091.00 40542.00 40.00 -0.317000E+02 0.563400E+01 1090 | 4091.00 40642.00 40.00 0.353400E+02 0.570680E+01 1091 | 4091.00 40742.00 40.00 0.907500E+02 0.681500E+01 1092 | 4091.00 40842.00 40.00 0.122920E+03 0.745840E+01 1093 | 4091.00 40942.00 40.00 0.131200E+03 0.762400E+01 1094 | 4091.00 41042.00 40.00 0.131730E+03 0.763460E+01 1095 | 4091.00 41142.00 40.00 0.128780E+03 0.757560E+01 1096 | 4091.00 41242.00 40.00 0.116710E+03 0.733420E+01 1097 | 4091.60 41342.00 40.00 0.966400E+02 0.693280E+01 1098 | 4092.00 41442.00 40.00 0.836200E+02 0.667240E+01 1099 | 4092.00 41542.00 40.00 0.589100E+02 0.617820E+01 1100 | 4092.00 41642.00 40.00 0.365000E+02 0.573000E+01 1101 | 4091.00 41742.00 40.00 0.171400E+02 0.534280E+01 1102 | 4090.00 41842.00 40.00 0.287000E+01 0.505740E+01 1103 | 4090.00 41942.00 40.00 -0.114800E+02 0.522960E+01 1104 | 4090.00 42042.00 40.00 -0.291200E+02 0.558240E+01 1105 | 4090.90 42142.00 40.00 -0.464300E+02 0.592860E+01 1106 | 4093.10 42242.00 40.00 -0.594900E+02 0.618980E+01 1107 | 4096.00 42342.00 40.00 -0.697800E+02 0.639560E+01 1108 | 4099.40 42441.00 40.00 -0.769200E+02 0.653840E+01 1109 | 4102.00 42541.00 40.00 -0.740400E+02 0.648080E+01 1110 | 4104.40 42641.00 40.00 -0.589200E+02 0.617840E+01 1111 | 4105.90 42741.00 40.00 -0.775000E+01 0.515500E+01 1112 | 4107.30 42841.00 40.00 0.841000E+02 0.668200E+01 1113 | 4108.90 42941.00 40.00 0.905500E+02 0.681100E+01 1114 | 4203.00 42959.00 40.00 0.850000E+02 0.670000E+01 1115 | 4204.30 42859.00 40.00 0.751200E+02 0.650240E+01 1116 | 4205.80 42759.00 40.00 0.204300E+02 0.540860E+01 1117 | 4207.00 42659.00 40.00 -0.289100E+02 0.557820E+01 1118 | 4208.00 42559.00 40.00 -0.482900E+02 0.596580E+01 1119 | 4209.00 42459.00 40.00 -0.498900E+02 0.599780E+01 1120 | 4207.60 42359.00 40.00 -0.445400E+02 0.589080E+01 1121 | 4205.20 42259.00 40.00 -0.337900E+02 0.567580E+01 1122 | 4202.80 42159.00 40.00 -0.223500E+02 0.544700E+01 1123 | 4200.00 42059.00 40.00 -0.649000E+01 0.512980E+01 1124 | 4199.00 41959.00 40.00 0.123700E+02 0.524740E+01 1125 | 4198.00 41859.00 40.00 0.326400E+02 0.565280E+01 1126 | 4196.80 41759.00 40.00 0.474000E+02 0.594800E+01 1127 | 4198.00 41659.00 40.00 0.629600E+02 0.625920E+01 1128 | 4200.10 41559.00 40.00 0.792000E+02 0.658400E+01 1129 | 4203.80 41459.00 40.00 0.936700E+02 0.687340E+01 1130 | 4207.30 41359.00 40.00 0.106030E+03 0.712060E+01 1131 | 4209.80 41259.00 40.00 0.115650E+03 0.731300E+01 1132 | 4212.00 41159.00 40.00 0.117340E+03 0.734680E+01 1133 | 4213.00 41059.00 40.00 0.110130E+03 0.720260E+01 1134 | 4214.10 40959.00 40.00 0.969500E+02 0.693900E+01 1135 | 4213.00 40859.00 40.00 0.771000E+02 0.654200E+01 1136 | 4210.00 40759.00 40.00 0.386000E+02 0.577200E+01 1137 | 4207.60 40659.00 40.00 -0.267500E+02 0.553500E+01 1138 | 4204.50 40559.00 40.00 -0.926200E+02 0.685240E+01 1139 | 4202.20 40459.00 40.00 -0.153640E+03 0.807280E+01 1140 | 4201.00 40359.00 40.00 -0.279870E+03 0.105974E+02 1141 | 4200.90 40259.00 40.00 -0.320550E+03 0.114110E+02 1142 | 4202.60 40159.00 40.00 0.716440E+03 0.193288E+02 1143 | 4204.00 40059.00 40.00 0.228653E+04 0.507306E+02 1144 | 4295.00 40030.00 40.00 0.395190E+04 0.840380E+02 1145 | 4293.00 40130.00 40.00 0.492219E+04 0.103444E+03 1146 | 4291.00 40230.00 40.00 0.302771E+04 0.655542E+02 1147 | 4289.40 40330.00 40.00 0.296320E+03 0.109264E+02 1148 | 4287.90 40430.00 40.00 -0.164830E+03 0.829660E+01 1149 | 4287.00 40530.00 40.00 -0.179250E+03 0.858500E+01 1150 | 4285.50 40630.00 40.00 -0.123930E+03 0.747860E+01 1151 | 4284.00 40730.00 40.00 -0.404400E+02 0.580880E+01 1152 | 4284.50 40830.00 40.00 0.343100E+02 0.568620E+01 1153 | 4286.10 40930.00 40.00 0.760400E+02 0.652080E+01 1154 | 4288.80 41030.00 40.00 0.978400E+02 0.695680E+01 1155 | 4292.40 41129.00 40.00 0.118730E+03 0.737460E+01 1156 | 4298.60 41229.00 40.00 0.130830E+03 0.761660E+01 1157 | 4302.00 41329.00 40.00 0.129740E+03 0.759480E+01 1158 | 4303.00 41429.00 40.00 0.133540E+03 0.767080E+01 1159 | 4302.50 41529.00 40.00 0.139440E+03 0.778880E+01 1160 | 4298.70 41629.00 40.00 0.111250E+03 0.722500E+01 1161 | 4294.30 41729.00 40.00 0.913200E+02 0.682640E+01 1162 | 4290.50 41829.00 40.00 0.700200E+02 0.640040E+01 1163 | 4286.90 41929.00 40.00 0.529700E+02 0.605940E+01 1164 | 4284.90 42029.00 40.00 0.342600E+02 0.568520E+01 1165 | 4285.10 42129.00 40.00 0.152500E+02 0.530500E+01 1166 | 4286.50 42229.00 40.00 -0.250000E+01 0.505000E+01 1167 | 4291.10 42329.00 40.00 -0.139300E+02 0.527860E+01 1168 | 4296.40 42429.00 40.00 -0.250400E+02 0.550080E+01 1169 | 4302.80 42528.00 40.00 -0.313000E+02 0.562600E+01 1170 | 4308.70 42628.00 40.00 -0.295900E+02 0.559180E+01 1171 | 4313.50 42728.00 40.00 -0.222100E+02 0.544420E+01 1172 | 4318.20 42828.00 40.00 0.343000E+01 0.506860E+01 1173 | 4320.60 42928.00 40.00 0.177200E+02 0.535440E+01 1174 | 4411.00 42966.00 40.00 0.970000E+01 0.519400E+01 1175 | 4414.50 42866.00 40.00 -0.160000E+00 0.500320E+01 1176 | 4413.70 42766.00 40.00 -0.107700E+02 0.521540E+01 1177 | 4409.30 42666.00 40.00 -0.135200E+02 0.527040E+01 1178 | 4405.20 42566.00 40.00 -0.928000E+01 0.518560E+01 1179 | 4401.00 42466.00 40.00 0.000000E+00 0.500000E+01 1180 | 4395.80 42366.00 40.00 0.788000E+01 0.515760E+01 1181 | 4392.00 42266.00 40.00 0.209000E+02 0.541800E+01 1182 | 4391.10 42166.00 40.00 0.384900E+02 0.576980E+01 1183 | 4392.00 42066.00 40.00 0.581800E+02 0.616360E+01 1184 | 4392.80 41966.00 40.00 0.798200E+02 0.659640E+01 1185 | 4393.10 41866.00 40.00 0.104000E+03 0.708000E+01 1186 | 4395.50 41766.00 40.00 0.129880E+03 0.759760E+01 1187 | 4398.10 41666.00 40.00 0.152920E+03 0.805840E+01 1188 | 4400.00 41567.00 40.00 0.180530E+03 0.861060E+01 1189 | 4400.00 41467.00 40.00 0.199070E+03 0.898140E+01 1190 | 4400.00 41367.00 40.00 0.194880E+03 0.889760E+01 1191 | 4400.00 41267.00 40.00 0.189940E+03 0.879880E+01 1192 | 4400.60 41167.00 40.00 0.183300E+03 0.866600E+01 1193 | 4402.00 41067.00 40.00 0.158790E+03 0.817580E+01 1194 | 4403.00 40967.00 40.00 0.133170E+03 0.766340E+01 1195 | 4403.10 40867.00 40.00 0.981800E+02 0.696360E+01 1196 | 4401.90 40767.00 40.00 0.276000E+01 0.505520E+01 1197 | 4400.60 40667.00 40.00 -0.111420E+03 0.722840E+01 1198 | 4399.10 40567.00 40.00 -0.176310E+03 0.852620E+01 1199 | 4399.00 40467.00 40.00 -0.519000E+02 0.603800E+01 1200 | 4399.00 40367.00 40.00 0.136289E+04 0.322578E+02 1201 | 4399.00 40267.00 40.00 0.147296E+04 0.344592E+02 1202 | 4399.50 40167.00 40.00 0.450201E+04 0.950402E+02 1203 | 4399.00 40067.00 40.00 0.402736E+04 0.855472E+02 1204 | 4481.00 40050.00 40.00 0.133900E+03 0.767800E+01 1205 | 4478.00 40150.00 40.00 0.522064E+04 0.109413E+03 1206 | 4475.90 40250.00 40.00 0.465382E+04 0.980764E+02 1207 | 4475.00 40350.00 40.00 0.360896E+04 0.771792E+02 1208 | 4476.00 40450.00 40.00 0.849790E+03 0.219958E+02 1209 | 4477.20 40550.00 40.00 -0.149590E+03 0.799180E+01 1210 | 4481.30 40650.00 40.00 -0.463700E+02 0.592740E+01 1211 | 4485.20 40749.00 40.00 -0.158400E+02 0.531680E+01 1212 | 4488.00 40849.00 40.00 0.816100E+02 0.663220E+01 1213 | 4489.90 40949.00 40.00 0.152230E+03 0.804460E+01 1214 | 4491.00 41049.00 40.00 0.187150E+03 0.874300E+01 1215 | 4491.00 41149.00 40.00 0.226210E+03 0.952420E+01 1216 | 4491.00 41249.00 40.00 0.242800E+03 0.985600E+01 1217 | 4491.00 41349.00 40.00 0.253700E+03 0.100740E+02 1218 | 4491.90 41449.00 40.00 0.270310E+03 0.104062E+02 1219 | 4494.60 41549.00 40.00 0.238800E+03 0.977600E+01 1220 | 4498.40 41649.00 40.00 0.229150E+03 0.958300E+01 1221 | 4501.10 41749.00 40.00 0.205790E+03 0.911580E+01 1222 | 4503.40 41849.00 40.00 0.178100E+03 0.856200E+01 1223 | 4504.60 41949.00 40.00 0.149380E+03 0.798760E+01 1224 | 4504.30 42049.00 40.00 0.120250E+03 0.740500E+01 1225 | 4502.00 42149.00 40.00 0.913300E+02 0.682660E+01 1226 | 4497.50 42249.00 40.00 0.632000E+02 0.626400E+01 1227 | 4492.50 42349.00 40.00 0.413500E+02 0.582700E+01 1228 | 4490.00 42449.00 40.00 0.281700E+02 0.556340E+01 1229 | 4490.00 42549.00 40.00 0.110500E+02 0.522100E+01 1230 | 4491.10 42649.00 40.00 -0.246000E+01 0.504920E+01 1231 | 4493.50 42749.00 40.00 -0.681000E+01 0.513620E+01 1232 | 4496.90 42849.00 40.00 -0.490000E+01 0.509800E+01 1233 | 507.00 39957.00 40.00 0.110830E+04 0.271660E+02 1234 | 509.30 39857.00 40.00 0.987840E+03 0.247568E+02 1235 | 511.80 39757.00 40.00 0.894130E+03 0.228826E+02 1236 | 513.00 39657.00 40.00 0.738730E+03 0.197746E+02 1237 | 513.00 39557.00 40.00 0.122473E+04 0.294946E+02 1238 | 513.00 39457.00 40.00 0.223945E+04 0.497890E+02 1239 | 511.00 39357.00 40.00 0.173811E+04 0.397622E+02 1240 | 509.00 39257.00 40.00 0.524200E+03 0.154840E+02 1241 | 507.00 39157.00 40.00 0.498460E+03 0.149692E+02 1242 | 507.00 39057.00 40.00 0.539920E+03 0.157984E+02 1243 | 604.00 39030.00 40.00 0.518100E+03 0.153620E+02 1244 | 607.00 39130.00 40.00 0.528380E+03 0.155676E+02 1245 | 605.60 39230.00 40.00 0.389740E+03 0.127948E+02 1246 | 603.00 39330.00 40.00 0.371470E+03 0.124294E+02 1247 | 601.00 39430.00 40.00 -0.308070E+03 0.111614E+02 1248 | 595.80 39529.00 40.00 0.183262E+04 0.416524E+02 1249 | 588.00 39629.00 40.00 0.883650E+03 0.226730E+02 1250 | 579.30 39729.00 40.00 0.851450E+03 0.220290E+02 1251 | 571.20 39828.00 40.00 0.826590E+03 0.215318E+02 1252 | 566.60 39928.00 40.00 0.104964E+04 0.259928E+02 1253 | 706.00 39949.00 40.00 0.100980E+04 0.251960E+02 1254 | 710.20 39849.00 40.00 0.904890E+03 0.230978E+02 1255 | 714.60 39749.00 40.00 0.839320E+03 0.217864E+02 1256 | 719.10 39649.00 40.00 0.179080E+04 0.408160E+02 1257 | 721.10 39549.00 40.00 0.136940E+04 0.323880E+02 1258 | 720.60 39449.00 40.00 0.250730E+03 0.100146E+02 1259 | 718.90 39349.00 40.00 0.194840E+03 0.889680E+01 1260 | 717.80 39249.00 40.00 0.277070E+03 0.105414E+02 1261 | 716.00 39149.00 40.00 0.395290E+03 0.129058E+02 1262 | 714.10 39049.00 40.00 0.434650E+03 0.136930E+02 1263 | 806.00 39028.00 40.00 0.330300E+03 0.116060E+02 1264 | 803.00 39128.00 40.00 0.390590E+03 0.128118E+02 1265 | 799.00 39228.00 40.00 0.265140E+03 0.103028E+02 1266 | 792.70 39327.00 40.00 0.191100E+03 0.882200E+01 1267 | 785.00 39427.00 40.00 0.179210E+03 0.858420E+01 1268 | 777.60 39527.00 40.00 0.454180E+03 0.140836E+02 1269 | 773.30 39627.00 40.00 0.185547E+04 0.421094E+02 1270 | 770.00 39727.00 40.00 0.163368E+04 0.376736E+02 1271 | 764.80 39827.00 40.00 0.788700E+03 0.207740E+02 1272 | 761.00 39926.00 40.00 0.853530E+03 0.220706E+02 1273 | 902.00 39959.00 40.00 0.716100E+03 0.193220E+02 1274 | 910.60 39859.00 40.00 0.705340E+03 0.191068E+02 1275 | 920.00 39760.00 40.00 0.107263E+04 0.264526E+02 1276 | 930.50 39660.00 40.00 0.300640E+03 0.110128E+02 1277 | 939.80 39561.00 40.00 0.279220E+03 0.105844E+02 1278 | 944.30 39461.00 40.00 0.186810E+03 0.873620E+01 1279 | 946.00 39361.00 40.00 0.177350E+03 0.854700E+01 1280 | 946.00 39261.00 40.00 0.190980E+03 0.881960E+01 1281 | 942.80 39161.00 40.00 0.220330E+03 0.940660E+01 1282 | 934.10 39061.00 40.00 0.230600E+03 0.961200E+01 1283 | 988.00 39029.00 40.00 0.165200E+03 0.830400E+01 1284 | 985.10 39129.00 40.00 0.210930E+03 0.921860E+01 1285 | 985.00 39229.00 40.00 0.189640E+03 0.879280E+01 1286 | 985.90 39329.00 40.00 0.177780E+03 0.855560E+01 1287 | 988.80 39429.00 40.00 0.177560E+03 0.855120E+01 1288 | 992.00 39529.00 40.00 0.164510E+03 0.829020E+01 1289 | 994.40 39629.00 40.00 0.178790E+03 0.857580E+01 1290 | 996.00 39729.00 40.00 0.499350E+03 0.149870E+02 1291 | 994.90 39829.00 40.00 0.764480E+03 0.202896E+02 1292 | 992.60 39928.00 40.00 0.622110E+03 0.174422E+02 1293 | 1117.00 39934.00 40.00 0.667100E+03 0.183420E+02 1294 | 1119.00 39834.00 40.00 0.890820E+03 0.228164E+02 1295 | 1119.00 39734.00 40.00 0.481020E+03 0.146204E+02 1296 | 1119.20 39634.00 40.00 0.189130E+03 0.878260E+01 1297 | 1120.00 39534.00 40.00 0.156030E+03 0.812060E+01 1298 | 1120.00 39434.00 40.00 0.163180E+03 0.826360E+01 1299 | 1119.00 39334.00 40.00 0.168170E+03 0.836340E+01 1300 | 1118.00 39234.00 40.00 0.193410E+03 0.886820E+01 1301 | 1114.50 39134.00 40.00 0.224360E+03 0.948720E+01 1302 | 1108.10 39034.00 40.00 0.159340E+03 0.818680E+01 1303 | 1191.00 39033.00 40.00 0.111900E+03 0.723800E+01 1304 | 1185.40 39133.00 40.00 0.174930E+03 0.849860E+01 1305 | 1180.00 39232.00 40.00 0.228660E+03 0.957320E+01 1306 | 1174.00 39332.00 40.00 0.181970E+03 0.863940E+01 1307 | 1167.60 39432.00 40.00 0.161610E+03 0.823220E+01 1308 | 1162.00 39532.00 40.00 0.159550E+03 0.819100E+01 1309 | 1162.60 39632.00 40.00 0.161300E+03 0.822600E+01 1310 | 1167.60 39732.00 40.00 0.271660E+03 0.104332E+02 1311 | 1173.10 39832.00 40.00 0.732350E+03 0.196470E+02 1312 | 1178.10 39931.00 40.00 0.674880E+03 0.184976E+02 1313 | 1303.00 39936.00 40.00 0.661100E+03 0.182220E+02 1314 | 1299.30 39836.00 40.00 0.423900E+03 0.134780E+02 1315 | 1295.50 39736.00 40.00 0.189930E+03 0.879860E+01 1316 | 1292.90 39636.00 40.00 0.152390E+03 0.804780E+01 1317 | 1291.80 39536.00 40.00 0.160250E+03 0.820500E+01 1318 | 1293.00 39436.00 40.00 0.173030E+03 0.846060E+01 1319 | 1294.00 39336.00 40.00 0.216990E+03 0.933980E+01 1320 | 1294.00 39236.00 40.00 0.248090E+03 0.996180E+01 1321 | 1294.00 39136.00 40.00 0.151040E+03 0.802080E+01 1322 | 1292.80 39036.00 40.00 0.117720E+03 0.735440E+01 1323 | 1388.00 39031.00 40.00 0.599000E+02 0.619800E+01 1324 | 1387.00 39131.00 40.00 0.127920E+03 0.755840E+01 1325 | 1387.00 39231.00 40.00 0.172920E+03 0.845840E+01 1326 | 1386.00 39331.00 40.00 0.266710E+03 0.103342E+02 1327 | 1384.00 39431.00 40.00 0.189050E+03 0.878100E+01 1328 | 1381.20 39531.00 40.00 0.152670E+03 0.805340E+01 1329 | 1378.10 39631.00 40.00 0.154020E+03 0.808040E+01 1330 | 1374.00 39731.00 40.00 0.159050E+03 0.818100E+01 1331 | 1373.00 39830.00 40.00 0.222670E+03 0.945340E+01 1332 | 1373.70 39930.00 40.00 0.581830E+03 0.166366E+02 1333 | 1522.00 39945.00 40.00 0.293400E+03 0.108680E+02 1334 | 1524.10 39845.00 40.00 0.189430E+03 0.878860E+01 1335 | 1526.40 39745.00 40.00 0.163070E+03 0.826140E+01 1336 | 1529.00 39645.00 40.00 0.170260E+03 0.840520E+01 1337 | 1532.50 39545.00 40.00 0.157150E+03 0.814300E+01 1338 | 1535.60 39445.00 40.00 0.199780E+03 0.899560E+01 1339 | 1537.00 39345.00 40.00 0.225270E+03 0.950540E+01 1340 | 1536.00 39245.00 40.00 0.101280E+03 0.702560E+01 1341 | 1534.30 39145.00 40.00 0.783000E+02 0.656600E+01 1342 | 1530.80 39045.00 40.00 0.466000E+02 0.593200E+01 1343 | 1591.00 39035.00 40.00 -0.490000E+01 0.509800E+01 1344 | 1587.50 39135.00 40.00 0.375700E+02 0.575140E+01 1345 | 1584.90 39235.00 40.00 0.702600E+02 0.640520E+01 1346 | 1583.00 39335.00 40.00 0.121900E+03 0.743800E+01 1347 | 1581.30 39435.00 40.00 0.220420E+03 0.940840E+01 1348 | 1579.80 39535.00 40.00 0.180980E+03 0.861960E+01 1349 | 1575.50 39634.00 40.00 0.167620E+03 0.835240E+01 1350 | 1570.70 39734.00 40.00 0.170860E+03 0.841720E+01 1351 | 1568.80 39834.00 40.00 0.161270E+03 0.822540E+01 1352 | 1569.30 39934.00 40.00 0.211640E+03 0.923280E+01 1353 | 1701.00 39934.00 40.00 0.199000E+03 0.898000E+01 1354 | 1703.30 39834.00 40.00 0.160710E+03 0.821420E+01 1355 | 1706.70 39734.00 40.00 0.173180E+03 0.846360E+01 1356 | 1710.40 39634.00 40.00 0.185290E+03 0.870580E+01 1357 | 1715.10 39534.00 40.00 0.233350E+03 0.966700E+01 1358 | 1719.70 39434.00 40.00 0.188280E+03 0.876560E+01 1359 | 1725.20 39334.00 40.00 0.112980E+03 0.725960E+01 1360 | 1729.30 39234.00 40.00 0.362900E+02 0.572580E+01 1361 | 1732.00 39134.00 40.00 -0.721000E+01 0.514420E+01 1362 | 1733.00 39034.00 40.00 -0.546400E+02 0.609280E+01 1363 | 1794.00 39029.00 40.00 -0.124300E+03 0.748600E+01 1364 | 1795.00 39129.00 40.00 -0.409600E+02 0.581920E+01 1365 | 1794.00 39229.00 40.00 0.562000E+01 0.511240E+01 1366 | 1793.20 39329.00 40.00 0.532100E+02 0.606420E+01 1367 | 1793.00 39429.00 40.00 0.159660E+03 0.819320E+01 1368 | 1794.00 39529.00 40.00 0.254630E+03 0.100926E+02 1369 | 1795.10 39629.00 40.00 0.217520E+03 0.935040E+01 1370 | 1794.00 39729.00 40.00 0.168990E+03 0.837980E+01 1371 | 1791.20 39829.00 40.00 0.154680E+03 0.809360E+01 1372 | 1790.00 39929.00 40.00 0.157090E+03 0.814180E+01 1373 | 1902.00 39938.00 40.00 0.144300E+03 0.788600E+01 1374 | 1902.50 39838.00 40.00 0.139240E+03 0.778480E+01 1375 | 1907.30 39738.00 40.00 0.163020E+03 0.826040E+01 1376 | 1914.10 39638.00 40.00 0.188750E+03 0.877500E+01 1377 | 1919.80 39538.00 40.00 0.230210E+03 0.960420E+01 1378 | 1926.80 39439.00 40.00 0.103480E+03 0.706960E+01 1379 | 1933.70 39339.00 40.00 0.375700E+02 0.575140E+01 1380 | 1939.30 39239.00 40.00 -0.913000E+01 0.518260E+01 1381 | 1941.90 39139.00 40.00 -0.105990E+03 0.711980E+01 1382 | 1942.60 39039.00 40.00 -0.247050E+03 0.994100E+01 1383 | 1993.00 39032.00 40.00 -0.396700E+03 0.129340E+02 1384 | 1994.80 39132.00 40.00 -0.187360E+03 0.874720E+01 1385 | 1995.30 39232.00 40.00 -0.482900E+02 0.596580E+01 1386 | 1995.90 39332.00 40.00 0.162400E+02 0.532480E+01 1387 | 1994.00 39432.00 40.00 0.595100E+02 0.619020E+01 1388 | 1992.20 39532.00 40.00 0.153240E+03 0.806480E+01 1389 | 1989.30 39632.00 40.00 0.198860E+03 0.897720E+01 1390 | 1984.60 39731.00 40.00 0.165840E+03 0.831680E+01 1391 | 1976.60 39831.00 40.00 0.149240E+03 0.798480E+01 1392 | 1971.10 39931.00 40.00 0.125660E+03 0.751320E+01 1393 | 2110.00 39954.00 40.00 0.107400E+03 0.714800E+01 1394 | 2118.00 39854.00 40.00 0.135280E+03 0.770560E+01 1395 | 2128.00 39755.00 40.00 0.160840E+03 0.821680E+01 1396 | 2137.90 39655.00 40.00 0.184380E+03 0.868760E+01 1397 | 2146.90 39556.00 40.00 0.161190E+03 0.822380E+01 1398 | 2153.50 39456.00 40.00 0.355200E+02 0.571040E+01 1399 | 2157.00 39356.00 40.00 -0.795000E+01 0.515900E+01 1400 | 2153.70 39256.00 40.00 -0.145400E+03 0.790800E+01 1401 | 2138.20 39157.00 40.00 -0.251750E+03 0.100350E+02 1402 | 2120.10 39059.00 40.00 -0.347240E+03 0.119448E+02 1403 | 2202.00 39032.00 40.00 0.525300E+03 0.155060E+02 1404 | 2200.20 39132.00 40.00 -0.173360E+03 0.846720E+01 1405 | 2198.50 39232.00 40.00 -0.181670E+03 0.863340E+01 1406 | 2196.80 39332.00 40.00 -0.876600E+02 0.675320E+01 1407 | 2193.00 39432.00 40.00 0.646000E+01 0.512920E+01 1408 | 2189.30 39532.00 40.00 0.612800E+02 0.622560E+01 1409 | 2185.10 39631.00 40.00 0.110160E+03 0.720320E+01 1410 | 2180.30 39731.00 40.00 0.186720E+03 0.873440E+01 1411 | 2173.90 39831.00 40.00 0.161380E+03 0.822760E+01 1412 | 2169.00 39931.00 40.00 0.115910E+03 0.731820E+01 1413 | 2310.00 39968.00 40.00 0.970000E+02 0.694000E+01 1414 | 2317.00 39868.00 40.00 0.114020E+03 0.728040E+01 1415 | 2326.10 39768.00 40.00 0.164780E+03 0.829560E+01 1416 | 2334.50 39669.00 40.00 0.126200E+03 0.752400E+01 1417 | 2342.00 39569.00 40.00 0.353100E+02 0.570620E+01 1418 | 2346.00 39469.00 40.00 -0.195300E+02 0.539060E+01 1419 | 2345.90 39369.00 40.00 -0.998200E+02 0.699640E+01 1420 | 2341.60 39269.00 40.00 -0.273080E+03 0.104616E+02 1421 | 2332.50 39170.00 40.00 -0.379180E+03 0.125836E+02 1422 | 2321.50 39070.00 40.00 -0.186250E+03 0.872500E+01 1423 | 2391.00 39028.00 40.00 0.478000E+02 0.595600E+01 1424 | 2398.00 39127.00 40.00 -0.180100E+03 0.860200E+01 1425 | 2403.80 39227.00 40.00 -0.398010E+03 0.129602E+02 1426 | 2407.80 39327.00 40.00 -0.235020E+03 0.970040E+01 1427 | 2410.60 39427.00 40.00 -0.907400E+02 0.681480E+01 1428 | 2411.00 39527.00 40.00 -0.357000E+01 0.507140E+01 1429 | 2410.00 39627.00 40.00 0.393400E+02 0.578680E+01 1430 | 2409.00 39727.00 40.00 0.132170E+03 0.764340E+01 1431 | 2408.30 39827.00 40.00 0.140950E+03 0.781900E+01 1432 | 2404.60 39927.00 40.00 0.964800E+02 0.692960E+01 1433 | 2485.00 39937.00 40.00 0.940000E+02 0.688000E+01 1434 | 2485.00 39837.00 40.00 0.103790E+03 0.707580E+01 1435 | 2488.20 39737.00 40.00 0.121500E+03 0.743000E+01 1436 | 2496.80 39637.00 40.00 0.304500E+02 0.560900E+01 1437 | 2505.50 39538.00 40.00 -0.268800E+02 0.553760E+01 1438 | 2511.80 39438.00 40.00 -0.106850E+03 0.713700E+01 1439 | 2518.40 39338.00 40.00 -0.184040E+03 0.868080E+01 1440 | 2525.00 39238.00 40.00 -0.328830E+03 0.115766E+02 1441 | 2529.30 39138.00 40.00 -0.166410E+03 0.832820E+01 1442 | 2532.00 39038.00 40.00 0.911500E+02 0.682300E+01 1443 | 2598.00 39034.00 40.00 0.924000E+03 0.234800E+02 1444 | 2595.10 39134.00 40.00 0.803800E+02 0.660760E+01 1445 | 2594.30 39234.00 40.00 -0.283950E+03 0.106790E+02 1446 | 2593.60 39334.00 40.00 -0.794900E+02 0.658980E+01 1447 | 2592.00 39434.00 40.00 -0.168320E+03 0.836640E+01 1448 | 2592.00 39534.00 40.00 -0.907200E+02 0.681440E+01 1449 | 2591.00 39634.00 40.00 -0.350500E+02 0.570100E+01 1450 | 2593.00 39734.00 40.00 0.480200E+02 0.596040E+01 1451 | 2592.80 39834.00 40.00 0.116540E+03 0.733080E+01 1452 | 2589.30 39934.00 40.00 0.757100E+02 0.651420E+01 1453 | 2688.00 39921.00 40.00 0.764000E+02 0.652800E+01 1454 | 2688.00 39821.00 40.00 0.915900E+02 0.683180E+01 1455 | 2688.00 39721.00 40.00 0.406000E+02 0.581200E+01 1456 | 2688.90 39621.00 40.00 -0.473400E+02 0.594680E+01 1457 | 2694.50 39521.00 40.00 -0.131850E+03 0.763700E+01 1458 | 2701.10 39421.00 40.00 -0.203380E+03 0.906760E+01 1459 | 2707.80 39321.00 40.00 -0.522600E+02 0.604520E+01 1460 | 2714.50 39222.00 40.00 -0.300300E+02 0.560060E+01 1461 | 2718.60 39122.00 40.00 0.804980E+03 0.210996E+02 1462 | 2798.00 39036.00 40.00 0.152390E+04 0.354780E+02 1463 | 2798.60 39136.00 40.00 0.116907E+04 0.283814E+02 1464 | 2800.00 39236.00 40.00 0.147810E+03 0.795620E+01 1465 | 2801.00 39336.00 40.00 -0.183000E+02 0.536600E+01 1466 | 2802.80 39436.00 40.00 -0.341200E+02 0.568240E+01 1467 | 2804.50 39536.00 40.00 -0.210260E+03 0.920520E+01 1468 | 2809.20 39636.00 40.00 -0.113110E+03 0.726220E+01 1469 | 2814.50 39735.00 40.00 -0.883000E+01 0.517660E+01 1470 | 2813.00 39835.00 40.00 0.752900E+02 0.650580E+01 1471 | 2806.90 39935.00 40.00 0.595000E+02 0.619000E+01 1472 | 2898.00 39933.00 40.00 0.505000E+02 0.601000E+01 1473 | 2898.10 39833.00 40.00 0.734800E+02 0.646960E+01 1474 | 2902.30 39733.00 40.00 -0.466100E+02 0.593220E+01 1475 | 2905.00 39633.00 40.00 -0.155910E+03 0.811820E+01 1476 | 2909.60 39533.00 40.00 -0.259340E+03 0.101868E+02 1477 | 2913.60 39433.00 40.00 0.526770E+03 0.155354E+02 1478 | 2918.10 39333.00 40.00 -0.154800E+02 0.530960E+01 1479 | 2921.80 39233.00 40.00 0.138880E+03 0.777760E+01 1480 | 2924.50 39133.00 40.00 0.836760E+03 0.217352E+02 1481 | 2927.10 39033.00 40.00 0.177186E+04 0.404372E+02 1482 | 3004.00 39030.00 40.00 0.177130E+04 0.404260E+02 1483 | 3002.00 39130.00 40.00 0.149970E+04 0.349940E+02 1484 | 3000.10 39230.00 40.00 0.639180E+03 0.177836E+02 1485 | 2999.00 39330.00 40.00 0.241260E+03 0.982520E+01 1486 | 2998.70 39430.00 40.00 0.114800E+03 0.729600E+01 1487 | 2998.00 39530.00 40.00 0.318210E+03 0.113642E+02 1488 | 2999.00 39630.00 40.00 -0.223370E+03 0.946740E+01 1489 | 3000.40 39730.00 40.00 -0.127330E+03 0.754660E+01 1490 | 3005.30 39829.00 40.00 -0.740000E+00 0.501480E+01 1491 | 3005.10 39929.00 40.00 -0.171500E+02 0.534300E+01 1492 | 3080.00 39964.00 40.00 0.489000E+02 0.597800E+01 1493 | 3082.40 39864.00 40.00 -0.141700E+02 0.528340E+01 1494 | 3084.00 39764.00 40.00 -0.129120E+03 0.758240E+01 1495 | 3085.60 39664.00 40.00 -0.106750E+03 0.713500E+01 1496 | 3088.50 39564.00 40.00 0.144300E+03 0.788600E+01 1497 | 3091.50 39464.00 40.00 0.104770E+03 0.709540E+01 1498 | 3096.30 39364.00 40.00 0.440750E+03 0.138150E+02 1499 | 3101.80 39264.00 40.00 0.103590E+04 0.257180E+02 1500 | 3106.80 39164.00 40.00 0.147216E+04 0.344432E+02 1501 | 3110.60 39064.00 40.00 0.251104E+04 0.552208E+02 1502 | 3173.00 39950.00 40.00 0.379900E+03 0.125980E+02 1503 | 3168.50 39850.00 40.00 0.255200E+02 0.551040E+01 1504 | 3165.20 39750.00 40.00 -0.172990E+03 0.845980E+01 1505 | 3162.80 39650.00 40.00 0.115770E+03 0.731540E+01 1506 | 3165.30 39550.00 40.00 0.167210E+03 0.834420E+01 1507 | 3171.50 39450.00 40.00 0.157420E+03 0.814840E+01 1508 | 3182.50 39351.00 40.00 0.493680E+03 0.148736E+02 1509 | 3191.70 39252.00 40.00 0.893670E+03 0.228734E+02 1510 | 3201.40 39152.00 40.00 0.120247E+04 0.290494E+02 1511 | 3210.60 39053.00 40.00 0.928200E+03 0.235640E+02 1512 | 3289.00 39032.00 40.00 -0.114000E+02 0.522800E+01 1513 | 3290.60 39132.00 40.00 0.108726E+04 0.267452E+02 1514 | 3291.00 39232.00 40.00 0.588460E+03 0.167692E+02 1515 | 3291.00 39332.00 40.00 0.142990E+03 0.785980E+01 1516 | 3292.00 39432.00 40.00 0.113600E+03 0.727200E+01 1517 | 3293.00 39532.00 40.00 0.372700E+02 0.574540E+01 1518 | 3294.00 39632.00 40.00 0.142500E+02 0.528500E+01 1519 | 3295.00 39732.00 40.00 -0.236080E+03 0.972160E+01 1520 | 3295.00 39832.00 40.00 -0.689300E+02 0.637860E+01 1521 | 3296.00 39932.00 40.00 0.273758E+04 0.597516E+02 1522 | 3399.00 39962.00 40.00 0.293250E+04 0.636500E+02 1523 | 3399.00 39862.00 40.00 0.147886E+04 0.345772E+02 1524 | 3398.10 39762.00 40.00 -0.863900E+02 0.672780E+01 1525 | 3397.90 39662.00 40.00 -0.186730E+03 0.873460E+01 1526 | 3397.60 39562.00 40.00 -0.745800E+02 0.649160E+01 1527 | 3401.10 39462.00 40.00 -0.420900E+02 0.584180E+01 1528 | 3406.30 39362.00 40.00 -0.859200E+02 0.671840E+01 1529 | 3413.90 39262.00 40.00 -0.595000E+02 0.619000E+01 1530 | 3420.10 39162.00 40.00 0.537430E+03 0.157486E+02 1531 | 3424.10 39063.00 40.00 -0.396500E+02 0.579300E+01 1532 | 3482.00 39042.00 40.00 -0.898000E+02 0.679600E+01 1533 | 3484.50 39142.00 40.00 -0.150500E+02 0.530100E+01 1534 | 3488.40 39242.00 40.00 0.107010E+03 0.714020E+01 1535 | 3492.20 39342.00 40.00 -0.118310E+03 0.736620E+01 1536 | 3496.90 39441.00 40.00 -0.962000E+02 0.692400E+01 1537 | 3502.60 39541.00 40.00 -0.790700E+02 0.658140E+01 1538 | 3504.50 39641.00 40.00 -0.138230E+03 0.776460E+01 1539 | 3502.00 39741.00 40.00 -0.211420E+03 0.922840E+01 1540 | 3502.50 39841.00 40.00 0.165293E+04 0.380586E+02 1541 | 3504.40 39941.00 40.00 0.146982E+04 0.343964E+02 1542 | 3611.00 39947.00 40.00 0.971100E+03 0.244220E+02 1543 | 3606.10 39847.00 40.00 0.953890E+03 0.240778E+02 1544 | 3602.00 39747.00 40.00 0.341400E+02 0.568280E+01 1545 | 3603.00 39647.00 40.00 -0.703100E+02 0.640620E+01 1546 | 3605.00 39547.00 40.00 -0.809800E+02 0.661960E+01 1547 | 3603.50 39447.00 40.00 -0.920100E+02 0.684020E+01 1548 | 3604.10 39347.00 40.00 -0.721900E+02 0.644380E+01 1549 | 3608.50 39247.00 40.00 0.978200E+02 0.695640E+01 1550 | 3613.20 39147.00 40.00 -0.269300E+02 0.553860E+01 1551 | 3618.00 39047.00 40.00 -0.520800E+02 0.604160E+01 1552 | 3697.00 39036.00 40.00 -0.448000E+02 0.589600E+01 1553 | 3702.00 39136.00 40.00 -0.539100E+02 0.607820E+01 1554 | 3706.00 39236.00 40.00 -0.349400E+02 0.569880E+01 1555 | 3708.80 39335.00 40.00 -0.798000E+01 0.515960E+01 1556 | 3709.00 39435.00 40.00 -0.858500E+02 0.671700E+01 1557 | 3704.50 39535.00 40.00 -0.431400E+02 0.586280E+01 1558 | 3700.40 39635.00 40.00 -0.758400E+02 0.651680E+01 1559 | 3694.70 39735.00 40.00 -0.274600E+02 0.554920E+01 1560 | 3688.30 39835.00 40.00 0.377770E+03 0.125554E+02 1561 | 3682.00 39935.00 40.00 0.106913E+04 0.263826E+02 1562 | 3810.00 39934.00 40.00 0.135050E+04 0.320100E+02 1563 | 3811.00 39834.00 40.00 0.405320E+03 0.131064E+02 1564 | 3811.00 39734.00 40.00 0.522800E+02 0.604560E+01 1565 | 3811.00 39634.00 40.00 -0.718000E+01 0.514360E+01 1566 | 3809.50 39534.00 40.00 -0.153900E+02 0.530780E+01 1567 | 3806.20 39434.00 40.00 -0.545400E+02 0.609080E+01 1568 | 3803.80 39334.00 40.00 -0.163000E+02 0.532600E+01 1569 | 3799.60 39234.00 40.00 -0.299400E+02 0.559880E+01 1570 | 3794.40 39134.00 40.00 -0.353700E+02 0.570740E+01 1571 | 3791.80 39034.00 40.00 -0.342100E+02 0.568420E+01 1572 | 3893.00 39034.00 40.00 -0.206000E+02 0.541200E+01 1573 | 3898.50 39134.00 40.00 -0.219100E+02 0.543820E+01 1574 | 3902.90 39233.00 40.00 -0.298100E+02 0.559620E+01 1575 | 3905.60 39333.00 40.00 0.173040E+03 0.846080E+01 1576 | 3905.80 39433.00 40.00 -0.120500E+02 0.524100E+01 1577 | 3905.00 39533.00 40.00 -0.219600E+02 0.543920E+01 1578 | 3903.30 39633.00 40.00 -0.160900E+02 0.532180E+01 1579 | 3897.90 39733.00 40.00 0.999000E+01 0.519980E+01 1580 | 3889.00 39833.00 40.00 0.134210E+03 0.768420E+01 1581 | 3881.80 39933.00 40.00 0.124345E+04 0.298690E+02 1582 | 4015.00 39955.00 40.00 0.135580E+04 0.321160E+02 1583 | 4009.60 39855.00 40.00 0.604420E+03 0.170884E+02 1584 | 4004.70 39755.00 40.00 0.155600E+02 0.531120E+01 1585 | 4000.30 39655.00 40.00 -0.352000E+01 0.507040E+01 1586 | 3988.50 39556.00 40.00 -0.674000E+01 0.513480E+01 1587 | 3979.60 39457.00 40.00 -0.562200E+02 0.612440E+01 1588 | 3976.00 39357.00 40.00 0.180250E+03 0.860500E+01 1589 | 3976.50 39257.00 40.00 -0.179300E+02 0.535860E+01 1590 | 3980.10 39157.00 40.00 -0.242100E+02 0.548420E+01 1591 | 3985.30 39057.00 40.00 -0.996000E+01 0.519920E+01 1592 | 4099.00 39035.00 40.00 0.171000E+02 0.534200E+01 1593 | 4104.00 39135.00 40.00 0.532000E+01 0.510640E+01 1594 | 4109.00 39234.00 40.00 -0.121100E+02 0.524220E+01 1595 | 4111.00 39334.00 40.00 -0.238300E+02 0.547660E+01 1596 | 4112.10 39434.00 40.00 0.490400E+02 0.598080E+01 1597 | 4114.00 39534.00 40.00 -0.589200E+02 0.617840E+01 1598 | 4110.20 39634.00 40.00 0.479300E+02 0.595860E+01 1599 | 4105.30 39734.00 40.00 0.647100E+02 0.629420E+01 1600 | 4101.50 39834.00 40.00 0.419900E+02 0.583980E+01 1601 | 4096.60 39934.00 40.00 0.854680E+03 0.220936E+02 1602 | 4205.00 39968.00 40.00 0.141780E+04 0.333560E+02 1603 | 4204.00 39868.00 40.00 0.407500E+02 0.581500E+01 1604 | 4203.90 39768.00 40.00 0.601000E+01 0.512020E+01 1605 | 4205.60 39668.00 40.00 0.185370E+03 0.870740E+01 1606 | 4206.00 39568.00 40.00 0.152970E+03 0.805940E+01 1607 | 4203.20 39468.00 40.00 -0.398000E+01 0.507960E+01 1608 | 4203.10 39368.00 40.00 0.160000E+01 0.503200E+01 1609 | 4204.80 39268.00 40.00 -0.245000E+02 0.549000E+01 1610 | 4206.00 39168.00 40.00 0.460000E+01 0.509200E+01 1611 | 4206.00 39068.00 40.00 0.270600E+02 0.554120E+01 1612 | 4320.00 39036.00 40.00 0.651000E+02 0.630200E+01 1613 | 4320.60 39136.00 40.00 0.458000E+02 0.591600E+01 1614 | 4321.00 39236.00 40.00 0.152000E+01 0.503040E+01 1615 | 4318.00 39336.00 40.00 -0.517700E+02 0.603540E+01 1616 | 4313.50 39436.00 40.00 0.389100E+02 0.577820E+01 1617 | 4308.50 39535.00 40.00 0.543650E+03 0.158730E+02 1618 | 4305.60 39635.00 40.00 0.996680E+03 0.249336E+02 1619 | 4304.00 39735.00 40.00 0.729500E+02 0.645900E+01 1620 | 4301.60 39835.00 40.00 -0.246000E+03 0.992000E+01 1621 | 4297.90 39935.00 40.00 -0.187810E+03 0.875620E+01 1622 | 4399.00 39950.00 40.00 0.220000E+02 0.544000E+01 1623 | 4399.00 39850.00 40.00 -0.290120E+03 0.108024E+02 1624 | 4399.00 39750.00 40.00 -0.203500E+03 0.907000E+01 1625 | 4399.80 39650.00 40.00 0.449500E+03 0.139900E+02 1626 | 4402.10 39550.00 40.00 0.137721E+04 0.325442E+02 1627 | 4405.20 39450.00 40.00 0.378650E+03 0.125730E+02 1628 | 4410.50 39350.00 40.00 -0.253200E+02 0.550640E+01 1629 | 4418.10 39250.00 40.00 0.260000E+01 0.505200E+01 1630 | 4422.80 39150.00 40.00 0.503100E+02 0.600620E+01 1631 | 4417.70 39051.00 40.00 0.736100E+02 0.647220E+01 1632 | 4491.00 39033.00 40.00 0.105800E+03 0.711600E+01 1633 | 4494.90 39133.00 40.00 0.794800E+02 0.658960E+01 1634 | 4498.90 39233.00 40.00 0.536800E+02 0.607360E+01 1635 | 4501.00 39333.00 40.00 0.129600E+02 0.525920E+01 1636 | 4500.00 39433.00 40.00 -0.465900E+02 0.593180E+01 1637 | 4497.10 39533.00 40.00 0.317880E+03 0.113576E+02 1638 | 4493.90 39632.00 40.00 0.528090E+03 0.155618E+02 1639 | 4491.10 39732.00 40.00 -0.109410E+03 0.718820E+01 1640 | 4488.10 39832.00 40.00 -0.222190E+03 0.944380E+01 1641 | 4484.80 39932.00 40.00 -0.233890E+03 0.967780E+01 1642 | -------------------------------------------------------------------------------- /environment.yml: -------------------------------------------------------------------------------- 1 | name: t21-tue-simpeg 2 | 3 | channels: 4 | - conda-forge 5 | - defaults 6 | 7 | dependencies: 8 | - python>=3.6 9 | - numpy 10 | - scipy 11 | - matplotlib 12 | - jupyter 13 | - ipywidgets>=0.6.0 14 | - mkl 15 | - simpeg 16 | - pyvista 17 | - dask 18 | - pip: 19 | - git+https://github.com/geoscixyz/geosci-labs.git 20 | -------------------------------------------------------------------------------- /environment_for_test.yml: -------------------------------------------------------------------------------- 1 | name: t21-tue-simpeg 2 | 3 | channels: 4 | - conda-forge 5 | - defaults 6 | 7 | dependencies: 8 | - python>=3.6 9 | - numpy 10 | - scipy 11 | - matplotlib 12 | - jupyter 13 | - ipywidgets>=0.6.0 14 | - mkl 15 | - simpeg 16 | - pyvista 17 | - dask 18 | -------------------------------------------------------------------------------- /figures/mag3disco.gif: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/simpeg/transform-2021-simpeg/b16c10094ee3b656bba5a631738bba28bc7d5a03/figures/mag3disco.gif -------------------------------------------------------------------------------- /figures/magobs.gif: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/simpeg/transform-2021-simpeg/b16c10094ee3b656bba5a631738bba28bc7d5a03/figures/magobs.gif -------------------------------------------------------------------------------- /figures/map.gif: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/simpeg/transform-2021-simpeg/b16c10094ee3b656bba5a631738bba28bc7d5a03/figures/map.gif -------------------------------------------------------------------------------- /figures/osmand-section.jpg: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/simpeg/transform-2021-simpeg/b16c10094ee3b656bba5a631738bba28bc7d5a03/figures/osmand-section.jpg -------------------------------------------------------------------------------- /pdfs/inversion-for-geologists-part1.pdf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/simpeg/transform-2021-simpeg/b16c10094ee3b656bba5a631738bba28bc7d5a03/pdfs/inversion-for-geologists-part1.pdf -------------------------------------------------------------------------------- /pdfs/inversion-for-geologists-part2.pdf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/simpeg/transform-2021-simpeg/b16c10094ee3b656bba5a631738bba28bc7d5a03/pdfs/inversion-for-geologists-part2.pdf -------------------------------------------------------------------------------- /requirements_dev.txt: -------------------------------------------------------------------------------- 1 | testipynb 2 | pytest -------------------------------------------------------------------------------- /tests/test_notebooks.py: -------------------------------------------------------------------------------- 1 | import os 2 | import numpy as np 3 | import testipynb 4 | import unittest 5 | 6 | NBDIR = os.path.sep.join( 7 | os.path.abspath(__file__).split(os.path.sep)[:-2] 8 | ) 9 | 10 | IGNORE = [ 11 | "1-magnetic-inversion-raglan-reproduce", 12 | "3-magnetic-inversion-raglan-lp", 13 | "2-linear-inversion-app" 14 | ] 15 | 16 | Test = testipynb.TestNotebooks(directory=NBDIR, timeout=2800) 17 | Test.ignore = IGNORE 18 | TestNotebooks = Test.get_tests() 19 | 20 | if __name__ == "__main__": 21 | unittest.main() 22 | --------------------------------------------------------------------------------