diff --git a/docs/source/pythonapi/examples/mgxs-part-iv.ipynb b/docs/source/pythonapi/examples/mgxs-part-iv.ipynb index bc85af5c2..823d67ae1 100644 --- a/docs/source/pythonapi/examples/mgxs-part-iv.ipynb +++ b/docs/source/pythonapi/examples/mgxs-part-iv.ipynb @@ -6,10 +6,10 @@ "source": [ "This Notebook illustrates the use of the openmc.mgxs.Library class specifically for application in OpenMC's multi-group mode. This example notebook follows the same process as was done in MGXS Part III, but instead uses OpenMC as the multi-group solver. This Notebook illustrates the following features:\n", "\n", - " Calculation of multi-group cross sections for a fuel assembly\n", - " Automated creation, manipulation and storage of MGXS with openmc.mgxs.Library\n", - " Validation of multi-group cross sections with OpenMC\n", - " Steady-state pin-by-pin fission rates comparison between Continuous-Energy mode and Multi-Group OpenMC.\n", + " - Calculation of multi-group cross sections for a fuel assembly\n", + " - Automated creation, manipulation and storage of MGXS with openmc.mgxs.Library\n", + " - Validation of multi-group cross sections with OpenMC\n", + " - Steady-state pin-by-pin fission rates comparison between Continuous-Energy mode and Multi-Group OpenMC.\n", "\n", "Note: This Notebook illustrates the use of Pandas DataFrames to containerize multi-group cross section data. We recommend using Pandas >v0.15.0 or later since OpenMC's Python API leverages the multi-indexing feature included in the most recent releases of Pandas.\n" ] @@ -170,22 +170,22 @@ "outputs": [], "source": [ "# Create a Universe to encapsulate a fuel pin\n", - "fuel_pin_universe = openmc.Universe(name='1.6% Fuel Pin')\n", + "fuel_pin_universe = openmc.Universe(name='1.6% Fuel Pin', universe_id=10)\n", "\n", "# Create fuel Cell\n", - "fuel_cell = openmc.Cell(name='1.6% Fuel')\n", + "fuel_cell = openmc.Cell(name='1.6% Fuel', cell_id=1)\n", "fuel_cell.fill = fuel\n", "fuel_cell.region = -fuel_outer_radius\n", "fuel_pin_universe.add_cell(fuel_cell)\n", "\n", "# Create a clad Cell\n", - "clad_cell = openmc.Cell(name='1.6% Clad')\n", + "clad_cell = openmc.Cell(name='1.6% Clad', cell_id=2)\n", "clad_cell.fill = zircaloy\n", "clad_cell.region = +fuel_outer_radius & -clad_outer_radius\n", "fuel_pin_universe.add_cell(clad_cell)\n", "\n", "# Create a moderator Cell\n", - "moderator_cell = openmc.Cell(name='1.6% Moderator')\n", + "moderator_cell = openmc.Cell(name='1.6% Moderator', cell_id=3)\n", "moderator_cell.fill = water\n", "moderator_cell.region = +clad_outer_radius\n", "fuel_pin_universe.add_cell(moderator_cell)" @@ -207,22 +207,22 @@ "outputs": [], "source": [ "# Create a Universe to encapsulate a control rod guide tube\n", - "guide_tube_universe = openmc.Universe(name='Guide Tube')\n", + "guide_tube_universe = openmc.Universe(name='Guide Tube', universe_id=20)\n", "\n", "# Create guide tube Cell\n", - "guide_tube_cell = openmc.Cell(name='Guide Tube Water')\n", + "guide_tube_cell = openmc.Cell(name='Guide Tube Water', cell_id=4)\n", "guide_tube_cell.fill = water\n", "guide_tube_cell.region = -fuel_outer_radius\n", "guide_tube_universe.add_cell(guide_tube_cell)\n", "\n", "# Create a clad Cell\n", - "clad_cell = openmc.Cell(name='Guide Clad')\n", + "clad_cell = openmc.Cell(name='Guide Clad', cell_id=5)\n", "clad_cell.fill = zircaloy\n", "clad_cell.region = +fuel_outer_radius & -clad_outer_radius\n", "guide_tube_universe.add_cell(clad_cell)\n", "\n", "# Create a moderator Cell\n", - "moderator_cell = openmc.Cell(name='Guide Tube Moderator')\n", + "moderator_cell = openmc.Cell(name='Guide Tube Moderator', cell_id=6)\n", "moderator_cell.fill = water\n", "moderator_cell.region = +clad_outer_radius\n", "guide_tube_universe.add_cell(moderator_cell)" @@ -239,12 +239,12 @@ "cell_type": "code", "execution_count": 8, "metadata": { - "collapsed": true + "collapsed": false }, "outputs": [], "source": [ "# Create fuel assembly Lattice\n", - "assembly = openmc.RectLattice(name='1.6% Fuel Assembly')\n", + "assembly = openmc.RectLattice(name='1.6% Fuel Assembly', lattice_id=100)\n", "assembly.dimension = (17, 17)\n", "assembly.pitch = (1.26, 1.26)\n", "assembly.lower_left = [-1.26 * 17. / 2.0] * 2" @@ -298,7 +298,7 @@ "outputs": [], "source": [ "# Create root Cell\n", - "root_cell = openmc.Cell(name='root cell')\n", + "root_cell = openmc.Cell(name='root cell', cell_id=0)\n", "root_cell.fill = assembly\n", "\n", "# Add boundary planes\n", @@ -347,7 +347,7 @@ "outputs": [], "source": [ "# OpenMC simulation parameters\n", - "batches = 200\n", + "batches = 500\n", "inactive = 10\n", "particles = 5000\n", "\n", @@ -434,7 +434,7 @@ "outputs": [ { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAPoAAAD6AgMAAAD1grKuAAAABGdBTUEAALGPC/xhBQAAACBjSFJN\nAAB6JgAAgIQAAPoAAACA6AAAdTAAAOpgAAA6mAAAF3CculE8AAAADFBMVEX///9yEhLpgJFNv8Tq\nQYT7AAAAAWJLR0QAiAUdSAAAAAd0SU1FB+AFBw4WAwCoz4wAAAWFSURBVGje7Zs7cttADIZ9CSvX\ncrP0iCxUqbBc8Ag6xR6BhV2EvYvwFD4CCx1ABT1jMdgndpegRQnOrCbjpPlGESISC4A/gd27e8H5\n83CX3b4+iKJrRHkS4vkghMPBonRYWGwtfgD2YN+dRDUOoh6lACw0Noi9w2fESuEoAR/uVuMolX03\n9oXGT7F3eFL2iEfhUX1f4cPdL/ishs+68ai+udE4xPhexbjX2FfjGNoPj/DPNX4Tsd+EODr8FvsV\ndf1Hd9P2VvCi4+s/aXvrf+upAD+1/9GV1mkOH5X9vV6THtfvACslcaUCbESL61drBPtdI8SrFMWr\nELsXCkuFDYW75gbiP7d9Cf7bAYI/aCwUShrBvh30+lWQkzVgZ/HD4OixNCgcQpJ3BxU/Ln91elKo\nM5VEE38QtJ+Yv6cQ9xjKNYayyl8TypP8DfJnQ2H/b/N3ye9P83cT33SQv/sQh9gV7zZ/0dNj5HQa\nC5vVzv9+/WFN2w8KVaZ2BwL1+pv4g0x1QRfjq0dB4Q3kT277oP6VNL6gKxNU9a8zK+WLbi/Wwpdi\nhbboKqyxFOulHMj6v4W/AXbmUeAxrv9J/CqEBXaRKsXaodD4nsYvkT/G6H1D4SR/iPy1Roj9JsQ5\ne18/7EUHv1+Fvx/Xj5V9Ugb5K8TW4TZEEdcvoz/up0VTe9qsVIppKVX6a7D6y9ZvwEKjrtQxPtv6\nfXII9vCxKOGaIeAIfEF8IvAG8ie3vRK9rRQl+PPpSctbhfpTUCpviH+kxsZgpT91+snoX1l49KK3\niUQvICRy5aUw6l8leoVwoo3Uv1rKreF/UFLY6d9QP4L9Wf2r7EP9GOSfcsjZ56f60kz+XmVPXv+R\nuP49ff0T/53Rv6n/7m2lvXT9Wqd/VUz8hvh5M/ED6ILmt4mfHYZSaePnTWpsf/SvqV9O6dLYYClL\nEetnoH/LBLFoBvrX189uTv8++kot5vTvQD4/9jP690g9P/4z/bvo/XVG/xYoZZx+8fr3MxAtsf7t\nUOkG2JqsTtCIpgCt/qX1226KqZS7gfzJbe+c9jLrtIZ8lXD+s4umlW6AKIVrlML2/cXjgPFjlJqI\nRC+Fj0bVJe+vSh56pSdR6YkQ1ygF10Wqf0FeLta/iKn9Mv1L24ti2e+7W4n1b3T/W+L+t9H9T/Sv\nVboUmqJJon1/hZq8LnzRDlDrX1u0xRT1+6vEpomMmyYkqi95vIH8yW1PN+122KkLcNLKi/WTF01z\n/cNASrWE/l3ev6T17zX909z9X27/euK/Rf3zWP+Waf9eEv37KkWJ+rfDl6ZglNDa+cEBhwYDvkoN\nP/rX69814NaI3imq0l7OYDy/qSdDGwr7r+Y3VbzoKZr6XX2lfxfOb87qXzr+b1j/Xlp/nP6dn98M\ncdH7cn7zjPObKsYWS3Eb9w8n85smHtqQuPuZ30T2dlIT6F9xFl+n8xslegL9a4c2KRr9W4rp/GYq\numiM9Nec/j2v/yj9u1h//hv9e93vc++f63/u+rPjL3f+5Lbn1j9m/eXWf+7zh/v8+2b9e/Hzn6s/\nuPqHrb8g71n6L3f+5Lbnvn8w33+4718/+5d47//c/gO7/5E7/nPbc/tv3P4fs//I7X9y+6/fqH+v\n6j9z+9/c/ju3/8+eP+TOn9z23PkXc/7Gnf9x5483q38Xzn+582fu/Js9fy8kb/6fO39y23P3n3S8\n/S/c/Tfc/T83uX/pgv1XE/9duP+Lu/+Mvf8td/znti8kb/8ld/9nx9t/Sjw/Ltr/yt1/+337f6/b\nf0zoB3nJ/ucVc/81d/83e/957vzJbc89/8A8f8E9/5HE78XnT/4H/cs5f8Q9/8Q9f8U+/5U7f3Lb\nc88fdrzzjyvm+cuf/Uu887/c88fs88954/8vO4SjPC+2QRIAAAAldEVYdGRhdGU6Y3JlYXRlADIw\nMTYtMDUtMDdUMTQ6MjI6MDMtMDQ6MDCiB/xLAAAAJXRFWHRkYXRlOm1vZGlmeQAyMDE2LTA1LTA3\nVDE0OjIyOjAzLTA0OjAw01pE9wAAAABJRU5ErkJggg==\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAPoAAAD6AgMAAAD1grKuAAAABGdBTUEAALGPC/xhBQAAACBjSFJN\nAAB6JgAAgIQAAPoAAACA6AAAdTAAAOpgAAA6mAAAF3CculE8AAAADFBMVEX///9yEhLpgJFNv8Tq\nQYT7AAAAAWJLR0QAiAUdSAAAAAd0SU1FB+AFDBUPAbLUzCgAAAWFSURBVGje7Zs7cttADIZ9CSvX\ncrP0iCxUqbBc8Ag6xR6BhV2EvYvwFD4CCx1ABT1jMdgndpegRQnOrCbjpPlGESISC4A/gd27e8H5\n83CX3b4+iKJrRHkS4vkghMPBonRYWGwtfgD2YN+dRDUOoh6lACw0Noi9w2fESuEoAR/uVuMolX03\n9oXGT7F3eFL2iEfhUX1f4cPdL/ishs+68ai+udE4xPhexbjX2FfjGNoPj/DPNX4Tsd+EODr8FvsV\ndf1Hd9P2VvCi4+s/aXvrf+upAD+1/9GV1mkOH5X9vV6THtfvACslcaUCbESL61drBPtdI8SrFMWr\nELsXCkuFDYW75gbiP7d9Cf7bAYI/aCwUShrBvh30+lWQkzVgZ/HD4OixNCgcQpJ3BxU/Ln91elKo\nM5VEE38QtJ+Yv6cQ9xjKNYayyl8TypP8DfJnQ2H/b/N3ye9P83cT33SQv/sQh9gV7zZ/0dNj5HQa\nC5vVzv9+/WFN2w8KVaZ2BwL1+pv4g0x1QRfjq0dB4Q3kT277oP6VNL6gKxNU9a8zK+WLbi/Wwpdi\nhbboKqyxFOulHMj6v4W/AXbmUeAxrv9J/CqEBXaRKsXaodD4nsYvkT/G6H1D4SR/iPy1Roj9JsQ5\ne18/7EUHv1+Fvx/Xj5V9Ugb5K8TW4TZEEdcvoz/up0VTe9qsVIppKVX6a7D6y9ZvwEKjrtQxPtv6\nfXII9vCxKOGaIeAIfEF8IvAG8ie3vRK9rRQl+PPpSctbhfpTUCpviH+kxsZgpT91+snoX1l49KK3\niUQvICRy5aUw6l8leoVwoo3Uv1rKreF/UFLY6d9QP4L9Wf2r7EP9GOSfcsjZ56f60kz+XmVPXv+R\nuP49ff0T/53Rv6n/7m2lvXT9Wqd/VUz8hvh5M/ED6ILmt4mfHYZSaePnTWpsf/SvqV9O6dLYYClL\nEetnoH/LBLFoBvrX189uTv8++kot5vTvQD4/9jP690g9P/4z/bvo/XVG/xYoZZx+8fr3MxAtsf7t\nUOkG2JqsTtCIpgCt/qX1226KqZS7gfzJbe+c9jLrtIZ8lXD+s4umlW6AKIVrlML2/cXjgPFjlJqI\nRC+Fj0bVJe+vSh56pSdR6YkQ1ygF10Wqf0FeLta/iKn9Mv1L24ti2e+7W4n1b3T/W+L+t9H9T/Sv\nVboUmqJJon1/hZq8LnzRDlDrX1u0xRT1+6vEpomMmyYkqi95vIH8yW1PN+122KkLcNLKi/WTF01z\n/cNASrWE/l3ev6T17zX909z9X27/euK/Rf3zWP+Waf9eEv37KkWJ+rfDl6ZglNDa+cEBhwYDvkoN\nP/rX69814NaI3imq0l7OYDy/qSdDGwr7r+Y3VbzoKZr6XX2lfxfOb87qXzr+b1j/Xlp/nP6dn98M\ncdH7cn7zjPObKsYWS3Eb9w8n85smHtqQuPuZ30T2dlIT6F9xFl+n8xslegL9a4c2KRr9W4rp/GYq\numiM9Nec/j2v/yj9u1h//hv9e93vc++f63/u+rPjL3f+5Lbn1j9m/eXWf+7zh/v8+2b9e/Hzn6s/\nuPqHrb8g71n6L3f+5Lbnvn8w33+4718/+5d47//c/gO7/5E7/nPbc/tv3P4fs//I7X9y+6/fqH+v\n6j9z+9/c/ju3/8+eP+TOn9z23PkXc/7Gnf9x5483q38Xzn+582fu/Js9fy8kb/6fO39y23P3n3S8\n/S/c/Tfc/T83uX/pgv1XE/9duP+Lu/+Mvf8td/znti8kb/8ld/9nx9t/Sjw/Ltr/yt1/+337f6/b\nf0zoB3nJ/ucVc/81d/83e/957vzJbc89/8A8f8E9/5HE78XnT/4H/cs5f8Q9/8Q9f8U+/5U7f3Lb\nc88fdrzzjyvm+cuf/Uu887/c88fs88954/8vO4SjPC+2QRIAAAAldEVYdGRhdGU6Y3JlYXRlADIw\nMTYtMDUtMTJUMjE6MTU6MDEtMDQ6MDDSFxCHAAAAJXRFWHRkYXRlOm1vZGlmeQAyMDE2LTA1LTEy\nVDIxOjE1OjAxLTA0OjAwo0qoOwAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] @@ -554,7 +554,7 @@ "cell_type": "code", "execution_count": 20, "metadata": { - "collapsed": true + "collapsed": false }, "outputs": [], "source": [ @@ -566,12 +566,31 @@ "cell_type": "markdown", "metadata": {}, "source": [ - "Lastly, we use the `Library` to construct the tallies needed to compute all of the requested multi-group cross sections in each domain and nuclide." + "Now that the `Library` has been setup, lets make sure it contains the types of cross sections which meet the needs of OpenMC's multi-group solver. Note that this step is done automatically when writing the Multi-Group Library file later in the process (as part of the `mgxs_lib.write_mg_library()`), but it is a good practice to also run this before spending all the time running OpenMC to generate the cross sections." ] }, { "cell_type": "code", "execution_count": 21, + "metadata": { + "collapsed": false + }, + "outputs": [], + "source": [ + "# Check the library - if no errors are raised, then the library is satisfactory.\n", + "mgxs_lib.check_library_for_openmc_mgxs()" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Lastly, we use the `Library` to construct the tallies needed to compute all of the requested multi-group cross sections in each domain and nuclide." + ] + }, + { + "cell_type": "code", + "execution_count": 22, "metadata": { "collapsed": true }, @@ -592,7 +611,7 @@ }, { "cell_type": "code", - "execution_count": 22, + "execution_count": 23, "metadata": { "collapsed": true }, @@ -612,7 +631,7 @@ }, { "cell_type": "code", - "execution_count": 23, + "execution_count": 24, "metadata": { "collapsed": true }, @@ -640,7 +659,7 @@ }, { "cell_type": "code", - "execution_count": 24, + "execution_count": 25, "metadata": { "collapsed": true }, @@ -652,7 +671,7 @@ }, { "cell_type": "code", - "execution_count": 25, + "execution_count": 26, "metadata": { "collapsed": false }, @@ -677,8 +696,8 @@ " Copyright: 2011-2016 Massachusetts Institute of Technology\n", " License: http://openmc.readthedocs.org/en/latest/license.html\n", " Version: 0.7.1\n", - " Git SHA1: 179e9ab147e505563d118ed58096b3d225160ffa\n", - " Date/Time: 2016-05-07 14:22:04\n", + " Git SHA1: c779ca42c41a062a6a813e03f2add2d182ca9190\n", + " Date/Time: 2016-05-12 21:15:02\n", " OpenMP Threads: 4\n", "\n", " ===========================================================================\n", @@ -906,7 +925,307 @@ " 198/1 1.01366 1.02596 +/- 0.00130\n", " 199/1 1.04471 1.02605 +/- 0.00130\n", " 200/1 1.02416 1.02604 +/- 0.00129\n", - " Creating state point statepoint.200.h5...\n", + " 201/1 1.01172 1.02597 +/- 0.00129\n", + " 202/1 1.01683 1.02592 +/- 0.00128\n", + " 203/1 1.01341 1.02586 +/- 0.00128\n", + " 204/1 1.01507 1.02580 +/- 0.00127\n", + " 205/1 1.02540 1.02580 +/- 0.00127\n", + " 206/1 1.00310 1.02568 +/- 0.00127\n", + " 207/1 1.02822 1.02570 +/- 0.00126\n", + " 208/1 1.01023 1.02562 +/- 0.00126\n", + " 209/1 1.04603 1.02572 +/- 0.00125\n", + " 210/1 1.00775 1.02563 +/- 0.00125\n", + " 211/1 1.01706 1.02559 +/- 0.00125\n", + " 212/1 0.99434 1.02543 +/- 0.00125\n", + " 213/1 1.03346 1.02547 +/- 0.00124\n", + " 214/1 1.05322 1.02561 +/- 0.00124\n", + " 215/1 1.03057 1.02563 +/- 0.00124\n", + " 216/1 1.00976 1.02556 +/- 0.00123\n", + " 217/1 1.02760 1.02557 +/- 0.00123\n", + " 218/1 1.01259 1.02550 +/- 0.00122\n", + " 219/1 1.02829 1.02552 +/- 0.00122\n", + " 220/1 1.02228 1.02550 +/- 0.00121\n", + " 221/1 1.06679 1.02570 +/- 0.00122\n", + " 222/1 1.03417 1.02574 +/- 0.00122\n", + " 223/1 1.04239 1.02582 +/- 0.00121\n", + " 224/1 1.02062 1.02579 +/- 0.00121\n", + " 225/1 1.00331 1.02569 +/- 0.00121\n", + " 226/1 1.00131 1.02557 +/- 0.00121\n", + " 227/1 1.01768 1.02554 +/- 0.00120\n", + " 228/1 1.00813 1.02546 +/- 0.00120\n", + " 229/1 1.05320 1.02558 +/- 0.00120\n", + " 230/1 1.03472 1.02563 +/- 0.00120\n", + " 231/1 1.01426 1.02557 +/- 0.00119\n", + " 232/1 1.00782 1.02549 +/- 0.00119\n", + " 233/1 1.02813 1.02551 +/- 0.00118\n", + " 234/1 1.01184 1.02545 +/- 0.00118\n", + " 235/1 1.02156 1.02543 +/- 0.00118\n", + " 236/1 0.99029 1.02527 +/- 0.00118\n", + " 237/1 1.04196 1.02535 +/- 0.00118\n", + " 238/1 1.01594 1.02531 +/- 0.00117\n", + " 239/1 1.02732 1.02531 +/- 0.00117\n", + " 240/1 0.98987 1.02516 +/- 0.00117\n", + " 241/1 1.03388 1.02520 +/- 0.00117\n", + " 242/1 1.01319 1.02515 +/- 0.00116\n", + " 243/1 1.02870 1.02516 +/- 0.00116\n", + " 244/1 1.01943 1.02514 +/- 0.00115\n", + " 245/1 1.04463 1.02522 +/- 0.00115\n", + " 246/1 1.03551 1.02526 +/- 0.00115\n", + " 247/1 1.00436 1.02517 +/- 0.00115\n", + " 248/1 1.03326 1.02521 +/- 0.00114\n", + " 249/1 1.05769 1.02534 +/- 0.00115\n", + " 250/1 1.01372 1.02530 +/- 0.00114\n", + " 251/1 1.02971 1.02531 +/- 0.00114\n", + " 252/1 1.01166 1.02526 +/- 0.00113\n", + " 253/1 1.03992 1.02532 +/- 0.00113\n", + " 254/1 1.01507 1.02528 +/- 0.00113\n", + " 255/1 1.03222 1.02530 +/- 0.00112\n", + " 256/1 1.03096 1.02533 +/- 0.00112\n", + " 257/1 1.01153 1.02527 +/- 0.00112\n", + " 258/1 1.03668 1.02532 +/- 0.00111\n", + " 259/1 1.03070 1.02534 +/- 0.00111\n", + " 260/1 1.01189 1.02529 +/- 0.00111\n", + " 261/1 1.00082 1.02519 +/- 0.00111\n", + " 262/1 1.03653 1.02523 +/- 0.00110\n", + " 263/1 1.02908 1.02525 +/- 0.00110\n", + " 264/1 1.00072 1.02515 +/- 0.00110\n", + " 265/1 1.00832 1.02509 +/- 0.00109\n", + " 266/1 1.04385 1.02516 +/- 0.00109\n", + " 267/1 1.00117 1.02507 +/- 0.00109\n", + " 268/1 1.02682 1.02507 +/- 0.00109\n", + " 269/1 1.03202 1.02510 +/- 0.00108\n", + " 270/1 1.01275 1.02505 +/- 0.00108\n", + " 271/1 1.02633 1.02506 +/- 0.00108\n", + " 272/1 1.04811 1.02514 +/- 0.00108\n", + " 273/1 1.02851 1.02516 +/- 0.00107\n", + " 274/1 1.01270 1.02511 +/- 0.00107\n", + " 275/1 1.06222 1.02525 +/- 0.00107\n", + " 276/1 1.02778 1.02526 +/- 0.00107\n", + " 277/1 1.02601 1.02526 +/- 0.00107\n", + " 278/1 1.02356 1.02526 +/- 0.00106\n", + " 279/1 1.00792 1.02519 +/- 0.00106\n", + " 280/1 1.02331 1.02518 +/- 0.00106\n", + " 281/1 1.00985 1.02513 +/- 0.00105\n", + " 282/1 1.02035 1.02511 +/- 0.00105\n", + " 283/1 0.98181 1.02495 +/- 0.00106\n", + " 284/1 1.01829 1.02493 +/- 0.00106\n", + " 285/1 1.02929 1.02494 +/- 0.00105\n", + " 286/1 1.03524 1.02498 +/- 0.00105\n", + " 287/1 1.01212 1.02493 +/- 0.00105\n", + " 288/1 1.03584 1.02497 +/- 0.00104\n", + " 289/1 1.02961 1.02499 +/- 0.00104\n", + " 290/1 0.99692 1.02489 +/- 0.00104\n", + " 291/1 1.03966 1.02494 +/- 0.00104\n", + " 292/1 1.00965 1.02489 +/- 0.00104\n", + " 293/1 1.02601 1.02489 +/- 0.00103\n", + " 294/1 1.03224 1.02492 +/- 0.00103\n", + " 295/1 1.01596 1.02489 +/- 0.00103\n", + " 296/1 1.06964 1.02504 +/- 0.00103\n", + " 297/1 1.03982 1.02509 +/- 0.00103\n", + " 298/1 0.99758 1.02500 +/- 0.00103\n", + " 299/1 1.01479 1.02496 +/- 0.00103\n", + " 300/1 1.04517 1.02503 +/- 0.00103\n", + " 301/1 0.99128 1.02492 +/- 0.00103\n", + " 302/1 1.01493 1.02488 +/- 0.00103\n", + " 303/1 1.00623 1.02482 +/- 0.00103\n", + " 304/1 1.02560 1.02482 +/- 0.00102\n", + " 305/1 1.00806 1.02477 +/- 0.00102\n", + " 306/1 1.03524 1.02480 +/- 0.00102\n", + " 307/1 0.99244 1.02469 +/- 0.00102\n", + " 308/1 0.98013 1.02454 +/- 0.00103\n", + " 309/1 1.00853 1.02449 +/- 0.00103\n", + " 310/1 1.00116 1.02441 +/- 0.00103\n", + " 311/1 1.01730 1.02439 +/- 0.00102\n", + " 312/1 1.01198 1.02435 +/- 0.00102\n", + " 313/1 1.02405 1.02435 +/- 0.00102\n", + " 314/1 1.01734 1.02432 +/- 0.00101\n", + " 315/1 1.02320 1.02432 +/- 0.00101\n", + " 316/1 1.03438 1.02435 +/- 0.00101\n", + " 317/1 1.00106 1.02428 +/- 0.00101\n", + " 318/1 1.03114 1.02430 +/- 0.00100\n", + " 319/1 1.04955 1.02438 +/- 0.00100\n", + " 320/1 1.03259 1.02441 +/- 0.00100\n", + " 321/1 1.00687 1.02435 +/- 0.00100\n", + " 322/1 1.05753 1.02446 +/- 0.00100\n", + " 323/1 1.03676 1.02450 +/- 0.00100\n", + " 324/1 0.99796 1.02441 +/- 0.00100\n", + " 325/1 1.03783 1.02445 +/- 0.00100\n", + " 326/1 1.02315 1.02445 +/- 0.00099\n", + " 327/1 1.04205 1.02451 +/- 0.00099\n", + " 328/1 1.01971 1.02449 +/- 0.00099\n", + " 329/1 1.02394 1.02449 +/- 0.00099\n", + " 330/1 1.03318 1.02452 +/- 0.00098\n", + " 331/1 1.01503 1.02449 +/- 0.00098\n", + " 332/1 1.07143 1.02463 +/- 0.00099\n", + " 333/1 1.00991 1.02459 +/- 0.00099\n", + " 334/1 1.03115 1.02461 +/- 0.00098\n", + " 335/1 1.04400 1.02467 +/- 0.00098\n", + " 336/1 1.03516 1.02470 +/- 0.00098\n", + " 337/1 1.02025 1.02468 +/- 0.00098\n", + " 338/1 1.03269 1.02471 +/- 0.00098\n", + " 339/1 1.03745 1.02475 +/- 0.00097\n", + " 340/1 1.03685 1.02478 +/- 0.00097\n", + " 341/1 1.01831 1.02476 +/- 0.00097\n", + " 342/1 1.01425 1.02473 +/- 0.00097\n", + " 343/1 1.02990 1.02475 +/- 0.00096\n", + " 344/1 1.02958 1.02476 +/- 0.00096\n", + " 345/1 1.03133 1.02478 +/- 0.00096\n", + " 346/1 1.02441 1.02478 +/- 0.00095\n", + " 347/1 1.07010 1.02492 +/- 0.00096\n", + " 348/1 1.02327 1.02491 +/- 0.00096\n", + " 349/1 1.03123 1.02493 +/- 0.00096\n", + " 350/1 1.03158 1.02495 +/- 0.00095\n", + " 351/1 1.03473 1.02498 +/- 0.00095\n", + " 352/1 1.04000 1.02502 +/- 0.00095\n", + " 353/1 1.01651 1.02500 +/- 0.00095\n", + " 354/1 1.03647 1.02503 +/- 0.00094\n", + " 355/1 1.04650 1.02509 +/- 0.00094\n", + " 356/1 1.04703 1.02516 +/- 0.00094\n", + " 357/1 1.00260 1.02509 +/- 0.00094\n", + " 358/1 1.00075 1.02502 +/- 0.00094\n", + " 359/1 1.04874 1.02509 +/- 0.00094\n", + " 360/1 1.03211 1.02511 +/- 0.00094\n", + " 361/1 1.02136 1.02510 +/- 0.00094\n", + " 362/1 1.00803 1.02505 +/- 0.00094\n", + " 363/1 1.00319 1.02499 +/- 0.00094\n", + " 364/1 1.01443 1.02496 +/- 0.00093\n", + " 365/1 1.02685 1.02496 +/- 0.00093\n", + " 366/1 1.02373 1.02496 +/- 0.00093\n", + " 367/1 1.02026 1.02495 +/- 0.00093\n", + " 368/1 1.01579 1.02492 +/- 0.00092\n", + " 369/1 1.08004 1.02508 +/- 0.00093\n", + " 370/1 1.01715 1.02505 +/- 0.00093\n", + " 371/1 0.98578 1.02494 +/- 0.00093\n", + " 372/1 1.03033 1.02496 +/- 0.00093\n", + " 373/1 1.03269 1.02498 +/- 0.00093\n", + " 374/1 1.04050 1.02502 +/- 0.00093\n", + " 375/1 1.00760 1.02498 +/- 0.00093\n", + " 376/1 1.04492 1.02503 +/- 0.00093\n", + " 377/1 1.04983 1.02510 +/- 0.00093\n", + " 378/1 1.06022 1.02519 +/- 0.00093\n", + " 379/1 1.02516 1.02519 +/- 0.00093\n", + " 380/1 1.01740 1.02517 +/- 0.00092\n", + " 381/1 1.02520 1.02517 +/- 0.00092\n", + " 382/1 1.02820 1.02518 +/- 0.00092\n", + " 383/1 1.00697 1.02513 +/- 0.00092\n", + " 384/1 1.03497 1.02516 +/- 0.00092\n", + " 385/1 0.98404 1.02505 +/- 0.00092\n", + " 386/1 1.05206 1.02512 +/- 0.00092\n", + " 387/1 1.01502 1.02509 +/- 0.00092\n", + " 388/1 1.02196 1.02508 +/- 0.00092\n", + " 389/1 1.02856 1.02509 +/- 0.00091\n", + " 390/1 1.01376 1.02506 +/- 0.00091\n", + " 391/1 1.01696 1.02504 +/- 0.00091\n", + " 392/1 1.03283 1.02506 +/- 0.00091\n", + " 393/1 1.00787 1.02502 +/- 0.00091\n", + " 394/1 1.02184 1.02501 +/- 0.00090\n", + " 395/1 1.03170 1.02503 +/- 0.00090\n", + " 396/1 1.04406 1.02508 +/- 0.00090\n", + " 397/1 1.03939 1.02511 +/- 0.00090\n", + " 398/1 1.00329 1.02506 +/- 0.00090\n", + " 399/1 1.04518 1.02511 +/- 0.00090\n", + " 400/1 1.03435 1.02513 +/- 0.00090\n", + " 401/1 1.00525 1.02508 +/- 0.00089\n", + " 402/1 1.03112 1.02510 +/- 0.00089\n", + " 403/1 1.00188 1.02504 +/- 0.00089\n", + " 404/1 1.01241 1.02501 +/- 0.00089\n", + " 405/1 1.01796 1.02499 +/- 0.00089\n", + " 406/1 1.02686 1.02499 +/- 0.00089\n", + " 407/1 1.01003 1.02496 +/- 0.00088\n", + " 408/1 1.02359 1.02495 +/- 0.00088\n", + " 409/1 1.01258 1.02492 +/- 0.00088\n", + " 410/1 1.04361 1.02497 +/- 0.00088\n", + " 411/1 1.00885 1.02493 +/- 0.00088\n", + " 412/1 1.00999 1.02489 +/- 0.00088\n", + " 413/1 0.97832 1.02477 +/- 0.00088\n", + " 414/1 1.04183 1.02482 +/- 0.00088\n", + " 415/1 1.02279 1.02481 +/- 0.00088\n", + " 416/1 1.04197 1.02485 +/- 0.00088\n", + " 417/1 1.04617 1.02491 +/- 0.00088\n", + " 418/1 1.01311 1.02488 +/- 0.00088\n", + " 419/1 1.03904 1.02491 +/- 0.00087\n", + " 420/1 1.00458 1.02486 +/- 0.00087\n", + " 421/1 0.98580 1.02477 +/- 0.00088\n", + " 422/1 1.01850 1.02475 +/- 0.00087\n", + " 423/1 1.03739 1.02478 +/- 0.00087\n", + " 424/1 1.02716 1.02479 +/- 0.00087\n", + " 425/1 1.00711 1.02475 +/- 0.00087\n", + " 426/1 1.01008 1.02471 +/- 0.00087\n", + " 427/1 1.03332 1.02473 +/- 0.00087\n", + " 428/1 1.00501 1.02468 +/- 0.00087\n", + " 429/1 1.04549 1.02473 +/- 0.00086\n", + " 430/1 1.00582 1.02469 +/- 0.00086\n", + " 431/1 1.00586 1.02464 +/- 0.00086\n", + " 432/1 1.00082 1.02459 +/- 0.00086\n", + " 433/1 1.00835 1.02455 +/- 0.00086\n", + " 434/1 1.03965 1.02458 +/- 0.00086\n", + " 435/1 1.02385 1.02458 +/- 0.00086\n", + " 436/1 1.01440 1.02456 +/- 0.00086\n", + " 437/1 1.03127 1.02458 +/- 0.00085\n", + " 438/1 1.02961 1.02459 +/- 0.00085\n", + " 439/1 0.99584 1.02452 +/- 0.00085\n", + " 440/1 1.04964 1.02458 +/- 0.00085\n", + " 441/1 0.99792 1.02452 +/- 0.00085\n", + " 442/1 1.04971 1.02457 +/- 0.00085\n", + " 443/1 1.01504 1.02455 +/- 0.00085\n", + " 444/1 1.04359 1.02460 +/- 0.00085\n", + " 445/1 1.01148 1.02457 +/- 0.00085\n", + " 446/1 1.01203 1.02454 +/- 0.00085\n", + " 447/1 1.02353 1.02454 +/- 0.00085\n", + " 448/1 1.06299 1.02462 +/- 0.00085\n", + " 449/1 1.00017 1.02457 +/- 0.00085\n", + " 450/1 1.01193 1.02454 +/- 0.00085\n", + " 451/1 1.00179 1.02449 +/- 0.00085\n", + " 452/1 1.02425 1.02449 +/- 0.00085\n", + " 453/1 1.03629 1.02451 +/- 0.00084\n", + " 454/1 1.01955 1.02450 +/- 0.00084\n", + " 455/1 1.00870 1.02447 +/- 0.00084\n", + " 456/1 1.04230 1.02451 +/- 0.00084\n", + " 457/1 1.05081 1.02457 +/- 0.00084\n", + " 458/1 1.00271 1.02452 +/- 0.00084\n", + " 459/1 1.01010 1.02448 +/- 0.00084\n", + " 460/1 1.04656 1.02453 +/- 0.00084\n", + " 461/1 1.00790 1.02450 +/- 0.00084\n", + " 462/1 1.02214 1.02449 +/- 0.00084\n", + " 463/1 1.04401 1.02453 +/- 0.00083\n", + " 464/1 1.02863 1.02454 +/- 0.00083\n", + " 465/1 0.99971 1.02449 +/- 0.00083\n", + " 466/1 1.00344 1.02444 +/- 0.00083\n", + " 467/1 1.02810 1.02445 +/- 0.00083\n", + " 468/1 1.02091 1.02444 +/- 0.00083\n", + " 469/1 1.00545 1.02440 +/- 0.00083\n", + " 470/1 1.01590 1.02438 +/- 0.00083\n", + " 471/1 1.04465 1.02443 +/- 0.00083\n", + " 472/1 1.02028 1.02442 +/- 0.00082\n", + " 473/1 1.01951 1.02441 +/- 0.00082\n", + " 474/1 1.03280 1.02443 +/- 0.00082\n", + " 475/1 1.04722 1.02447 +/- 0.00082\n", + " 476/1 1.03587 1.02450 +/- 0.00082\n", + " 477/1 1.02234 1.02449 +/- 0.00082\n", + " 478/1 1.07848 1.02461 +/- 0.00082\n", + " 479/1 1.04759 1.02466 +/- 0.00082\n", + " 480/1 1.07189 1.02476 +/- 0.00083\n", + " 481/1 1.05811 1.02483 +/- 0.00083\n", + " 482/1 1.04554 1.02487 +/- 0.00083\n", + " 483/1 1.01956 1.02486 +/- 0.00083\n", + " 484/1 1.01055 1.02483 +/- 0.00083\n", + " 485/1 1.00845 1.02480 +/- 0.00082\n", + " 486/1 1.04607 1.02484 +/- 0.00082\n", + " 487/1 1.05955 1.02492 +/- 0.00083\n", + " 488/1 1.02245 1.02491 +/- 0.00082\n", + " 489/1 0.98206 1.02482 +/- 0.00083\n", + " 490/1 1.03786 1.02485 +/- 0.00083\n", + " 491/1 1.02973 1.02486 +/- 0.00082\n", + " 492/1 1.02890 1.02487 +/- 0.00082\n", + " 493/1 1.02086 1.02486 +/- 0.00082\n", + " 494/1 1.01194 1.02483 +/- 0.00082\n", + " 495/1 1.01902 1.02482 +/- 0.00082\n", + " 496/1 1.01783 1.02481 +/- 0.00082\n", + " 497/1 1.02129 1.02480 +/- 0.00081\n", + " 498/1 1.02407 1.02480 +/- 0.00081\n", + " 499/1 1.02873 1.02480 +/- 0.00081\n", + " 500/1 1.00998 1.02477 +/- 0.00081\n", + " Creating state point statepoint.500.h5...\n", "\n", " ===========================================================================\n", " ======================> SIMULATION FINISHED <======================\n", @@ -915,27 +1234,27 @@ "\n", " =======================> TIMING STATISTICS <=======================\n", "\n", - " Total time for initialization = 1.4810E+00 seconds\n", - " Reading cross sections = 1.1600E+00 seconds\n", - " Total time in simulation = 9.8823E+01 seconds\n", - " Time in transport only = 9.8622E+01 seconds\n", - " Time in inactive batches = 2.1290E+00 seconds\n", - " Time in active batches = 9.6694E+01 seconds\n", - " Time synchronizing fission bank = 1.4000E-02 seconds\n", - " Sampling source sites = 1.1000E-02 seconds\n", - " SEND/RECV source sites = 3.0000E-03 seconds\n", + " Total time for initialization = 1.5880E+00 seconds\n", + " Reading cross sections = 1.2650E+00 seconds\n", + " Total time in simulation = 2.6051E+02 seconds\n", + " Time in transport only = 2.6013E+02 seconds\n", + " Time in inactive batches = 2.0990E+00 seconds\n", + " Time in active batches = 2.5841E+02 seconds\n", + " Time synchronizing fission bank = 6.5000E-02 seconds\n", + " Sampling source sites = 4.4000E-02 seconds\n", + " SEND/RECV source sites = 2.1000E-02 seconds\n", " Time accumulating tallies = 2.0000E-03 seconds\n", " Total time for finalization = 0.0000E+00 seconds\n", - " Total time elapsed = 1.0031E+02 seconds\n", - " Calculation Rate (inactive) = 23485.2 neutrons/second\n", - " Calculation Rate (active) = 9824.81 neutrons/second\n", + " Total time elapsed = 2.6211E+02 seconds\n", + " Calculation Rate (inactive) = 23820.9 neutrons/second\n", + " Calculation Rate (active) = 9480.98 neutrons/second\n", "\n", " ============================> RESULTS <============================\n", "\n", - " k-effective (Collision) = 1.02505 +/- 0.00122\n", - " k-effective (Track-length) = 1.02604 +/- 0.00129\n", - " k-effective (Absorption) = 1.02501 +/- 0.00111\n", - " Combined k-effective = 1.02544 +/- 0.00091\n", + " k-effective (Collision) = 1.02480 +/- 0.00073\n", + " k-effective (Track-length) = 1.02477 +/- 0.00081\n", + " k-effective (Absorption) = 1.02552 +/- 0.00068\n", + " Combined k-effective = 1.02519 +/- 0.00055\n", " Leakage Fraction = 0.00000 +/- 0.00000\n", "\n" ] @@ -946,7 +1265,7 @@ "0" ] }, - "execution_count": 25, + "execution_count": 26, "metadata": {}, "output_type": "execute_result" } @@ -965,14 +1284,14 @@ }, { "cell_type": "code", - "execution_count": 26, + "execution_count": 27, "metadata": { - "collapsed": true + "collapsed": false }, "outputs": [], "source": [ "# Move the StatePoint File\n", - "ce_spfile = './ce.h5'\n", + "ce_spfile = './ce_statepoint.h5'\n", "os.rename('statepoint.' + str(batches) + '.h5', ce_spfile)\n", "# Move the Summary file\n", "ce_sumfile = './ce_summary.h5'\n", @@ -985,44 +1304,26 @@ "source": [ "# Tally Data Processing\n", "\n", - "Our simulation ran successfully and created statepoint and summary output files. We begin our analysis by instantiating a `StatePoint` object." - ] - }, - { - "cell_type": "code", - "execution_count": 27, - "metadata": { - "collapsed": true - }, - "outputs": [], - "source": [ - "# Load the statepoint file\n", - "sp = openmc.StatePoint(ce_spfile)" - ] - }, - { - "cell_type": "markdown", - "metadata": {}, - "source": [ - "Next we will save the value of keff from the continuous-energy calculation for later comparison" + "Our simulation ran successfully and created statepoint and summary output files. Let's begin by loading the StatePoint file, but not automatically linking the summary file." ] }, { "cell_type": "code", "execution_count": 28, "metadata": { - "collapsed": true + "collapsed": false }, "outputs": [], "source": [ - "ce_keff = sp.k_combined" + "# Load the statepoint file, but not the summary file, as it is a different filename than expected.\n", + "sp = openmc.StatePoint(ce_spfile, autolink=False)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ - "In addition to the statepoint file, our simulation also created a summary file which encapsulates information about the materials and geometry. This is necessary for the `openmc.mgxs` module to properly process the tally data. We first create a `Summary` object and link it with the statepoint." + "In addition to the statepoint file, our simulation also created a summary file which encapsulates information about the materials and geometry. This is necessary for the `openmc.mgxs` module to properly process the tally data. We first create a `Summary` object and link it with the statepoint. Normally this would not need to be performed, but since we have renamed our summary file to avoid conflicts with the Multi-Group calculation's summary file, we will load this in explicitly." ] }, { @@ -1037,32 +1338,6 @@ "sp.link_with_summary(su)" ] }, - { - "cell_type": "markdown", - "metadata": {}, - "source": [ - "Next we will extract our fission distribution results from the statepoint for later comparison." - ] - }, - { - "cell_type": "code", - "execution_count": 30, - "metadata": { - "collapsed": false - }, - "outputs": [], - "source": [ - "# Get the OpenMC fission rate mesh tally data\n", - "mesh_tally = sp.get_tally(name='mesh tally')\n", - "openmc_fission_rates = mesh_tally.get_values(scores=['fission'])\n", - "\n", - "# Reshape array to 2D for plotting\n", - "openmc_fission_rates.shape = (17,17)\n", - "\n", - "# Normalize to the average pin power\n", - "openmc_fission_rates /= np.mean(openmc_fission_rates)" - ] - }, { "cell_type": "markdown", "metadata": {}, @@ -1072,7 +1347,7 @@ }, { "cell_type": "code", - "execution_count": 31, + "execution_count": 30, "metadata": { "collapsed": false }, @@ -1105,7 +1380,7 @@ }, { "cell_type": "code", - "execution_count": 32, + "execution_count": 31, "metadata": { "collapsed": false }, @@ -1114,53 +1389,38 @@ "name": "stderr", "output_type": "stream", "text": [ - "/home/nelsonag/git/openmc/openmc/tallies.py:1996: RuntimeWarning: invalid value encountered in true_divide\n", + "/home/nelsonag/git/openmc/openmc/tallies.py:1988: RuntimeWarning: invalid value encountered in true_divide\n", " self_rel_err = data['self']['std. dev.'] / data['self']['mean']\n", - "/home/nelsonag/git/openmc/openmc/tallies.py:1997: RuntimeWarning: invalid value encountered in true_divide\n", + "/home/nelsonag/git/openmc/openmc/tallies.py:1989: RuntimeWarning: invalid value encountered in true_divide\n", " other_rel_err = data['other']['std. dev.'] / data['other']['mean']\n", - "/home/nelsonag/git/openmc/openmc/tallies.py:1998: RuntimeWarning: invalid value encountered in true_divide\n", + "/home/nelsonag/git/openmc/openmc/tallies.py:1990: RuntimeWarning: invalid value encountered in true_divide\n", " new_tally._mean = data['self']['mean'] / data['other']['mean']\n" ] - }, - { - "data": { - "text/plain": [ - "{10000: 'fuel.2g',\n", - " 10001: 'fuel_clad.2g',\n", - " 10002: 'fuel_mod.2g',\n", - " 10003: 'gt_inmod.2g',\n", - " 10004: 'gt_clad.2g',\n", - " 10005: 'gt_outmod.2g'}" - ] - }, - "execution_count": 32, - "metadata": {}, - "output_type": "execute_result" } ], "source": [ "mgxs_lib.write_mg_library(filename='mgxs', xs_type='macro',\n", " domain_names=['fuel', 'fuel_clad', 'fuel_mod',\n", " 'gt_inmod', 'gt_clad', 'gt_outmod'],\n", - " xs_ids='2g')" + " xs_ids='2m')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ - "Now we will need to recreate similar xml files from above, beginning with materials.xml" + "Now we will need to recreate similar xml files from above, beginning with materials.xml. Similar to how continuous-energy cross section libraries are named, the `openmc.Macroscopic` quantities below can either have their `xs_id` included (i.e., `'.2m'`), or this can be left off but the `default_xs` parameter of the materials file be used instead to be set to the `'xs_id'` of interest (which is `'.2m'` in this case as defined in the previous cell)." ] }, { "cell_type": "code", - "execution_count": 33, + "execution_count": 32, "metadata": { "collapsed": false }, "outputs": [], "source": [ - "# Instantiate our Macroscopic Data\n", + "# Instantiate our Macroscopic Data using mat_names for the name\n", "fuel_macro = openmc.Macroscopic('fuel')\n", "fuel_clad_macro = openmc.Macroscopic('fuel_clad')\n", "fuel_mod_macro = openmc.Macroscopic('fuel_mod')\n", @@ -1171,39 +1431,39 @@ "# Now define the materials\n", "\n", "# 1.6 enriched fuel UO2\n", - "fuel = openmc.Material(name='1.6% Fuel UO2')\n", + "fuel = openmc.Material(name='1.6% Fuel UO2', material_id=1)\n", "fuel.set_density('macro', 1.0)\n", "fuel.add_macroscopic(fuel_macro)\n", "\n", "# 1.6 enriched fuel cladding\n", - "fuel_clad = openmc.Material(name='1.6% Fuel Clad')\n", + "fuel_clad = openmc.Material(name='1.6% Fuel Clad', material_id=2)\n", "fuel_clad.set_density('macro', 1.0)\n", "fuel_clad.add_macroscopic(fuel_clad_macro)\n", "\n", "# 1.6 enriched fuel moderator\n", - "fuel_mod = openmc.Material(name='1.6% Fuel Water')\n", + "fuel_mod = openmc.Material(name='1.6% Fuel Water', material_id=3)\n", "fuel_mod.set_density('macro', 1.0)\n", "fuel_mod.add_macroscopic(fuel_mod_macro)\n", "\n", "# Guide Tube Inner Moderator\n", - "gt_inmod = openmc.Material(name='GT Inner Water')\n", + "gt_inmod = openmc.Material(name='GT Inner Water', material_id=4)\n", "gt_inmod.set_density('macro', 1.0)\n", "gt_inmod.add_macroscopic(gt_inmod_macro)\n", "\n", "# Guide Tube Cladding\n", - "gt_clad = openmc.Material(name='GT Clad')\n", + "gt_clad = openmc.Material(name='GT Clad', material_id=5)\n", "gt_clad.set_density('macro', 1.0)\n", "gt_clad.add_macroscopic(gt_clad_macro)\n", "\n", "# Guide Tube Outer Moderator\n", - "gt_outmod = openmc.Material(name='GT Outer Water')\n", + "gt_outmod = openmc.Material(name='GT Outer Water', material_id=6)\n", "gt_outmod.set_density('macro', 1.0)\n", "gt_outmod.add_macroscopic(gt_outmod_macro)\n", "\n", "# Finally, instantiate our Materials object\n", "materials_file = openmc.Materials((fuel, fuel_clad, fuel_mod,\n", " gt_inmod, gt_clad, gt_outmod))\n", - "materials_file.default_xs = '2g'\n", + "materials_file.default_xs = '2m'\n", "\n", "# Export to \"materials.xml\"\n", "materials_file.export_to_xml()\n" @@ -1213,61 +1473,62 @@ "cell_type": "markdown", "metadata": {}, "source": [ - "For our geometry files we will simply repeat what as done for continuous-energy mode, except change the cell fill (i.e., the material) to use our newly defined materials." + "\n", + "For our geometry files we will do the same as before but now we will be pointing at our newly created materials instead." ] }, { "cell_type": "code", - "execution_count": 34, + "execution_count": 33, "metadata": { "collapsed": false }, "outputs": [], "source": [ "# Create a Universe to encapsulate a fuel pin\n", - "fuel_pin_universe = openmc.Universe(name='1.6% Fuel Pin')\n", + "fuel_pin_universe = openmc.Universe(name='1.6% Fuel Pin', universe_id=10)\n", "\n", "# Create fuel Cell\n", - "fuel_cell = openmc.Cell(name='1.6% Fuel')\n", + "fuel_cell = openmc.Cell(name='1.6% Fuel', cell_id=1)\n", "fuel_cell.fill = fuel\n", "fuel_cell.region = -fuel_outer_radius\n", "fuel_pin_universe.add_cell(fuel_cell)\n", "\n", "# Create a clad Cell\n", - "clad_cell = openmc.Cell(name='1.6% Clad')\n", + "clad_cell = openmc.Cell(name='1.6% Clad', cell_id=2)\n", "clad_cell.fill = fuel_clad\n", "clad_cell.region = +fuel_outer_radius & -clad_outer_radius\n", "fuel_pin_universe.add_cell(clad_cell)\n", "\n", "# Create a moderator Cell\n", - "moderator_cell = openmc.Cell(name='1.6% Moderator')\n", + "moderator_cell = openmc.Cell(name='1.6% Moderator', cell_id=3)\n", "moderator_cell.fill = fuel_mod\n", "moderator_cell.region = +clad_outer_radius\n", "fuel_pin_universe.add_cell(moderator_cell)\n", "\n", "# Create a Universe to encapsulate a control rod guide tube\n", - "guide_tube_universe = openmc.Universe(name='Guide Tube')\n", + "guide_tube_universe = openmc.Universe(name='Guide Tube', universe_id=20)\n", "\n", "# Create guide tube Cell\n", - "guide_tube_cell = openmc.Cell(name='Guide Tube Water')\n", + "guide_tube_cell = openmc.Cell(name='Guide Tube Water', cell_id=4)\n", "guide_tube_cell.fill = gt_inmod\n", "guide_tube_cell.region = -fuel_outer_radius\n", "guide_tube_universe.add_cell(guide_tube_cell)\n", "\n", "# Create a clad Cell\n", - "clad_cell = openmc.Cell(name='Guide Clad')\n", + "clad_cell = openmc.Cell(name='Guide Clad', cell_id=5)\n", "clad_cell.fill = gt_clad\n", "clad_cell.region = +fuel_outer_radius & -clad_outer_radius\n", "guide_tube_universe.add_cell(clad_cell)\n", "\n", "# Create a moderator Cell\n", - "moderator_cell = openmc.Cell(name='Guide Tube Moderator')\n", + "moderator_cell = openmc.Cell(name='Guide Tube Moderator', cell_id=6)\n", "moderator_cell.fill = gt_outmod\n", "moderator_cell.region = +clad_outer_radius\n", "guide_tube_universe.add_cell(moderator_cell)\n", "\n", "# Create fuel assembly Lattice\n", - "assembly = openmc.RectLattice(name='1.6% Fuel Assembly')\n", + "assembly = openmc.RectLattice(name='1.6% Fuel Assembly', lattice_id=100)\n", "assembly.dimension = (17, 17)\n", "assembly.pitch = (1.26, 1.26)\n", "assembly.lower_left = [-1.26 * 17. / 2.0] * 2\n", @@ -1289,7 +1550,7 @@ "assembly.universes = universes\n", "\n", "# Create root Cell\n", - "root_cell = openmc.Cell(name='root cell')\n", + "root_cell = openmc.Cell(name='root cell', cell_id=0)\n", "root_cell.fill = assembly\n", "\n", "# Add boundary planes\n", @@ -1316,7 +1577,7 @@ }, { "cell_type": "code", - "execution_count": 35, + "execution_count": 34, "metadata": { "collapsed": true }, @@ -1334,50 +1595,14 @@ "cell_type": "markdown", "metadata": {}, "source": [ - "Finally, lets tell OpenMC we want to tally fissions over a mesh for comparison. " + "Finally, since we want similar tally data in the end, we will leave our pre-existing `tallies.xml` file for this calculation.\n", + "\n", + "At this point, the problem is set up and we can run the multi-group calculation." ] }, { "cell_type": "code", - "execution_count": 36, - "metadata": { - "collapsed": true - }, - "outputs": [], - "source": [ - "# Instantiate a tally Mesh\n", - "mesh = openmc.Mesh(mesh_id=1)\n", - "mesh.type = 'regular'\n", - "mesh.dimension = [17, 17]\n", - "mesh.lower_left = [-10.71, -10.71]\n", - "mesh.upper_right = [+10.71, +10.71]\n", - "\n", - "# Instantiate tally Filter\n", - "mesh_filter = openmc.Filter()\n", - "mesh_filter.mesh = mesh\n", - "\n", - "# Instantiate the Tally\n", - "tally = openmc.Tally(name='mesh tally')\n", - "tally.filters = [mesh_filter]\n", - "tally.scores = ['fission']\n", - "\n", - "# Add tally to collection\n", - "tallies_file.append(tally)\n", - "\n", - "# Export all tallies to a \"tallies.xml\" file\n", - "tallies_file.export_to_xml()" - ] - }, - { - "cell_type": "markdown", - "metadata": {}, - "source": [ - "Before we run the calculation we will close the StatePoint file (as we are about to over-write it), and then we can run the multi-group calculation." - ] - }, - { - "cell_type": "code", - "execution_count": 37, + "execution_count": 35, "metadata": { "collapsed": false }, @@ -1402,8 +1627,8 @@ " Copyright: 2011-2016 Massachusetts Institute of Technology\n", " License: http://openmc.readthedocs.org/en/latest/license.html\n", " Version: 0.7.1\n", - " Git SHA1: 179e9ab147e505563d118ed58096b3d225160ffa\n", - " Date/Time: 2016-05-07 14:23:45\n", + " Git SHA1: c779ca42c41a062a6a813e03f2add2d182ca9190\n", + " Date/Time: 2016-05-12 21:19:25\n", " OpenMP Threads: 4\n", "\n", " ===========================================================================\n", @@ -1417,12 +1642,12 @@ " Reading tallies XML file...\n", " Building neighboring cells lists for each surface...\n", " Loading Cross Section Data...\n", - " Loading fuel.2g Data...\n", - " Loading fuel_clad.2g Data...\n", - " Loading fuel_mod.2g Data...\n", - " Loading gt_inmod.2g Data...\n", - " Loading gt_clad.2g Data...\n", - " Loading gt_outmod.2g Data...\n", + " Loading fuel.2m Data...\n", + " Loading fuel_clad.2m Data...\n", + " Loading fuel_mod.2m Data...\n", + " Loading gt_inmod.2m Data...\n", + " Loading gt_clad.2m Data...\n", + " Loading gt_outmod.2m Data...\n", " Initializing source particles...\n", "\n", " ===========================================================================\n", @@ -1431,207 +1656,507 @@ "\n", " Bat./Gen. k Average k \n", " ========= ======== ==================== \n", - " 1/1 1.01863 \n", - " 2/1 1.02630 \n", - " 3/1 1.03077 \n", - " 4/1 0.99715 \n", - " 5/1 1.02328 \n", - " 6/1 1.02283 \n", - " 7/1 1.00540 \n", - " 8/1 1.02232 \n", - " 9/1 0.99782 \n", - " 10/1 1.00838 \n", - " 11/1 1.01803 \n", - " 12/1 1.02530 1.02167 +/- 0.00363\n", - " 13/1 1.00514 1.01616 +/- 0.00589\n", - " 14/1 0.98994 1.00960 +/- 0.00777\n", - " 15/1 1.01028 1.00974 +/- 0.00602\n", - " 16/1 1.04607 1.01580 +/- 0.00780\n", - " 17/1 1.03300 1.01825 +/- 0.00703\n", - " 18/1 1.03149 1.01991 +/- 0.00631\n", - " 19/1 0.98692 1.01624 +/- 0.00667\n", - " 20/1 1.05205 1.01982 +/- 0.00695\n", - " 21/1 1.01572 1.01945 +/- 0.00630\n", - " 22/1 1.02517 1.01993 +/- 0.00577\n", - " 23/1 1.00274 1.01861 +/- 0.00547\n", - " 24/1 1.04739 1.02066 +/- 0.00547\n", - " 25/1 1.01883 1.02054 +/- 0.00509\n", - " 26/1 1.02021 1.02052 +/- 0.00476\n", - " 27/1 1.04696 1.02207 +/- 0.00474\n", - " 28/1 1.02751 1.02238 +/- 0.00448\n", - " 29/1 1.09537 1.02622 +/- 0.00572\n", - " 30/1 1.03685 1.02675 +/- 0.00545\n", - " 31/1 0.99812 1.02539 +/- 0.00536\n", - " 32/1 1.02526 1.02538 +/- 0.00511\n", - " 33/1 1.05466 1.02665 +/- 0.00505\n", - " 34/1 1.04816 1.02755 +/- 0.00491\n", - " 35/1 1.00148 1.02651 +/- 0.00483\n", - " 36/1 1.02315 1.02638 +/- 0.00464\n", - " 37/1 1.05771 1.02754 +/- 0.00461\n", - " 38/1 1.01675 1.02715 +/- 0.00446\n", - " 39/1 1.03707 1.02749 +/- 0.00432\n", - " 40/1 1.01903 1.02721 +/- 0.00418\n", - " 41/1 1.00332 1.02644 +/- 0.00412\n", - " 42/1 1.02533 1.02641 +/- 0.00399\n", - " 43/1 0.98531 1.02516 +/- 0.00406\n", - " 44/1 1.00406 1.02454 +/- 0.00399\n", - " 45/1 1.01057 1.02414 +/- 0.00389\n", - " 46/1 1.02755 1.02424 +/- 0.00378\n", - " 47/1 1.02783 1.02433 +/- 0.00368\n", - " 48/1 1.00003 1.02369 +/- 0.00364\n", - " 49/1 1.00442 1.02320 +/- 0.00358\n", - " 50/1 1.03215 1.02342 +/- 0.00350\n", - " 51/1 1.01672 1.02326 +/- 0.00341\n", - " 52/1 1.03702 1.02359 +/- 0.00335\n", - " 53/1 1.02063 1.02352 +/- 0.00327\n", - " 54/1 1.04596 1.02403 +/- 0.00323\n", - " 55/1 1.01926 1.02392 +/- 0.00316\n", - " 56/1 1.03058 1.02407 +/- 0.00310\n", - " 57/1 1.06126 1.02486 +/- 0.00313\n", - " 58/1 1.06411 1.02568 +/- 0.00317\n", - " 59/1 1.03278 1.02582 +/- 0.00311\n", - " 60/1 1.04472 1.02620 +/- 0.00307\n", - " 61/1 1.00186 1.02572 +/- 0.00305\n", - " 62/1 1.01133 1.02545 +/- 0.00300\n", - " 63/1 1.03713 1.02567 +/- 0.00295\n", - " 64/1 1.01363 1.02544 +/- 0.00291\n", - " 65/1 0.98126 1.02464 +/- 0.00296\n", - " 66/1 1.01500 1.02447 +/- 0.00292\n", - " 67/1 1.02437 1.02447 +/- 0.00286\n", - " 68/1 1.05057 1.02492 +/- 0.00285\n", - " 69/1 1.04903 1.02533 +/- 0.00283\n", - " 70/1 1.02199 1.02527 +/- 0.00278\n", - " 71/1 1.00536 1.02494 +/- 0.00276\n", - " 72/1 1.01658 1.02481 +/- 0.00272\n", - " 73/1 1.00866 1.02455 +/- 0.00268\n", - " 74/1 1.01800 1.02445 +/- 0.00264\n", - " 75/1 0.99176 1.02395 +/- 0.00265\n", - " 76/1 1.03336 1.02409 +/- 0.00262\n", - " 77/1 1.02699 1.02413 +/- 0.00258\n", - " 78/1 1.01596 1.02401 +/- 0.00254\n", - " 79/1 1.02292 1.02400 +/- 0.00250\n", - " 80/1 1.04804 1.02434 +/- 0.00249\n", - " 81/1 0.99494 1.02393 +/- 0.00249\n", - " 82/1 1.02646 1.02396 +/- 0.00246\n", - " 83/1 1.01223 1.02380 +/- 0.00243\n", - " 84/1 1.02572 1.02383 +/- 0.00239\n", - " 85/1 1.02709 1.02387 +/- 0.00236\n", - " 86/1 1.00315 1.02360 +/- 0.00235\n", - " 87/1 1.01809 1.02353 +/- 0.00232\n", - " 88/1 1.01566 1.02342 +/- 0.00229\n", - " 89/1 1.01093 1.02327 +/- 0.00227\n", - " 90/1 1.02812 1.02333 +/- 0.00224\n", - " 91/1 1.02288 1.02332 +/- 0.00221\n", - " 92/1 1.04070 1.02353 +/- 0.00219\n", - " 93/1 1.03697 1.02370 +/- 0.00217\n", - " 94/1 1.03486 1.02383 +/- 0.00215\n", - " 95/1 1.06359 1.02430 +/- 0.00218\n", - " 96/1 1.04811 1.02457 +/- 0.00217\n", - " 97/1 1.01303 1.02444 +/- 0.00215\n", - " 98/1 1.01243 1.02430 +/- 0.00213\n", - " 99/1 1.03238 1.02439 +/- 0.00211\n", - " 100/1 1.02054 1.02435 +/- 0.00208\n", - " 101/1 1.00402 1.02413 +/- 0.00207\n", - " 102/1 1.03800 1.02428 +/- 0.00206\n", - " 103/1 1.02541 1.02429 +/- 0.00203\n", - " 104/1 1.06867 1.02476 +/- 0.00207\n", - " 105/1 1.03192 1.02484 +/- 0.00205\n", - " 106/1 1.00100 1.02459 +/- 0.00204\n", - " 107/1 1.01098 1.02445 +/- 0.00202\n", - " 108/1 1.02930 1.02450 +/- 0.00200\n", - " 109/1 1.02173 1.02447 +/- 0.00198\n", - " 110/1 1.01411 1.02437 +/- 0.00197\n", - " 111/1 1.03920 1.02452 +/- 0.00195\n", - " 112/1 1.01984 1.02447 +/- 0.00193\n", - " 113/1 1.03912 1.02461 +/- 0.00192\n", - " 114/1 1.04124 1.02477 +/- 0.00191\n", - " 115/1 1.04802 1.02499 +/- 0.00190\n", - " 116/1 1.04129 1.02515 +/- 0.00189\n", - " 117/1 1.03072 1.02520 +/- 0.00187\n", - " 118/1 1.05167 1.02544 +/- 0.00187\n", - " 119/1 0.99954 1.02521 +/- 0.00187\n", - " 120/1 1.00093 1.02499 +/- 0.00187\n", - " 121/1 1.04929 1.02520 +/- 0.00186\n", - " 122/1 1.04556 1.02539 +/- 0.00185\n", - " 123/1 1.03298 1.02545 +/- 0.00184\n", - " 124/1 1.01603 1.02537 +/- 0.00182\n", - " 125/1 1.03522 1.02546 +/- 0.00181\n", - " 126/1 1.05644 1.02572 +/- 0.00181\n", - " 127/1 1.03754 1.02582 +/- 0.00180\n", - " 128/1 1.01524 1.02573 +/- 0.00179\n", - " 129/1 1.01263 1.02562 +/- 0.00178\n", - " 130/1 0.99835 1.02540 +/- 0.00178\n", - " 131/1 1.01268 1.02529 +/- 0.00177\n", - " 132/1 1.03975 1.02541 +/- 0.00175\n", - " 133/1 1.00702 1.02526 +/- 0.00175\n", - " 134/1 1.02335 1.02525 +/- 0.00173\n", - " 135/1 1.04378 1.02539 +/- 0.00173\n", - " 136/1 1.04610 1.02556 +/- 0.00172\n", - " 137/1 1.02284 1.02554 +/- 0.00171\n", - " 138/1 1.05720 1.02578 +/- 0.00171\n", - " 139/1 1.00965 1.02566 +/- 0.00170\n", - " 140/1 1.03719 1.02575 +/- 0.00169\n", - " 141/1 1.02413 1.02574 +/- 0.00168\n", - " 142/1 1.03125 1.02578 +/- 0.00167\n", - " 143/1 1.03641 1.02586 +/- 0.00166\n", - " 144/1 1.02137 1.02582 +/- 0.00164\n", - " 145/1 1.01522 1.02575 +/- 0.00163\n", - " 146/1 1.05163 1.02594 +/- 0.00163\n", - " 147/1 1.03612 1.02601 +/- 0.00162\n", - " 148/1 1.03346 1.02606 +/- 0.00161\n", - " 149/1 1.02306 1.02604 +/- 0.00160\n", - " 150/1 1.01764 1.02598 +/- 0.00159\n", - " 151/1 1.01787 1.02592 +/- 0.00158\n", - " 152/1 1.03263 1.02597 +/- 0.00157\n", - " 153/1 1.01877 1.02592 +/- 0.00156\n", - " 154/1 1.02870 1.02594 +/- 0.00155\n", - " 155/1 1.03071 1.02597 +/- 0.00154\n", - " 156/1 1.04229 1.02609 +/- 0.00153\n", - " 157/1 1.03973 1.02618 +/- 0.00152\n", - " 158/1 1.02180 1.02615 +/- 0.00151\n", - " 159/1 1.01067 1.02604 +/- 0.00151\n", - " 160/1 1.02888 1.02606 +/- 0.00150\n", - " 161/1 1.01711 1.02600 +/- 0.00149\n", - " 162/1 1.01087 1.02590 +/- 0.00148\n", - " 163/1 1.01886 1.02586 +/- 0.00147\n", - " 164/1 1.02210 1.02583 +/- 0.00146\n", - " 165/1 1.04020 1.02593 +/- 0.00146\n", - " 166/1 1.03658 1.02600 +/- 0.00145\n", - " 167/1 1.03222 1.02603 +/- 0.00144\n", - " 168/1 1.03247 1.02608 +/- 0.00143\n", - " 169/1 0.99739 1.02590 +/- 0.00143\n", - " 170/1 1.02464 1.02589 +/- 0.00142\n", - " 171/1 1.04623 1.02601 +/- 0.00142\n", - " 172/1 1.04328 1.02612 +/- 0.00142\n", - " 173/1 1.00812 1.02601 +/- 0.00141\n", - " 174/1 1.01224 1.02593 +/- 0.00141\n", - " 175/1 1.00882 1.02582 +/- 0.00140\n", - " 176/1 1.01286 1.02574 +/- 0.00140\n", - " 177/1 1.02048 1.02571 +/- 0.00139\n", - " 178/1 1.04269 1.02581 +/- 0.00138\n", - " 179/1 1.05862 1.02601 +/- 0.00139\n", - " 180/1 1.02924 1.02603 +/- 0.00138\n", - " 181/1 1.01491 1.02596 +/- 0.00137\n", - " 182/1 1.04255 1.02606 +/- 0.00137\n", - " 183/1 0.99191 1.02586 +/- 0.00137\n", - " 184/1 1.00392 1.02573 +/- 0.00137\n", - " 185/1 1.02982 1.02576 +/- 0.00137\n", - " 186/1 1.02682 1.02576 +/- 0.00136\n", - " 187/1 1.01484 1.02570 +/- 0.00135\n", - " 188/1 1.02825 1.02572 +/- 0.00134\n", - " 189/1 0.98954 1.02551 +/- 0.00135\n", - " 190/1 1.00522 1.02540 +/- 0.00135\n", - " 191/1 1.03762 1.02547 +/- 0.00134\n", - " 192/1 1.02091 1.02544 +/- 0.00134\n", - " 193/1 1.04549 1.02555 +/- 0.00133\n", - " 194/1 1.05531 1.02572 +/- 0.00134\n", - " 195/1 1.01479 1.02566 +/- 0.00133\n", - " 196/1 1.01337 1.02559 +/- 0.00132\n", - " 197/1 0.99187 1.02541 +/- 0.00133\n", - " 198/1 1.01280 1.02534 +/- 0.00132\n", - " 199/1 1.00049 1.02521 +/- 0.00132\n", - " 200/1 1.01879 1.02518 +/- 0.00132\n", - " Creating state point statepoint.200.h5...\n", + " 1/1 1.01702 \n", + " 2/1 0.99463 \n", + " 3/1 1.02321 \n", + " 4/1 0.98628 \n", + " 5/1 1.03122 \n", + " 6/1 1.00774 \n", + " 7/1 1.05616 \n", + " 8/1 1.03051 \n", + " 9/1 1.02321 \n", + " 10/1 1.04380 \n", + " 11/1 1.05837 \n", + " 12/1 1.01514 1.03676 +/- 0.02161\n", + " 13/1 1.06720 1.04690 +/- 0.01608\n", + " 14/1 1.01696 1.03942 +/- 0.01361\n", + " 15/1 1.03549 1.03863 +/- 0.01057\n", + " 16/1 1.01599 1.03486 +/- 0.00942\n", + " 17/1 1.03070 1.03427 +/- 0.00799\n", + " 18/1 1.03778 1.03470 +/- 0.00693\n", + " 19/1 1.03042 1.03423 +/- 0.00613\n", + " 20/1 1.01047 1.03185 +/- 0.00598\n", + " 21/1 1.03251 1.03191 +/- 0.00541\n", + " 22/1 1.02047 1.03096 +/- 0.00503\n", + " 23/1 1.01729 1.02991 +/- 0.00474\n", + " 24/1 1.02948 1.02988 +/- 0.00439\n", + " 25/1 1.01963 1.02919 +/- 0.00414\n", + " 26/1 1.00626 1.02776 +/- 0.00413\n", + " 27/1 1.04531 1.02879 +/- 0.00402\n", + " 28/1 0.99936 1.02716 +/- 0.00412\n", + " 29/1 1.04497 1.02809 +/- 0.00401\n", + " 30/1 1.02429 1.02790 +/- 0.00381\n", + " 31/1 1.05112 1.02901 +/- 0.00379\n", + " 32/1 1.01843 1.02853 +/- 0.00365\n", + " 33/1 1.04478 1.02924 +/- 0.00355\n", + " 34/1 1.01719 1.02873 +/- 0.00344\n", + " 35/1 0.99873 1.02753 +/- 0.00351\n", + " 36/1 1.00054 1.02649 +/- 0.00353\n", + " 37/1 1.03986 1.02699 +/- 0.00343\n", + " 38/1 1.02243 1.02683 +/- 0.00331\n", + " 39/1 1.02744 1.02685 +/- 0.00319\n", + " 40/1 1.01174 1.02634 +/- 0.00313\n", + " 41/1 1.04973 1.02710 +/- 0.00312\n", + " 42/1 0.99564 1.02612 +/- 0.00317\n", + " 43/1 1.03022 1.02624 +/- 0.00308\n", + " 44/1 1.03526 1.02650 +/- 0.00300\n", + " 45/1 1.02143 1.02636 +/- 0.00292\n", + " 46/1 1.03264 1.02653 +/- 0.00284\n", + " 47/1 1.03868 1.02686 +/- 0.00278\n", + " 48/1 1.02385 1.02678 +/- 0.00271\n", + " 49/1 1.03897 1.02710 +/- 0.00266\n", + " 50/1 1.01267 1.02674 +/- 0.00261\n", + " 51/1 0.99683 1.02601 +/- 0.00265\n", + " 52/1 1.04189 1.02638 +/- 0.00261\n", + " 53/1 1.02871 1.02644 +/- 0.00255\n", + " 54/1 1.02564 1.02642 +/- 0.00250\n", + " 55/1 1.02955 1.02649 +/- 0.00244\n", + " 56/1 1.02390 1.02643 +/- 0.00239\n", + " 57/1 1.03342 1.02658 +/- 0.00234\n", + " 58/1 1.01430 1.02633 +/- 0.00231\n", + " 59/1 0.99242 1.02563 +/- 0.00236\n", + " 60/1 1.00442 1.02521 +/- 0.00235\n", + " 61/1 1.03870 1.02547 +/- 0.00232\n", + " 62/1 1.02146 1.02540 +/- 0.00228\n", + " 63/1 1.04782 1.02582 +/- 0.00227\n", + " 64/1 1.02872 1.02587 +/- 0.00223\n", + " 65/1 1.02420 1.02584 +/- 0.00219\n", + " 66/1 1.01974 1.02573 +/- 0.00215\n", + " 67/1 1.00774 1.02542 +/- 0.00214\n", + " 68/1 1.01323 1.02521 +/- 0.00211\n", + " 69/1 1.01468 1.02503 +/- 0.00208\n", + " 70/1 1.02869 1.02509 +/- 0.00205\n", + " 71/1 1.02284 1.02505 +/- 0.00202\n", + " 72/1 1.04815 1.02543 +/- 0.00202\n", + " 73/1 1.01119 1.02520 +/- 0.00200\n", + " 74/1 1.03314 1.02533 +/- 0.00197\n", + " 75/1 1.02333 1.02529 +/- 0.00194\n", + " 76/1 1.04030 1.02552 +/- 0.00193\n", + " 77/1 1.02537 1.02552 +/- 0.00190\n", + " 78/1 1.02875 1.02557 +/- 0.00187\n", + " 79/1 1.03588 1.02572 +/- 0.00185\n", + " 80/1 1.05250 1.02610 +/- 0.00186\n", + " 81/1 1.00477 1.02580 +/- 0.00186\n", + " 82/1 1.03903 1.02598 +/- 0.00184\n", + " 83/1 1.02378 1.02595 +/- 0.00182\n", + " 84/1 1.01107 1.02575 +/- 0.00180\n", + " 85/1 1.01550 1.02561 +/- 0.00178\n", + " 86/1 1.00540 1.02535 +/- 0.00178\n", + " 87/1 1.03056 1.02542 +/- 0.00176\n", + " 88/1 1.01742 1.02531 +/- 0.00174\n", + " 89/1 0.99730 1.02496 +/- 0.00175\n", + " 90/1 1.03569 1.02509 +/- 0.00174\n", + " 91/1 1.04514 1.02534 +/- 0.00173\n", + " 92/1 1.02757 1.02537 +/- 0.00171\n", + " 93/1 1.00610 1.02514 +/- 0.00171\n", + " 94/1 1.03576 1.02526 +/- 0.00169\n", + " 95/1 1.03732 1.02540 +/- 0.00168\n", + " 96/1 1.04784 1.02567 +/- 0.00168\n", + " 97/1 1.06507 1.02612 +/- 0.00172\n", + " 98/1 1.03673 1.02624 +/- 0.00170\n", + " 99/1 1.01270 1.02609 +/- 0.00169\n", + " 100/1 1.01980 1.02602 +/- 0.00167\n", + " 101/1 1.01357 1.02588 +/- 0.00166\n", + " 102/1 1.03125 1.02594 +/- 0.00164\n", + " 103/1 1.01527 1.02582 +/- 0.00163\n", + " 104/1 1.02403 1.02580 +/- 0.00161\n", + " 105/1 1.03435 1.02589 +/- 0.00160\n", + " 106/1 1.04113 1.02605 +/- 0.00159\n", + " 107/1 1.03291 1.02612 +/- 0.00157\n", + " 108/1 1.02478 1.02611 +/- 0.00156\n", + " 109/1 1.05814 1.02643 +/- 0.00158\n", + " 110/1 1.02647 1.02643 +/- 0.00156\n", + " 111/1 0.98951 1.02607 +/- 0.00159\n", + " 112/1 1.00739 1.02589 +/- 0.00158\n", + " 113/1 1.04165 1.02604 +/- 0.00157\n", + " 114/1 1.00047 1.02579 +/- 0.00158\n", + " 115/1 1.02550 1.02579 +/- 0.00156\n", + " 116/1 1.02408 1.02577 +/- 0.00155\n", + " 117/1 1.03110 1.02582 +/- 0.00153\n", + " 118/1 1.02874 1.02585 +/- 0.00152\n", + " 119/1 1.02348 1.02583 +/- 0.00151\n", + " 120/1 1.01969 1.02577 +/- 0.00149\n", + " 121/1 1.02312 1.02575 +/- 0.00148\n", + " 122/1 1.03261 1.02581 +/- 0.00147\n", + " 123/1 0.98394 1.02544 +/- 0.00150\n", + " 124/1 1.03771 1.02555 +/- 0.00149\n", + " 125/1 1.01857 1.02549 +/- 0.00148\n", + " 126/1 1.00066 1.02527 +/- 0.00148\n", + " 127/1 1.02372 1.02526 +/- 0.00147\n", + " 128/1 1.03307 1.02533 +/- 0.00146\n", + " 129/1 1.00889 1.02519 +/- 0.00145\n", + " 130/1 1.02053 1.02515 +/- 0.00144\n", + " 131/1 1.00943 1.02502 +/- 0.00144\n", + " 132/1 1.07225 1.02541 +/- 0.00148\n", + " 133/1 1.04068 1.02553 +/- 0.00147\n", + " 134/1 1.03509 1.02561 +/- 0.00146\n", + " 135/1 1.01250 1.02550 +/- 0.00145\n", + " 136/1 1.02179 1.02547 +/- 0.00144\n", + " 137/1 1.05685 1.02572 +/- 0.00145\n", + " 138/1 1.04217 1.02585 +/- 0.00144\n", + " 139/1 1.02793 1.02586 +/- 0.00143\n", + " 140/1 1.01207 1.02576 +/- 0.00143\n", + " 141/1 1.03445 1.02582 +/- 0.00142\n", + " 142/1 1.03579 1.02590 +/- 0.00141\n", + " 143/1 1.00786 1.02576 +/- 0.00140\n", + " 144/1 0.99089 1.02550 +/- 0.00142\n", + " 145/1 1.02617 1.02551 +/- 0.00141\n", + " 146/1 1.01691 1.02545 +/- 0.00140\n", + " 147/1 1.00692 1.02531 +/- 0.00139\n", + " 148/1 0.97702 1.02496 +/- 0.00143\n", + " 149/1 1.04002 1.02507 +/- 0.00142\n", + " 150/1 1.01262 1.02498 +/- 0.00141\n", + " 151/1 1.03613 1.02506 +/- 0.00141\n", + " 152/1 1.02920 1.02509 +/- 0.00140\n", + " 153/1 1.02199 1.02507 +/- 0.00139\n", + " 154/1 1.03421 1.02513 +/- 0.00138\n", + " 155/1 1.05882 1.02536 +/- 0.00139\n", + " 156/1 1.02649 1.02537 +/- 0.00138\n", + " 157/1 1.01933 1.02533 +/- 0.00137\n", + " 158/1 1.04269 1.02545 +/- 0.00137\n", + " 159/1 0.99604 1.02525 +/- 0.00137\n", + " 160/1 1.04748 1.02540 +/- 0.00137\n", + " 161/1 1.00501 1.02526 +/- 0.00137\n", + " 162/1 1.00550 1.02513 +/- 0.00137\n", + " 163/1 1.00115 1.02498 +/- 0.00137\n", + " 164/1 1.02283 1.02496 +/- 0.00136\n", + " 165/1 1.01964 1.02493 +/- 0.00135\n", + " 166/1 1.02287 1.02491 +/- 0.00134\n", + " 167/1 1.05498 1.02511 +/- 0.00134\n", + " 168/1 1.05267 1.02528 +/- 0.00135\n", + " 169/1 1.00474 1.02515 +/- 0.00135\n", + " 170/1 1.03469 1.02521 +/- 0.00134\n", + " 171/1 1.02499 1.02521 +/- 0.00133\n", + " 172/1 1.03961 1.02530 +/- 0.00132\n", + " 173/1 1.01240 1.02522 +/- 0.00132\n", + " 174/1 1.00762 1.02511 +/- 0.00132\n", + " 175/1 1.00200 1.02497 +/- 0.00131\n", + " 176/1 1.01449 1.02491 +/- 0.00131\n", + " 177/1 1.01111 1.02483 +/- 0.00130\n", + " 178/1 1.01208 1.02475 +/- 0.00130\n", + " 179/1 1.03304 1.02480 +/- 0.00129\n", + " 180/1 1.04504 1.02492 +/- 0.00129\n", + " 181/1 1.03476 1.02498 +/- 0.00128\n", + " 182/1 1.02124 1.02495 +/- 0.00128\n", + " 183/1 0.98855 1.02474 +/- 0.00128\n", + " 184/1 1.04689 1.02487 +/- 0.00128\n", + " 185/1 1.00618 1.02476 +/- 0.00128\n", + " 186/1 1.02012 1.02474 +/- 0.00127\n", + " 187/1 1.00162 1.02461 +/- 0.00127\n", + " 188/1 1.03269 1.02465 +/- 0.00127\n", + " 189/1 1.04772 1.02478 +/- 0.00127\n", + " 190/1 1.01132 1.02471 +/- 0.00126\n", + " 191/1 1.02669 1.02472 +/- 0.00125\n", + " 192/1 1.01154 1.02464 +/- 0.00125\n", + " 193/1 1.05795 1.02483 +/- 0.00126\n", + " 194/1 1.01615 1.02478 +/- 0.00125\n", + " 195/1 1.03828 1.02485 +/- 0.00125\n", + " 196/1 1.00695 1.02476 +/- 0.00124\n", + " 197/1 1.04126 1.02484 +/- 0.00124\n", + " 198/1 1.02834 1.02486 +/- 0.00123\n", + " 199/1 1.01000 1.02478 +/- 0.00123\n", + " 200/1 0.99294 1.02462 +/- 0.00123\n", + " 201/1 1.00248 1.02450 +/- 0.00123\n", + " 202/1 1.03461 1.02455 +/- 0.00123\n", + " 203/1 1.06289 1.02475 +/- 0.00124\n", + " 204/1 1.03010 1.02478 +/- 0.00123\n", + " 205/1 1.04636 1.02489 +/- 0.00123\n", + " 206/1 1.05434 1.02504 +/- 0.00123\n", + " 207/1 1.03993 1.02512 +/- 0.00123\n", + " 208/1 1.02672 1.02512 +/- 0.00122\n", + " 209/1 1.04958 1.02525 +/- 0.00122\n", + " 210/1 0.99194 1.02508 +/- 0.00123\n", + " 211/1 1.01570 1.02503 +/- 0.00122\n", + " 212/1 1.04079 1.02511 +/- 0.00122\n", + " 213/1 1.02961 1.02513 +/- 0.00121\n", + " 214/1 1.03797 1.02520 +/- 0.00121\n", + " 215/1 1.03714 1.02526 +/- 0.00120\n", + " 216/1 1.03299 1.02529 +/- 0.00120\n", + " 217/1 1.00461 1.02519 +/- 0.00120\n", + " 218/1 1.02386 1.02519 +/- 0.00119\n", + " 219/1 1.01955 1.02516 +/- 0.00119\n", + " 220/1 1.04372 1.02525 +/- 0.00118\n", + " 221/1 1.01694 1.02521 +/- 0.00118\n", + " 222/1 0.99642 1.02507 +/- 0.00118\n", + " 223/1 1.00999 1.02500 +/- 0.00118\n", + " 224/1 1.02703 1.02501 +/- 0.00117\n", + " 225/1 1.00236 1.02491 +/- 0.00117\n", + " 226/1 1.02825 1.02492 +/- 0.00117\n", + " 227/1 1.04535 1.02502 +/- 0.00116\n", + " 228/1 1.01779 1.02498 +/- 0.00116\n", + " 229/1 1.01058 1.02492 +/- 0.00116\n", + " 230/1 1.00391 1.02482 +/- 0.00115\n", + " 231/1 1.05990 1.02498 +/- 0.00116\n", + " 232/1 1.01885 1.02495 +/- 0.00116\n", + " 233/1 1.03204 1.02498 +/- 0.00115\n", + " 234/1 0.99396 1.02485 +/- 0.00115\n", + " 235/1 1.01828 1.02482 +/- 0.00115\n", + " 236/1 1.08225 1.02507 +/- 0.00117\n", + " 237/1 1.00335 1.02498 +/- 0.00117\n", + " 238/1 1.03097 1.02500 +/- 0.00117\n", + " 239/1 1.01738 1.02497 +/- 0.00116\n", + " 240/1 1.02261 1.02496 +/- 0.00116\n", + " 241/1 1.02814 1.02497 +/- 0.00115\n", + " 242/1 1.01158 1.02491 +/- 0.00115\n", + " 243/1 1.03507 1.02496 +/- 0.00114\n", + " 244/1 1.01914 1.02493 +/- 0.00114\n", + " 245/1 1.04555 1.02502 +/- 0.00114\n", + " 246/1 1.02459 1.02502 +/- 0.00113\n", + " 247/1 1.05827 1.02516 +/- 0.00114\n", + " 248/1 1.02549 1.02516 +/- 0.00113\n", + " 249/1 1.03354 1.02520 +/- 0.00113\n", + " 250/1 1.04186 1.02526 +/- 0.00113\n", + " 251/1 1.00466 1.02518 +/- 0.00112\n", + " 252/1 0.99065 1.02504 +/- 0.00113\n", + " 253/1 1.03065 1.02506 +/- 0.00112\n", + " 254/1 1.02167 1.02505 +/- 0.00112\n", + " 255/1 1.01700 1.02501 +/- 0.00112\n", + " 256/1 1.03619 1.02506 +/- 0.00111\n", + " 257/1 1.01833 1.02503 +/- 0.00111\n", + " 258/1 1.02211 1.02502 +/- 0.00110\n", + " 259/1 1.04348 1.02509 +/- 0.00110\n", + " 260/1 1.03444 1.02513 +/- 0.00110\n", + " 261/1 1.05597 1.02525 +/- 0.00110\n", + " 262/1 1.02085 1.02524 +/- 0.00110\n", + " 263/1 1.00552 1.02516 +/- 0.00109\n", + " 264/1 1.03976 1.02522 +/- 0.00109\n", + " 265/1 1.02810 1.02523 +/- 0.00109\n", + " 266/1 1.00911 1.02516 +/- 0.00108\n", + " 267/1 1.01963 1.02514 +/- 0.00108\n", + " 268/1 1.03732 1.02519 +/- 0.00108\n", + " 269/1 1.02422 1.02519 +/- 0.00107\n", + " 270/1 1.01546 1.02515 +/- 0.00107\n", + " 271/1 1.05488 1.02526 +/- 0.00107\n", + " 272/1 1.01709 1.02523 +/- 0.00107\n", + " 273/1 1.05629 1.02535 +/- 0.00107\n", + " 274/1 1.03864 1.02540 +/- 0.00107\n", + " 275/1 1.01472 1.02536 +/- 0.00106\n", + " 276/1 1.03425 1.02539 +/- 0.00106\n", + " 277/1 1.00663 1.02532 +/- 0.00106\n", + " 278/1 1.03326 1.02535 +/- 0.00106\n", + " 279/1 1.02571 1.02535 +/- 0.00105\n", + " 280/1 1.00525 1.02528 +/- 0.00105\n", + " 281/1 1.00451 1.02520 +/- 0.00105\n", + " 282/1 1.04016 1.02526 +/- 0.00105\n", + " 283/1 0.98343 1.02510 +/- 0.00105\n", + " 284/1 1.04843 1.02519 +/- 0.00105\n", + " 285/1 1.01807 1.02516 +/- 0.00105\n", + " 286/1 1.02393 1.02516 +/- 0.00105\n", + " 287/1 1.01851 1.02514 +/- 0.00104\n", + " 288/1 1.03976 1.02519 +/- 0.00104\n", + " 289/1 1.03153 1.02521 +/- 0.00104\n", + " 290/1 1.00416 1.02514 +/- 0.00104\n", + " 291/1 1.01426 1.02510 +/- 0.00103\n", + " 292/1 1.02583 1.02510 +/- 0.00103\n", + " 293/1 1.01680 1.02507 +/- 0.00103\n", + " 294/1 1.04578 1.02514 +/- 0.00103\n", + " 295/1 1.03162 1.02517 +/- 0.00102\n", + " 296/1 1.01682 1.02514 +/- 0.00102\n", + " 297/1 1.00488 1.02507 +/- 0.00102\n", + " 298/1 1.03057 1.02508 +/- 0.00101\n", + " 299/1 1.01126 1.02504 +/- 0.00101\n", + " 300/1 1.03528 1.02507 +/- 0.00101\n", + " 301/1 1.05548 1.02518 +/- 0.00101\n", + " 302/1 1.02994 1.02519 +/- 0.00101\n", + " 303/1 1.03010 1.02521 +/- 0.00100\n", + " 304/1 1.04031 1.02526 +/- 0.00100\n", + " 305/1 1.05866 1.02537 +/- 0.00101\n", + " 306/1 1.03602 1.02541 +/- 0.00100\n", + " 307/1 1.01362 1.02537 +/- 0.00100\n", + " 308/1 1.01318 1.02533 +/- 0.00100\n", + " 309/1 1.04262 1.02539 +/- 0.00100\n", + " 310/1 1.01626 1.02536 +/- 0.00099\n", + " 311/1 1.00285 1.02528 +/- 0.00099\n", + " 312/1 0.98155 1.02514 +/- 0.00100\n", + " 313/1 1.05649 1.02524 +/- 0.00100\n", + " 314/1 1.00960 1.02519 +/- 0.00100\n", + " 315/1 1.05350 1.02528 +/- 0.00100\n", + " 316/1 1.03842 1.02533 +/- 0.00100\n", + " 317/1 1.01394 1.02529 +/- 0.00100\n", + " 318/1 1.01830 1.02527 +/- 0.00099\n", + " 319/1 1.02050 1.02525 +/- 0.00099\n", + " 320/1 1.03402 1.02528 +/- 0.00099\n", + " 321/1 1.04547 1.02534 +/- 0.00099\n", + " 322/1 1.02579 1.02534 +/- 0.00098\n", + " 323/1 1.01922 1.02533 +/- 0.00098\n", + " 324/1 1.01050 1.02528 +/- 0.00098\n", + " 325/1 1.01426 1.02524 +/- 0.00098\n", + " 326/1 1.03283 1.02527 +/- 0.00097\n", + " 327/1 1.03859 1.02531 +/- 0.00097\n", + " 328/1 1.01536 1.02528 +/- 0.00097\n", + " 329/1 1.03149 1.02530 +/- 0.00097\n", + " 330/1 1.04328 1.02535 +/- 0.00096\n", + " 331/1 1.01949 1.02534 +/- 0.00096\n", + " 332/1 1.02319 1.02533 +/- 0.00096\n", + " 333/1 1.01704 1.02530 +/- 0.00096\n", + " 334/1 1.02691 1.02531 +/- 0.00095\n", + " 335/1 1.03188 1.02533 +/- 0.00095\n", + " 336/1 1.03107 1.02535 +/- 0.00095\n", + " 337/1 1.02410 1.02534 +/- 0.00094\n", + " 338/1 0.99917 1.02526 +/- 0.00094\n", + " 339/1 1.03593 1.02529 +/- 0.00094\n", + " 340/1 1.02286 1.02529 +/- 0.00094\n", + " 341/1 1.04154 1.02534 +/- 0.00094\n", + " 342/1 1.01664 1.02531 +/- 0.00094\n", + " 343/1 1.01041 1.02527 +/- 0.00093\n", + " 344/1 1.02033 1.02525 +/- 0.00093\n", + " 345/1 1.03137 1.02527 +/- 0.00093\n", + " 346/1 1.02162 1.02526 +/- 0.00093\n", + " 347/1 1.00835 1.02521 +/- 0.00092\n", + " 348/1 1.01168 1.02517 +/- 0.00092\n", + " 349/1 1.01168 1.02513 +/- 0.00092\n", + " 350/1 1.03509 1.02516 +/- 0.00092\n", + " 351/1 1.01883 1.02514 +/- 0.00092\n", + " 352/1 1.04314 1.02519 +/- 0.00091\n", + " 353/1 0.99067 1.02509 +/- 0.00092\n", + " 354/1 1.03100 1.02511 +/- 0.00091\n", + " 355/1 1.01664 1.02508 +/- 0.00091\n", + " 356/1 1.02193 1.02507 +/- 0.00091\n", + " 357/1 1.03213 1.02509 +/- 0.00091\n", + " 358/1 1.00555 1.02504 +/- 0.00091\n", + " 359/1 1.04849 1.02511 +/- 0.00091\n", + " 360/1 1.02174 1.02510 +/- 0.00090\n", + " 361/1 1.05064 1.02517 +/- 0.00090\n", + " 362/1 1.05274 1.02525 +/- 0.00091\n", + " 363/1 1.00932 1.02520 +/- 0.00090\n", + " 364/1 1.03400 1.02523 +/- 0.00090\n", + " 365/1 1.00149 1.02516 +/- 0.00090\n", + " 366/1 1.01631 1.02514 +/- 0.00090\n", + " 367/1 1.03928 1.02517 +/- 0.00090\n", + " 368/1 1.01318 1.02514 +/- 0.00090\n", + " 369/1 1.04610 1.02520 +/- 0.00090\n", + " 370/1 1.04338 1.02525 +/- 0.00089\n", + " 371/1 1.01638 1.02523 +/- 0.00089\n", + " 372/1 1.04056 1.02527 +/- 0.00089\n", + " 373/1 1.00090 1.02520 +/- 0.00089\n", + " 374/1 1.01261 1.02517 +/- 0.00089\n", + " 375/1 1.03919 1.02520 +/- 0.00089\n", + " 376/1 0.99900 1.02513 +/- 0.00089\n", + " 377/1 1.00168 1.02507 +/- 0.00089\n", + " 378/1 0.99476 1.02499 +/- 0.00089\n", + " 379/1 1.04960 1.02505 +/- 0.00089\n", + " 380/1 0.99797 1.02498 +/- 0.00089\n", + " 381/1 1.04956 1.02505 +/- 0.00089\n", + " 382/1 1.02803 1.02505 +/- 0.00089\n", + " 383/1 0.99388 1.02497 +/- 0.00089\n", + " 384/1 1.00767 1.02492 +/- 0.00089\n", + " 385/1 1.00856 1.02488 +/- 0.00089\n", + " 386/1 1.02997 1.02489 +/- 0.00088\n", + " 387/1 0.97841 1.02477 +/- 0.00089\n", + " 388/1 0.99712 1.02470 +/- 0.00089\n", + " 389/1 0.99072 1.02461 +/- 0.00089\n", + " 390/1 1.02439 1.02461 +/- 0.00089\n", + " 391/1 1.02769 1.02462 +/- 0.00089\n", + " 392/1 1.02205 1.02461 +/- 0.00089\n", + " 393/1 1.03702 1.02464 +/- 0.00088\n", + " 394/1 1.00274 1.02458 +/- 0.00088\n", + " 395/1 1.00131 1.02452 +/- 0.00088\n", + " 396/1 1.00130 1.02446 +/- 0.00088\n", + " 397/1 1.00472 1.02441 +/- 0.00088\n", + " 398/1 1.00724 1.02437 +/- 0.00088\n", + " 399/1 1.03061 1.02438 +/- 0.00088\n", + " 400/1 0.99651 1.02431 +/- 0.00088\n", + " 401/1 0.99290 1.02423 +/- 0.00088\n", + " 402/1 1.02166 1.02423 +/- 0.00088\n", + " 403/1 1.01691 1.02421 +/- 0.00088\n", + " 404/1 1.00492 1.02416 +/- 0.00088\n", + " 405/1 1.00663 1.02411 +/- 0.00088\n", + " 406/1 1.01865 1.02410 +/- 0.00087\n", + " 407/1 1.02717 1.02411 +/- 0.00087\n", + " 408/1 1.01793 1.02409 +/- 0.00087\n", + " 409/1 1.02606 1.02410 +/- 0.00087\n", + " 410/1 1.03809 1.02413 +/- 0.00087\n", + " 411/1 1.03780 1.02417 +/- 0.00086\n", + " 412/1 1.02782 1.02418 +/- 0.00086\n", + " 413/1 1.03077 1.02419 +/- 0.00086\n", + " 414/1 1.00651 1.02415 +/- 0.00086\n", + " 415/1 1.05594 1.02423 +/- 0.00086\n", + " 416/1 0.99558 1.02416 +/- 0.00086\n", + " 417/1 1.00689 1.02411 +/- 0.00086\n", + " 418/1 1.02932 1.02413 +/- 0.00086\n", + " 419/1 1.03552 1.02415 +/- 0.00086\n", + " 420/1 1.03735 1.02419 +/- 0.00085\n", + " 421/1 1.02402 1.02419 +/- 0.00085\n", + " 422/1 1.04227 1.02423 +/- 0.00085\n", + " 423/1 1.03087 1.02425 +/- 0.00085\n", + " 424/1 1.04363 1.02429 +/- 0.00085\n", + " 425/1 1.02676 1.02430 +/- 0.00085\n", + " 426/1 1.03739 1.02433 +/- 0.00085\n", + " 427/1 1.02977 1.02434 +/- 0.00084\n", + " 428/1 1.02547 1.02435 +/- 0.00084\n", + " 429/1 1.03552 1.02437 +/- 0.00084\n", + " 430/1 1.04282 1.02442 +/- 0.00084\n", + " 431/1 1.03171 1.02443 +/- 0.00084\n", + " 432/1 1.01030 1.02440 +/- 0.00084\n", + " 433/1 1.04168 1.02444 +/- 0.00084\n", + " 434/1 0.98994 1.02436 +/- 0.00084\n", + " 435/1 0.98166 1.02426 +/- 0.00084\n", + " 436/1 1.00178 1.02421 +/- 0.00084\n", + " 437/1 1.03801 1.02424 +/- 0.00084\n", + " 438/1 1.02099 1.02423 +/- 0.00084\n", + " 439/1 1.01305 1.02421 +/- 0.00084\n", + " 440/1 1.02286 1.02420 +/- 0.00083\n", + " 441/1 1.03697 1.02423 +/- 0.00083\n", + " 442/1 0.99050 1.02415 +/- 0.00083\n", + " 443/1 1.02238 1.02415 +/- 0.00083\n", + " 444/1 1.05188 1.02421 +/- 0.00083\n", + " 445/1 1.03150 1.02423 +/- 0.00083\n", + " 446/1 1.01071 1.02420 +/- 0.00083\n", + " 447/1 1.03713 1.02423 +/- 0.00083\n", + " 448/1 1.03631 1.02426 +/- 0.00083\n", + " 449/1 1.02968 1.02427 +/- 0.00083\n", + " 450/1 1.03031 1.02428 +/- 0.00082\n", + " 451/1 1.02161 1.02428 +/- 0.00082\n", + " 452/1 0.99036 1.02420 +/- 0.00082\n", + " 453/1 1.02581 1.02420 +/- 0.00082\n", + " 454/1 1.03140 1.02422 +/- 0.00082\n", + " 455/1 1.01962 1.02421 +/- 0.00082\n", + " 456/1 1.00680 1.02417 +/- 0.00082\n", + " 457/1 1.00178 1.02412 +/- 0.00082\n", + " 458/1 1.02306 1.02412 +/- 0.00082\n", + " 459/1 1.02653 1.02412 +/- 0.00081\n", + " 460/1 1.02934 1.02413 +/- 0.00081\n", + " 461/1 1.00872 1.02410 +/- 0.00081\n", + " 462/1 1.00012 1.02405 +/- 0.00081\n", + " 463/1 0.99057 1.02397 +/- 0.00081\n", + " 464/1 1.02353 1.02397 +/- 0.00081\n", + " 465/1 1.01402 1.02395 +/- 0.00081\n", + " 466/1 1.01651 1.02393 +/- 0.00081\n", + " 467/1 1.01024 1.02390 +/- 0.00081\n", + " 468/1 1.02504 1.02391 +/- 0.00080\n", + " 469/1 1.00891 1.02387 +/- 0.00080\n", + " 470/1 1.04038 1.02391 +/- 0.00080\n", + " 471/1 1.04346 1.02395 +/- 0.00080\n", + " 472/1 1.02634 1.02396 +/- 0.00080\n", + " 473/1 1.01207 1.02393 +/- 0.00080\n", + " 474/1 1.00787 1.02390 +/- 0.00080\n", + " 475/1 1.03591 1.02392 +/- 0.00080\n", + " 476/1 1.04257 1.02396 +/- 0.00080\n", + " 477/1 1.00536 1.02392 +/- 0.00079\n", + " 478/1 1.07545 1.02403 +/- 0.00080\n", + " 479/1 1.02306 1.02403 +/- 0.00080\n", + " 480/1 1.02733 1.02404 +/- 0.00080\n", + " 481/1 1.00990 1.02401 +/- 0.00080\n", + " 482/1 0.99031 1.02394 +/- 0.00080\n", + " 483/1 0.98006 1.02384 +/- 0.00080\n", + " 484/1 1.05635 1.02391 +/- 0.00080\n", + " 485/1 1.02410 1.02391 +/- 0.00080\n", + " 486/1 1.01227 1.02389 +/- 0.00080\n", + " 487/1 1.00614 1.02385 +/- 0.00080\n", + " 488/1 1.01837 1.02384 +/- 0.00080\n", + " 489/1 1.02565 1.02384 +/- 0.00080\n", + " 490/1 1.00530 1.02381 +/- 0.00079\n", + " 491/1 1.01958 1.02380 +/- 0.00079\n", + " 492/1 1.04490 1.02384 +/- 0.00079\n", + " 493/1 1.02567 1.02384 +/- 0.00079\n", + " 494/1 1.03865 1.02387 +/- 0.00079\n", + " 495/1 1.03990 1.02391 +/- 0.00079\n", + " 496/1 0.98352 1.02382 +/- 0.00079\n", + " 497/1 1.00909 1.02379 +/- 0.00079\n", + " 498/1 1.03661 1.02382 +/- 0.00079\n", + " 499/1 1.04423 1.02386 +/- 0.00079\n", + " 500/1 1.06406 1.02394 +/- 0.00079\n", + " Creating state point statepoint.500.h5...\n", "\n", " ===========================================================================\n", " ======================> SIMULATION FINISHED <======================\n", @@ -1640,27 +2165,27 @@ "\n", " =======================> TIMING STATISTICS <=======================\n", "\n", - " Total time for initialization = 6.3000E-02 seconds\n", + " Total time for initialization = 5.3000E-02 seconds\n", " Reading cross sections = 5.0000E-03 seconds\n", - " Total time in simulation = 7.3280E+01 seconds\n", - " Time in transport only = 7.3104E+01 seconds\n", - " Time in inactive batches = 1.1200E+00 seconds\n", - " Time in active batches = 7.2160E+01 seconds\n", - " Time synchronizing fission bank = 2.5000E-02 seconds\n", - " Sampling source sites = 1.8000E-02 seconds\n", - " SEND/RECV source sites = 7.0000E-03 seconds\n", - " Time accumulating tallies = 0.0000E+00 seconds\n", + " Total time in simulation = 1.8631E+02 seconds\n", + " Time in transport only = 1.8590E+02 seconds\n", + " Time in inactive batches = 1.1710E+00 seconds\n", + " Time in active batches = 1.8514E+02 seconds\n", + " Time synchronizing fission bank = 7.3000E-02 seconds\n", + " Sampling source sites = 5.1000E-02 seconds\n", + " SEND/RECV source sites = 2.2000E-02 seconds\n", + " Time accumulating tallies = 4.0000E-03 seconds\n", " Total time for finalization = 0.0000E+00 seconds\n", - " Total time elapsed = 7.3353E+01 seconds\n", - " Calculation Rate (inactive) = 44642.9 neutrons/second\n", - " Calculation Rate (active) = 13165.2 neutrons/second\n", + " Total time elapsed = 1.8637E+02 seconds\n", + " Calculation Rate (inactive) = 42698.5 neutrons/second\n", + " Calculation Rate (active) = 13233.2 neutrons/second\n", "\n", " ============================> RESULTS <============================\n", "\n", - " k-effective (Collision) = 1.02597 +/- 0.00117\n", - " k-effective (Track-length) = 1.02518 +/- 0.00132\n", - " k-effective (Absorption) = 1.02581 +/- 0.00070\n", - " Combined k-effective = 1.02562 +/- 0.00068\n", + " k-effective (Collision) = 1.02403 +/- 0.00071\n", + " k-effective (Track-length) = 1.02394 +/- 0.00079\n", + " k-effective (Absorption) = 1.02539 +/- 0.00044\n", + " Combined k-effective = 1.02518 +/- 0.00042\n", " Leakage Fraction = 0.00000 +/- 0.00000\n", "\n" ] @@ -1671,15 +2196,12 @@ "0" ] }, - "execution_count": 37, + "execution_count": 35, "metadata": {}, "output_type": "execute_result" } ], "source": [ - "# Close the StatePoint File\n", - "sp._f.close()\n", - "\n", "# Run the Multi-Group OpenMC Simulation\n", "openmc.run()" ] @@ -1691,12 +2213,13 @@ "# Results Comparison\n", "Now we can compare the multi-group and continuous-energy results.\n", "\n", - "We will begin by loading the multi-group statepoint file we just finished writing and extracting the calculated keff." + "We will begin by loading the multi-group statepoint file we just finished writing and extracting the calculated keff.\n", + "Since we did not rename the summary file, we do not need to load it separately this time." ] }, { "cell_type": "code", - "execution_count": 38, + "execution_count": 36, "metadata": { "collapsed": false }, @@ -1704,11 +2227,27 @@ "source": [ "# Load the last statepoint file and keff value\n", "mgsp = openmc.StatePoint('statepoint.' + str(batches) + '.h5')\n", - "mgsu = openmc.Summary('summary.h5')\n", - "mgsp.link_with_summary(mgsu)\n", "mg_keff = mgsp.k_combined" ] }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Next, we can load the continuous-energy eigenvalue for comparison." + ] + }, + { + "cell_type": "code", + "execution_count": 37, + "metadata": { + "collapsed": true + }, + "outputs": [], + "source": [ + "ce_keff = sp.k_combined" + ] + }, { "cell_type": "markdown", "metadata": {}, @@ -1718,7 +2257,7 @@ }, { "cell_type": "code", - "execution_count": 39, + "execution_count": 38, "metadata": { "collapsed": false }, @@ -1727,9 +2266,9 @@ "name": "stdout", "output_type": "stream", "text": [ - "Continuous-Energy keff = 1.025440\n", - "Multi-Group keff = 1.025621\n", - "bias [pcm]: -18.1\n" + "Continuous-Energy keff = 1.025194\n", + "Multi-Group keff = 1.025183\n", + "bias [pcm]: 1.1\n" ] } ], @@ -1745,7 +2284,7 @@ "cell_type": "markdown", "metadata": {}, "source": [ - "We see quite good agreement with only an 18pcm difference between the two." + "We see quite good agreement with only an 1 pcm difference between the two. While these results are quite favorable, due to the high degree of approximations inherent in practical application of multi-group theory, one should not expect results of such fidelity always for multi-group Monte Carlo calculations." ] }, { @@ -1764,16 +2303,9 @@ "First, we extract volume-integrated fission rates from the Multi-Group calculation's mesh fission rate tally for each pin cell in the fuel assembly." ] }, - { - "cell_type": "markdown", - "metadata": {}, - "source": [ - "Now we can do the same for the Multi-Group results." - ] - }, { "cell_type": "code", - "execution_count": 40, + "execution_count": 39, "metadata": { "collapsed": false }, @@ -1790,6 +2322,32 @@ "mgopenmc_fission_rates /= np.mean(mgopenmc_fission_rates)" ] }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Now we can do the same for the Multi-Group results." + ] + }, + { + "cell_type": "code", + "execution_count": 40, + "metadata": { + "collapsed": true + }, + "outputs": [], + "source": [ + "# Get the OpenMC fission rate mesh tally data\n", + "mesh_tally = sp.get_tally(name='mesh tally')\n", + "openmc_fission_rates = mesh_tally.get_values(scores=['fission'])\n", + "\n", + "# Reshape array to 2D for plotting\n", + "openmc_fission_rates.shape = (17,17)\n", + "\n", + "# Normalize to the average pin power\n", + "openmc_fission_rates /= np.mean(openmc_fission_rates)" + ] + }, { "cell_type": "markdown", "metadata": {}, @@ -1807,7 +2365,7 @@ { "data": { "text/plain": [ - "" + "" ] }, "execution_count": 41, @@ -1816,9 +2374,9 @@ }, { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAADDCAYAAACS2+oqAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAHExJREFUeJzt3Xu8HGV9x/HP75B7SCBcEkwiQUwFaosQAbGixdvxUimX\nokKqoqGIr0KttaKiFBAvQKUYK6H1JUhBuaitEbxUoiIKRQUFvLRcDUmIxxwgCUm4BpJf/3hmzWSz\nu/Pbs7tndyff9+u1r7Nn59l5np35zW9nZ+aZx9wdERHpfwPdboCIiLSHErqISEkooYuIlIQSuohI\nSSihi4iUhBK6iEhJ9FxCN7PfmNkrut2O7ZmZfcfM3t7C+//NzD7azjZtj8xss5nt3WB6abcVxeAI\nuXvhA5gP3AZsAH4HfBt4WeS9BfO9DDin1fl085F9hqeB9dljA3BHt9sVaPdZwMZcm9cDH+h2u5po\n8xrgZuDQJt7/Q2DBKLRzGfAUsEvV63cCm4E9g/PZBOydi7OmthVgLHAmcHe2jh/Mtt3Xdntd1lif\nisE2PAr30M3s/cCFwCeA6cCewMXAXxa9dztyvrtPzR5T3P3AdldgZju0e57ANbk2T3X3CzpQR7td\n4+5Tgd2AG4Gvdbc5NTnwAHB85QUz+xNgQjYtylpsx38BRwBvA6YBzwM+C7yxZmWdibEiisF2Kvg2\nmUr65jymQZlxwELSnvtK4DPA2Gzan5P2Ct4PDGdl3plNO4n0TfcU6dvu2uz1B4BX5b4NvwJcnpX5\nNTAvV/dmsj2Y7P+t9mKyOu4DHgG+ATwne31O9t6BWt+cwPNJK+pR4CHg6gafv+6eU66edwDLs3l9\nJDfdgA8D9wMPA9cAO1e9d0H23huz199B2gN8GDijsryAGcDjwLTc/F+c1blDnT2NK4r2Ihoti2xd\nD2fT7gT+uJn1kFuHJwP3AquBiwr2jq7I/b8faS921+z/nYFvZu1cnT2fmU37BPAs8EQWS/+avb4v\nsCQrfxfw5tz83wj8b1b+QeD9wb2wB4CPALfmXvs0cHrW3j1r7a0BJwA3Vcc3gW2lRhtek8XDcwJt\n/SDwS+BJ0mHY/bK2rSVtc0fUio0Gbf474LfZevjn6PpUDLYeg0V76C8FxmcLoJ4zgEOA/YEXZc/P\nyE3fA5gCzAT+BlhkZju5+xeAK0krfKq7H1ln/kcAVwE7ZQtnUW5a3b0dM3sV8CngWOA5wApSwix8\nL/Bx4Hp33xmYDXyuQdmIlwF/RNrIzjSzfbLX/570S+flpOWzlvTrJ+8VpBX+OjPbj/T5jyd9pp2y\n9+Huw6SN4C259/41Kfg3tdD2msvCzAaBw4C52bS3kgJyK4H1APAXpC+fA4C3ZPNuyMzGkZLJatJy\ng5SMvgg8l/RL8gmyeHH3M4CbgFOzeHuvmU0ibUhfJu1tHQ9cnC1ngEuAkzztjf0JcENRu3J+Ckwx\ns33MbIC0Xr5M8V73NnHZxLaS92rgZ+7++0DZ44A3kJLRAHAd8F1gd+C9wJVm9kdNtPkoYF72ONLM\nFgTa0IhiMBiDRQl9V+ARd9/coMx84GPuvtrdVwMfA/InMzYCH3f3Te7+38BjwD415lPPze5+vaev\nqy+RvjgqGm0c84FL3f2X7v4Mae/opWa2Z6DOZ4A5ZjbL3Te6+y0F5U8zszVmtjb7e1lumgNnZ/P5\nFWlP6EXZtHcDH3X332dtPAc4NksAlfee5e5PuvvTpIC8zt1/4u7Pko6P5l1BtuyzeRxPWmb1vLWq\n3Xs0sSyeIX1R/7GZmbvfk32pVIush3PdfYO7P0j6UjqgqM2kDeVE4NhKfLr7Gndf7O5Pu/vjwLmk\nL8R63gQ84O5XeHIn6TDFsdn0jcALzWyKu6/LpjfjS6QN/rWk49hDTb6/FbsBqyr/mNm0bD0/amZP\nVpX9rLsPZTF2KDDZ3c9392fd/YfAt8gdPgo4L1teK0m/3hu9VzHYxhgsSuirgd1yCaaWmaRvvIrl\n2Wt/mEfVF8ITwI4F9eatyj1/AphQ0J58u5ZX/skW7mpgVuC9p5GWza1m9mszexeAmZ1uZhvMbL2Z\n5fekP+3uu7j7tOzvu6rmlw+y/OefAyzOAnkN8H+kIJ2RK7+y6jM9mPtMT7L1Hsm1wH5mthcwCDzq\n7j9v8Dm/UtXuVTXK1FwW2YZ+EWnvY5WZ/buZ1VqvkfVQb/nUbTPpfM5vgIMqE8xsopl93syWmdmj\nwI+Anc2s3hf/HODQyvI3s7Wkjb+y/P+KtOe23Mx+aGaHNmhXLV/O5vdO0pdtx2RxWYnN2aRl/JzK\ndHdf6+7TSHuh46reXjfGMsuJbTe15ledD6opBtsYg0WJ8Sek43ZHNSjzu6xR+QZG90SaOUFUyxPA\npNz/+W/3oXy7zGwy6RfHStKxReq9190fcvd3u/ss4D2kn0B7u/u5vuXkzd+22HZIX4RvyAK5EtST\nq34m55fR70k/OSufaWL2mSrtfhr4Kukk2NtovHceUm9ZZNMucveDgBeSfnWdVmMWjdZDK+1ak7Xn\nbDOrBP8/kg5tHZz9BK/sGVU2pup4e5B0biK//Ke6+6lZHb9w96NIhx6uJS3bZtq4gnSM+g3A12sU\neZz68bvN7ArqmpKLzZXAD4CDzaxWMq1OLvl5D5EOF+TtSdrOo23Ov39PWvxlohiMx2DDhO7u60kn\nARaZ2ZHZt88YM3uDmZ2XFbsGOMPMdjOz3YB/Ip5IhkknfZqRD8Y7gPlmNmBmryedhK24CniXme1v\nZuNJx9B+6u4PuvsjpAB9W/beBaQTL6kCs2PNrPLt/SjppMlIj0M3Oiz0eeBTlZ9+Zra7meWvHqp+\n738CR5jZoWY2lnR4q9qXSHuER5D2EFtSb1mY2UFmdoiZjSGdTHuK2suo7npotW3ufg/pWO+Hspem\nZG1Zb2a7AGdXvaU63r4FvMDM3pbF9djsc+2bPZ9vZlM9nYPYQDqh1awFpBOX1Yc5IJ3EOybbruaS\nfr7X09S24u7fIx06+Ea2nsZm6+qlNP5y+BnwuJl9MFsmh5MOC1zdRJtPM7Odzey5pPNE1cerm6IY\njMdg4aELd/8M6SqVM0hnblcAf8uWE6WfAH4OVI4P/xz4ZKNZ5p5fSjo+tMbMvl5jetH730c6qbiW\ndJxuca7dN5C+XL5OSt7PI538qTiJdHb/EdKZ6v/JTTsY+JmZrc8+53vdfTn1fTD7qbs++9n7UJ32\nVv//WdK37hIzWwfcQjqpXPO97v5/pCsIvkLa61hHWidP58rcQgr427M9xJHI11tvWUwFvkC6FvcB\n0nLc5pKzwHpotHwiLgBOynYmFpL2Hh8hLcvvVJX9LPBmM1ttZgvd/THSoanjSMtzCDiPLYck3g48\nkP10fjfpJHPEHz6Duz/g7rfXmka6QuMZ0mHFy9j2C7jVbeUYUsL4MmkbWUraTl5Xpw6yY8x/Sbq6\n4hHSIY23u/t9wTZDiulfALeTLmT4YkE7a1EMJk3FoLm3etRDuiX76fgo6Sz/8tzrPwCudPeRbEgi\nI2Zmm0nxuLTbbdke9VzXf2nMzN6U/dydDPwL8KuqZH4wcCBpL15EtiNK6P3nSNLPspWk4/5/+Olo\nZv9Buqb177Mz+SKjTT/5u0iHXERESkJ76CIiJTGmEzPNLiFcSPrCuNTdz69RRj8NpKPcvdWbW21D\nsS29oF5st/2Qi6VenPeS7iUxRLrt7nHufndVOd88e8v/Z6+Ds3eqanTgKPCzgavD1z9WXGZDo5sb\n1LGQdN1ksyYGymwomD4hMI9aFz5/jnTdY0XkwupIe8cGyszavbjMhvXbvnbus3B6k7seOz3d/oTe\nTGw/lut688mN8NFc38zlTxTXNTXQnsgiCVTFQzVeu4R046WKSJxMD5R5JlDmqUCZWvG/CDgl9//a\nGmWqzSkuEvrsk4qLMKbqXpbnb4YPVR0jmViwYQ+8epCJ1y2pG9udOORyCHCfuy/Prmm9hnQiT6Tf\nKbalp3Uioc9i63tBrKS5+0CI9CrFtvS0ThxDr/VToOZxnbPXbXm+c9uPdnZes3dq6gWHFBfpOYcF\ndjtu2gw3j+CwWZPCsf3JjVue79SHsT2v2w0YgYO73YAmvSwYFz/eBDdlh5bt7vsblu1EQl9JuiFP\nxWzq3Jyn+ph5v+nHhP6SbjdgBF4eSOgvH9i63Hmt3AG+vnBs54+Z96N+TOj9trNyWDChv2KH9AAY\n2Hcun7y3fifcThxyuQ2Ya2ZzLN0A/jjSDfNF+p1iW3pa2/fQ3X2TmZ1K6rFYubTrrnbXIzLaFNvS\n6zpyHbq7f5fmRiUS6QuKbellHUnoUWsKRjucMrl4HnfXuG65WqAIUwJliq4Nh8ZDs1RsM+hhDUX3\nvY20N3L9bGTZ7FdcJMQDDXry6eIy0wPXs/NwoEwHrW1wAfj0wIHOZwMneGtdP14tsn4j14avCZSJ\nXIf+20CZSHtmFxcJxX9km46I1DU1cF6n+lr1akVjtanrv4hISSihi4iUhBK6iEhJKKGLiJSEErqI\nSEkooYuIlIQSuohISSihi4iURFc7FhWNrbEm0CsiMjxH5EP+KlDmBM4sLPN5ziksE2nzewrquiJQ\nT+Rzzw98pm8E6op0PPnTwIgDuwTmsyZSWZf9rtHEQKehWoOTVIt0Gop0YjsxEAPfCsTAjYG6iuIa\nYttQJLaPDtT1tUBdkXWxV6BMKGwLVuq4ghFLtIcuIlISSugiIiWhhC4iUhJK6CIiJaGELiJSEkro\nIiIloYQuIlIS5kUXg3eqYjO/t6DMxoLpAIsDZc4MXI+6MHA9auTG+5MCZZ4fKBO5xrhI5PrZyLWx\nkZv3nxZYxlcHlvGrA7sYFhhcd7dN4O7BYXjby8z8zgbTI+vlvkCZdvWLmBCoKzLmdWQ+7Yq3fQNl\nImMDjg2UiVynH1nOkUGsiz771MFBXrBkSd3Y1h66iEhJKKGLiJSEErqISEkooYuIlIQSuohISSih\ni4iUhBK6iEhJKKGLiJREVzsW/bKgzA6B+dwaKNNwsIHMvECZhwJlpgfKRAYdWFEw/cDAPIYDZQJj\nTjAlUGZWoEykw8jLAit9TKDMzhu727FoSYPpkeUQ6VQXievIeonEwMxAmcjnisR+ZKVFYjKyvUZE\nOhTOCJSZHShTNHCHOhaJiGwnlNBFREpCCV1EpCSU0EVESkIJXUSkJJTQRURKQgldRKQklNBFREqi\n6Dr2ETGzZcA6YDPwjLvXHKyjqONLpMtTZDSRMwOjiUQ6PCxo08glkY4K/1RQV+QzRTo5fSDwmS4K\n1BUZYecfAnXdtqm4rombApV1SDS2G43esyFQz9FtGmUrMvJVJK4XBeqaFqgrMsrSuYG6IqN+nRyo\n64I2fa43B+q6IVBXUUIe3+L7R2ozcLi7RzqhifQTxbb0rE4dcrEOzlukmxTb0rM6FZgOXG9mt5nZ\nSR2qQ6QbFNvSszp1yOXP3H2Vme0OfM/M7nL3mztUl8hoUmxLz+pIQnf3Vdnfh81sMXAIsE3QX5F7\n/qLsITIStwK3jUI90dj+Yu75gcTujilSyx3ZA2DC/fc3LNv2hG5mk4ABd3/MzCYDg8DHapV9R7sr\nl+3WIdmj4uIO1NFMbC/oQP2yfcrvEOw8dy7/tnRp3bKd2EOfASw2M8/mf6V7w9tDi/QLxbb0tLYn\ndHd/ADig3fMV6TbFtvS6ro5Y9JuCMpERUG4PlIl8a00MlGnXCDFjA2WWFUyPtDdST8QugTJPBsqs\nCZSJjOoSWcYH0d0Ri+5sMD3SsSwyEtd+wfYUWRYos1eb5hOJgUi8RUYsiuSPyKhGkdGIGnUkq9gj\nUGb3gukTBweZqRGLRETKTwldRKQklNBFREpCCV1EpCSU0EVESkIJXUSkJJTQRURKQgldRKQkOnW3\nxZCiji+RjjF7BcpEOnJEOiFEOrRERj2IdAoaVzA9smwiK/epQJmpgTKRDiORZRzpoDS3aOEAbAyU\n6aBGcRCJx70CZSId3SLrblKgTGTko4hIT69Ip6FI/K8IlNkzUCayjUTiNtJprqiuHQqmaw9dRKQk\nlNBFREpCCV1EpCSU0EVESkIJXUSkJJTQRURKQgldRKQklNBFREqiqx2LitwVKHM0ZxaWuZpzCssU\nXbAP8JZAXYsCdUU6lpxSUNfCNtVzWps+U6Qfzz8E6vp2oK4NXe40FNFoHLAnAu+PxPVVgWUVqWt+\noK5vBeqKdJg7MVDXBYG6ikb2AfhAoK7LA3VFOnB9JFDXDYG6ijpMbSqYrj10EZGSUEIXESkJJXQR\nkZJQQhcRKQkldBGRklBCFxEpCSV0EZGSUEIXESkJc2/UBaKDFZv5DwvKTAvMZ1mgzEOBMnMCZYYD\nZdo1IktkBKAiGwJlIj3L5gbKRDpfRD7TawNl9g0MwzN+Pbh7ZHW0nZn5j1ucR2Td7RIoExlJJxLX\nMwJlIus3Em+R0ZHaNYpWZBlGREZQiowMNbNoHoODPHfJkrqxrT10EZGSUEIXESkJJXQRkZJQQhcR\nKQkldBGRklBCFxEpCSV0EZGSUEIXESmJEY9YZGaXAm8Cht19/+y1acBXSP10lgFvcfd19eZRdBF9\npINBROSi/0gnhLWBMpFOOBFFPWIio8PMCpT5baBMZPlF6no2UCYyytKGxwOFWtCO2G60zCIb3fJA\nmUhHlYin2lQm8rkinfymB8pERNoc6cAV6QgYmU9kOyrKMZsLpreyh34Z8Lqq1z4MfN/d9wFuAE5v\nYf4i3aLYlr404oTu7jez7RfKkcDl2fPLgaNGOn+RblFsS79q9zH06e4+DODuq4iN5SrSDxTb0vN0\nUlREpCRGfFK0jmEzm+Huw2a2BwXnQD6Xe34I8JI2N0a2Hzc7/E9nbxzaVGxfkns+L3uIjMQvgNuz\n5xPuv79h2VYTurH1BRnXAe8EzgdOAK5t9Oa/a7FykYrDLD0qPr2p5Vm2FNt/03L1IsmLswfATnPn\nsmjp0rplR3zIxcyuAm4BXmBmK8zsXcB5wGvN7B7gNdn/In1FsS39asR76O4+v86k14x0niK9QLEt\n/ardx9DbWnnkYv2jObOwzELOKSwTOfz6vkBdXwzUtTpQ12kFdX01UE+ks9Qpgc90QaCu5wfqOjlQ\n108CdT3b+uGUjmvUOSYyutM7AstqUZvi+tRAXZ8P1BXpCFgU1xDbhnYM1BWJ7Uhdke3oxEBdVwfq\niozE1IiuchERKQkldBGRklBCFxEpCSV0EZGSUEIXESkJJXQRkZJQQhcRKQkldBGRkjD3zt7RqG7F\nZn5PQZnIBf3LAmUio6RERiUZDpTZL1AmMipPkV0DZSIdPSJlIp0dIqM5RTqDRG7QFhn5aG/A3YsG\nfuoIM/OfNpg+LTCPnwXKPBYos2+gTCSuIzHwRKBMpLPgjECZiEgHvkgsRUaPmhMoExnNrGg57zg4\nyNwlS+rGtvbQRURKQgldRKQklNBFREpCCV1EpCSU0EVESkIJXUSkJJTQRURKQgldRKQkujpi0axJ\njadPDfRUiHSKiIw0dHlgNJFZgboiHWzasdAjo95EOktFBv+JdD5aEFjGNwSWcaTjyS6BMt02s8G0\nuwPvnx0o8+rAMr8isMwjcfJkoEzB5gzEto9Ih8LItjg2UCayHZ0TWM7fCCznSNxObHG69tBFREpC\nCV1EpCSU0EVESkIJXUSkJJTQRURKQgldRKQklNBFREpCCV1EpCS6OmLRuvGNy0SatnpjcZnIiCOR\nUYReH+hgsCjQwWCHQF3vKahrcaCeSKeronqidRV1eIDYiC27FsQEwMRAmfHruzti0YoG0yMdZyLu\nC5SJjMgzv00xEOmAdkKgrq8G6op0Gjo6UNelbYrtPw2UicynqKPfpMFBZmvEIhGR8lNCFxEpCSV0\nEZGSUEIXESkJJXQRkZJQQhcRKQkldBGRklBCFxEpiRF3LDKzS4E3AcPuvn/22lnAScBDWbGPuPt3\n67zf/7egjr0mF7djTKSXTsBwoFfE7YH5tGs0naLOJ5F62jX6z/RAmciINjMCuw8W6Aq0S6BBA78f\neceidsT20gbzj4xYFenoFhFZLz8KlNkrUCYwwFioo1O7RiOKdKyLdPaZESgTCbTIyFBFsTFhcJDp\nHepYdBnwuhqvX+ju87JHzYAX6XGKbelLI07o7n4ztYcI7Ep3a5F2UWxLv+rEMfRTzOxOM7vEzHbq\nwPxFukWxLT2tHQPQ510MnOPubmafAC4ETqxXeFHu+cHAIW1ujGw/bnwabgzcqK0FTcX2wtzzQ7OH\nyEj8JHsAjLn//oZl25rQ3f3h3L9fAL7ZqPwp7axctmuHj0+PinMea+/8m43t97W3etmOvTR7AEyY\nO5cLltY/5d7qIRcjd1zRzPbITTsG+E2L8xfpFsW29J0R76Gb2VXA4cCuZrYCOAt4pZkdAGwGlgEn\nt6GNIqNKsS39asQJ3d3n13j5shbaItITFNvSr9p9UrQp++3eePr6dcXz2PB4e9oyMXDw6SWbi8tE\nOldEOkXsFShTZOq44jJrAicSp08qLvNY4IPvOrO4zDOBDl5WNKxLD2jUYWViYL20y3Bg/UZG24l0\nCIp0nIl05Kl1vWi1yDYUucY0ENohkVGoIsun6HMV9aNU138RkZJQQhcRKQkldBGRklBCFxEpiZ5J\n6B3u5dcRt3S7ASNwU+DEbq/5UeSMXA/rxzi5o9sNGIHI3VB7yU87MM/eSejtul/oKOrHDfXmPkzo\nP1ZCH3VK6J1X6oQuIiKt6ep16Ow/b8vzpUOw99YXKg8E7scxpk179gORr7aqvduBoSHGzNy6zZFL\njEdroQ/UuKjVVg4xMHtLm8cFlt/AhOIyY54KNCgwUoDVWufLh7A5ueU8vkaZat/v7v7amHlbYrs6\nTmqtl04ZH1i/O9Z4bdzQEDvm2hz5kRSJ68iqi9RVa+ybsUNDTM61OXI7zMAYOqFr5yNpqDo37DA0\nxLiq/FG0DMfMnQtLltSdPuIRi1plZt2pWLYbIx2xqFWKbem0erHdtYQuIiLtpWPoIiIloYQuIlIS\nPZHQzez1Zna3md1rZh/qdnsizGyZmf3SzO4ws1u73Z5azOxSMxs2s1/lXptmZkvM7B4zu76XhlKr\n096zzGylmd2ePV7fzTY2q99iW3HdGaMV211P6GY2AFxEGmX9hcDxZrZvd1sVshk43N0PdPdeHT2v\n1uj1Hwa+7+77ADcAp496q+qr1V6AC919Xvb47mg3aqT6NLYV150xKrHd9YROGkr0Pndf7u7PANcA\nR3a5TRFGbyy/uuqMXn8kcHn2/HLgqFFtVAN12guxO6H2on6MbcV1B4xWbPfCipsFPJj7f2X2Wq9z\n4Hozu83MTup2Y5ow3d2HAdx9FVBwV/qecIqZ3Wlml/TaT+kC/RjbiuvR1dbY7oWEXusbqh+upfwz\ndz8IeCNppRzW7QaV1MXA8939AGAVcGGX29OMfoxtxfXoaXts90JCXwnsmft/NjDUpbaEZXsBldHg\nF5N+XveDYTObAX8Y+PihLrenIXd/2Ld0lvgCcHA329OkvottxfXo6URs90JCvw2Ya2ZzzGwccBxw\nXZfb1JCZTTKzHbPnk4FBencU+K1Gryct23dmz08Arh3tBhXYqr3ZxllxDL27nGvpq9hWXHdcx2O7\nu/dyAdx9k5mdCiwhfcFc6u53dblZRWYAi7Mu3mOAK929/g0WuqTO6PXnAV8zswXACuDN3Wvh1uq0\n95VmdgDp6otlwMlda2CT+jC2FdcdMlqxra7/IiIl0QuHXEREpA2U0EVESkIJXUSkJJTQRURKQgld\nRKQklNBFREpCCV1EpCSU0EVESuL/ASY96jLsHTMbAAAAAElFTkSuQmCC\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAADDCAYAAACS2+oqAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAG9ZJREFUeJzt3Xm4HFWZx/HvG0iAkATCEpYAUYkCLohIAuig6OhFHRGG\nQQXGhUXUEcRlRAUxIKDAyEBQ0GEkIgiC44yI20hUdB6QxSDizhJDAuGaACEQkCUxeeePUw2Vprvr\n7b7Vt7srv8/z9HP7dp2uc7rqrberq+rUMXdHREQG35heN0BERMqhhC4iUhFK6CIiFaGELiJSEUro\nIiIVoYQuIlIRfZfQzez3ZvaqXrdjXWZmPzSzd47g/V82s0+V2aZ1kZmtMbPntZhe2W1FMdghdy98\nAIcB84BHgfuAHwCvjLy3YL4XA6eOdD69fGSf4SlgRfZ4FPh1r9sVaPfJwMpcm1cAH+t1u9po80PA\n9cBebbz/Z8CRo9DOhcCTwGZ1r98GrAF2CM5nNfC8XJy1ta0AY4FZwO3ZOr4323Zf3+t12WB9KgZL\neBTuoZvZR4FzgNOBKcAOwJeAtxS9dx1ylrtPyh4T3f1lZVdgZuuVPU/gylybJ7n72V2oo2xXuvsk\nYAvg58C3etuchhy4Gzi09oKZvRjYMJsWZSNsx/8A+wPvACYDzwXOA97UsLLuxFgRxWCZCr5NJpG+\nOQ9qUWYcMJu0574YOBcYm017NWmv4KPA0qzM4dm0o0nfdE+Svu2uzl6/G3ht7tvwm8AlWZnfAbvn\n6l5DtgeT/b/WXkxWx13Ag8B3gG2y16dl7x3T6JsT2JG0oh4G7geuaPH5m+455ep5F7Aom9eJuekG\nfBKYDzwAXAlsWvfeI7P3/jx7/V2kPcAHgJNqywvYCvgrMDk3/5dnda7XZE/j0qK9iFbLIlvXS7Np\ntwEvbGc95Nbh+4A7gWXA+QV7R5fm/t+FtBe7efb/psD3snYuy55vm007Hfgb8HgWS1/IXt8ZmJuV\n/xPw1tz83wT8ISt/L/DR4F7Y3cCJwC9zr30eOCFr7w6N9taAdwPX1cc3gW2lQRtel8XDNoG2fhz4\nDfAE6TDsLlnblpO2uf0bxUaLNn8Q+HO2Hv4tuj4VgyOPwaI99L2BDbIF0MxJwExgV+Cl2fOTctO3\nBiYC2wLvAS4ws03c/SvA5aQVPsndD2gy//2BbwCbZAvngty0pns7ZvZa4HPAwcA2wD2khFn4XuA0\n4Bp33xTYDvhii7IRrwSeT9rIZpnZTtnrHyL90tmHtHyWk3795L2KtML3M7NdSJ//UNJn2iR7H+6+\nlLQRvC333n8mBf/qEbS94bIwsyHg74Dp2bS3kwJyLYH1APAPpC+f3YC3ZfNuyczGkZLJMtJyg5SM\nvgpsT/ol+ThZvLj7ScB1wLFZvB1nZuNJG9JlpL2tQ4EvZcsZ4CLgaE97Yy8Gri1qV85NwEQz28nM\nxpDWy2UU73U/Ky7b2Fby/h642d3/Eih7CPBGUjIaA3wX+BGwJXAccLmZPb+NNh8I7J49DjCzIwNt\naEUxGIzBooS+OfCgu69pUeYw4DPuvszdlwGfAfInM1YCp7n7anf/X+AxYKcG82nmene/xtPX1ddJ\nXxw1rTaOw4A57v4bd19F2jva28x2CNS5CphmZlPdfaW731BQ/ngze8jMlmd/L85Nc+CUbD6/Je0J\nvTSb9l7gU+7+l6yNpwIHZwmg9t6T3f0Jd3+KFJDfdfcb3f1vpOOjeZeSLftsHoeSllkzb69r99Zt\nLItVpC/qF5qZufsd2ZdKvch6OMPdH3X3e0lfSrsVtZm0oRwFHFyLT3d/yN2vcven3P2vwBmkL8Rm\n3gzc7e6XenIb6TDFwdn0lcCLzGyiuz+STW/H10kb/OtJx7GH23z/SGwBLKn9Y2aTs/X8sJk9UVf2\nPHcfzmJsL2Bjdz/L3f/m7j8Dvk/u8FHAmdnyWkz69d7qvYrBEmOwKKEvA7bIJZhGtiV949Usyl57\neh51XwiPAxMK6s1bknv+OLBhQXvy7VpU+ydbuMuAqYH3Hk9aNr80s9+Z2REAZnaCmT1qZivMLL8n\n/Xl338zdJ2d/j6ibXz7I8p9/GnBVFsgPAX8kBelWufKL6z7TvbnP9ARr75FcDexiZs8BhoCH3f2W\nFp/zm3XtXtKgTMNlkW3o55P2PpaY2X+YWaP1GlkPzZZP0zaTzuf8HtijNsHMNjKzC81soZk9DPwf\nsKmZNfvinwbsVVv+ZractPHXlv8/kfbcFpnZz8xsrxbtauSybH6Hk75suyaLy1psbkdaxtvUprv7\ncnefTNoLHVf39qYxlllEbLtpNL/6fFBPMVhiDBYlxhtJx+0ObFHmvqxR+QZG90TaOUHUyOPA+Nz/\n+W/34Xy7zGxj0i+OxaRjizR7r7vf7+7vdfepwPtJP4Ge5+5n+DMnbz4wwrZD+iJ8YxbItaDeuO5n\ncn4Z/YX0k7P2mTbKPlOt3U8B/0U6CfYOWu+dhzRbFtm08919D+BFpF9dxzeYRav1MJJ2PZS15xQz\nqwX/v5IObc3IfoLX9oxqG1N9vN1LOjeRX/6T3P3YrI5fufuBpEMPV5OWbTttvId0jPqNwLcbFPkr\nzeP3WbMrqGtiLjYXAz8FZphZo2Ran1zy8x4mHS7I24G0nUfbnH//Dozwl4liMB6DLRO6u68gnQS4\nwMwOyL591jezN5rZmVmxK4GTzGwLM9sC+DTxRLKUdNKnHflg/DVwmJmNMbM3kE7C1nwDOMLMdjWz\nDUjH0G5y93vd/UFSgL4je++RpBMvqQKzg82s9u39MOmkSafHoVsdFroQ+Fztp5+ZbWlm+auH6t/7\n38D+ZraXmY0lHd6q93XSHuH+pD3EEWm2LMxsDzObaWbrk06mPUnjZdR0PYy0be5+B+lY7yeylyZm\nbVlhZpsBp9S9pT7evg+8wMzekcX12Oxz7Zw9P8zMJnk6B/Eo6YRWu44knbisP8wB6STeQdl2NZ30\n872ZtrYVd/8x6dDBd7L1NDZbV3vT+svhZuCvZvbxbJnsSzoscEUbbT7ezDY1s+1J54nqj1e3RTEY\nj8HCQxfufi7pKpWTSGdu7wE+wDMnSk8HbgFqx4dvAT7bapa553NIx4ceMrNvN5he9P4Pk04qLicd\np7sq1+5rSV8u3yYl7+eSTv7UHE06u/8g6Uz1L3LTZgA3m9mK7HMe5+6LaO7j2U/dFdnP3vubtLf+\n//NI37pzzewR4AbSSeWG73X3P5KuIPgmaa/jEdI6eSpX5gZSwN+a7SF2Il9vs2UxCfgK6Vrcu0nL\n8VmXnAXWQ6vlE3E2cHS2MzGbtPf4IGlZ/rCu7HnAW81smZnNdvfHSIemDiEtz2HgTJ45JPFO4O7s\np/N7SSeZI57+DO5+t7vf2mga6QqNVaTDihfz7C/gkW4rB5ESxmWkbWQBaTvZr0kdZMeY30K6uuJB\n0iGNd7r7XcE2Q4rpXwG3ki5k+GpBOxtRDCZtxaC5j/Soh/RK9tPxYdJZ/kW5138KXO7unWxIIh0z\nszWkeFzQ67asi/qu67+0ZmZvzn7ubgz8O/DbumQ+A3gZaS9eRNYhSuiD5wDSz7LFpOP+T/90NLOv\nka5p/VB2Jl9ktOknfw/pkIuISEVoD11EpCLW78ZMs0sIZ5O+MOa4+1kNyuingXSVu4/05lbPotiW\nftAstks/5GKpF+edpHtJDJNuu3uIu99eV87XbPPM/6c8CqdMXHteqwJHgVcFrgx+4sniMita3dwg\ns7zu/wtJd/TJi3xDrgqUKbrofXKH9VwAHJP7v/4zNbJZoMykQJltA41+4qlnv/bZlfCpXN/G1YEe\nAZs8VX5Cbye285d4zCZdX9uOjQJlVpRUptEmdBHpxks1jwbmE4n9wKbIlECZRlnrP0nX9bVT19hA\nmYgtA2Xqw/YLpBvltGP80BDbzZ3bNLa7cchlJnCXuy/Krmm9knQiT2TQKbalr3UjoU9l7XtBLKa9\n+0CI9CvFtvS1bhxDb/RToOFxnVNyv+M2Lf1oZ/e9vNcN6MCMXjegA/sEhl24bg1cHzhsNkLh2J6d\nex45HNVvdu91AzowaNvjnsFyN2cPgLHz57cs242Evph0Q56a7Whyc576Y+aDZo/iIn1nZnGRvvOq\nQELfZ0x61Jw5kjvANxeO7XaPmfcbJfTuiyb0PXNlx0+fzrkLmnfC7cYhl3nAdDObZukG8IeQbpgv\nMugU29LXSt9Dd/fVZnYsqcdi7dKuP5Vdj8hoU2xLv+vKdeju/iPaG5VIZCAotqWfdSWhR/3t8ZHP\nYzgwj4cC84lcG35fcRGeEygT+dhFbW50c+16kZUbub44ck105Hre5ZGLogMmR84yNriefTS1iqfI\nTdUj6yWyOCMxcH9xkVC8RU6JtboHdU2kzZ3cmL6RyHYfCbfQ9fWBA9xPFJzYrx9qqp66/ouIVIQS\nuohIRSihi4hUhBK6iEhFKKGLiFSEErqISEUooYuIVIQSuohIRfS0Y9HjBVfjrwh0Dol09nkkUGZh\noMxHmFVY5kJOLSwT6VhUVNdXA/VEhi45qqTPFOlY9JLADbO2CsxnbAkd0rqtVWhH1kukI0+kY9HS\nQJmyYmBloK5jAnWNZmxfGqgr0vkoUmajwN1AizrxFd2nTnvoIiIVoYQuIlIRSugiIhWhhC4iUhFK\n6CIiFaGELiJSEUroIiIVYe6RKzq7ULGZ/6GgTOQ625uLi4SuH58VuB41ctH+5oEy2wbKlDEWRORa\n5shACpH5zCrp+uJXBuraLDBo9Barwd0tMLvSmZnPazE9sjwXBsq8O7DMZ5d0nfVmgTKR7SMy2MzW\ngTIbBspErsGPDJRxbEmxvUugrqLlPH5oiO3mzm0a29pDFxGpCCV0EZGKUEIXEakIJXQRkYpQQhcR\nqQgldBGRilBCFxGpCCV0EZGK6GnHojsKykQ6IfwiUKasjhMF43EAsU5DkY4lfy6YvnNgHpHltzxQ\nZstAmUgHjUinq0jni6JBAABeQG87Fv2sxfTIQBC/D5SJDNwyNVDmsUCZiSWViXRkmxwoE4m34UCZ\nSAelSP4YFyjznECZKQXTJw0N8QJ1LBIRqT4ldBGRilBCFxGpCCV0EZGKUEIXEakIJXQRkYpQQhcR\nqQgldBGRiogMMtI2M1tI6vewBljl7jMblSvqYLMwUNfxJY1GFOnsExn56OxAXZGuXKcW1HVaoJ4J\ngXo+HfhMZ5Q06k1k5JcfBOqKdALrlmhst1rHkdGoPhZYVueWFNeRui4oKQbeV9LoP5FtqKzPNTZQ\n1zElxXZRQi7aA+9KQicF+77uHumIKDJIFNvSt7p1yMW6OG+RXlJsS9/qVmA6cI2ZzTOzo7tUh0gv\nKLalb3XrkMsr3H2JmW0J/NjM/uTu13epLpHRpNiWvtWVhO7uS7K/D5jZVcBM4FlB/+Xc8z2AGd1o\njKwTfgXcOgr1RGP7a7nnu2UPkU7MA27Jno+bP79l2dITupmNB8a4+2NmtjEwBHymUdl/KbtyWWe9\nPHvUzOlCHe3E9uFdqF/WTTN4Zmd3wvTpfHHBgqZlu7GHvhVwlZl5Nv/L3X1uF+oRGW2KbelrpSd0\nd78b/cKUClJsS7/r6YhFRaPyREbc+V2gTKTDQ2QUnEWBMtMCZSKfq2hElqKRTSA2OkxZo7FERpmJ\ndKiJjGoUKfMaejti0dUtpj8emEekQ1BkBK3IelkYKLNJoExkZKtI7Ee2xcjyicR/ZESnyDYSsUWg\nTNEe9uZDQ7xcIxaJiFSfErqISEUooYuIVIQSuohIRSihi4hUhBK6iEhFKKGLiFSEErqISEV0626L\nIUUjBBR1rgEYHyizYaBMpMNDpK5IZ49JgTJFy2Z1YB6RkX3uD5SJdE6JdOJ4LFBm60CZSAeWXmu1\nYZWx/ovqqCmrg1Kkc1lkPpGeXpERgiK5IVJXZD6R9kTqKqODUtF2rz10EZGKUEIXEakIJXQRkYpQ\nQhcRqQgldBGRilBCFxGpCCV0EZGKUEIXEamInnYsKuqwEhmN6FBmFZa5gFMLy6wM1PWxQF1nB+qK\njBF1fEFdkc/0SKCej5T0mSKdtyLL74pAXZGOHv0sMnLTMYFlNTuwrCKjI50aqOuMkup6f0mfK5K4\nIvEWqSuyvUa2owsDdc0smD6hYLr20EVEKkIJXUSkIpTQRUQqQgldRKQilNBFRCpCCV1EpCKU0EVE\nKkIJXUSkIsw9ctl8Fyo285sKyiwLzGdxpK5AmUhHhUiHkMioRpERgG4vmD41MI/IaDWRUVQiI+xE\nREaFenWgTGQkph0Bd4+s+tKZmf+gxfRIB5zIsop0sIrE9X2BMpEYiNQVGSFoSqBM5LPPD5SJbK+R\nEbIiI3ZFttlNiqYPDbHL3LlNY1t76CIiFaGELiJSEUroIiIVoYQuIlIRSugiIhWhhC4iUhFK6CIi\nFaGELiJSER2PWGRmc4A3A0vdfdfstcnAN4FpwELgbe7edOCcMro0RToGREbTiXTkeDJQZlygTKQT\nQtHnivSYiXTAWRQos1WgTESkc0pkGa8eaUMKlBHbrTpsRTp8LQ+UiSzPohFuILY8NwqUiYjUFYmB\nyDKMdD6KfK5IZ6hIXWXkqqJ6RrKHfjGwX91rnwR+4u47AdcCJ4xg/iK9otiWgdRxQnf363n2jsQB\nwCXZ80uAAzudv0ivKLZlUJV9DH2Kuy8FcPclxG6DIDIIFNvS93RSVESkIjo+KdrEUjPbyt2XmtnW\nwP2tCl+Ue7579hDpxM3Zo4vaiu0rcs9fDLykq02TKpsH3JI9Hze/9T0kR5rQjbUvuPgucDhwFvBu\n4OpWb37PCCsXqdkze9ScP/JZjii2Dx15/SIAzMgeABOmT+eLCxY0LdvxIRcz+wZwA/ACM7vHzI4A\nzgReb2Z3AK/L/hcZKIptGVQd76G7+2FNJr2u03mK9APFtgyqso+ht6VotJylgXkcxazCMqdxamGZ\nyChCJwTqOiNQV2QkplkFdUXqiXTcOjHwmWYF6poYqOv4QF3fCtQ1LVBXr7XqsBIZISgSa5H10rTn\nU86nS4rrSAe+jwTqOjtQ13ol1TU7UFdkhKljA3VdEajrxQXTiz63rnIREakIJXQRkYpQQhcRqQgl\ndBGRilBCFxGpCCV0EZGKUEIXEakIJXQRkYow9zLGDeqgYjOfV1AmMorQwkCZxwJlIp0iVgTKREb3\niYy2UlRXpJ5Ir7GWd5jKRDoNRUbYicxnZqBMZFSolwLuHhnYqXRFsb0wMI9I56PIMo90mIvEY2T0\nq4jI6D+RUYQiZSIdEyMBEqkrkj+2D5TZtmD6+KEhtp87t2lsaw9dRKQilNBFRCpCCV1EpCKU0EVE\nKkIJXUSkIpTQRUQqQgldRKQilNBFRCqipyMW7TK+9fSFgaFCFgbq+VhgNJHzA6OJFDQ3rIzOFU8G\n5jGhpLY8GigTGWEnMmLL6kBdUwJleq3Vct088P5IZ5/IiDyRkYYinWsmBcpE4mRsoEykPZH4j7Q5\n0nnxwyWNtBXpoFT0uYoStvbQRUQqQgldRKQilNBFRCpCCV1EpCKU0EVEKkIJXUSkIpTQRUQqQgld\nRKQiejpi0b0FZSK9nhYHyiwMlIl0ePjHQAeD2YEOBqsCdR1fUNecQD2REZYinVMuCdQV6cQxNVAm\nMjJOpOPJ8+jtiEW/bzE90pkl0rEoMiJPpLPPsYEYODsQAysDdZ0YqOvCQF2R3HBUSdtrZNSn5wbK\nTAuUKVrvGrFIRGQdoYQuIlIRSugiIhWhhC4iUhFK6CIiFaGELiJSEUroIiIVoYQuIlIRHXcsMrM5\nwJuBpe6+a/baycDRwP1ZsRPd/UdN3u93FtQxuaSvmxVrisssCswn0pEjMipJRFEHg7LqiXSoinQI\ninSWinQaCo1YNK64zKYrO+9YVEZs39Fi/pFOQ5HOR2WJdM6LxFtk9KvIZ4+MDDYxUCayvUZiO/LZ\nI22OxH9Rh6kNh4aY0qWORRcD+zV4/Rx33z17NAx4kT6n2JaB1HFCd/frgeUNJvWku7VIWRTbMqi6\ncQz9GDO7zcwuMrNNujB/kV5RbEtfi9zjph1fAk51dzez04FzgKOaFf5C7vme2UOkE9etgesD50pG\noK3Y/mLu+UwU29K5G7MHwPrz57csW2pCd/cHcv9+Bfheq/LHlVm5rNP2GZMeNWdFbv3XhnZj+4Pl\nVi/rsL2zB8CG06dz9oIFTcuO9JCLkTuuaGZb56YdBLS6i6hIP1Nsy8DpeA/dzL4B7Atsbmb3ACcD\nrzGz3YA1pNuQv6+ENoqMKsW2DKqOE7q7H9bg5YtH0BaRvqDYlkFV9knRthRdsD820LqHAsdKIx0e\npgTKRDoYRDo8RDpXFNUV+UyR9t5fXKS0jhWbBzoErQr0LLIBuHiwVRPLiqNIDDwZKPP8QJlG13DW\nmxAoE2lzpJNapEPccwJlIusiUiYyYtdmgfgvsn5BTlTXfxGRilBCFxGpCCV0EZGKUEIXEamIvkno\nNxYX6Tu39roBHbi51w3owC86uyFo3xjEZT6IsX1brxvQpuu60LNZCX0EBjHof9nrBnRACX30DWJs\nD1pC78atKvomoYuIyMj09Dr0sbvv/vTzMcPDjN1227WmjwlcbDoucNFq5FtrvUCZ+i/UscPDbFzX\n5g0D84lcRr1BwfTIQBCN5rHe8DAb5Nq8cQltgdhniqzP9Rrstdh9w6w39Zk2j4msrJt6u4+5QS62\n169b5pEfHJHr+iMxENnAG63f+tiO1BVpc+Q69EiZRp9r3PAwE+q2xyKR7TUS/5FLzOvj3xYPM2a7\n9to7ZsfpwNym0zsesWikzGzAf0hLv+t0xKKRUmxLtzWL7Z4ldBERKZeOoYuIVIQSuohIRfRFQjez\nN5jZ7WZ2p5l9otftiTCzhWb2GzP7tZn15dWAZjbHzJaa2W9zr002s7lmdoeZXdNPQ6k1ae/JZrbY\nzG7NHm/oZRvbNWixrbjujtGK7Z4ndDMbA5xPGmX9RcChZrZzb1sVsgbY191f5u4ze92YJhqNXv9J\n4CfuvhNwLXDCqLequUbtBTjH3XfPHj8a7UZ1akBjW3HdHaMS2z1P6KQhF+9y90Xuvgq4Ejigx22K\nMPpj+TXVZPT6A4BLsueXAAeOaqNaaNJeiF0V2Y8GMbYV110wWrHdDytuKnBv7v/F2Wv9zoFrzGye\nmR3d68a0YYq7LwVw9yXAlj1uT8QxZnabmV3Ubz+lCwxibCuuR1epsd0PCb3RN9QgXEv5CnffA3gT\naaX8Xa8bVFFfAnZ0992AJcA5PW5POwYxthXXo6f02O6HhL4Y2CH3/3bAcI/aEpbtBdRGg7+K9PN6\nECw1s63g6YGPI4MW9Yy7P+DPdJb4CjCjl+1p08DFtuJ69HQjtvshoc8DppvZNDMbBxwCfLfHbWrJ\nzMab2YTs+cbAEP07Cvxao9eTlu3h2fN3A1ePdoMKrNXebOOsOYj+Xc6NDFRsK667ruux3dN7uQC4\n+2ozO5Z0g4IxwBx3/1OPm1VkK+CqrIv3+sDl7t78Bgs90mT0+jOBb5nZkcA9wFt718K1NWnva8xs\nN9LVFwuB9/WsgW0awNhWXHfJaMW2uv6LiFREPxxyERGREiihi4hUhBK6iEhFKKGLiFSEErqISEUo\noYuIVIQSuohIRSihi4hUxP8DxMTTRn5AfRMAAAAASUVORK5CYII=\n", "text/plain": [ - "" + "" ] }, "metadata": {}, @@ -1845,6 +2403,15 @@ "source": [ "We also see very good agreement between the fission rate distributions." ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], + "source": [] } ], "metadata": { diff --git a/openmc/mgxs/library.py b/openmc/mgxs/library.py index 82b761673..334a0fa9c 100644 --- a/openmc/mgxs/library.py +++ b/openmc/mgxs/library.py @@ -724,7 +724,7 @@ class Library(object): def write_mg_library(self, xs_type='macro', domain_names=None, xs_ids=None, filename='mg_cross_sections', directory='./', - return_names=True): + return_names=False): """Creates a cross-section data library file for the Multi-Group mode of OpenMC. @@ -749,7 +749,7 @@ class Library(object): return_names : bool Flag to indicate if the user would like the names of the materials generated by this function returned with completion. - Defaults to True. + Defaults to False, indicating that no names will be returned. Returns -------