Comments and feedback are welcome! (see AUTHOR below)
Input file must be in SNAP205 ASCII or TIPSY ASCII format. Output files are in ASCII.
x,y,z: coordinates of
the subsystem center
radius: cutoff radius for particles to be examined (i.e. the radius if the individual system under study
nmin best value is strongly dependent from the steepness of the density profile, for homogeneous systems is recommanded to use nmin near the total number of particles, otherwise the center of density will be chosen in an almost random way!!
An example paramter file:
king input file name 1 input file name type king root file name 1 4096 first and last 0 multi-system option 256,0.95,3 nmin, facr & radius 128 particles per bin 64 number of resamplings 1234567890 random seed
xxxxrho r r^1/4 Log(r) rho(r) errrho(r) m(r) errm(r) xxxxrms r vrmsr err Vrmsth err Vrmsph err xxxxvel r <vr> err <vtheta> err <vphi> err xxxxsigma r sigmar err sigmat err sigmaph err beta err r = mean radius of particles belonging to each shell rho(r) = density in each shell errrho(r) = 1-sigma error on density of each shell m(r) = cumulative mass function errm(r) = 1-sigma error on m(r) vrmsr = root mean square of radial velocity at radius r vrmsth = root mean square of theta velocity at radius r vrmsph = root mean square of phi velocity at radius r err = 1-sigma error of previous estimated quantity <vr> = mean radial velocity at radius r <vtheta> = mean theta velocity at radius r <vphi> = mean phi velocity at radius r sigmar = standard deviation (sigma) for radial velocity at radius r sigmat = standard deviation (sigma) for theta velocity at radius r sigmaph = standard deviation (sigma) for phi velocity at radius r beta = anisotropy parameter as in Jeans Equation; beta=1-(sigmath**2/sigmar**2)In source file, expressions for alternative formulations of anisotropy parameters are available; since a bootstrap is made, beta is not simply the result of a calculation on other columns in file; ( f(<q>) not always equal to <f(q)>!!).
$NEMO/usr/governato/ (location in NEMO) rhobocop.f the source code rhobocop.doc this document king King model W0=7 N=4096 (SNAP205 format) rbpar parameters for analysis of king
10-dec-1994 V1.1 released version FL&FG
Table of Contents