Hi, I have a problem when using yt to obtain the X-ray luminosity of a simulated region. I found that yt 4.3 and yt 4.4 give very different results, i.e., about 5 vs. 10 (* 1e39 erg/s). Do you know what causes this difference? Here is the code. Thanks! import yt ds = yt.load('snapshot_155.hdf5') xray_fields = yt.add_xray_emissivity_field(ds, 0.5, 2.0, table_type='apec', metallicity=1.0) sp = ds.sphere("c", (100.0, "kpc")) shpere= sp.quantities.total_quantity("xray_luminosity_0.5_2.0_keV") print(shpere) Best, Dr. Xiaoxia Zhang (张小霞) Department of Astronomy Xiamen University
Hi, What is the code used by snapshot_155.hdf5? And can you paste the result of the log output when running yt.load for yt 4.3 and yt 4.4? Best, John
On Dec 10, 2023, at 9:10 PM, 张小霞 <zhangxx@xmu.edu.cn> wrote:
Hi,
I have a problem when using yt to obtain the X-ray luminosity of a simulated region. I found that yt 4.3 and yt 4.4 give very different results, i.e., about 5 vs. 10 (* 1e39 erg/s). Do you know what causes this difference? Here is the code. Thanks!
import yt ds = yt.load('snapshot_155.hdf5') xray_fields = yt.add_xray_emissivity_field(ds, 0.5, 2.0, table_type='apec', metallicity=1.0) sp = ds.sphere("c", (100.0, "kpc"))
shpere= sp.quantities.total_quantity("xray_luminosity_0.5_2.0_keV")
print(shpere)
Best, Dr. Xiaoxia Zhang (张小霞) Department of Astronomy Xiamen University
_______________________________________________ yt-users mailing list -- yt-users@python.org <mailto:yt-users@python.org> To unsubscribe send an email to yt-users-leave@python.org <mailto:yt-users-leave@python.org> https://mail.python.org/mailman3/lists/yt-users.python.org/ Member address: jzuhone@gmail.com <mailto:jzuhone@gmail.com>
participants (2)
-
John ZuHone -
张小霞