function fh = open(fname)
% fh=open(fname) Opens a Gauss data set for reading and stores info in
% (fh)'th row of the global matrix _gaussfh
% USE:  fh = open(fname) where fname is a string
global gaussfh; % gauss file handle info
fnamehd=strrep(fname,'.dat','.dht');
faux=fopen(fnamehd,'r');
if faux==-1;
   'Could not open Gauss dataset header file';
   'The following file was not found or is already open:';
   fnamehd;
end;
[headdat,headlen]=fread(faux,inf);
fclose(faux);
fh=fopen(fname,'r');
if fh==-1;
   'Could not open Gauss dataset ';
   'The following file was not found or is already open:';
   fname;
end;
gaussfh(fh,1)=headdat(7);
gaussfh(fh,2)=headdat(13);
return;