fea_import

function fea_import

[Filename Path]=uigetfile(['*.fea']);
if strcmp('0',num2str(Filename)) | strcmp('0',num2str(Path))
   return
else
    global zeng2 
    zeng2.FEA.xyzv=dlmread([Path Filename]);
    %Find yz---------------------------------------------
    minmax=[min(zeng2.FEA.xyzv(:,1)) max(zeng2.FEA.xyzv(:,1))];
    xmax=find(zeng2.FEA.xyzv(:,1)==minmax(2));
    xmin=find(zeng2.FEA.xyzv(:,1)==minmax(1));
    xmean=mean(minmax);
    temp=min(abs(zeng2.FEA.xyzv(:,1)-xmean));
    xmean=find(abs(zeng2.FEA.xyzv(:,1)-xmean)==temp(1));
    if length(xmin)==length(xmax)==length(xmean)
        mychoice=xmin;
    else
        lengthmax=max([length(xmin) length(xmax) length(xmean)]);
        if length(xmin)==lengthmax
            mychoice=xmin;
        elseif length(xmax)==lengthmax
            mychoice=xmax;
        else
            mychoice=xmean;
        end
    end
    %/Find yz---------------------------------------------
    zeng2.FEA.yz=zeng2.FEA.xyzv(mychoice,[2 3]);
    zeng2.FEA.number_sites=length(zeng2.FEA.yz(:,1));
    zeng2.FEA.xminmax=[min(zeng2.FEA.xyzv(:,1)) max(zeng2.FEA.xyzv(:,1))];
    hd=fea_yz;
end