function y=lddat(nam)
  d = load(nam);
  y.x = d(:,1);
  y.dat = d;
  y.y = d(:,3) ./ d(:,5);
  y.shift = d(:,8);
  y.pol = y.shift(1)/abs(y.shift(1));
  y.B = 0.0;
  y.sample = "";
end