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

rr [first-terminus]
?first terminus?

rr [last-terminus]
?last terminus?

md
results off

rr [label-dim]
?temporary dimension of labelling?

; --- Input Parameters ---
; [node-num]        node number
; [label-dim]       minimum width for labeling
; [first-terminus]  first terminus
; [last-terminus]   last terminus
; -------- Input --------
; [particles]       particle template
; [node_doc]        list of particle in each node (modified/labelled)
; ------- Outputs -------
; [node_avg]        output average

; label node doc (at bottom)
if([node-num].eq.1) then
    ; is a copy of the dendrogram
    SD /    PARTICLE_NUM DENDRO_HEIGHT
    [node_doc][node-num]
else
    SD /    PARTICLE_NUM  DENDRO_KEY
    [node_doc][node-num]
endif

sd e  ; close doc
[node_doc][node-num]

; average images
as r
[particles]
[node_doc][node-num]
A                     ; _A_ll
[node_avg][node-num]  ; OUTPUT
_2                    ; dummy file: variance


; CHECK IF AVERAGE NEEDS TO BE RESIZED FOR LABELING

; get image dimension
fi [orig-img-dim]
_2    ; INPUT: variance
(12)  ; header position for x-dim

; if smaller than minimum, interpolate
if([orig-img-dim].lt.[label-dim]) then
    ip
    [node_avg][node-num]
    _1
    [label-dim],[label-dim]
else
    cp
    [node_avg][node-num]
    _1
endif


md
results on

; LABEL

; get number of particles
ud n,[num-parts-node]
[node_doc][node-num]

; get #particles
[num-particles] = [last-terminus] - [first-terminus] + 1
; should be the same as [num-parts-node], but...

; label average w/#particles
la b
_1  ; INPUT: average
_4
{****[node-num]},n={***[num-particles]}
;{***[node-num]},n={*****[num-particles]}

; if average was resized, then re-resize
if([orig-img-dim].lt.[label-dim]) then
    ; get new y-dimension (should be +36)
    fi [labeled-ydim]
    _4
    (2)  ; header position for y-dim

    ; calculate new y-dimension
    [new-ydim] = [labeled-ydim] * [orig-img-dim] / [label-dim] 

    ; interpolate back to original width
    ip
    _4  ; INPUT: labelled average
    [labeled_avg][node-num]
    [orig-img-dim],[new-ydim]
else
    cp
    _4
    [labeled_avg][node-num]
endif

re

; Modified 2009-06-02
;    2009-06-02 (trs) -- improved annotation