function [a,b,c,d,e,f] = read_domain_legacy(file)
  % [a,b,c,d,e,f] = read_domain_legacy(file)
  % Reads domain boundaries from 'current_domain.tmp'
  % Input
  %  file         path to 'current_domain.tmp'
  % Output
  % a,b,c,d,e,f   domain boundaries
  
  border=load(file);
  a=border(1);
  b=border(2);
  c=border(3);
  d=border(4);
  e=border(5);
  f=border(6);
end