#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 {
    {BIC2 kf2 {Boron Boron Interstitial}}
    {BIC2 kr2 BIC2}
    {Interstitial Tkf2 {Boron Interstitial Boron}}
    {Interstitial Tkr2 BIC2}
    {Boron Tkf2 {Boron Boron Interstitial}}
    {Boron Tkf2 {Boron Boron Interstitial}}
    }}}
#define the rate constants for the names we used above!
array set $Base {kf2 {Double {[pdbGetDouble Si Interstitial Capture] * [pdbGetDouble Si B Interstitial D0] * [pdbGetDouble Si B I Binding]}}}

array set $Base {kr2 {Double {-[pdbGetDouble Si BIC2 kf2] * [Arrhenius 5.0e22 2.5] / [pdbGetDouble Si B I Binding]}}}


array set $Base {tkf2 {Double {- [pdbGetDouble Si BIC2 kf2]}}}
array set $Base {tkr2 {Double {- [pdbGetDouble Si BIC2 kr2]}}}
 
return {Model Reaction kf2 kr2 tkf2 tkr2}
