This script generates the data file for an idealized EARth equilibrium. By J. Decker <jodecker@alum.mit.edu> and Y. Peysson <yves.peysson@cea.fr>
0001 % 0002 % This script generates the data file for an idealized EARth equilibrium. 0003 % 0004 % By J. Decker <jodecker@alum.mit.edu> and Y. Peysson <yves.peysson@cea.fr> 0005 % 0006 clear all 0007 % 0008 id = 'EARTH';% scenario identification 0009 % 0010 RE = 6371e3;% earth radius (m) 0011 % 0012 [qe,me,mp,mn,e0,mu0,re,mc2,clum] = pc_dke_yp; 0013 % 0014 npsi = 101;% number of points in psi grid 0015 ntheta = 161;%65;% number of points in theta grid 0016 % 0017 % Equilibrium parameters 0018 % 0019 RL = 5*RE;% Average limit of the plasmasphere from earth center (m) 0020 % 0021 M = -7.84e22;% 1990 measurements 0022 % 0023 equil_prof.zZi = [1,1,1];% Ion types: (1) H/D/T, (2) He, ..., (6) C [1,p] (WARNING: Zi must be [1,1,1,imp1,imp2] for hydrogen plasmas) 0024 equil_prof.zmi = [1,2,3];% Ion mass (uma) [1,p] (WARNING: Zi must be [1,2,3,mimp1,mimp2] for hydrogen plasmas) 0025 zfi = [1,0,0];% Hydrogen isotopic fraction (H/D/T) [1,3] (WARNING: only used when hydrogen plasmas are considered) 0026 % 0027 Te0 = 0.2e-3;% Core electron temperature (keV) 0028 Tea = 1.2e-3;% Edge electron temperature (keV) 0029 % 0030 ne0 = 1e10;% Core electron density (m-3) 0031 nea = 1e9;% Edge electron density (m-3) 0032 % 0033 Ti0 = Te0;% Core ion temperature (keV) 0034 Tia = Tea;% Edge ion temperature (keV) 0035 % 0036 % Plasma parameters 0037 % 0038 prho = linspace(0,1,npsi); 0039 % 0040 equil_prof.pTe = Te0 + prho*(Tea-Te0);%linear approximation 0041 pTi = Ti0 + prho*(Tia-Ti0);%linear approximation 0042 % 0043 equil_prof.pne = ne0*exp(log(nea/ne0)*prho); 0044 pni = equil_prof.pne; 0045 % 0046 equil_prof.pzTi = ones(3,1)*pTi; 0047 equil_prof.pzni = zfi.'*pni; 0048 % 0049 %Magnetic equilibrium 0050 % 0051 equil_magnetic = equildipole_fxd(RL-RE,RE,M,prho,ntheta); 0052 % 0053 equil = conc_struct_jd(equil_magnetic,equil_prof); 0054 equil.id = id; 0055 % 0056 filename = ['EQUIL_',equil.id,'.mat']; 0057 % 0058 save(filename,'equil'); 0059 %