[yt-users] PHOX mock observation errors
John ZuHone
jzuhone at gmail.com
Tue Sep 15 09:14:10 PDT 2015
Hi Barbara,
Sorry to hear you’re having trouble at the moment. I have a few other questions.
1) Can you post the entire script? Use http://paste.yt-project.org and send a link.
2) Can you give me a sense of what the mass and temperature of your cluster is? It sounds like it is fairly large. And the radius?
3) What code did the simulation come from?
After knowing these, I think we can get somewhere. I have some ideas.
Best,
John
> On Sep 15, 2015, at 10:59 AM, Barbara Ramirez <Barbara.Ramirez-Mosquera at uibk.ac.at> wrote:
>
> Hello!
>
> I have written before about trying to implement the xray photon generator.
> This time I've updated all the libraries, YT, Heasoft... etc.
>
> I'm just trying to generate a photon list for the moment but I'm getting two sorts of errors:
> 1)
>
> ---------------------------------------------------------------------------
> RuntimeError Traceback (most recent call last)
> /home/barbara/hires_final/photon_trial.py in <module>()
> 90 cosmo = Cosmology()
> 91
> ---> 92 photons = PhotonList.from_scratch(Roi, redshift, A, exp_time,thermal_model,center="center",cosmology=cosmo)
> 93 photons.write_h5_file("my_photons.h5")
> 94
>
> /home/barbara/yt-x86_64/src/yt-hg/yt/analysis_modules/photon_simulator/photon_simulator.pyc in from_scratch(cls, data_source, redshift, area, exp_time, photon_model, parameters, center, dist, cosmology)
> 329 parameters["Width"] = 2.*width.in_units("kpc")
> 330
> --> 331 photons = photon_model(data_source, parameters)
> 332
> 333 mylog.info("Finished generating photons.")
>
> /home/barbara/yt-x86_64/src/yt-hg/yt/analysis_modules/photon_simulator/photon_models.pyc in __call__(self, data_source, parameters)
> 188 raise RuntimeError("Number of photons generated for this chunk "+
> 189 "exceeds photons_per_chunk (%d)! " % self.photons_per_chunk +
> --> 190 "Increase photons_per_chunk!")
> 191
> 192 if self.method == "invert_cdf":
>
> RuntimeError: Number of photons generated for this chunk exceeds photons_per_chunk (1000000000)! Increase photons_per_chunk!
>
>
> Ok, so I increase the number of photons per chunk and I get:
> 2)
> ---------------------------------------------------------------------------
> MemoryError Traceback (most recent call last)
> /home/barbara/hires_final/photon_trial.py in <module>()
> 90 cosmo = Cosmology()
> 91
> ---> 92 photons = PhotonList.from_scratch(Roi, redshift, A, exp_time,thermal_model,center="center",cosmology=cosmo)
> 93 photons.write_h5_file("my_photons.h5")
> 94
>
> /home/barbara/yt-x86_64/src/yt-hg/yt/analysis_modules/photon_simulator/photon_simulator.pyc in from_scratch(cls, data_source, redshift, area, exp_time, photon_model, parameters, center, dist, cosmology)
> 329 parameters["Width"] = 2.*width.in_units("kpc")
> 330
> --> 331 photons = photon_model(data_source, parameters)
> 332
> 333 mylog.info("Finished generating photons.")
>
> /home/barbara/yt-x86_64/src/yt-hg/yt/analysis_modules/photon_simulator/photon_models.pyc in __call__(self, data_source, parameters)
> 156
> 157 number_of_photons = np.zeros(num_cells, dtype="uint64")
> --> 158 energies = np.zeros(self.photons_per_chunk)
> 159
> 160 start_e = 0
>
> MemoryError:
>
>
> This is just a trial, so I'm using almost a copy of the code that's on the documentation:
>
> Roi = ds.sphere("center",radius=width_half)
> apec_model = TableApecModel("/home/barbara/yt-x86_64/atomdb_v2.0.2",0.5,10.0,100,thermal_broad=False,apec_vers="2.0.2")
> thermal_model = ThermalPhotonModel(apec_model, X_H=0.75,photons_per_chunk=1000000000, Zmet=0.3,method="invert_cdf")
>
> A = 6000.
> exp_time = 100.0e3
> redshift = 0.05
> cosmo = Cosmology()
>
> photons = PhotonList.from_scratch(Roi, redshift, A, exp_time,thermal_model,center="center",cosmology=cosmo)
> photons.write_h5_file("my_photons.h5")
>
> My snapshot information:
>
> length unit 0.999752515846 cm
> time unit 1.0 s
> mass unit 1.0 g
> velocity unit 1.0 cm/s
> Enter total Width of plot in Mpccm: 0.5
> smallest cell: 14.6469091426 kpccm/h
> minima and maxima 1.04339659244e-29 g/cm**3 1.5187912088e-24 g/cm**3
> The simulation is [ 89.40546647 89.40546647 89.40546647] Mpccm
> Centering region in:, [ 44.58268195 44.53902694 44.51719943] Mpccm
> yt : [INFO ] 2015-09-15 16:49:43,219 Cosmology: h = 0.71, omega_matter = 0.27, omega_lambda = 0.73
>
>
> Am I just using wrong parameters to try to simulate the Xray emission of this simulated cluster?
>
> Thanks beforehand for the help!
>
> --
> Barbara L. Ramirez MSc.
> Institut für Astro- und Teilchenphysik
> University of Innsbruck
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20150915/cc231076/attachment.html>
More information about the yt-users
mailing list