#set the interstitial model information
array set $Base {Solve {Boolean 0}}

#List the reaction terms
#each reaction has three components
#  1) a list of equations to enter it into
#  2) the rate constant name
#  3) a list of terms in the reaction

array set $Base {Reaction {String {
    {B2I kf1 {Boron Boron Interstitial}}
    {B2I kr1 B2I}
    {Interstitial Tkf1 {Boron Boron Interstitial}}
    {Interstitial Tkr1 B2I}
    {Boron Tkf1 {Boron Boron Interstitial}}
    {Boron Tkf1 {Boron Boron Interstitial}}
    }}}
#define the rate constants for the names we used above!
array set $Base {kf1 {Double {[pdbGetDouble Si I Capture] * [pdbGetDouble Si I D0] * [pdbGetDouble Si B Interstitial Binding]}}}

array set $Base {kr1 {Double {-[pdbGetDouble Si B2I kf1] * [Arrhenius 5.0e22 2.5] / [pdbGetDouble Si B Interstitial Binding]}}}

array set $Base {Tkf1 {Double {- [pdbGetDouble Si B2I kf1]}}}
array set $Base {Tkr1 {Double {- [pdbGetDouble Si B2I kr1]}}}

return {Model Reaction kf1 kr1 Tkf1 Tkr1}
