Up: No Title
Previous Page: Future
Next Page: References
To illustrates some of the features in NEMO we show the following C-shell script that analyzes the results of a contrived encounter between two spherical Plummer models:
# make two plummer spheres, the second one not in virial equilibrium mkplummer out=p1 nbody=1000 mkplummer out=- nbody=1000 | snapvirial in=- out=p2 'virial=1/sqrt(3)' snapstack in1=p1 in2=p2 out=p.in deltar=10,0,0 deltav=-1.5,0,0 # integrate hackcode1 in=p.in out=p.out tstop=10 options=phi,acc freqout=1 > p.log # analyze last 3 dynamical times foreach t (8 9 10) snaptrim in=p.out out=- times=$t | snapcenter in=- out=- | unbind in=- out=- | snaprect in=- out=p.out.$t weight="-phi*phi*phi" > p.rect.$t # freeze potential, and classify 50 orbits in the rectified potential foreach i (`nemoinp 1:1000:20`) stoo in=p.out.$t out=- ibody=$i | orbint in=- out=- nsteps=10000 dt=0.01 potname=hackcode1 potfile=p.out.$t | orbname in=- >> p.orb.$t end end
Acknowledgements: NEMO was initially developed by Joshua Barnes, Piet Hut and myself during the 1986-87 academic year in Princeton, I wish to thank them for the stimulating environment in those initial years. I would also like to thank the Monbusho International Scientific Research Program : Joint Research 03044032, for a visit to Tokyo, where some of the work described here has been done.