Masetti Mobility Model

proc mobility.masetti {} {
	set mu_max_e 1417.0;		set mu_max_h 470.5
	set mu_e_min1 52.2;			set mu_h_min1 44.9
	set mu_e_min2 52.2;			set mu_h_min2 0.0
	set mu_e1 43.4;				set mu_h1 29.0
	set Pc_e 0.0;				set Pc_h 9.23e16
	set Cr_e 9.68e16;			set Cr_h 2.23e17
	set Cs_e 3.34e20;			set Cs_h 6.1e20
	set a_e 0.68;				set a_h 0.719
	set b_e 2.0;				set b_h 2.0

	term add name=mob1e eqn = "$mu_e_min1*exp(-$Pc_e/(ND+NA))"
	term add name=mob2e eqn = "($mu_max_e-$mu_e_min2)/(1+(((ND+NA)/$Cr_e)^($a_e)))"
	term add name=mob3e eqn = "$mu_e1/(1+(($Cs_e/(ND+NA))^($b_e)))"
	term add name=Emob_Mas eqn = "mob1e+mob2e-mob3e"

	term add name=mob1h eqn = "$mu_h_min1*exp(-$Pc_h/(ND+NA))"
	term add name=mob2h eqn = "($mu_max_h-$mu_h_min2)/(1+(((ND+NA)/$Cr_h)^($a_h)))"
	term add name=mob3h eqn = "$mu_h1/(1+(($Cs_h/(ND+NA))^($b_h)))"
	term add name=Hmob_Mas eqn = "mob1h+mob2h-mob3h"
}

NOTE: "term" will eventually be written out of the FLOOXS code in a future release, but you can change to "solution" at that point. "solution" is more flexible.