-
slevy authored
Examples of new CM-node-ring-marker and particle tracking. Replace nclump.cmap with one for use with the new nclump codes, where nclump is positive for stars and negative for CM nodes and rings. The new cmap uses the same colors for both, but it could be changed.
slevy authoredExamples of new CM-node-ring-marker and particle tracking. Replace nclump.cmap with one for use with the new nclump codes, where nclump is positive for stars and negative for CM nodes and rings. The new cmap uses the same colors for both, but it could be changed.
data2.cf 1.96 KiB
#! /usr/bin/env partiview
# The above trick searches for "partiview" on $PATH,
# making this .cf file be an executable script:
# invoke with either
# partiview whatever.cf
# or just
# whatever.cf
# These are the data fields filled in by the kira reader as
# compiled into partiview (kira_parti.cc) as of Dec 30 2000.
# If the code changes, these datavar definitions should change too.
# The field names given here are arbitrary, but indices are hard-wired.
# id: >0 positive small integer index for leaf nodes (stars),
# <0 negative of small integer index for non-leaves (CM nodes)
# mass: mass
# nclump: number of stars in clump. 1=singles, 2=binaries, ...
# topnode: small int index, equal for all stars/CM nodes in clump.
# topnode == abs(id) for root node of a clump.
# treeaddr: bit-encoded, giving location in a clump tree.
# 0 for single stars (nclump=1). Otherwise,
# 1 for root of a clump
# 2*(parent treeaddr) for left-child of parent
# 2*(parent treeaddr)+1 for right-child of parent
# ringsize: zero for stars, nonzero world-unit radius for (some) CM nodes
datavar 0 id
datavar 1 mass
datavar 2 nclump
datavar 3 topnode
datavar 4 treeaddr
datavar 5 ringsize
# ``read kira data from file "data2.out"''
kira data2.out
# ``scale particle brightness by mass (over range 0 - 0.01)''
eval lum mass 0 0.01
# ``brighten them up'' Intrinsic brightness = <lum function> * psize * slum
eval psize 100
# use bright-colored colormap, and
# color by number-of-stars-in-clump. colorindex = nclump + (-1)
# so e.g. nclump=1 => colorindex=0
eval cmap nclump.cmap
eval color nclump exact 8
# Turn center-of-mass nodes off,
# draw ring markers for multiple stars,
# choose marker radius from abs(semimajor axis) * 1.5
eval kiractl nodes off
eval kiractl rings on
eval kiractl ringsize semimajor
eval kiractl ringscale 1.5
# Could track a particle, with:
## eval kiractl track 12
# ``load data for time 0 and report to user how many particles are in it''
eval step 0
eval every