NBC is dependent on pH_i and pH_e, 
J_kinetic is scaled by (1+10^(ni(pH_i-pKi)))*(1+10^(-ne(pH_e-pKe)))
This scale is rearranged to 
scale = (1+Kmi^ni*[H_i]^-ni)*(1+Kme^-ne*[H_e]^ne)

2 extra reactions are added for each of the 6 transitions between channel states.
-ni and ne are the stoichiometries of H_i and H_e (only applied in the forward reaction rates)

In final bond graph, each state Sx is the sum of 3 velocities.