%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%  spcore.m - determines d,h & l for core transformer given N
%             and h/w ratio. See data1cor.m for req'd inputs.
%             Also, calculates approx magnetizing current and
%             weights of conductor and core material.
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear;
global Ggammam GSF Gkl GV1 GI1 Gf Gdelc Grw Grhoc 
global GBm GHm Gpm GPm Gkc GN
data1cor,
x0=[GV1*GI1/1000 5 2]'; % Solution is x0 dependent.
delx=zeros(3,1); f=delx; x=x0;
% Newton-Raphson solution, x1=d x2=w x3=l
for i=1:100
   [f,J]=core(x0); delx=-inv(J)*f; x=x0+delx;
      if all(abs(delx)<1e-12);break; else; end
      x0=x;
end
l=x(3); d=x(1); w=x(2); h=Grw*x(2); N=GN; 
disp(' '); disp(' ');
disp('   CORE SIZING for CORE-TYPE TRANSFORMER');
if any(x<0); disp('EXTRANEOUS SOLUTION'); pause; else; end
format compact
disp(' '); disp(' ')
disp(['    l (in)' '    d (in)' '    w (in)  '...
      '  h (in) ' '  N (turns)  ']);
disp([ l d w h N ]);
% Magnetizing current - A & pu ( Joint/Cross grain neglected )
Im=2*GHm*(h+w+l)/N/sqrt(2); puIm=Im/GI1; 
disp('  '); disp('  ') 
disp(['    Im (A)   ' '  puIm  ']);
disp([ Im puIm ]);
% Conductor & core weight
Wfe=2*Ggammam*(2*l+(Grw+1)*w)*l*d*GSF;
Wc=2*Gkc*0.322*h*w*(l+d+2*w);
disp('  '); disp('  ')
disp(['   Wfe (lb) ' ' Wc (lb)     ']);
disp([ Wfe Wc ]);