disp_ebwES_jd

0001 function [xNperp] = disp_ebwES_jd(xB,xTe,xne,omega,guess_lambda,xguess_Nperp);
0002 %
0003 if nargin >= 6
0004     isNperp = 1;
0005 else
0006     isNperp = 0;
0007 end    
0008 if nargin < 5
0009     guess_lambda = 10;
0010 end
0011 %
0012 N = 100;
0013 %
0014 [qe,me,mp,mn,e0,mu0,re,mc2] = pc_dke_yp;
0015 %
0016 xwpe = sqrt(xne*qe^2/e0/me);
0017 xwce = qe*xB/me;
0018 xbetath = sqrt(xTe/mc2);
0019 %
0020 xabr = xwpe./xwce;
0021 %
0022 xq = omega./xwce;
0023 %
0024 if isNperp == 1
0025     xguess_lambda = (xguess_Nperp.*xq.*xbetath).^2;
0026 else
0027     xguess_lambda = guess_lambda*ones(size(xB));
0028 end
0029 %
0030 [lambda,Yval,Yflag] = inv_ebwES_jd(xq,N,xabr,xguess_lambda);
0031 %
0032 kprho = sqrt(lambda);
0033 xNperp = kprho./xq./xbetath;