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@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@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org