[FLASH-USERS] Visualizing 3D Cylindrical grid

Rahul Kashyap rkashyap at umassd.edu
Sat Jan 28 19:47:29 EST 2017


Hi Nathan,

Thanks. I think that was the problem. I was annotating grids. I have yt
3.3.2; most recent is 3.3.3.

Is it possible to do so now for cylindrical grids? It might be required for
my analysis at some point so, may be you could point me towards tweaking yt
to get it working.

Best,
-Rahul


%matplotlib inline
import yt

path='super3d_hdf5_chk_0000'
ds = yt.load(path)

s = yt.SlicePlot(ds, "r", ["density","temperature"], center=[0,0,0])
s.zoom(20.0)
#s.annotate_grids()
s.show()
s.save('temp_data/')
#"""

On Sat, Jan 28, 2017 at 7:10 PM, Nathan Goldbaum <nathan12343 at gmail.com>
wrote:

> Hi Rahul,
>
> I'm not able to reproduce the behavior you're seeing. The particular error
> you showed indicates that you are using some of the plot annotation
> functions that are attached to yt plot objects. Not all of these annotation
> functions support curvilinear geometries. Note that in principle many of
> these *could* be implemented, that just hasn't happened yet. Unfortunately
> the error doesn't say which of these functions is triggering the issue, so
> I can't tell that myself based on the information you've given.
>
> Is there any chance you can share a minimal working example that
> demonstrates the issue you're having?
>
> -Nathan
>
> On Sat, Jan 28, 2017 at 6:00 PM, Rahul Kashyap <rkashyap at umassd.edu>
> wrote:
>
>> Hi John,
>>
>> Thanks. Below is the error along with the call back. I'm using SlicePlot
>> function here.
>>
>> s = yt.SlicePlot(ds, "z", ["density","temperature",
>> "velocity_magnitude","pressure"], center="max")
>>
>> It gives me grid informations but, perhaps treats it as cartesian. The
>> angles are in radians.
>> I have provided link to a small FLASH file where data is trivial --
>> https://www.dropbox.com/s/4e01kp3fsd6e3fh/super3d_hdf5_chk_0000?dl=0
>>
>> Appreciated,
>> -Rahul
>>
>> ---------------------------------------------------------------------------YTDataTypeUnsupported                     Traceback (most recent call last)<ipython-input-9-5ab982893667> in <module>()    102 #                weight_field = ("CellMass")).save()    103 #s.show()--> 104 s.save('temp_data/')    105 #"""    106
>> /home/rkashyap/sw/miniconda2/lib/python2.7/site-packages/yt/visualization/plot_container.pyc in newfunc(*args, **kwargs)     76             # it is the responsibility of _setup_plots to     77             # call args[0].run_callbacks()---> 78             args[0]._setup_plots()     79         rv = f(*args, **kwargs)     80         return rv
>> /home/rkashyap/sw/miniconda2/lib/python2.7/site-packages/yt/visualization/plot_window.pyc in _setup_plots(self)    954     955         self._set_font_properties()--> 956         self.run_callbacks()    957         self._plot_valid = True    958
>> /home/rkashyap/sw/miniconda2/lib/python2.7/site-packages/yt/visualization/plot_window.pyc in run_callbacks(self)   1007                     callback(cbw)   1008                 except YTDataTypeUnsupported as e:-> 1009                     six.reraise(YTDataTypeUnsupported, e)   1010                 except Exception as e:   1011                     six.reraise(YTPlotCallbackError,
>> /home/rkashyap/sw/miniconda2/lib/python2.7/site-packages/six.pyc in reraise(tp, value, tb)
>> YTDataTypeUnsupported: This operation is not supported for data of geometry cylindrical; It supports data of geometries ('cartesian', 'spectral_cube')
>>
>>
>>
>> On Sat, Jan 28, 2017 at 6:18 PM, John ZuHone <jzuhone at gmail.com> wrote:
>>
>>> Hi Rahul,
>>>
>>> Could you detail exactly what the issue with yt was? We'd like to help
>>> out.
>>>
>>> John
>>>
>>> Sent from my iPhone
>>>
>>> On Jan 28, 2017, at 3:12 PM, Rahul Kashyap <rkashyap at umassd.edu> wrote:
>>>
>>> Dear FLASH users,
>>>
>>> Has anyone come up with visualizing 3D cylindrical data set by FLASH?
>>>
>>> I tried VisIt by transforming coordinate from cylindrical to cartesian.
>>> This sets up the grid with z-coordinate from 0 to 6.28 i.e. 2*pi i.e. does
>>> not consider it as angle.
>>> It might be that it's due to the order of coordinates used by visit
>>> which is r,θ,z​ as opposed to r,z,θ​ in FLASH.
>>> yt has also not worked.
>>>
>>> Any help is highly appreciated.
>>>
>>> Thanks,
>>> -Rahul
>>>
>>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20170128/2c83f974/attachment-0001.htm>


More information about the flash-users mailing list