[yt-users] Make projections for uniform grid FLASH4 data
John ZuHone
jzuhone at gmail.com
Tue Aug 7 12:42:48 PDT 2012
Suoqing,
Just a quick question... I notice that your dataset is 2D. Have you tried using a SlicePlot instead? If your dataset is 2D, projections and slices should be identical.
I'm not saying that the projection should fail like it does, but I suspect that you probably would not have trouble with a slice and you would get the same result.
Best,
John Z
On Aug 7, 2012, at 3:40 PM, Ji Suoqing wrote:
> Ah, thanks a lot! With your suggestion I have axis ticks now.
>
> And the ug projection is still a problem.
>
> Best wishes,
> Suoqing
>
> 在 2012-8-7，下午3:34，Ji Suoqing <jisuoqing at gmail.com> 写道：
>
>> Hi Matt,
>>
>> Thanks for you response. I've tried:
>>
>>>>> from yt.mods import *
>>>>> pf=load('relax_hdf5_chk_0000')
>> yt : [INFO ] 2012-08-07 23:39:54,570 Parameters: current_time = 0.0
>> yt : [INFO ] 2012-08-07 23:39:54,571 Parameters: domain_dimensions = [512 256 1]
>> yt : [INFO ] 2012-08-07 23:39:54,571 Parameters: domain_left_edge = [ 1.50000000e+07 -6.55360000e+09 0.00000000e+00]
>> yt : [INFO ] 2012-08-07 23:39:54,571 Parameters: domain_right_edge = [ 2.62144000e+10 6.55360000e+09 1.74532925e-02]
>> yt : [INFO ] 2012-08-07 23:39:54,572 Parameters: cosmological_simulation = 0.0
>>>>> p = ProjectionPlot(pf, "z", "Density")
>> yt : [INFO ] 2012-08-07 23:40:02,771 Adding unknown field fspd to list of fields
>> yt : [INFO ] 2012-08-07 23:40:02,772 Adding unknown field spec to list of fields
>> Traceback (most recent call last):
>> File "<stdin>", line 1, in <module>
>> File "/aphelion/jisuoqing/yt-x86_64/src/yt-hg/yt/visualization/plot_window.py", line 875, in __init__
>> proj = pf.h.proj(axis,fields,weight_field=weight_field,max_level=max_level,center=center)
>> File "/aphelion/jisuoqing/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 1697, in __init__
>> self._initialize_source(source)
>> File "/aphelion/jisuoqing/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 1724, in _initialize_source
>> source = self.pf.h.all_data()
>> File "/aphelion/jisuoqing/yt-x86_64/src/yt-hg/yt/data_objects/hierarchy.py", line 188, in all_data
>> pf.domain_left_edge, pf.domain_right_edge)
>> File "/aphelion/jisuoqing/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 3221, in __init__
>> self._refresh_data()
>> File "/aphelion/jisuoqing/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 315, in _refresh_data
>> self.get_data()
>> File "/aphelion/jisuoqing/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 2508, in get_data
>> raise YTNoDataInObjectError(self)
>> yt.utilities.exceptions.YTNoDataInObjectError: The object requested has no data included in it.
>>
>> All best,
>> Suoqing
>>
>> 在 2012-8-7，下午3:24，Matthew Turk <matthewturk at gmail.com> 写道：
>>
>>> Hi Suoqing,
>>>
>>> Could you try it with the ProjectionPlot call? It should look like this:
>>>
>>> p = ProjectionPlot(pf, "z", "Density")
>>> p.save()
>>>
>>> -Matt
>>>
>>> On Tue, Aug 7, 2012 at 3:14 PM, Ji Suoqing <jisuoqing at gmail.com> wrote:
>>>> Hi all,
>>>>
>>>> I'm trying to make projections for uniform grid FLASH4 data with the most updated YT, but there're errors when adding projections:
>>>>
>>>> yt : [WARNING ] 2012-08-07 23:13:56,239 integer runtime parameter checkpointfilenumber overwrites a simulation scalar of the same name
>>>> yt : [INFO ] 2012-08-07 23:13:56,244 Parameters: current_time = 0.0
>>>> yt : [INFO ] 2012-08-07 23:13:56,245 Parameters: domain_dimensions = [512 256 1]
>>>> yt : [INFO ] 2012-08-07 23:13:56,245 Parameters: domain_left_edge = [ 1.50000000e+07 -6.55360000e+09 0.00000000e+00]
>>>> yt : [INFO ] 2012-08-07 23:13:56,245 Parameters: domain_right_edge = [ 2.62144000e+10 6.55360000e+09 1.74532925e-02]
>>>> yt : [INFO ] 2012-08-07 23:13:56,245 Parameters: cosmological_simulation = 0.0
>>>> yt : [INFO ] 2012-08-07 23:13:56,587 Max Value is 2.46340e+06 at 38378027.3437500000000000 6400000.0000000000000000 0.5000000000000000 in grid FLASHGrid_0033 ([128 128 1]) at level 0 (1, 0, 0)
>>>> yt : [INFO ] 2012-08-07 23:13:56,587 Created plot collection with default plot-center = [13114700000.0, 0.0, 0.0087266462599716477]
>>>> Traceback (most recent call last):
>>>> File "proj.py", line 18, in <module>
>>>> p1 = pc.add_projection(my_var, 2)
>>>> File "/aphelion/jisuoqing/yt-x86_64/src/yt-hg/yt/visualization/plot_collection.py", line 784, in add_projection
>>>> **field_parameters)
>>>> File "/aphelion/jisuoqing/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 1697, in __init__
>>>> self._initialize_source(source)
>>>> File "/aphelion/jisuoqing/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 1724, in _initialize_source
>>>> source = self.pf.h.all_data()
>>>> File "/aphelion/jisuoqing/yt-x86_64/src/yt-hg/yt/data_objects/hierarchy.py", line 188, in all_data
>>>> pf.domain_left_edge, pf.domain_right_edge)
>>>> File "/aphelion/jisuoqing/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 3221, in __init__
>>>> self._refresh_data()
>>>> File "/aphelion/jisuoqing/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 315, in _refresh_data
>>>> self.get_data()
>>>> File "/aphelion/jisuoqing/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 2508, in get_data
>>>> raise YTNoDataInObjectError(self)
>>>> yt.utilities.exceptions.YTNoDataInObjectError: The object requested has no data included in it.
>>>>
>>>>
>>>> And another separated question: how can I add axis ticks onto projection or slice plot to show the spacial scale? I suspect there should exist some easy way, but I failed to find it.
>>>>
>>>> Thanks for your care!
>>>>
>>>> All best,
>>>> Suoqing
>>>> _______________________________________________
>>>> yt-users mailing list
>>>> yt-users at lists.spacepope.org
>>>> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>>> _______________________________________________
>>> yt-users mailing list
>>> yt-users at lists.spacepope.org
>>> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>>
>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
More information about the yt-users
mailing list