Quasi-Fermi Contact

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 "Doping - $ni*exp((DevPsi-Qfn)/$vt) + $ni*exp((Qfp-DevPsi)/$vt)"
pdbSetString $Contact Qfn Equation "Qfn-$Contact"
pdbSetString $Contact Qfp Equation "Qfp-$Contact"
}