function [hpol, fluxo] = hpol_espira(bobines, Rw, Zw)
% Calculo do campo H poloidal provocado por uma espira na coordeenada R=Rw Z= Zw
% nos pontos incluidos da estrutura bobines
% bobines.R 	--coord. R das bobines
% bobines.Z	--coord. R das bobines
% bobines.eRZ	--coord. em R do versor em  (R,Z) da direcção poloidal
% bobines.epZ	--coord. em Z do versor em  (R,Z) da direcção poloidal
% Normalizado a uma corrente = 1 Amp

epR = bobines.epR;
epZ = bobines.epZ;

[hR, hZ, fluxo] = hRZ_espira(bobines.R, bobines.Z, Rw, Zw);
%disp(['hR = ' num2str(hR(1))]);
%disp(['hZ = ' num2str(hZ(1))]);
hpol=hR .* epR + hZ .* epZ; 	% produto interno de campo pelo versor