<p>Hi Patrick, </p>
<p>Is this a FLASH dataset, by any chance? If so, I think I know the problem for why it is not adding the units. I will take a look and get back to you. </p>
<p>Matt</p>
<div class="gmail_quote">On Nov 13, 2012 9:00 AM, "Patrick Rieser" <<a href="mailto:patrick.rieser@uibk.ac.at">patrick.rieser@uibk.ac.at</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
mpccm doesn't seem to be in the output. Do I have to set a switch somewhere so that it is created?<br>
<br>
<br>
level # grids # cells # cells^3<br>
------------------------------<u></u>----------------<br>
0 1 4096 15<br>
1 8 32768 31<br>
2 64 262144 63<br>
3 512 2097152 127<br>
4 512 2097152 127<br>
5 512 2097152 127<br>
6 512 2097152 127<br>
7 512 2097152 127<br>
------------------------------<u></u>----------------<br>
2633 10784768<br>
<br>
<br>
t = 3.03585620e+17 = 3.03585620e+17 s = 9.62004779e+09 years<br>
<br>
Smallest Cell:<br>
Width: 4.883e-04 unitary<br>
Width: 1.562e-02 mpc<br>
Width: 1.562e+01 kpc<br>
Width: 1.562e+04 pc<br>
Width: 3.223e+09 au<br>
Width: 6.934e+11 rsun<br>
Width: 2.995e+17 miles<br>
Width: 4.821e+22 cm<br>
Width: 4.821e+22 1<br>
<br>
<br>
Best wishes,<br>
Patrick<br>
<br>
Am 13.11.2012 14:55, schrieb Matthew Turk:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi Patrick,<br>
<br>
It gets autocreated by concatenating "mpc" with "cm", so it won't show<br>
up directly in the source. If you run pf.h.print_stats() it will<br>
print out all the available length units -- what shows up in there?<br>
<br>
-Matt<br>
<br>
On Tue, Nov 13, 2012 at 7:55 AM, Patrick Rieser<br>
<<a href="mailto:patrick.rieser@uibk.ac.at" target="_blank">patrick.rieser@uibk.ac.at</a>> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
When I try<br>
<br>
pc.set_width((28,'mpccm'))<br>
<br>
it raises a key error:<br>
<br>
<br>
Traceback (most recent call last):<br>
File "projections.py", line 89, in <module><br>
pc.set_width((28,'mpccm'))<br>
File<br>
"/home/patrick/Programme/yt-<u></u>x86_64/src/yt-hg/yt/<u></u>visualization/plot_window.py",<br>
line 81, in newfunc<br>
rv = f(*args, **kwargs)<br>
File<br>
"/home/patrick/Programme/yt-<u></u>x86_64/src/yt-hg/yt/<u></u>visualization/plot_window.py",<br>
line 458, in set_width<br>
width = [w / <a href="http://self.pf" target="_blank">self.pf</a>[unit] for w in width]<br>
File<br>
"/home/patrick/Programme/yt-<u></u>x86_64/src/yt-hg/yt/data_<u></u>objects/static_output.py",<br>
line 143, in __getitem__<br>
raise KeyError(key)<br>
KeyError: 'mpccm'<br>
<br>
<br>
Also when I search through the yt source code, "mpccm" is only found in<br>
universal_fields.py (in the WeakLensingConvergence field) and in<br>
light_ray.py.<br>
<br>
Best wishes,<br>
Patrick<br>
<br>
<br>
Am 12.11.2012 20:32, schrieb Britton Smith:<br>
<br>
Patrick,<br>
<br>
Don't worry about it. I do that stuff all the time.<br>
<br>
Britton<br>
<br>
<br>
On Mon, Nov 12, 2012 at 2:26 PM, Patrick Rieser <<a href="mailto:patrick.rieser@uibk.ac.at" target="_blank">patrick.rieser@uibk.ac.at</a>><br>
wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Ah damn, now I feel stupid. I should have realized that. I will try it<br>
tomorrow, but yes, that should fix it!<br>
<br>
Thanks a lot!<br>
<br>
Best wishes,<br>
Patrick<br>
<br>
Am 12.11.2012 19:38, schrieb Nathan Goldbaum:<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
If Britton is correct (and looking closely at your script, I suspect he<br>
is) you'll likely fix the issue by specifying the width of the image in<br>
comoving units:<br>
<br>
pc.set_width((28,'mpccm'))<br>
<br>
Right now the plot axes will still choose to annotate the x and y axes<br>
using physical units, even though the width of the image will be the<br>
comoving width you requested. Sorry for that bit of nonintuitiveness. Sam<br>
Skillman just filed a bug about this issue which I will fix soon:<br>
<a href="https://bitbucket.org/yt_analysis/yt/issue/463/plotwindow-axes-unit-defaults" target="_blank">https://bitbucket.org/yt_<u></u>analysis/yt/issue/463/<u></u>plotwindow-axes-unit-defaults</a><br>
<br>
-Nathan<br>
<br>
On 11/12/12 10:25 AM, Britton Smith wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi Patrick,<br>
<br>
Is this a cosmology simulation? If so, I think the issue is that the<br>
window is in physical coordinates, which are increasing with time.<br>
<br>
Britton<br>
<br>
<br>
On Mon, Nov 12, 2012 at 1:21 PM, Nathan Goldbaum <<a href="mailto:nathan12343@gmail.com" target="_blank">nathan12343@gmail.com</a><br>
<mailto:<a href="mailto:nathan12343@gmail.com" target="_blank">nathan12343@gmail.com</a>><u></u>> wrote:<br>
<br>
Hi Patrick,<br>
<br>
I'm unable to reproduce the issue you're seeing. I'm running the<br>
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".<br>
I did this following Britton's recipe from the workshop:<br>
<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<br>
<a href="http://yt-project.org/data" target="_blank">yt-project.org/data</a> <<a href="http://yt-project.org/data" target="_blank">http://yt-project.org/data</a>><br>
<br>
It would help me track down what's going wrong if you could come<br>
up with a somewhat simpler script that reproduces the error,<br>
preferably using one of the datasets on <a href="http://yt-project.org" target="_blank">yt-project.org</a><br>
<<a href="http://yt-project.org" target="_blank">http://yt-project.org</a>>. Please feel free to e-mail me in private<br>
or join us on irc so we can iterate on this.<br>
<br>
Cheers,<br>
<br>
Nathan<br>
<br>
<br>
On 11/12/12 9:50 AM, Patrick Rieser wrote:<br>
<br>
Hey all,<br>
<br>
So I still get some strange behaviour with the ProjectionPlot.<br>
If I use the projection.set_width() with the clumps callback<br>
method I get the real image shrinked inside the the plot<br>
window. Strange thing is, that it grows with each step! Here<br>
is a demonstration how it looks (you have to look closely as<br>
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<br>
fine, but the scale on the axis grows with each step. So it<br>
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<br>
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 =<br>
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" %<br>
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'<br>
% (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>
<mailto:<a href="mailto:yt-users@lists.spacepope.org" target="_blank">yt-users@lists.<u></u>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>
<br>
<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> <mailto:<a href="mailto:yt-users@lists.spacepope.org" target="_blank">yt-users@lists.<u></u>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>
<br>
<br>
<br>
<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>
</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>
</blockquote>
<br>
<br>
<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>
<br>
<br>
<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>
<br>
</blockquote>
______________________________<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>
</blockquote></div>