* Diode Vth

.lib ../../../libraries/sky130_fd_pr/PDKVERSION/models/sky130.lib.spice CORNER

X1 3 0 DEVICENAME area=1.0 M=1
Rd 3 4 100

* DC source for current measure
Vid 4 5
Vdd 5 0 DC 0V

.control
* Sweep Vdd from 0 to 2.0V
dc Vdd 0 2.0 0.05
* NOTE:  Internally accessed names MUST be in lowercase
wrdata DEVICENAME.data -Vid#branch V(3)
* Find threshold
let ih=-Vid#branch[38]
let il=-Vid#branch[25]
let vh=V(3)[38]
let vl=V(3)[25]
let vth=((vl - vh) / (ih - il)) * ih + vh
echo threshold voltage
print vth
quit
.endc
.end