%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%  joint.m - determines the mmf per joint for a transformer
%            core. Calls Hm7.m where B-H curve is in lines/
%            sq.in. & Amp-turns/in. Material is M-7 grain 
%            oriented ESS.
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%clear; clf;  % Disable for use with tranmag.m
BFe=linspace(0,200000, 200); n=length(BFe);
del=0.010;  % Avg joint length (in)
for i=1:n
HFe=Hm7(BFe(i));
FT(i)=HFe*del;
BT(i)=0.5*(BFe(i)/64.52e03+pi*4e-07*HFe*39.37)*64.52;
end
for i=1:n; if BT(i)>=50; m=i; break; end; end
% Suppress plot for use with tranmag.m
%plot(0,0,FT(m:n),BT(m:n)); grid
%title('Joint mmf drop');
%xlabel('Amp-turns per joint'); 
%ylabel('Flux density - kilolines/sq.in.');
BT=1000*BT; % BT calculated in kilolines/in^2