[hd]=fea_yz

function [hd]=fea_yz
global zeng2
if isempty(zeng2.FEA)
    hd=[];
    zexst('err','No FEA file');

else
    hd=figure(1000);
    delete(get(hd,'children'));
    axes;
    for i=1:zeng2.FEA.number_sites
        patch(zeng2.FEA.yz(i,1),zeng2.FEA.yz(i,2),'b',...
            'marker','o',...
            'LineWidth',1.25,...   
            'MarkerEdgeColor','b',...
            'MarkerFaceColor','none',...
            'MarkerSize',8,...         %'tag','patch axon');
            'tag','axon positions');
    end
end