conductivity_model_yp

  Conductivity calculated by various models for test only (circular concentric plasma cross-section)
  
 by Y.Peysson (CEA-DRFC) <yves.peysson@cea.fr> and J. Decker (CEA-DRFC) <joan.decker@cea.fr>

0001 function [iar,sigman] = conductivity_model_yp(equil,model,xrho)
0002 %
0003 %  Conductivity calculated by various models for test only (circular concentric plasma cross-section)
0004 %
0005 % by Y.Peysson (CEA-DRFC) <yves.peysson@cea.fr> and J. Decker (CEA-DRFC) <joan.decker@cea.fr>
0006 %
0007 if nargin < 2,
0008     error('Not enough input arguments in conductivity_model_yp.m !');
0009 end
0010 %
0011 if nargin < 3,
0012     %
0013     [equilDKE] = equilibrium_jd(equil);%Magnetic equilibrium for DKE
0014     %
0015 else
0016     %
0017     xpsin_f = rho2psi_jd(equil,xrho);
0018     %
0019     radialDKE.xpsin_S_dke = NaN;
0020     radialDKE.xpsin_f = xpsin_f;
0021     %
0022     [equilDKE] = equilibrium_jd(equil,radialDKE);%Magnetic equilibrium for DKE
0023     %
0024 end
0025 %
0026 iar = equilDKE.xrho*equilDKE.ap/equilDKE.Rp;
0027 %
0028 if model == 1,
0029     sigman = 1 - 1.95*sqrt(iar) + 0.95*iar;%Sigmar-Coppi-Kim model (low inverse aspect ratio limit)
0030 elseif model == 2, 
0031     Iho= hinton_oberman_integral_yp(equilDKE);%Hinton-Oberman Integral (Nucl. Fusion, 9 (1969) 319
0032     sigman = Iho.*(1 + 0.039.*(1 - Iho))./(1 + 0.471.*(1 - Iho));
0033 end