[yt-users] Blocky Projections with RAMSES data when log=True

Butler,Michael J butler85 at astro.ufl.edu
Fri Oct 4 07:09:29 PDT 2013


  

Hi Matt, 

Thanks. The problem is (which was going to be my next
email, and the reason I was trying out 2.5), I can't seem to make plots
in the 3.0 branch anymore since updating it. I don't recall which
version it was when I could, but now I get the error at the end of the
email when I try and make a slice/projection (it happens both with the
PlotCollection method as well as SlicePlot).  

-Mike 

In [5]:
pc.add_slice("Density",2)
---------------------------------------------------------------------------
AssertionError
Traceback (most recent call
last)
/home/itp/butler/yt-x86_64/src/yt-hg/scripts/iyt in ()
----> 1
pc.add_slice("Density",2)

/home/itp/butler/yt-x86_64/src/yt-hg/yt/visualization/plot_collection.py
in add_slice(self, field, axis, coord, center, use_colorbar, figure,
axes, fig_size, periodic, obj, field_parameters)
 427 if obj is None:

428 if field_parameters is None: field_parameters = {}
--> 429 obj =
self.pf.hierarchy.slice(axis, coord, field,
 430 center=center,
**field_parameters)
 431 p =
self._add_plot(PCSlicePlot(

/home/itp/butler/yt-x86_64/src/yt-hg/yt/data_objects/static_output.py
in hierarchy(self)
 223 raise RuntimeError("You should not instantiate
StaticOutput.")
 224 self._instantiated_hierarchy =
self._hierarchy_class(
--> 225 self, data_style=self.data_style)
 226
return self._instantiated_hierarchy
 227 h = hierarchy #
alias

/home/itp/butler/yt-x86_64/src/yt-hg/yt/frontends/ramses/data_structures.py
in __init__(self, pf, data_style)
 301 
 302 self.float_type =
np.float64
--> 303 super(RAMSESGeometryHandler, self).__init__(pf,
data_style)
 304 
 305 def
_initialize_oct_handler(self):

/home/itp/butler/yt-x86_64/src/yt-hg/yt/geometry/geometry_handler.py
in __init__(self, pf, data_style)
 66 
 67 mylog.debug("Setting up
domain geometry.")
---> 68 self._setup_geometry()
 69 
 70
mylog.debug("Initializing data grid data
IO")

/home/itp/butler/yt-x86_64/src/yt-hg/yt/geometry/oct_geometry_handler.py
in _setup_geometry(self)
 48 def _setup_geometry(self):
 49
mylog.debug("Initializing Octree Geometry Handler.")
---> 50
self._initialize_oct_handler()
 51 
 52 def
get_smallest_dx(self):

/home/itp/butler/yt-x86_64/src/yt-hg/yt/frontends/ramses/data_structures.py
in _initialize_oct_handler(self)
 306 nv = len(self.fluid_field_list)

307 self.domains = [RAMSESDomainFile(self.parameter_file, i + 1, nv)
-->
308 for i in range(self.parameter_file['ncpu'])]
 309 total_octs =
sum(dom.local_oct_count #+ dom.ngridbound.sum()
 310 for dom in
self.domains)

/home/itp/butler/yt-x86_64/src/yt-hg/yt/frontends/ramses/data_structures.py
in __init__(self, pf, domain_id, nvar)
 70 self._read_amr_header()
 71
self._read_particle_header()
---> 72 self._read_amr()
 73 
 74
_hydro_offset =
None

/home/itp/butler/yt-x86_64/src/yt-hg/yt/frontends/ramses/data_structures.py
in _read_amr(self)
 248 assert(pos.shape[0] == ng)
 249 n =
self.oct_handler.add(cpu + 1, level - min_level, pos)
--> 250 assert(n
== ng)
 251 self.oct_handler.finalize()
 252 

AssertionError: 

On Fri,
4 Oct 2013 09:11:36 -0400, Matthew Turk wrote: 

> Hi Mike, 
> 
> The
stable branch of yt probably doesn't work with RAMSES anymore... I think
we even disabled it recently. 
> 
> However, the 3.0 branch (listed as
"bleeding edge" on the website) full supports Ramses and is considerably
faster and more memory conservative. I'd suggest you give it a shot - at
this point it's quite stable and production ready. The docs are still
lagging but that's next up; for Ramses it should be mostly identical in
usage. If you do run into any trouble, drop a line here? 
> 
> Matt 
>
On Oct 4, 2013 9:04 AM, "Butler,Michael J" wrote:
> 
>> Hi everyone,
>>

>> I have been working with the stable branch of yt to look at RAMSES
data and have run into a problem. When doing simple projections of e.g.
Density with take_log=True, the projection becomes blocky (see attached
image). This doesn't happen if I take out the take_log line. Here is
what I am doing :
>> 
>> pf = load("output_00427/info_00427.txt")
>>
pf.h
>> pf.field_info['Density'].take_log = True
>>
p=SlicePlot(pf,2,"Density")
>> 
>> Any ideas what is going on here?
>>

>> Thanks,
>> 
>> Mike
>>
_______________________________________________
>> yt-users mailing
list
>> yt-users at lists.spacepope.org [1]
>>
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org [2]

 


Links:
------
[1] mailto:yt-users at lists.spacepope.org
[2]
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
[3]
mailto:butler85 at astro.ufl.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20131004/e7982303/attachment.htm>


More information about the yt-users mailing list