[yt-users] projection depth

[Luigi Iapichino]

Dear all,

by comparing volume renderings and projections in my simulation I realised 
that I don't have control on the projection depth, in the latter case. 
This is the relevant part of my script:

...

L = 0.1 # projection width

pf = load(data_dir+fn) # load data

reg = pf.h.region([0.5, 0.5, 0.5],[0.45,0.45,0.45],[0.55,0.55,0.55]) # a small 
region
#reg = pf.h.region([0.5, 0.5, 0.5],[0,0,0],[1,1,1]) # the whole box

pc = PlotCollection(reg, center=[0.5,0.5,0.5])

p = pc.add_projection("Density", 0, weight_field="Density",fig_size=(7.5,6))
    
p.modify["contour"]("Density", ncont=5,take_log=True)

pc.set_width(L, '1')

pc.save("%s_test-proj" % pf,override=True)


If I set L = 0.2, I don't see in the central region of the projection more 
structures than with L=0.1 (only the plot width is, of course, larger). 
Moreover, the result do not depend on the region that I select (either third 
or fourth line). It seems to me that, in all cases, the projection is 
actually performed along the whole simulation box, and not on the region that 
I select. In other words, I can control the plot width, but not the 
projection depth. In the volume rendering, on the other hand, the plot width 
sets the depth, too. Is anything wrong or missing in my script?

If needed, I can send you a few explanatory plots and the whole sample script 
for reproducing my issue.
Looking forward to your suggestions,

  Luigi

-- 

---------------------------------------------------------------

Luigi Iapichino
Zentrum fuer Astronomie der Universitaet Heidelberg
Institut fuer Theoretische Astrophysik
Albert-Ueberle-Str. 2, D-69120 Heidelberg, Germany
Tel: +49 6221 548983, Fax: +49 6221 544221
e-mail: luigi at ita.uni-heidelberg.de
URL: http://www.ita.uni-heidelberg.de/~luigi/