FUNCTION ra_stimpos, im stimindex = findgen(40) + 1 stim = float(im[1:40, 22]) stimpos = total(stim * stimindex)/total(stim) return, stimpos END