<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Thanks, it works now! :)<div>-Ryan</div><div><br></div><div><div>On Oct 18, 2012, at 3:22 PM, Nathan Goldbaum wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div>Hi Ryan,</div><div><br></div><div>Thanks for the bug report. I've fixed the issue and we've merged the changes into the tip of the development branch.</div><div><br></div><div>Matt has also added tests to help ensure that this sort of bug won't crop up again in the future.</div><div><br></div><div>Please let us know if you have further trouble.</div><div><br></div><div>Cheers,</div><div><br></div><div>Nathan</div><br><div><div>On Oct 18, 2012, at 8:44 AM, M Ryan Joung wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div>Hello,</div><div><br>I am trying to produce an all-sky map projection of a field. Maybe I did something wrong here, but even after the fix described in </div><div> <a href="https://bitbucket.org/yt_analysis/yt/issue/451/radialvelocity-now-seems-to-require-a" target="_blank">https://bitbucket.org/yt_analysis/yt/issue/451/radialvelocity-now-seems-to-require-a</a></div><div>(which sets 'normal' to [0,0,1]), there seems to be an array dimension mismatch between Jy and theta, phi, as you can see below. Can anyone help me with this problem? I attach a logfile (output with --detail) with this error.<br><br>Thank you,</div><div>-Ryan</div><div><br></div><div><br></div><div>The part of the script where the error occurs:</div><div>-----------------------------</div><div>[snipped]</div><div><div>for _fields in grouper(3, field_names):</div><div> f = h5py.File("hi_maps.h5", "a")</div><div> fields = [fn for fn in _fields if fn not in f["/"]]</div><div> allsky = camera.allsky_projection(pf, c,</div><div> 300.0/pf['kpc'], NP, fields,</div><div> inner_radius = rr, rotation = rotation) <-- line 60</div></div><div>[snipped]</div><div>-----------------------------</div><div>-----------------------------<br>Traceback (most recent call last):<br> File "sky_maps.py", line 60, in <module><br> inner_radius = rr, rotation = rotation)<br> File "/Users/mkjoung/work/yt/yt-i386/src/yt-hg/yt/visualization/volume_rendering/camera.py", line 1536, in allsky_projection<br> for field in fields]<br> File "/Users/mkjoung/work/yt/yt-i386/src/yt-hg/yt/data_objects/grid_patch.py", line 157, in __getitem__<br> self.get_data(key)<br> File "/Users/mkjoung/work/yt/yt-i386/src/yt-hg/yt/data_objects/grid_patch.py", line 200, in get_data<br> self._generate_field(field)<br> File "/Users/mkjoung/work/yt/yt-i386/src/yt-hg/yt/data_objects/grid_patch.py", line 145, in _generate_field<br> self[field] = self.pf.field_info[field](self)<br> File "/Users/mkjoung/work/yt/yt-i386/src/yt-hg/yt/data_objects/field_info_container.py", line 316, in __call__<br> dd = self._function(self, data)<br> File "sky_maps.py", line 50, in VelocityTransmittedHI<br> (((vel * 1e5) < data["RadialVelocity"])<br> File "/Users/mkjoung/work/yt/yt-i386/src/yt-hg/yt/data_objects/grid_patch.py", line 157, in __getitem__<br> self.get_data(key)<br> File "/Users/mkjoung/work/yt/yt-i386/src/yt-hg/yt/data_objects/grid_patch.py", line 200, in get_data<br> self._generate_field(field)<br> File "/Users/mkjoung/work/yt/yt-i386/src/yt-hg/yt/data_objects/grid_patch.py", line 145, in _generate_field<br> self[field] = self.pf.field_info[field](self)<br> File "/Users/mkjoung/work/yt/yt-i386/src/yt-hg/yt/data_objects/field_info_container.py", line 316, in __call__<br> dd = self._function(self, data)<br> File "/Users/mkjoung/work/yt/yt-i386/src/yt-hg/yt/data_objects/universal_fields.py", line 885, in _RadialVelocity<br> return get_sph_r_component(velocities, theta, phi, normal)<br> File "/Users/mkjoung/work/yt/yt-i386/src/yt-hg/yt/utilities/math_utils.py", line 797, in get_sph_r_component<br> Jy*np.sin(theta)*np.sin(phi) + \<br>ValueError: operands could not be broadcast together with shapes (128,128,128,3) (128,384,128)<br>-----------------------------</div><div><br></div><div></div></div><span><log3></span><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div></div></div>_______________________________________________<br>yt-users mailing list<br><a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br><a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br></blockquote></div><br></div>_______________________________________________<br>yt-users mailing list<br><a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org<br></blockquote></div><br></body></html>