EbyP	 abs(Vx)/Pgas	                            
Eave	 FEave(EbyP)	                             
Te	   if((ne/ne_center)>0.001,abs(ee/ne),0)	   
Ri	   Pgas*FAi(Te)*FVe(Te)*ne	                 
Lion	 Ri*ETHion	                               
Lexc	 Ri*ETHexc	                               
