puts $Base

#set the boron model information
array set $Base {DiffModel {Switch 1 {Constant Fermi Pair React}}}
array set $Base {Charge {Switch 1 {Neutral Acceptor Donor}}}
array set $Base {ActiveModel {Switch 1 {None Solid Precipitation}}}
array set $Base {Discrete {Switch 2 {linear log s2 s4}}}

#set the error bounds for the dopant
array set $Base {Abs.Error {Double 1.0e10}}

array set $Base {Defects {String "Int Vac"}}

array set $Base {Fi {Double { ([pdbGetDouble Si B I D0] + [pdbGetDouble Si B I Dp]) / [pdbGetDouble Si B Dstar] }}}

array set $Base {Dstar {Double { [pdbGetDouble Si B I D0] + [pdbGetDouble Si B I Dp] + [pdbGetDouble Si B V D0] + [pdbGetDouble Si B V Dp] }}}

array set $Base {D0 {Double { [pdbGetDouble Si B I D0] + [pdbGetDouble Si B V D0] }}}

array set $Base {Dp {Double { [pdbGetDouble Si B I Dp] + [pdbGetDouble Si B V Dp] }}}

#solubility data as a function of temperature ... this needs to be sorted!
array set $Base {Solubility {Double {[Arrhenius 7.68e22 0.7086]}}}

#array parameters for defect interactions
array set $Base {Interstitial Interstitial}
array set $Base {interstitial Interstitial}
array set $Base {Interst Interstitial}
array set $Base {interst Interstitial}
array set $Base {inter Interstitial}
array set $Base {Inter Interstitial}
array set $Base {I Interstitial}
array set $Base {Int Interstitial}
array set $Base {Vacancy Vacancy}
array set $Base {vacancy Vacancy}
array set $Base {Vacan Vacancy}
array set $Base {vacan Vacancy}
array set $Base {V Vacancy}
array set $Base {Vac Vacancy}

proc InitBoron {} {
	if { [pdbGetBoolean Si BIC Solve] } {
		sel z=3e8 name=BIC store
	}
	if { [pdbGetBoolean Si BIC1 Solve] } {
		sel z=1e8 name=BIC1 store
		sel z=Boron*.3 name=BoronSub store
	}
	if { [pdbGetBoolean Si BIC2 Solve] } {
		sel z=1e8 name=BIC2 store
		sel z=Boron*.3 name=BoronSub store
	}
}


return { DiffModel Charge ActiveModel Discrete Abs.Error Defects Fi Dstar D0 Dp Solubility Interstitial interstitial Interst interst inter Inter I Vacancy vacancy Vacan vacan V Int Vac }

