<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">Hi Aris,<div><br></div><div>I assume you are dealing with FLASH 2-d cylindrical data. What do you mean by saying you got wrong results? Have you used the correct definition of cell volume (2 * pi * r * dr * dz for 2d cylindrical) when you compute volume-weighted or cellmass-weighted quantities?</div><div><br></div><div>Best wishes,</div><div>
<div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><div>--</div><div>Suoqing JI</div><div>Ph.D Student</div><div>Department of Physics</div><div>University of California, Santa Barbara</div><div>CA 93106, USA</div></div></div></div>
</div>
<br><div><div>On May 9, 2014, at 10:11 AM, Τρίτσης Άρης- Ευάγγελος <<a href="mailto:tritsis@physics.uoc.gr">tritsis@physics.uoc.gr</a>> wrote:</div><blockquote type="cite"><br>Hi Matt,<br><br>Thanks a million for you immediate responses!! Since my main objective is to produce column density plots<br>and then projection plots can you also tell me if such a feature is implemented for cylindrical coordinates?<br>When I try to produce such plots with the ProjectionPlot and ProfilePlot routines I get wrong results. I<br>have attached the plots I get from yt plus a dump from my model so you can check it if you like. Again I would really<br>appreciate any help :).<br><br>Aris<span><FlatCollIso_hdf5_plt_cnt_0157_Projection_theta_dens.png></span><span><FlatCollIso_hdf5_plt_cnt_0157_1d-Profile_r_density.png></span><span><FlatCollIso_hdf5_plt_cnt_0157></span>_______________________________________________<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</blockquote><blockquote type="cite">On 2014-05-09 19:04, Matthew Turk wrote:<br><blockquote type="cite">Hi Aris,<br>Ah!  I see the problem -- the off-axis slice plot probably isn't<br>implemented.  I'm sorry to say that might take a bit more work and<br>won't be a readily added improvement; on-axis should work, however.<br>-Matt<br>On Fri, May 9, 2014 at 12:02 PM, Τρίτσης Άρης- Ευάγγελος<br><<a href="mailto:tritsis@physics.uoc.gr">tritsis@physics.uoc.gr</a>> wrote:<br><blockquote type="cite">On 2014-05-09 18:43, Matthew Turk wrote:<br><blockquote type="cite">Hi Aris,<br>On Fri, May 9, 2014 at 9:21 AM, Τρίτσης Άρης- Ευάγγελος<br><<a href="mailto:tritsis@physics.uoc.gr">tritsis@physics.uoc.gr</a>> wrote:<br><blockquote type="cite">Hi yt users!<br>I am new to yt and would like to ask if anyone has come across issues<br>similar to mine.<br>Is it possible to make slice,projection and 1D-profle plots in<br>cylindrical<br>coordinates in a Flash output file<br>in any orientation and if not does anyone know when such a feature will<br>be<br>implemented to yt?<br>Secondly, when I try to update all dependencies I get the following<br>error:<br></blockquote>In terms of cylindrical coordinates, yt-3.0 ("bleeding edge" on the<br>website) should fully support them.  I'm not entirely sure I<br>understand the question of orientation, though.<br><blockquote type="cite">/home/*/yt-x86_64/bin/python2.7: can't open file '/distribute_setup.py':<br>[Errno 2] No such file or directory<br></blockquote>Ah, I think this might mean you are missing the libssl-dev<br>dependencies?  I'm not sure, though.  Kacper?<br>-Matt<br><blockquote type="cite">I would appreciate any help :)<br>Aris<br>_______________________________________________<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>_______________________________________________<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>Hi Matt,<br>I do use the bleeding edge version which I installed from the installation<br>script. However, when I try to produce an off axis slice plot I get the<br>following:<br>YTCoordinateNotImplemented                Traceback (most recent call last)<br><ipython-input-67-440c2559450b> in <module>()<br>----> 1 cut=OffAxisSlicePlot(pf, L, 'density', width=(25, 'pc'),<br>north_vector=north_vector)<br>/home/tritsis/yt-x86_64/src/yt-hg/yt/visualization/plot_window.pyc in<br>__init__(self, pf, normal, fields, center, width, axes_unit, north_vector,<br>fontsize, field_parameters)<br>   1246         PWViewerMPL.__init__(self, cutting, bounds, fields=fields,<br>   1247                              origin='center-window',periodic=False,<br>-> 1248                              oblique=True, fontsize=fontsize)<br>   1249         if axes_unit is None:<br>   1250             axes_unit = get_axes_unit(width, pf)<br>/home/tritsis/yt-x86_64/src/yt-hg/yt/visualization/plot_window.pyc in<br>__init__(self, *args, **kwargs)<br>    675         if self._plot_type is None:<br>    676             self._plot_type = kwargs.pop("plot_type")<br>--> 677         PlotWindow.__init__(self, *args, **kwargs)<br>    678<br>    679     def _setup_origin(self):<br>/home/tritsis/yt-x86_64/src/yt-hg/yt/visualization/plot_window.pyc in<br>__init__(self, data_source, bounds, buff_size, antialias, periodic, origin,<br>oblique, window_size, fields, fontsize, aspect, setup)<br>    300         self.fields = fields<br>    301         super(PlotWindow, self).__init__(data_source, window_size,<br>fontsize)<br>--> 302         self._set_window(bounds) # this automatically updates the<br>data and plot<br>    303         self.origin = origin<br>    304         if self.data_source.center is not None and oblique == False:<br>/home/tritsis/yt-x86_64/src/yt-hg/yt/visualization/plot_container.pyc in<br>newfunc(*args, **kwargs)<br>     29         args[0]._plot_valid = False<br>     30         if hasattr(args[0], '_recreate_frb'):<br>---> 31             args[0]._recreate_frb()<br>     32         if args[0]._initfinished:<br>     33             args[0]._setup_plots()<br>/home/tritsis/yt-x86_64/src/yt-hg/yt/visualization/plot_window.pyc in<br>_recreate_frb(self)<br>    349                                         periodic=self._periodic)<br>    350         if old_fields is None:<br>--> 351             self.frb._get_data_source_fields()<br>    352         else:<br>    353             for key, unit in zip(old_fields, old_units):<br>/home/tritsis/yt-x86_64/src/yt-hg/yt/visualization/fixed_resolution.pyc in<br>_get_data_source_fields(self)<br>    140         for f in fields:<br>    141             if f not in exclude and f[0] not in<br>self.data_source.pf.particle_types:<br>--> 142                 self[f]<br>    143<br>    144<br>/home/tritsis/yt-x86_64/src/yt-hg/yt/visualization/fixed_resolution.pyc in<br>__getitem__(self, item)<br>    370     def __getitem__(self, item):<br>    371         if item in self.data: return self.data[item]<br>--> 372         indices = np.argsort(self.data_source['dx'])[::-1]<br>    373         bounds = []<br>    374         for b in self.bounds:<br>/home/tritsis/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.pyc in<br>__getitem__(self, key)<br>    210                 return self.field_data[f]<br>    211             else:<br>--> 212                 self.get_data(f)<br>    213         # fi.units is the unit expression string. We depend on the<br>registry<br>    214         # hanging off the dataset to define this unit object.<br>/home/tritsis/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.pyc in<br>get_data(self, fields)<br>    619<br>    620         fields_to_generate += gen_fluids + gen_particles<br>--> 621         self._generate_fields(fields_to_generate)<br>    622<br>    623     def _generate_fields(self, fields_to_generate):<br>/home/tritsis/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.pyc in<br>_generate_fields(self, fields_to_generate)<br>    636                 fi = self.pf._get_field_info(*field)<br>    637                 try:<br>--> 638                     fd = self._generate_field(field)<br>    639                     if type(fd) == np.ndarray:<br>    640                         fd = self.pf.arr(fd, fi.units)<br>/home/tritsis/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.pyc in<br>_generate_field(self, field)<br>    247                 tr = self._generate_particle_field(field)<br>    248             else:<br>--> 249                 tr = self._generate_fluid_field(field)<br>    250             if tr is None:<br>    251                 raise YTCouldNotGenerateField(field, self.pf)<br>/home/tritsis/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.pyc in<br>_generate_fluid_field(self, field)<br>    267             rv = self._generate_spatial_fluid(field,<br>ngt_exception.ghost_zones)<br>    268         else:<br>--> 269             rv = finfo(gen_obj)<br>    270         return rv<br>    271<br>/home/tritsis/yt-x86_64/src/yt-hg/yt/fields/derived_field.pyc in<br>__call__(self, data)<br>    176                 "for %s" % (self.name,))<br>    177         with self.unit_registry(data):<br>--> 178             dd = self._function(self, data)<br>    179         for field_name in data.keys():<br>    180             if field_name not in original_fields:<br>/home/tritsis/yt-x86_64/src/yt-hg/yt/geometry/coordinate_handler.pyc in<br>_unknown_coord(field, data)<br>     31<br>     32 def _unknown_coord(field, data):<br>---> 33     raise YTCoordinateNotImplemented<br>     34<br>     35 def _get_coord_fields(axi, units = "code_length"):<br>YTCoordinateNotImplemented: This coordinate is not implemented for this<br>geometry type<br>Plus I checked and I have all libssl-dev dependencies. Thanks again.<br>_______________________________________________<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>_______________________________________________<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><br></blockquote></div><br></body></html>