proc qf.ohmic.contact {Contact} {
set vt 0.02558357
set ni 1.1e10
pdbSetBoolean $Contact Qfn Flux 1
pdbSetBoolean $Contact Qfp Flux 1
pdbSetBoolean $Contact DevPsi Flux 1
pdbSetBoolean $Contact Qfn Fixed 1
pdbSetBoolean $Contact Qfp Fixed 1
pdbSetBoolean $Contact DevPsi Fixed 1
pdbSetDouble $Contact Qfp Flux.Scale 1.619e-19
pdbSetDouble $Contact Qfn Flux.Scale 1.619e-19
pdbSetString $Contact DevPsi Equation "ND - NA - $ni*exp((DevPsi-Qfn)/$vt) + $ni*exp((Qfp-DevPsi)/$vt)"
pdbSetString $Contact Qfn Equation "Qfn-$Contact"
pdbSetString $Contact Qfp Equation "Qfp-$Contact"
pdbSetString $Contact Equation "1.619e-19 * (Flux_Qfp - Flux_Qfn)"
}