plt_Jel_Jt_phi_fea_int.m

function plt_Jel_Jt_phi_fea_int(result)

% Get data from Result to plot easier
phi_fea = result.fea.Phi;
Jt = result.fea.Jtotal_Avg;
Jel = result.fea.Jel_EPFM;
num_steps = result.fea.num_steps;
J_tear = result.fea.J_tear_phi;
Jel_tear = result.fea.Jel_tear_phi;
%CMOD_tear = result.tear_CMOD;
%tear_CMOD_analysis  = result.fea.tear_CMOD_analysis;
if ~isempty(result.tear_angle)
    phi_tear = result.tear_angle;
else
    phi_tear = 30.0; %set equal to 30 deg. for generic plot without test data
end
%find phi result closest to phi_tear
%[C,I] = min(abs(phi_fea-phi_tear));
%phi_str = num2str(phi_fea(I), '%5.2f');
phi_str = num2str(phi_tear, '%5.2f');

plot(phi_fea,Jel(:,num_steps),'--ob')
ylabel(result.J_str);
xlabel('\phi (deg.)');
%    title('J_{total} and J_{elastic} vs. \phi at max. load - FEA');
hold on

plot(phi_fea,Jt(:,num_steps),'-sb','MarkerFaceColor','b')
str1 = 'J_{el} Final Step';
str2 = 'J_{total} Final Step';
str3 = 'J_{el} @ Tear Point';
str4 = 'J_{total} @ Tear Point';
str5 = strcat('tear location at \phi =', phi_str);
% create line at phi = phi tear
x = [phi_tear phi_tear];
y = [0 max(Jt(:,num_steps))*1.02];


if J_tear(1) ~= 0
    plot(phi_fea,Jel_tear,'--dr')
    plot(phi_fea,J_tear,'-^r','MarkerFaceColor','r')
    plot(x,y,'--k');
    legend_text = {str1, str2, str3, str4, str5};
elseif isempty(result.tear_angle)
    %plot(x,y,'--k');
    legend_text = {str1, str2};
else
    plot(x,y,'--k');
    legend_text = {str1, str2, str5};
end

ylim([0 inf])
legend(legend_text, 'Location', 'best');
hold off

%print figure to emf file and save fig file
%     filename = strcat(FilePrefix, '_Jk');
%     print_file = strcat(pathname, '\', filename);
%     print(gcf,'-dmeta','-r300',print_file);
%     FigFile = strcat(print_file, '.fig');
%     saveas(gcf, FigFile);