FUNCTION STRPOSMULTI, s, subs, cpt
;+ 
; NAME: STRPOSMULTI
; PURPOSE: 
;   extract substring successive positions within a string, using strpos
; CATEGORY: I-4-a
; CALLING SEQUENCE: 
;   index=STRPOSMULTI(s, subs [, cpt] )
; INPUTS: 
;   s      -- string : initial string
;   subs   -- string : substring, if '', return 0 and cpt eq 1
; OPTIONAL INPUT PARAMETERS: 
;   none
; KEYED INPUTS: 
;   none
; OUTPUTS: 
;   index  -- long array : positions (if none -1)
; OPTIONAL OUTPUT PARAMETERS: 
;   cpt    -- long : number of psoitions
; EXAMPLE: 
; ALGORITHM: 
;   straightforward loop using strpos
; DEPENDENCIES: 
;   none
; COMMON BLOCKS: 
;    SESSION_BLOCK, SESSION_MODE, ERROR_CURRENT, STATUS_BOOL 
; SIDE EFFECTS: 
;   none
; RESTRICTIONS: 
;       UNTESTED
; CALLED PROCEDURES AND FUNCTIONS: 
;   none
; MODIFICATION HISTORY: 
;    2-Nov-1995  written with template_gen                  FV IAS
;-
 
;------------------------------------------------------------
; common blocks 
;------------------------------------------------------------
 
; environment parameters 
 COMMON SESSION_BLOCK, SESSION_MODE, ERROR_CURRENT, STATUS_BOOL
 
;------------------------------------------------------------
; on error conditions
;------------------------------------------------------------
 ON_ERROR,  ERROR_CURRENT
 
;------------------------------------------------------------
; initialization
;------------------------------------------------------------
 
 ROUTINE_NAME = 'STRPOSMULTI'
 VERSION = '1.0' 
 CATEGORY = 'I-4-a'
 STATUS = ['SUCCESS', 'S', ROUTINE_NAME+ ' V.' + VERSION, CATEGORY]
 VAR_NAMES= ['s', 'subs', 'cpt', 'output']
  s_s = CONV_STRING(s)
  s_subs = CONV_STRING(subs)
  s_cpt = CONV_STRING(cpt)
 CALL_VAL = [s_s, s_subs, s_cpt, '']
 output= -1
 
;------------------------------------------------------------
; parameters check
;------------------------------------------------------------
 
 IF N_PARAMS() LT 2 THEN BEGIN
   PRINT, 'CALLING SEQUENCE: ', $ 
    'output=STRPOSMULTI(s, subs  [, cpt] )'
   STATUS(0) = ['PARAMETER MISSING', 'E']
   GOTO, CLOSING
 ENDIF

 IF subs eq '' THEN BEGIN & output=[0l] & cpt=0 & GOTO, CLOSING & END
 
;------------------------------------------------------------
; function body
;------------------------------------------------------------

 pos = strpos(s, subs)
 WHILE pos ne -1 DO BEGIN
   output = [output, pos]
   pos = strpos(s, subs, pos+1)
 ENDWHILE

 index = where(output ne -1, cpt)
 IF cpt GT 0 THEN output=output(index) ELSE output = -1
  
;------------------------------------------------------------
; closing
;------------------------------------------------------------
 
 CLOSING:
  s_s = CONV_STRING(s)
  s_subs = CONV_STRING(subs)
  s_cpt = CONV_STRING(cpt)
  s_output = CONV_STRING(output)
  ACTL_VAL = [s_s, s_subs, s_cpt, s_output]
 
  IF (STRMID(STATUS(1),0,1) NE 'S') THEN STATUS_BOOL=1 
  RECORD_LOGFILE, STATUS, CALL_VAL, VAR_NAMES, ACTL_VAL
 
  RETURN, output
 
 END