function read_net2(file)
%     function read_net2(file)
%
% Autor: Alexander Pruell
% Datum: 08.06.00
%
% liest aus der Datei file die Netzdaten und schreibt diese
% in globale Variablen:
%
%    N_nodes  ... Anzahl der Knoten
%    node(i)  ... Datensruktur mit:
%      node(i).koo ... Koordinaten des i-ten Knotens
%      node(i).T0  ... Anfangstemperatur am i-ten Knoten
%
%    N_ele    ... Anzahl der Vierecke
%    ele(i,:) ... 4d-Vektor mit den Knotennummern des
%                 i-ten Vierecks
%
%    N_neum   ... Anzahl der Randstuecke mit Neumann-
%                 Randbedingungen
%    neum(i,:)... 2d-Vektor mit den Knotennummern des
%                 i-ten Randstuecks mit Neumann-
%                 Randbedingungen
%
%    N_diri   ... Anzahl der Knoten mit Dirichlet-
%                 Randbedingungen
%    diri(i)  ... Knotennummer des i-ten Knotens mit 
%                 Dirichlet-Randbedingungen


global N_nodes N_ele N_neum N_diri node ele neum diri


fid = fopen(file);

N_nodes = fscanf(fid,'%f',1);
N_ele   = fscanf(fid,'%f',1);
N_neum  = fscanf(fid,'%f',1);
N_diri  = 0;

for i=1:N_nodes
  node(i).koo       = fscanf(fid,'%f',2);
  dirichlet         = fscanf(fid,'%s',1)=='j';
  node(i).T0        = fscanf(fid,'%f',1);
  if dirichlet
    N_diri          = N_diri + 1;
    diri(N_diri)    = i;
  end
end

ele = zeros(N_ele,4);
for i=1:N_ele
  ele(i,:)          = fscanf(fid,'%f',4)';
end

neum = zeros(N_neum,2);
for i=1:N_neum
  neum(i,:)         = fscanf(fid,'%f',2)';
end

fclose(fid);
