{ "metadata": { "name": "" }, "nbformat": 3, "nbformat_minor": 0, "worksheets": [ { "cells": [ { "cell_type": "code", "collapsed": false, "input": [ "from pyne.xs.channels import sigma_t\n", "from pyne.material import Material, from_atom_frac\n", "from pyne.mesh import Mesh, IMeshTag, MetadataTag, ComputedTag\n", "from yt.config import ytcfg; ytcfg[\"yt\",\"suppressStreamLogging\"] = \"True\"\n", "from yt.visualization.plot_window import SlicePlot\n", "from yt.frontends.moab.api import PyneMoabHex8StaticOutput" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 1 }, { "cell_type": "code", "collapsed": false, "input": [ "def cantor(n):\n", " return [0.] + cant(0., 1., n) + [1.]\n", "\n", "def cant(x, y, n):\n", " if n == 0:\n", " return []\n", " new_pts = [2.*x/3. + y/3., x/3. + 2.*y/3.]\n", " return cant(x, new_pts[0], n-1) + new_pts + cant(new_pts[1], y, n-1)" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 2 }, { "cell_type": "code", "collapsed": false, "input": [ "#c5 = cantor(5)\n", "c5x = [0.] + cant(0., 1500., 5) + [1500.]\n", "c5y = [-850.] + cant(-850., 850., 5) + [850.]\n", "coords = [c5x, c5y, [-850.0, 850.0]]" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 22 }, { "cell_type": "code", "collapsed": false, "input": [ "m = Mesh(structured_coords=coords, structured=True)\n", "fuel = from_atom_frac({'U235': 0.045, 'U238': 0.955, 'O16': 2.0}, density=10.7)\n", "cool = from_atom_frac({'H1': 2.0, 'O16': 1.0}, density=1.0)\n", "for i in range(len(m)):\n", " m.mats[i] = cool\n", "m.mats[len(m)/2] = fuel\n", "m.mats[len(m)/4] = fuel\n", "for i, c in enumerate(c5[:-1]):\n", " m.mats[i*len(c5)] = fuel" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 23 }, { "cell_type": "code", "collapsed": false, "input": [ "pf = PyneMoabHex8StaticOutput(m)\n", "s = SlicePlot(pf, 'z', 'density', origin='native')\n", "s.annotate_triangle_facets(\"bllite30matls_zip.h5m\", 1, 0.0, plot_args={\"colors\": 'w'})\n", "s.display()" ], "language": "python", "metadata": {}, "outputs": [ { "html": [ "
" ], "metadata": {}, "output_type": "display_data", "text": [ "" ] } ], "prompt_number": 26 }, { "cell_type": "code", "collapsed": false, "input": [], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 5 } ], "metadata": {} } ] }