{
"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": {}
}
]
}