!
! Compute mean of subscans
!
pen 1 /wei 4
pen 2 /wei 4
def int jj j1 j2 k n
let k &1
def real sum
let jj  0
def rea a[2] b[2] dummy[nrecord,nchan]

define inte num0 num1 num2
define real d e 
let num0 1
def integer c[4] ic
if (onoff_mode.eq.1) then
  let c 1 1 2 2
else
  let c 2 1 2 1
endif
def character tel*3
let tel 'telescope'
!
! Pico Veleta
! -----------
if (tel.ne."CSO") then
 for i 1 to nsubscan
  let num1 nsubrec[i]
  let num2 num1+num0-1
  define real dx[num1] ey[num1]
  let dx signal[num0:num2,k]
  let ey blanking[1] /where dx.eq.blanking[1]
  compute d mean dx
  compute e mean ey
  del /var dx ey
  let d d-e
  if (e.eq.blanking[1]) then
    let d blanking[1]
  endif
  let ic i-int(i|4)*4+1
  pen c[ic]
  draw r lst[num0] d /user
  draw l lst[num2] d /user
  let num0 num2+1 
 next
else 
! CSO
! ---
! First level ON
  def real lst1 lst2
  def int ii
  pen 1
  let lst1 lst[1]-(lst[2]-lst[1])|2.
  let lst2 lst[1]
  let d signal[1,k]
  draw r lst1 d /user
  draw l lst2 d /user
! Following "OFF OFF ON ON" levels
  for i 2 to nrecord-1 by 2
    let ic i-int(i|4)*4+1
    pen c[ic]
    let ii i+1
    let lst1 lst[i]
    let lst2 lst[ii]
    let d (signal[i,k]+signal[ii,k])|2.
    draw r lst1 d /user
    draw l lst2 d /user
  next
! Last level ON
  pen 1
  let lst1 lst[nrecord]
  let ii nrecord-1
  let lst2 lst1+(lst[nrecord]-lst[ii])|2.
  let d signal[nrecord,k]
  draw r lst1 d /user
  draw l lst2 d /user
endif
pen 2 /def
pen 1 /def
pen 0