gspdisk

gspdisk                           Make exponential disk in a gsp halo

  gsp=                            Input gsp for halo mass profile

  out=                            Output N-body model of disk

  mdisk=0.1875                    Total disk mass

  alpha=12.0                      Inverse exponential scale length

  z0=0.01                         Vertical scale height

  epsilon=-1                      Plummer smoothing parameter.
				  No smoothing (cf. v1.1) if < 0.

  mu=2.0                          Ratio of sig_r/sig_z.
				  Constant unless r_mu > 0.

  r_mu=-1                         Scale radius for mu(R) function.
				  If > 0 then mu -> 1 as R -> 0.

  eta=-1                          Velocity distribution parameter.
				  Use Gaussian distrib. if <= 0.

  rcut=1.0                        Outer disk cutoff radius

  ndisk=12288                     Number of disk particles

  seed=54321                      Seed for random number generator

  zerocm=true                     If TRUE, translate CM to origin

  VERSION=1.2                     Josh Barnes  27 March 2002

Joshua E. Barnes      (barnes at ifa.hawaii.edu) Created: 9 June 2008 http://www.ifa.hawaii.edu/~barnes/zeno/gspdisk.html