function bobines = bobines_isttok_Port03_p15d0
%BOBINES = BOBINES_ISTTOK
%Calculo da posi??o das bobines no ISTTOK
% .N   --  numero de bobines de teste
% .RM  --  Raio Maior (R do centro da camara)
% .rb  --  Raio das bobines contado a partir de R=Rm
% .Rp  --  Raio-- das coordenadas toroidais
% .R   --  Coord. R das bobines
% .Z   --  Coord. Z das bobines
% .epR --  Versores da direc??o poloidal em cada bobine no sentido directo de Iplasma
% .epZ
% .sinh_eta -- Coordenadas Toroidas da bobines sinh(eta)
% .cosh_eta =  RM/rb -- "aspect ratio" 
% .vartheta 
% .labels_ab -- rotulos da bobines 

% 

N=12;				%  N -- numero de bobines de teste
RM=.46; 			%  RM-- Raio maior - Centro do vessel 
rb= .187/2;			%  rb-- raio das bobines
a = 0.085;          %  a -- raio menor
Rp   = sqrt(RM^2-rb^2);		% 

hn = floor(N/2);
n=0:N-1;
theta0 = 15.0;
% modificado pela posi??o inicial das bobines
R  = rb * cos(theta0 / 180 * pi + n/6 * pi) + RM; 	% ordem b1, b2 ...b6, a6,  a2, a1
Z  = rb * sin(theta0 / 180 * pi + n/6 * pi);
epR =  sin(theta0 / 180 * pi - n*pi/6);
epZ = -cos(theta0 / 180 * pi - n*pi/6);

cosh_eta_b = RM/rb;		%  a_rt-- "aspect ratio" =  cosh(eta_b)
sinh_eta_b=sqrt(cosh_eta_b^2 -  1);
%lame = R / sinh_eta_b ;
sin_varth_b = sinh_eta_b * Z ./R;
cos_varth_b = cosh_eta_b - Rp * sinh_eta_b ./R ; % = - Rp/lame + cosh(eta)
varth_b     = atan2(sin_varth_b, cos_varth_b);


bobines.N   = N;
bobines.RM  = RM;
bobines.rb  = rb;
bobines.a   = a;
bobines.Rp  = Rp;
bobines.R   = R;
bobines.R2   = R .^2;
bobines.Z   = Z;
bobines.epR = epR;
bobines.epZ = epZ;
bobines.sinh_eta =  sinh_eta_b;
bobines.cosh_eta =  cosh_eta_b;
bobines.vartheta =  varth_b;
bobines.labels_ab= 'b1|b2|b3|b4|b5|b6|b7|b8|b9|b10|b11|b12';

% plot(bobines.R,bobines.Z)
% text(bobines.R,bobines.Z,bobines.labels_ab)