K_junc	  2*Rg*T/Far*(1+dlnfdlnc)*(1-t_plus)	 Junction potantial coefficient  
i1	      -K1eff*phi1x	                       Solid phase current density     
i2	      -K2eff*(phi2x-K_junc*1/c2*c2x)	     Electrolyte current density     
soc_pos	 c1s_pos/c1max_pos	                  State of charge Positive        
soc_neg	 c1s_neg/c1max_neg	                  State of charge Negative