%**************************************************************************
%This function plots the multipole field strengths of int_bx_dz and int_by_dz.
%
%Input:
%plt_file, plot file for results
%al_n, skew multipole field strength at the reference radius (Tm)
%bl_n, normal multipole field strength at the reference radius (Tm)
%r0, reference radius (m)
%
%Zachary Wolf, 10/5/10

function [] = magic_fingers_plot_har_intbxydz(plt_file, al_n, bl_n, r0)
  
%Check
if length(al_n) ~= length(bl_n)
    return;
end
if length(bl_n) == 1
    return;
end

%Initialize
figure;

%Plot the field harmonics
nn = 1 : length(bl_n);
%stem(nn, ann, 'linewidth', 0.5, 'markerfacecolor', 'red', 'marker', 'o', 'markersize', '4');
semilogy(nn, abs(al_n), 'r*', 'linewidth', 0.5);
hold on;
%stem(nn, bn, 'linewidth', 0.5, 'markerfacecolor', 'blue', 'marker', 's', 'markersize', '4');
semilogy(nn, abs(bl_n), 'b*', 'linewidth', 0.5);

%Labels
legend('|ALn|', '|BLn|');
title(['Multipole Field Strengths At R0 = ', num2str(r0), ' (m)']);
xlabel('Harmonic Number');
ylabel('Absolute Value Of ALn, BLn (Tm)');

%Finalize
grid on;
zoom on;
hold off;

%Print
orient landscape;
print('-append', '-dpsc', plt_file);

%Done
