How to visualize a tipsy file with glnemo2
Convert to gadget2 file format using '''tipsy2gadget''' program
So far glnemo2 does not support natively tipsy file format, therefore you have to convert it in a known file format like gadget2
tipsy2gadget in=??? out=??? scale=100 tipsy2gadget help=h in : tipsy input file [???] out : gadget2 output file) [???] scale : scaling factor (pos,hsml) [100]
example :
tipsy2gadget in=planck.64_64.00047 out=planck.64_64.00047.g2
planck.64_64.00047.g2 is the gadget2 file created
get some info from gadget file using '''uns_info''' program
uns_info planck.64_64.00047.g2 ************************************************** File name : planck.64_64.00047.g2 File type : Gadget2 Nbody selected = 262144 Time=0 ================================================== halo : 262144 mass[1] = 1.211165e-06 1.211165e-06 1.211165e-06 1.211165e-06 1.211165e-06 1.211165e-06 . . . pos [3] = -4.950922e-01 -4.887503e-01 -4.955622e-01 -4.788181e-01 -4.887079e-01 -4.943947e-01 . . . vel [3] = -2.016661e-01 2.457134e-01 -2.518843e-01 -1.523429e-01 2.570927e-01 -1.726928e-01 . . . id [1] = 0 1 2 3 4 5 . . .
We can see that there is only one component here : dark matter halo
- Visualize gadget file with glnemo2
glnemo2 planck.64_64.00047.g2
Compute the density with '''uns_density''' program
uns_density help=h in : UNS input file (gadget|nemo) [???] out : output file [???] select : component selected (disk,stars,halo,gas,range) [???] m : 0: Ferrers method, 1:hackdens method [0] K : number of neighbours [32] N : order of Ferrers kernel [1] ncrit : ncrit value, if empty=max(1,K/4) [] time : selected time [all] first : add a trailing numbering to the first output file [f] offset : +/- time offset [0.01]
uns_density planck.64_64.00047.g2 planck.64_64.00047.g2.dens halo
We compute density on halo component
- visualize with glnemo2