Skip to content

Instantly share code, notes, and snippets.

@jzuhone

jzuhone/offaxis_bug.ipynb

Created April 22, 2014 20:27
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save jzuhone/11193110 to your computer and use it in GitHub Desktop.
Save jzuhone/11193110 to your computer and use it in GitHub Desktop.
Download ZIP
offaxis_bug
Display the source blob
Display the rendered blob
Raw
offaxis_bug.ipynb
{
"metadata": {
"name": "",
"signature": "sha256:e3be5a8f54fdecf49009c353813a324fa9d2288984a6c366644aa513cd4267e0"
},
"nbformat": 3,
"nbformat_minor": 0,
"worksheets": [
{
"cells": [
{
"cell_type": "code",
"collapsed": false,
"input": [
"%matplotlib inline\n",
"import yt"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 1
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"ds = yt.load(\"GasSloshingLowRes/sloshing_low_res_hdf5_plt_cnt_0300\")"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stderr",
"text": [
"yt : [WARNING ] 2014-04-22 16:23:22,072 integer runtime parameter checkpointfilenumber overwrites a simulation scalar of the same name\n"
]
},
{
"output_type": "stream",
"stream": "stderr",
"text": [
"yt : [WARNING ] 2014-04-22 16:23:22,073 integer runtime parameter forcedplotfilenumber overwrites a simulation scalar of the same name\n"
]
},
{
"output_type": "stream",
"stream": "stderr",
"text": [
"yt : [WARNING ] 2014-04-22 16:23:22,073 integer runtime parameter nbegin overwrites a simulation scalar of the same name\n"
]
},
{
"output_type": "stream",
"stream": "stderr",
"text": [
"yt : [WARNING ] 2014-04-22 16:23:22,073 integer runtime parameter plotfilenumber overwrites a simulation scalar of the same name\n"
]
},
{
"output_type": "stream",
"stream": "stderr",
"text": [
"yt : [INFO ] 2014-04-22 16:23:22,089 Parameters: current_time = 9.46712365603e+16\n"
]
},
{
"output_type": "stream",
"stream": "stderr",
"text": [
"yt : [INFO ] 2014-04-22 16:23:22,090 Parameters: domain_dimensions = [8 8 8]\n"
]
},
{
"output_type": "stream",
"stream": "stderr",
"text": [
"yt : [INFO ] 2014-04-22 16:23:22,091 Parameters: domain_left_edge = [ -3.08560000e+24 -3.08560000e+24 -3.08560000e+24]\n"
]
},
{
"output_type": "stream",
"stream": "stderr",
"text": [
"yt : [INFO ] 2014-04-22 16:23:22,091 Parameters: domain_right_edge = [ 3.08560000e+24 3.08560000e+24 3.08560000e+24]\n"
]
},
{
"output_type": "stream",
"stream": "stderr",
"text": [
"yt : [INFO ] 2014-04-22 16:23:22,092 Parameters: cosmological_simulation = 0.0\n"
]
}
],
"prompt_number": 2
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"slc = yt.SlicePlot(ds, [1.0,1.0,1.0], [\"density\"])\n",
"slc.set_axes_unit(\"code_length\")\n",
"slc.show()"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stderr",
"text": [
"yt : [INFO ] 2014-04-22 16:23:22,647 Loading field plugins.\n"
]
},
{
"output_type": "stream",
"stream": "stderr",
"text": [
"yt : [INFO ] 2014-04-22 16:23:22,647 Loaded angular_momentum (8 new fields)\n"
]
},
{
"output_type": "stream",
"stream": "stderr",
"text": [
"yt : [INFO ] 2014-04-22 16:23:22,648 Loaded astro (14 new fields)\n"
]
},
{
"output_type": "stream",
"stream": "stderr",
"text": [
"yt : [INFO ] 2014-04-22 16:23:22,648 Loaded cosmology (20 new fields)\n"
]
},
{
"output_type": "stream",
"stream": "stderr",
"text": [
"yt : [INFO ] 2014-04-22 16:23:22,649 Loaded fluid (54 new fields)\n"
]
},
{
"output_type": "stream",
"stream": "stderr",
"text": [
"yt : [INFO ] 2014-04-22 16:23:22,651 Loaded fluid_vector (86 new fields)\n"
]
},
{
"output_type": "stream",
"stream": "stderr",
"text": [
"yt : [INFO ] 2014-04-22 16:23:22,652 Loaded geometric (101 new fields)\n"
]
},
{
"output_type": "stream",
"stream": "stderr",
"text": [
"yt : [INFO ] 2014-04-22 16:23:22,652 Loaded local (101 new fields)\n"
]
},
{
"output_type": "stream",
"stream": "stderr",
"text": [
"yt : [INFO ] 2014-04-22 16:23:22,652 Loaded magnetic_field (107 new fields)\n"
]
},
{
"output_type": "stream",
"stream": "stderr",
"text": [
"yt : [INFO ] 2014-04-22 16:23:22,653 Loaded species (107 new fields)\n"
]
},
{
"output_type": "stream",
"stream": "stderr",
"text": [
"yt : [INFO ] 2014-04-22 16:23:23,371 xlim = -3085599999999999991611392.000000 3085599999999999991611392.000000\n"
]
},
{
"output_type": "stream",
"stream": "stderr",
"text": [
"yt : [INFO ] 2014-04-22 16:23:23,371 ylim = -3085599999999999991611392.000000 3085599999999999991611392.000000\n"
]
},
{
"ename": "UnitParseError",
"evalue": "Could not find unit symbol 'code_length' in the provided symbols.",
"output_type": "pyerr",
"traceback": [
"\u001b[0;31m---------------------------------------------------------------------------\u001b[0m\n\u001b[0;31mUnitParseError\u001b[0m Traceback (most recent call last)",
"\u001b[0;32m<ipython-input-3-0fc19cdcfd4d>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m()\u001b[0m\n\u001b[1;32m 1\u001b[0m \u001b[0mslc\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0myt\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mSlicePlot\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mds\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m[\u001b[0m\u001b[0;36m1.0\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;36m1.0\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;36m1.0\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m[\u001b[0m\u001b[0;34m\"density\"\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m----> 2\u001b[0;31m \u001b[0mslc\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mset_axes_unit\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m\"code_length\"\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 3\u001b[0m \u001b[0mslc\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mshow\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m/Users/jzuhone/Source/yt-3.x/yt/visualization/plot_container.pyc\u001b[0m in \u001b[0;36mnewfunc\u001b[0;34m(*args, **kwargs)\u001b[0m\n\u001b[1;32m 54\u001b[0m \u001b[0mrv\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mf\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m*\u001b[0m\u001b[0margs\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 55\u001b[0m \u001b[0margs\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_plot_valid\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mFalse\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 56\u001b[0;31m \u001b[0margs\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_setup_plots\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 57\u001b[0m \u001b[0;32mreturn\u001b[0m \u001b[0mrv\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 58\u001b[0m \u001b[0;32mreturn\u001b[0m \u001b[0mnewfunc\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m/Users/jzuhone/Source/yt-3.x/yt/visualization/plot_window.pyc\u001b[0m in \u001b[0;36m_setup_plots\u001b[0;34m(self)\u001b[0m\n\u001b[1;32m 713\u001b[0m \u001b[0maspect\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mnp\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mfloat64\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mpf\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mquan\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;36m1.0\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0munit_y\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m/\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mpf\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mquan\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;36m1.0\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0munit_x\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 714\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 715\u001b[0;31m \u001b[0mextentx\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0;34m[\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mxlim\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mi\u001b[0m\u001b[0;34m]\u001b[0m \u001b[0;34m-\u001b[0m \u001b[0mxc\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0min_units\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0munit_x\u001b[0m\u001b[0;34m)\u001b[0m \u001b[0;32mfor\u001b[0m \u001b[0mi\u001b[0m \u001b[0;32min\u001b[0m \u001b[0;34m(\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;36m1\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 716\u001b[0m \u001b[0mextenty\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0;34m[\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mylim\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mi\u001b[0m\u001b[0;34m]\u001b[0m \u001b[0;34m-\u001b[0m \u001b[0myc\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0min_units\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0munit_y\u001b[0m\u001b[0;34m)\u001b[0m \u001b[0;32mfor\u001b[0m \u001b[0mi\u001b[0m \u001b[0;32min\u001b[0m \u001b[0;34m(\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;36m1\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 717\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m/Users/jzuhone/Source/yt-3.x/yt/units/yt_array.pyc\u001b[0m in \u001b[0;36min_units\u001b[0;34m(self, units)\u001b[0m\n\u001b[1;32m 374\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 375\u001b[0m \"\"\"\n\u001b[0;32m--> 376\u001b[0;31m \u001b[0mnew_units\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_unit_repr_check_same\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0munits\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 377\u001b[0m \u001b[0mconversion_factor\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0munits\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mget_conversion_factor\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mnew_units\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 378\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m/Users/jzuhone/Source/yt-3.x/yt/units/yt_array.pyc\u001b[0m in \u001b[0;36m_unit_repr_check_same\u001b[0;34m(self, units)\u001b[0m\n\u001b[1;32m 320\u001b[0m \u001b[0;31m# let Unit() handle units arg if it's not already a Unit obj.\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 321\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0;32mnot\u001b[0m \u001b[0misinstance\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0munits\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mUnit\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 322\u001b[0;31m \u001b[0munits\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mUnit\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0munits\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mregistry\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0munits\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mregistry\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 323\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 324\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0;32mnot\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0munits\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0msame_dimensions_as\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0munits\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m/Users/jzuhone/Source/yt-3.x/yt/units/unit_object.pyc\u001b[0m in \u001b[0;36m__new__\u001b[0;34m(cls, unit_expr, cgs_value, dimensions, registry, **assumptions)\u001b[0m\n\u001b[1;32m 186\u001b[0m \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 187\u001b[0m \u001b[0;31m# lookup the unit symbols\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 188\u001b[0;31m \u001b[0mcgs_value\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdimensions\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0m_get_unit_data_from_expr\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0munit_expr\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mregistry\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mlut\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 189\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 190\u001b[0m \u001b[0;31m# Create obj with superclass construct.\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m/Users/jzuhone/Source/yt-3.x/yt/units/unit_object.pyc\u001b[0m in \u001b[0;36m_get_unit_data_from_expr\u001b[0;34m(unit_expr, unit_symbol_lut)\u001b[0m\n\u001b[1;32m 439\u001b[0m \u001b[0;31m# Now for the sympy possibilities\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 440\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0misinstance\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0munit_expr\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mSymbol\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 441\u001b[0;31m \u001b[0;32mreturn\u001b[0m \u001b[0m_lookup_unit_symbol\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mstr\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0munit_expr\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0munit_symbol_lut\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 442\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 443\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0misinstance\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0munit_expr\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mNumber\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;32m/Users/jzuhone/Source/yt-3.x/yt/units/unit_object.pyc\u001b[0m in \u001b[0;36m_lookup_unit_symbol\u001b[0;34m(symbol_str, unit_symbol_lut)\u001b[0m\n\u001b[1;32m 506\u001b[0m \u001b[0;31m# no dice\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 507\u001b[0m raise UnitParseError(\"Could not find unit symbol '%s' in the provided \" \\\n\u001b[0;32m--> 508\u001b[0;31m \"symbols.\" % symbol_str)\n\u001b[0m\u001b[1;32m 509\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 510\u001b[0m \u001b[0;32mdef\u001b[0m \u001b[0mvalidate_dimensions\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mdimensions\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;31mUnitParseError\u001b[0m: Could not find unit symbol 'code_length' in the provided symbols."
]
}
],
"prompt_number": 3
},
{
"cell_type": "code",
"collapsed": false,
"input": [],
"language": "python",
"metadata": {},
"outputs": []
}
],
"metadata": {}
}
]
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment