Previous: Future
Up: No Title
Next: References
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 deltar=10,0,0 deltav=-1.5,0,0
  # integrate 
  hackcode1 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

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.