<div dir="ltr">Hi all,<div><br></div><div>I'm trying to create two derived fields (one for stellar density, and one for dark matter density) based <a href="https://bitbucket.org/yt_analysis/yt/src/cee2a160c5351c376afea65bea3462792406154e/yt/frontends/enzo/fields.py?at=yt-2.x&fileviewer=file-view-default#fields.py-350">on</a> <a href="https://bitbucket.org/yt_analysis/yt/src/cee2a160c5351c376afea65bea3462792406154e/yt/frontends/enzo/fields.py?at=yt-2.x&fileviewer=file-view-default#fields.py-394">these</a> from yt2. I'm running yt 3.2.3, and the star_density and dm_density fields seem to have been removed. So far, the yt2 code has transferred fine, with minimal changes:</div><ul><li>amrutils.CICDeposit_3() -> amrutils.CICDeposit.CICDeposit_3()<br></li><li>np.float64(data['dx']) -> np.float64(np.ravel(data['dx'])[0]</li><ul><li>(data['dx'] appears to be a multidimensional array of all the same values of dx, so I'm just taking the first one -- is there a better way to do this?)</li></ul></ul><div>Besides those two quirks, the code runs (I'm able to add a field using ds.add_field), but I get an error when I try to access the data using ds.all_data()['star_density'] (it seems now is when the field is actually computed).</div><div><br></div><div>I get the following error:</div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><span class=""><ipython-input-5-1af8fbcb718a></span><span class=""> in </span><span class="">_spdensity</span><span class="">(field, data)<br></span><span class=""> 1</span><span class=""> </span><span class="">def</span><span class=""> _spdensity</span><span class="">(</span><span class="">field</span><span class="">,</span><span class=""> data</span><span class="">):<br></span><span class="">----> 2</span><span class=""> </span><span class="">blank </span><span class="">=</span><span class=""> np</span><span class="">.</span><span class="">zeros</span><span class="">(</span><span class="">data</span><span class="">.</span><span class="">ActiveDimensions</span><span class="">,</span><span class=""> dtype</span><span class="">=</span><span class="">'float64'</span><span class="">)<br></span><span class=""> 3</span><span class=""> </span><span class="">if</span><span class=""> data</span><span class="">[</span><span class="">'particle_position_x'</span><span class="">].</span><span class="">size </span><span class="">==</span><span class=""> </span><span class="">0</span><span class="">:</span><span class=""> </span><span class="">return</span><span class=""> blank<br></span><span class=""> 4</span><span class=""> filter </span><span class="">=</span><span class=""> data</span><span class="">[</span><span class="">'creation_time'</span><span class="">]</span><span class=""> </span><span class="">></span><span class=""> </span><span class="">0.0<br></span><span class=""> 5</span><span class=""> </span><span class="">if</span><span class=""> </span><span class="">not</span><span class=""> filter</span><span class="">.</span><span class="">any</span><span class="">():</span><span class=""> </span><span class="">return</span><span class=""> blank</span><br><span class="">AttributeError</span><span class="">: 'YTRegion' object has no attribute 'ActiveDimensions'</span></blockquote>
<p class=""><span class=""></span></p>
<div>
<p class=""><span class=""></span></p>
</div><div><br></div><div>All the examples of ActiveDimensions that I saw in the docs and source accessed ActiveDimensions as a property of some grid (i.e. g.ActiveDimensions). Is there a simple way to access that from the data/YTRegion object that I'm missing? I can access ds.index.grids and take the first grid there, but ActiveDimensions differs across the grids (And when I try to do this from within _spdensity, data doesn't have an index.grids object, just fake_index (which can't access the grids)). Is there something I'm missing?</div><div><br></div><div><br></div><div>Thanks,</div><div><br></div><div>Rasmi</div></div>