test_Mr.py
View Options

			#!/usr/bin/env python


			from pNbody import *
			import sys
			import iclib

			n = 10000

			Rmax = 100.
			rs = arange(0,Rmax,1)
			Mr = rs
			pos = iclib.generic_Mr(n,50,rs.astype(float32),Mr.astype(float32),-1)

			pos = pos.astype(float32)


			nb = Nbody('Mr.dat',ftype='gadget',pos=pos,status='new')
			nb.npart = array([0,0,n,0,0,0])
			nb.npart_tot = nb.npart
			nb.write()

pNbody/src/iclib/test_Mr.py4fbc3430dc96master