(x80,x61,x70,x28)

; --- Input Parameters ---
; x80          node number
; x61          first terminus
; x70          last terminus
; x28          minimum width for labeling
; -------- Input --------
; [particles]  
; ------- Outputs -------
; [node_doc]   
; [node_avg]   

; label node doc (at bottom)
SD /     PARTICLE#    DENDRO_KEY
[node_doc]
sd e
[node_doc]

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


; CHECK IF AVERAGE NEEDS TO BE RESIZED FOR LABELING

; get image dimension
fi x65
_2    ; INPUT: variance
(12)  ; header position for x-dim

; if smaller than minimum, interpolate
if(x65.lt.x28) then
    ip
    [node_avg]
    _1
    x28,x28  ; new dimensions
else
    cp
    [node_avg]
    _1
endif


; LABEL

; get number of particles
ud n,x66
[node_doc]

; get #particles
x60=x70-x61+1
; should be the same as x66, but...

; label average w/#particles
la b
_1  ; INPUT: average
_4  ; WAS [node_avg]
{***x80},n={*****x60}

; if average was resized, then re-resize
if(x65.lt.x28) then
    ; get new y-dimension (should be +36)
    fi x22
    _4
    (2)  ; header position for y-dim

    ; calculate new y-dimension
    x51=x65/x28*x22
    ; x65==original x-dim, x28==new x-dim, x22==interpolated y-dim

    ip
    _4  ; INPUT: labelled average
    [labeled_avg]
    x65,x51
else
    cp
    _4
    [labeled_avg]
endif

re