([node-num])
;([node-num],[first-terminus],[last-terminus])

md
vb off

rr [first-terminus]
?first terminus?

rr [last-terminus]
?last terminus?

; --- Input Parameters ---
; [node-num]        node number
; [first-terminus]  first terminus
; [last-terminus]   last terminus
; -------- Input --------
; [node_avg]        output average
; ------- Output -------
; [node_lut]        node-to-particle lookup table

if([first-terminus].gt.0) then
    ; get #particles
    [num-particles] = [last-terminus] - [first-terminus] + 1

    ; save to lookup-table
    sd [node-num],[first-terminus],[last-terminus],[num-particles]
    [temp_lut]
else
    [num-particles] = 0

    ; save to lookup-table
    sd [node-num],[first-terminus],[last-terminus],[num-particles]
    [temp_lut]
endif

sd e  ; close doc
[temp_lut]

; increment node#
[node-num]=[node-num]+1

; delete pre-existing files for next node#
de
[node_doc][node-num]
de
[node_avg][node-num]
de
[labeled_avg][node-num]
; Since there are gaps now, DE A will give up after a gap of 10.

md
vb on

re

; Modified 2009-06-02
;    2009-06-02 (trs) -- first particle# of 0 creates dummy entries
;    2009-06-02 (trs) -- improved annotation