#!/usr/bin/octave
# load data
if (length(who("calibdata")) == 0) 
  calibdata = load("calibrate-50a.log");
endif

# split it into each keithley value
if (length(who("a")) == 0)
  for i = 1 : (length(calibdata) / 4000)
    block=1:4000;
    a(i, block, 1) = calibdata((i-1) * 4000 + block, 2);
    a(i, block, 2) = calibdata((i-1) * 4000 + block, 3);
    a(i, block, 3) = calibdata((i-1) * 4000 + block, 1);
  endfor
endif

# plot them all
#clf;
#figure(1);
#hold on

out = [];
for i = 1:2:21
  range = 1:4000;
  x = a(i, range, 1);
  y = a(i, range, 2);
  out = [ out x' y' ];
endfor

save("-ascii","calib.dat","out");