Hi Patrick,<br><br>Is this a cosmology simulation? If so, I think the issue is that the window is in physical coordinates, which are increasing with time.<br><br>Britton<br><div class="gmail_extra"><br><br><div class="gmail_quote">
On Mon, Nov 12, 2012 at 1:21 PM, Nathan Goldbaum <span dir="ltr"><<a href="mailto:nathan12343@gmail.com" target="_blank">nathan12343@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi Patrick,<br>
<br>
I'm unable to reproduce the issue you're seeing. I'm running the following script:<br>
<br>
from yt.mods import * # set up our namespace<br>
from yt.analysis_modules.level_<u></u>sets.api import *<br>
<br>
pf = load("galaxy0030/galaxy0030")<br>
<br>
master_clump = pf.h.load_object('My_Clumps')<br>
<br>
all_clumps = get_lowest_clumps(master_<u></u>clump)<br>
<br>
for i in range(1,3):<br>
<br>
prj = ProjectionPlot(pf,2,'Density',<u></u>center='c')<br>
<br>
prj.set_width((20,'kpc'))<br>
<br>
slc.annotate_text([0.0,1.05], "Clump %s" % i)<br>
<br>
slc.annotate_clumps([all_<u></u>clumps[i]])<br>
<br>
slc.save(str(i))<br>
<br>
I've previously saved the clump objects in the file "My_clumps". I did this following Britton's recipe from the workshop:<br>
<br>
<a href="https://bitbucket.org/brittonsmith/yt.workshop2012.clump-finding/src/1e7af99cec95fb307bb79055f908d39a200b7091/scripts/find_clumps_and_save.py?at=default" target="_blank">https://bitbucket.org/<u></u>brittonsmith/yt.workshop2012.<u></u>clump-finding/src/<u></u>1e7af99cec95fb307bb79055f908d3<u></u>9a200b7091/scripts/find_<u></u>clumps_and_save.py?at=default</a><br>
<br>
When I run the script, I get the following two images:<br>
<br>
<a href="http://i.imgur.com/MhkaQ.png" target="_blank">http://i.imgur.com/MhkaQ.png</a><br>
<a href="http://i.imgur.com/ZGzTl.png" target="_blank">http://i.imgur.com/ZGzTl.png</a><br>
<br>
This script uses the IsolatedGalaxy dataset stored at <a href="http://yt-project.org/data" target="_blank">yt-project.org/data</a><br>
<br>
It would help me track down what's going wrong if you could come up with a somewhat simpler script that reproduces the error, preferably using one of the datasets on <a href="http://yt-project.org" target="_blank">yt-project.org</a>. Please feel free to e-mail me in private or join us on irc so we can iterate on this.<br>
<br>
Cheers,<br>
<br>
Nathan<div class="HOEnZb"><div class="h5"><br>
<br>
On 11/12/12 9:50 AM, Patrick Rieser wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hey all,<br>
<br>
So I still get some strange behaviour with the ProjectionPlot. If I use the projection.set_width() with the clumps callback method I get the real image shrinked inside the the plot window. Strange thing is, that it grows with each step! Here is a demonstration how it looks (you have to look closely as it doesn't grow that much each step):<br>
<br>
<a href="http://www.flickr.com/photos/87912862@N05/sets/72157631990474435/" target="_blank">http://www.flickr.com/photos/<u></u>87912862@N05/sets/<u></u>72157631990474435/</a><br>
<br>
If I don't use the set_width() then the image itself looks fine, but the scale on the axis grows with each step. So it seems to be the same error than with fixed width.<br>
<br>
<br>
Best wishes,<br>
Patrick<br>
<br>
<br>
<br>
Here is the code that I am using:<br>
<br>
# ID is just a list of tuples with clumps/their position in an array<br>
# sim_files and clump_files is a list containing the filenames<br>
# get_myclump() simply returns a specific clump<br>
<br>
for j, series in enumerate(ID):<br>
for i in series:<br>
pf = load(sim_files[i[0]])<br>
myclump = get_myclump(clump_files[i[0]], i[1])<br>
<br>
pc = ProjectionPlot(pf, axis, my_field, center = pf.domain_center, weight_field = my_wfield)<br>
pc.set_width(28, "mpc")<br>
pc.annotate_text([0.0,1.05], "Clump %s" % j)<br>
<br>
if redshift:<br>
pc.annotate_text([0.917,1.05], "z = %0.4f" % pf.current_redshift)<br>
if grid:<br>
pc.annotate_grids()<br>
if fix_scale:<br>
pc.set_zlim(my_field,<u></u>boundaries[0],boundaries[1])<br>
<br>
pc.set_cmap("Density", "idl01")<br>
pc.annotate_clumps(myclump)<br>
<br>
pc.save('%s/Clump_%s_File_%s_<u></u>Nr_%s_%s_Projection.png' % (myBASEDIR, j, i[0], i[1], axis))<br>
______________________________<u></u>_________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org" target="_blank">yt-users@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" target="_blank">http://lists.spacepope.org/<u></u>listinfo.cgi/yt-users-<u></u>spacepope.org</a><br>
</blockquote>
<br>
______________________________<u></u>_________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org" target="_blank">yt-users@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" target="_blank">http://lists.spacepope.org/<u></u>listinfo.cgi/yt-users-<u></u>spacepope.org</a><br>
</div></div></blockquote></div><br></div>