<div dir="ltr">Hi all,<div><br></div><div>I am trying to write a bunch of information from a gas-only single enzo output to a file for colleagues to read.  I don't think they are familiar with HDF5, so I am hoping to write to a fits file.  What I want to put in the file is the x,y,z coordinates, the x,y,z velocity, and the total energy of the gas--which is a user-defined function in my case.  Ideally I also only want to include gas that has a tracer fraction of some value.  </div><div><br></div><div>so basically I have something like:</div><div><br></div><div>ds = yt.load("file")</div><div>galgas = ds.cut_region(["obj['Metal_Density']/obj['density'] > 0.1"])</div><div><br></div><div>can I just create an array of the fields I want from galgas?</div><div><br></div><div>Thanks,</div><div>Stephanie</div><div><br></div><div><div><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div><div><font face="garamond,serif">--<br>Dr. Stephanie Tonnesen<br></font></div><font face="garamond,serif">Alvin E. Nashman Postdoctoral Fellow<br></font></div><font face="garamond,serif">Carnegie Observatories, Pasadena, CA<br></font></div><font face="garamond,serif"><a href="mailto:stonnes@gmail.com" target="_blank">stonnes@gmail.com</a><br></font></div></div></div>
</div></div>