question on the halo profiler
Dear all, I have the following (rather standard, I guess) problem: I want to produce radial profiles of selected quantities, from a number of halos in my cosmological box. I am using the yt halo profiler, in the unstable branch. So far, I am able to add filters in such a way to select as many halos as I want, but I have problems with setting the outer radius of my profiles. Here are the relevant parts of my script: ... import yt.analysis_modules.halo_profiler.api as HP ... hp = HP.HaloProfiler(data_dir, halos="multiple",halo_radius=500,radius_units="kpc",n_profile_bins=20) ... In particular, I am not able to change the outer radius by setting halo_radius. I read in the manual: " halo_radius (float): if no halo radii are provided in the halo list file, this parameter ..." etc. Now, in the halo list file provided by the HOP tool (within yt) there is an entry called max_r. Is it a problem if I want to override it and set the outer radius differently? As a side note, however max_r is >not< the current outer radius that I found in my in Halo_0000_profile.dat . I look forward to your suggestions, Cheers, Luigi -- --------------------------------------------------------------- Luigi Iapichino Zentrum fuer Astronomie der Universitaet Heidelberg Institut fuer Theoretische Astrophysik Albert-Ueberle-Str. 2, D-69120 Heidelberg, Germany Tel: +49 6221 548983, Fax: +49 6221 544221 e-mail: luigi@ita.uni-heidelberg.de URL: http://www.ita.uni-heidelberg.de/~luigi/
Hi Luigi, The original behavior of the HaloProfiler was to only use the halo_radius keyword if the radius was not in the halo list. I just pushed a change such that if you specify the radius with the halo_radius command, it will use that even if another radius is present in the halo list. This should not change the behavior for anyone using the halo profiler with radii from the halo list, but someone please let me know if any issues crop up. If you update your yt source and reinstall you should be all set. Britton On Thu, May 19, 2011 at 4:29 AM, Luigi Iapichino < luigi@ita.uni-heidelberg.de> wrote:
Dear all,
I have the following (rather standard, I guess) problem: I want to produce radial profiles of selected quantities, from a number of halos in my cosmological box. I am using the yt halo profiler, in the unstable branch. So far, I am able to add filters in such a way to select as many halos as I want, but I have problems with setting the outer radius of my profiles.
Here are the relevant parts of my script:
... import yt.analysis_modules.halo_profiler.api as HP
...
hp = HP.HaloProfiler(data_dir, halos="multiple",halo_radius=500,radius_units="kpc",n_profile_bins=20) ...
In particular, I am not able to change the outer radius by setting halo_radius. I read in the manual: " halo_radius (float): if no halo radii are provided in the halo list file, this parameter ..." etc.
Now, in the halo list file provided by the HOP tool (within yt) there is an entry called max_r. Is it a problem if I want to override it and set the outer radius differently? As a side note, however max_r is >not< the current outer radius that I found in my in Halo_0000_profile.dat . I look forward to your suggestions, Cheers,
Luigi
--
---------------------------------------------------------------
Luigi Iapichino Zentrum fuer Astronomie der Universitaet Heidelberg Institut fuer Theoretische Astrophysik Albert-Ueberle-Str. 2, D-69120 Heidelberg, Germany Tel: +49 6221 548983, Fax: +49 6221 544221 e-mail: luigi@ita.uni-heidelberg.de URL: http://www.ita.uni-heidelberg.de/~luigi/ _______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
Thanks Britton, now it works in the way I wished! Cheers, Luigi On Thursday 19 May 2011, Britton Smith wrote:
Hi Luigi,
The original behavior of the HaloProfiler was to only use the halo_radius keyword if the radius was not in the halo list. I just pushed a change such that if you specify the radius with the halo_radius command, it will use that even if another radius is present in the halo list. This should not change the behavior for anyone using the halo profiler with radii from the halo list, but someone please let me know if any issues crop up.
If you update your yt source and reinstall you should be all set.
Britton
On Thu, May 19, 2011 at 4:29 AM, Luigi Iapichino <
luigi@ita.uni-heidelberg.de> wrote:
Dear all,
I have the following (rather standard, I guess) problem: I want to produce radial profiles of selected quantities, from a number of halos in my cosmological box. I am using the yt halo profiler, in the unstable branch. So far, I am able to add filters in such a way to select as many halos as I want, but I have problems with setting the outer radius of my profiles.
Here are the relevant parts of my script:
... import yt.analysis_modules.halo_profiler.api as HP
...
hp = HP.HaloProfiler(data_dir, halos="multiple",halo_radius=500,radius_units="kpc",n_profile_bins=20) ...
In particular, I am not able to change the outer radius by setting halo_radius. I read in the manual: " halo_radius (float): if no halo radii are provided in the halo list file, this parameter ..." etc.
Now, in the halo list file provided by the HOP tool (within yt) there is an entry called max_r. Is it a problem if I want to override it and set the outer radius differently? As a side note, however max_r is >not< the current outer radius that I found in my in Halo_0000_profile.dat . I look forward to your suggestions, Cheers,
Luigi
--
---------------------------------------------------------------
Luigi Iapichino Zentrum fuer Astronomie der Universitaet Heidelberg Institut fuer Theoretische Astrophysik Albert-Ueberle-Str. 2, D-69120 Heidelberg, Germany Tel: +49 6221 548983, Fax: +49 6221 544221 e-mail: luigi@ita.uni-heidelberg.de URL: http://www.ita.uni-heidelberg.de/~luigi/ _______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
-- --------------------------------------------------------------- Luigi Iapichino Zentrum fuer Astronomie der Universitaet Heidelberg Institut fuer Theoretische Astrophysik Albert-Ueberle-Str. 2, D-69120 Heidelberg, Germany Tel: +49 6221 548983, Fax: +49 6221 544221 e-mail: luigi@ita.uni-heidelberg.de URL: http://www.ita.uni-heidelberg.de/~luigi/
participants (2)
-
Britton Smith
-
Luigi Iapichino