deltacalmj3

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
 [delta,ddelta_dR,ddelta_dz,rho]=delta_fit(R,z);
 Approximation du ripple sur Tore Supra
 Formule Basiuk
 [R,z] scalaires, vecteurs ou matrices, unit�m
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 Derni�e modif. par LC le 03/08/00
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

0001 function ripple=deltacalmj3(R,z);
0002 
0003 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0004 % [delta,ddelta_dR,ddelta_dz,rho]=delta_fit(R,z);
0005 % Approximation du ripple sur Tore Supra
0006 % Formule Basiuk
0007 % [R,z] scalaires, vecteurs ou matrices, unit�m
0008 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0009 % Derni�e modif. par LC le 03/08/00
0010 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0011 
0012 x=R-2.04;
0013 racine=sqrt(1+1.04*x-0.2704*z.*z);
0014 rho2=1+0.52*x-racine;
0015 rho1=sqrt(rho2);
0016 ripple=2.2e-4*exp(13.6*rho1+11.83*rho2);