Skip to content
Snippets Groups Projects
Commit 98955481 authored by slevy's avatar slevy
Browse files

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.
parent 45248e49
No related branches found
No related tags found
No related merge requests found
......@@ -8,7 +8,7 @@
# whatever.cf
# These are the data fields filled in by the kira reader as
# compiled into partiview (kira_parti.cc) as of Dec 21 2000.
# 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.
......@@ -23,6 +23,7 @@
# 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
......@@ -34,9 +35,6 @@ datavar 5 ringsize
# ``read kira data from file "data2.out"''
kira data2.out
eval step 0
eval every
# ``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
......@@ -46,7 +44,7 @@ eval psize 100
# 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 -1
eval color nclump exact 8
# Turn center-of-mass nodes off,
# draw ring markers for multiple stars,
......@@ -55,6 +53,8 @@ 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
......
8
0.740741 0.709877 0.370370 0.988426
0.319273 0.907407 0.384621 0.987269
0.976852 0.085927 0.407650 0.984375
0.976852 0.896955 0.018090 0.997685
# This colormap is expected to be used as in
# eval cmap nclump.cmap
# eval color nclump exact 8
17
0.918210 0.972222 0.972222 1.000000 # (#0) CM nodes/rings for systems with >=8 stars
0.913446 0.114133 0.986111 1.000000
0.972415 0.493056 0.986111 1.000000
0.013439 0.967593 0.888080 0.995370 # ...
0.976852 0.896955 0.018090 0.997685 # (#4) CM nodes/rings for quadruple systems
0.976852 0.085927 0.407650 0.984375 # (#5) CM nodes/rings for triple systems
0.319273 0.907407 0.384621 0.987269 # (#6) CM nodes/rings for double systems
0.740741 0.709877 0.370370 0.988426 # (#7) unused (single CM nodes)
1 1 1 1 # (#8) unused
0.740741 0.709877 0.370370 0.988426 # (#9) single stars
0.319273 0.907407 0.384621 0.987269 # (#10) doubles
0.976852 0.085927 0.407650 0.984375 # (#11) triples
0.976852 0.896955 0.018090 0.997685 # ...
0.013439 0.967593 0.888080 0.995370
0.972415 0.493056 0.986111 1.000000
0.913446 0.114133 0.986111 1.000000
0.918210 0.972222 0.972222 1.000000
0.918210 0.972222 0.972222 1.000000 # (#16) stars for systems with >=8 stars
#! /usr/bin/env partiview
filepath +:.:../data
# 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 toplevelnode
datavar 4 treeaddr # (0=single, 1=root)
datavar 5 ringsize
kira primbin16.out
eval lum mass 0 0.01
eval psize 100
## eval cment 1 1 .7 .3
## eval color nclump exact
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 instantaneous separation * 1.5
eval kiractl nodes off
eval kiractl rings on
eval kiractl ringsize sep
eval kiractl ringscale 1.5
## Could track a particle, with:
## eval kiractl track 12
eval every
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment