this function calculates the critical field
0001 function xE_c = Ecritical_jd(equil,xrho) 0002 % 0003 % this function calculates the critical field 0004 % 0005 [qe,me,~,~,e0,~,~,~,clum] = pc_dke_yp;%Universal physics constants 0006 % 0007 plnc_e = 31.3 - 0.5*log(equil.pne) + log(equil.pTe*1000);%Coulomb logarithm (Sauter et al. Phys. Plasmas, 6 (1999) 2834) 0008 % 0009 pE_c = equil.pne*qe^3.*plnc_e/(4*pi*e0^2*me*clum^2); 0010 % 0011 if nargin < 2 0012 xE_c = pE_c; 0013 else 0014 prho = equil.ptx(:,1)/equil.ptx(end,1); 0015 xE_c = interp1(prho,pE_c,xrho,'spline'); 0016 end