example_mkmap2dn.py
6c0821250383master

example_mkmap2dn.py
View Options


			from pNbody import *
			from pNbody import io
			import mapping

			nb = Nbody('../../examples/snap.dat',ftype='gadget')



			mx = max(nb.rxyz())

			kx = 256
			ky = 256
			kz = 256

			pos = ( nb.pos + [mx,mx,mx] )/(2*mx) # between 0,1
			#pos = pos*[kx,ky,kz]
			pos = pos.astype(float32)

			mat = mapping.mkmap2dn(pos,nb.mass,nb.mass,(kx,ky))



			#io.WriteFits(mat,'qq.fits')
			mplot(mat)