source: issm/branches/trunk-jpl-damage/src/m/utils/Exp/readwrite/expread.m@ 11427

function Struct=expread(filename,whole);
%EXPREAD - read a file exp and build a Structure
%
%   This routine reads a file .exp and build a Structure containing the 
%   fields x and y corresponding to the coordinates, one for the filename of
%   the exp file, for the density, for the nodes, and a field closed to 
%   indicate if the domain is closed. 
%   The first argument is the .exp file to be read and the second one (optional) 
%   indicate if the last point shall be read (1 to read it, 0 not to).
%
%   Usage:
%      Struct=expread(filename,whole)
%  
%   Example:
%      Struct=expread('domainoutline.exp')
%      Struct=expread('domainoutline.exp',1)
%
%   See also EXPDOC, EXPWRITEASVERTICES

%some checks
if ~exist(filename),
        error(['expread error message: file ' filename ' not found!']);
end
if nargin<2,
        whole=1;
end

%initialize number of profile
count=0;

%open file
fid=fopen(filename,'r');

%loop over the number of profiles
while (~feof(fid)),

        %update number of profiles
   count=count+1;

   %Get file name
        A=fscanf(fid,'%s %s',2);
        if ~strncmp(A,'##Name:',7), break; end
        if length(A)>7, 
                Struct(count).name=A(8:end);
        else
                Struct(count).name='';
        end

        %Get Icon
        A=fscanf(fid,'%s %s',2);
        if ~strncmp(A,'##Icon:',6), break; end

        %Get Info
        A=fscanf(fid,'%s %s %s %s',4);
        if ~strncmp(A,'#Points',7), break; end

        %Get number of nods and density
   A=fscanf(fid,'%f %f',[1 2]);
   Struct(count).nods=A(1);
   Struct(count).density=A(2);

        %Get Info
        A=fscanf(fid,'%s %s %s %s',5);
        if ~strncmp(A,'#XposYpos',9), break; end

        %Get Coordinates
        A=fscanf(fid,'%f %f',[2 Struct(count).nods]);
        Struct(count).x=A(1,:)';
        Struct(count).y=A(2,:)';

        if(Struct(count).nods~=length(Struct(count).x))error(['Profile ' num2str(count) ' reports incorrect length']); end;

        %Check if closed
        if (Struct(count).nods > 1) && ...
           (Struct(count).x(end) == Struct(count).x(1)) && ...
           (Struct(count).y(end) == Struct(count).y(1))
                Struct(count).closed=true;
                %skip last coordinate if whole=0,
                if whole==0
                        Struct(count).nods=Struct(count).nods-1;
                        Struct(count).x   =Struct(count).x(1:end-1,1);
                        Struct(count).y   =Struct(count).y(1:end-1,1);
                end
        else
                Struct(count).closed=false;
        end

end

%close file
fclose(fid);