%%% Exemples pour un calcule de précession
%
% $Revision: 15 $
% $Author: lenertz $
% $Date: 2015-07-30 17:36:07 +0200 (jeu., 30 juil. 2015) $
% contact: lenertz@cea.fr


%% Déclaration des paramètres pour une modification plus simple de l'exemple
% voir doc de computePolarization pour plus d'informations sur les
% arguments

% paramètres optionnelles
<<<<<<< .mine

H = [0.00,0,0.008]; % vecteur Hxyz en T
M = 0.6; % en T (peut aussi être un vecteur dans une direction différente de H)
MagMod = @wireVortexField_vector_V2; % déclaration de la fonction calculant le champ magnétique

=======
ang = 0.1;
H = 1e-4*[50, 000, 1]; % vecteur Hxyz en T
M =  0.6*[sind(ang), 0, cosd(ang)];%0.6*[sind(ang)*0.6,0, cosd(ang)]; % en T (peut aussi être un vecteur dans une direction différente de H)
%MagMod = @wireVortexField_vector_V4; % déclaration de la fonction calculant le champ magnétique
MagMod = @wirePlanField_vector;
>>>>>>> .r16
% Paramètres optionnelles
longueur = 200e-6; %en m /!\ attention à choisir en accord avec MagMod
x_step = 0.1e-6; % en m
largeur = 350e-6; % en m
y_step = 2e-6; % en m
r = 100e-6; % en m
wavelength = 4e-10; % en m
wl_sigma = wavelength*0.07; % en m
wl_nbStep = 3; % en m
wl_mesure = 1; % en unité de wl_sigma
beam_width = 30e-6; % en m 
BO_save = false; % booléen
<<<<<<< .mine
startPolar = 0.9; % entre 1 et -1
    %angPhase = 0:2*pi/10+0.001:2*pi; % en rad
    angPhase = 0; % en rad
=======
supMagParams = [0.95*r, 0.1, 0.9];
>>>>>>> .r16

startPolar = 1; % entre 1 et -1
angPhase = 0.00:2*pi/5:2*pi-0.001;%0:2*pi/6+0.001:2*pi; % en rad


Params = [longueur, x_step, largeur, y_step, r, wavelength, wl_sigma,...
          wl_nbStep, wl_mesure, beam_width, BO_save];


%% Execution
% Calcul
<<<<<<< .mine
p = computePolarization(H,M,MagMod, Params, startPolar, angPhase,1,2);
=======
p = computePolarization(H,M,MagMod, Params, supMagParams, startPolar, angPhase);
>>>>>>> .r16



% % Affichage du résultat en polarisation sur une même figure pour plusieurs
% % largeurs de faisceau:
% fig_polar = p.plotPolar(10e-6);
% p.plotPolar(20e-6, fig_polar);
% p.plotPolar(40e-6, fig_polar);
% p.plotPolar(100e-6, fig_polar);
