<div dir="ltr"><div><div><div><div>Hello,<br><br></div>I am starting upon an adventure of actually using yt and I am already getting the impression I may be using a chainsaw to cut kindling.<br><br></div>My data is uniform grid, not that large (128^3 for now) and I need to access the fields, both on disk and derived in 3D array format. I have tried using a covering_grid thus:<br>

<br>```<br>ts = yt.load("/archive/io_testing/3D_tube128_Slog_p90_A20r2-3_B005/3D_tube128_Slog_p90_A20r2-3_B005_*.gdf")<br><br>ds = ts[10]<br><br>cg = ds.h.covering_grid(0, ds.domain_left_edge, ds.domain_dimensions)<br>

<br>cg['density_pert']<br>```<br><br></div>which soemtimes seems to give me an array full of -999 and sometimes gives me this lovely traceback:<br><font size="1">```<br><br>---------------------------------------------------------------------------<br>

DataCubeError                             Traceback (most recent call last)<br><ipython-input-55-62e42995fe53> in <module>()<br>----> 1 cg['density_pert']<br><br>/home/stuart/BitBucket/yt/yt/data_objects/data_containers.pyc in __getitem__(self, key)<br>

    316             if key not in self.fields:<br>    317                 self.fields.append(key)<br>--> 318             self.get_data(key)<br>    319         return self.field_data[key]<br>    320 <br><br>/home/stuart/BitBucket/yt/yt/data_objects/data_containers.pyc in get_data(self, fields)<br>

   3757         for i, grid in enumerate(self._grids):<br>   3758             if self._use_pbar: pbar.update(i)<br>-> 3759             count -= self._get_data_from_grid(grid, obtain_fields)<br>   3760             if count <= 0: break<br>

   3761         if self._use_pbar: pbar.finish()<br><br>/home/stuart/BitBucket/yt/yt/data_objects/data_containers.pyc in save_state(self, grid, field, *args, **kwargs)<br>     79         old_keys = grid.field_data.keys()<br>

     80         grid.field_parameters = self.field_parameters<br>---> 81         tr = func(self, grid, field, *args, **kwargs)<br>     82         grid.field_parameters = old_params<br>     83         grid.field_data = YTFieldData( [(k, grid.field_data[k]) for k in old_keys] )<br>

<br>/home/stuart/BitBucket/yt/yt/data_objects/data_containers.pyc in _get_data_from_grid(self, grid, fields)<br>   3801             c_fields, g_fields,<br>   3802             self.ActiveDimensions, grid.ActiveDimensions,<br>

-> 3803             grid.child_mask, self.domain_width, ll, 0)<br>   3804         return count<br>   3805 <br><br>DataCubeError: DataCubeGeneric: Invalid parameters.<br>```</font><br><br></div>I only want my arrays :'(<br>

<div><br>On another, related, note:<br>I have a bit of code that executes in serial at the top of a otherwise embarrassingly parallel piece of vtk code, I need the most optimal io possible for this, as it is all io limited, how well is yt suited to getting one slice out of my gdf files as fast as humanly possible?<br>

<br>Stuart<br></div></div>