<div dir="ltr"><div><div><div><div><div><div>Hello Everyone!<br><br></div>So, I have a question about how does covering grid actually works. I have a Flash simulation output of a 28*28 Mpc box containing a cluster. I am trying to simulate the spectra from a region of size (0.08 Mpc x, 0.08 Mpc y, 28 Mpc z) around the centre, (I will assume z is my line of sight) so I use this line to define the region:<br>
<br>#The shape of my region <br>cube = pf.h.covering_grid(level_max, left_edge=newlep, right_edge=newrep,dims=[6, 6, 128]) <br><br></div>where level = 5 (max level of refinement for my Flash data), "newlep" is centre coordinates (Mpc)  minus 0.08 Mpc and newrep is centre coordinates plus 0.08 Mpc, and then into percentage units.<br>
<br> I put them into an array of dimensions [6,6,128] because I want 36 spectra (one for each pixel of the ccd) so I will fix an (x,y) position and write an XSPEC file for all z elements. (I need basically temperature and density of all the elements).<br>
<br></div>My question is: how does covering grid makes a single field value (temperature, for example)  inside the (x,y,z) position? because I am arbitrarily saying to divide the region into 6 pixel, but the data file might have  many simulation pixels inside , so how are the field values calculated? is it just an average? or are just certain points taken into consideration?<br>
<br></div>Thanks for any hints/tips/help!<br></div>Barbara Ramirez<br></div>University of Innsbruck, Austria<br></div>