; pnums.spi - gets the range of particle numbers associated with
;            each micrograph.
;
;  Output file lists total particles and 1st & last particle numbers.

; -------------- Parameter --------------
[first-part] = 1               ; first particle number

; ----------- Input files -------------- 
FR G
[FILENUMS]../sel_micrograph    ; file numbers
FR G
[sndc]coords/sndc{****x55}     ;  coordinate doc file with particle numbers

; ----------- Output files --------------
FR G
[order]order_picked               ; output doc file 

; -------------- END BATCH HEADER -------------

DE
[order]

; get the number of files
UD n,x20
[FILENUMS]


SD /     MICROGRAPH  # PARTICLES  CUM.TOTAL   FIRST       LAST
[order]

x91= [first-part] - 1 

DO LB1 x11=1,x20

   UD x11,x55
   [FILENUMS]
   ; x55 is now the micrograph number

   UD N,x77      ; number of particles for this micrograph
   [sndc]

   UD 1,x33,x34,x78     ; x78 = first particle number
   [sndc]

   UD x77,x33,x34,x79   ; x79 = last particle number
   [sndc]

   x91=x91+x77        ; cumulative total

   SD x11,x55,x77,x91,x78,x79
   [order]

LB1

EN D

; Modified 2008-09-17
;    2008-09-17 (trs) -- include first particle# as input