Re: [yt-users] adding quantities to halo catalogs (Britton Smith)

Hi Britton, Thanks for the help. I'm getting an error trying to use halo.data_container. 14 hc.add_callback("sphere",factor=1.0) 15 def _gas_mass(halo): ---> 16 return halo.data_container.quantities.total_quantity("cell_mass") 17 18 add_quantity("gas_mass",_gas_mass) AttributeError: 'Halo' object has no attribute 'data_container' I've made my code very simple and am using the example data to debug this. I'm just trying to get the total gas mass for now data_ds = yt.load('Enzo_64/RD0006/RedshiftOutput0006') halos_ds = yt.load('rockstar_halos/halos_0.0.bin') def _gas_mass(halo): return halo.data_container.quantities.total_quantity("cell_mass") add_quantity("gas_mass",_gas_mass) # Instantiate a catalog using those two paramter files hc = HaloCatalog(data_ds=data_ds, halos_ds=halos_ds) hc.add_callback("sphere",factor=1.0) hc.add_quantity("gas_mass") hc.create() Any help would be appreciated. Thanks, Ari

On Wed, Jun 24, 2015 at 3:22 PM, Professor Maller <amaller@citytech.cuny.edu
wrote:
Hi Britton,
Thanks for the help. I'm getting an error trying to use halo.data_container. 14 hc.add_callback("sphere",factor=1.0) 15 def _gas_mass(halo): ---> 16 return halo.data_container.quantities.total_quantity("cell_mass") 17 18 add_quantity("gas_mass",_gas_mass)
AttributeError: 'Halo' object has no attribute 'data_container'
I've made my code very simple and am using the example data to debug this. I'm just trying to get the total gas mass for now
data_ds = yt.load('Enzo_64/RD0006/RedshiftOutput0006') halos_ds = yt.load('rockstar_halos/halos_0.0.bin')
def _gas_mass(halo): return halo.data_container.quantities.total_quantity("cell_mass")
Looking at the code, I think this should be "halo.data_object". You can see the code for the halo_sphere callback in yt/analysis_modules/halo_analysis/halo_callbacks.py, which might be helpful to look at when debugging issues like this.
add_quantity("gas_mass",_gas_mass)
# Instantiate a catalog using those two paramter files hc = HaloCatalog(data_ds=data_ds, halos_ds=halos_ds)
hc.add_callback("sphere",factor=1.0) hc.add_quantity("gas_mass") hc.create()
Any help would be appreciated.
Thanks,
Ari
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
participants (2)
-
Nathan Goldbaum
-
Professor Maller