sp 1250*(1-exp(-0.091*p0))+400*exp(-0.337*p0) "Auxiliary expression for Th"
Th (300+sp*(log10(Pabs/40)/log10(40))) "Gas temperature"
Lambda_diff ((pi/L)^2+(2.405/Rad)^2)^-0.5 "Effective diffusion length"
alpha_s_b "abs(a1*erf(a2*rho+a3)*exp(-a4/rho^1.35) / ( exp( (gamma-1)/(2*gamma) - 0.49 ) ))" "Ratio of Sheath electronegativity and bulk electronegativity"
rho "abs(alpha+a5*(   exp(a6*(gamma-50))-1)  )" "Auxiliary expression for alpha_s_b"
gamma (e_const*plas.Te/k_B_const)/Th "Ratio of Te and Th"
alpha plas.n_wCl_1m/plas.ne "Electronegativity"
alpha_s abs(alpha_s_b)*alpha "Sheath electronegativity"
a1 0.607 "Data for alpha_s_b and rho"
a2 5.555 "Data for alpha_s_b and rho"
a3 -11.16 "Data for alpha_s_b and rho"
a4 1.634 "Data for alpha_s_b and rho"
a5 12e-3 "Data for alpha_s_b and rho"
a6 -107e-3 "Data for alpha_s_b and rho"
Epsilon_p "0.5*plas.Te*(1+alpha_s) / ( 1+alpha_s*gamma )" "Plasma potential"
Epsilon_s "-plas.Te*log(  abs(4*uB_avg*(1+alpha_s) / (veth*(1+alpha_s*(vi/veth)^2)) ) )" "Sheath potential"
Epsilon_e 2*plas.Te "Mean kinetic energy lost per electron lost"
veth sqrt(8*e_const*plas.Te/(pi*me)) "Electron mean thermal velocity"
vi sqrt(8*k_B_const*Th/(pi*mCl)) "Cl- mean thermal velocity"
hR 0.8*(4+Rad/Lambda_ion)^(-0.5)/(1+alpha0) "Auxiliary expression for the correction factors"
alpha0 alpha*3/2 "Auxiliary expression for the correction factors"
hL 0.86*(3+L/2/Lambda_ion)^(-0.5)/(1+alpha0) "Auxiliary expression for the correction factors"
hcCl2 "1/( gamma^0.5  + gamma^0.5*NCl2_s^0.5*plas.n_wCl2_1p*(plas.n_wCl_1m)^(-3/2)  )" "Auxiliary expression for the correction factors"
NCl2_s 15*vpCl2/(56*Krec*Lambda_ion) "Auxiliary expression for the correction factors"
Krec 5e-14*(300/Th)^0.5 "Rate coefficient for mutual neutralization"
hcCl "1/( gamma^0.5  + gamma^0.5*NCl_s^0.5*plas.n_wCl_1p*(plas.n_wCl_1m)^(-3/2)  )" "Auxiliary expression for the correction factors"
NCl_s 15*vpCl/(56*Krec*Lambda_ion) "Auxiliary expression for the correction factors"
hRCl2 "sqrt( hR^2 + hcCl2^2  )" "Radial correction factor for Cl2+"
hRCl "sqrt( hR^2 + hcCl^2  )" "Radial correction factor for Cl+"
hLCl2 "sqrt( hL^2 + hcCl2^2  )" "Longitudinal correction factor for Cl2+"
hLCl "sqrt( hL^2 + hcCl^2  )" "Longitudinal correction factor for Cl+"
gammaCl_steel "10^( -1.22 -1.34*exp(-1.48*plas.n_wCl/plas.n_wCl2) )" "Cl recombination coefficient"
mCl plas.M_wCl/N_A_const "Cl mass"
mCl2 plas.M_wCl2/N_A_const "Cl2 mass"
vpCl "sqrt( 8*k_B_const*Th / (pi*mCl) )" "Cl mean thermal velocity"
vpCl2 "sqrt( 8*k_B_const*Th / (pi*mCl2) )" "Cl2 mean thermal velocity"
me 9.1e-31[kg] "Electron mass"
Lambda_ion k_B_const*Th/(p0*sigma_ion) "Mean-free-path"
sigma_ion 7e-19[m^2] "Effective scattering cross-section"
uB_avg "(plas.n_wCl2_1p*plas.vBohm_wCl2_1p + plas.n_wCl_1p*plas.vBohm_wCl_1p) / (plas.n_wCl_1p+plas.n_wCl2_1p)" "Bohm velocity weighted average"