#NMOS IDVG curve
go atlas
set Vd=1
set Vg=2
# IDVG 
MESH infile=NMOS.str  space.mult=1.0

# Specify models
#MODELS   MOS PRINT  SRH  BGN  BTBT  TEMP=300
models srh cvt print BTBT temperature=300 

# Set the workfunction of the Poly gate
CONTACT NAME=gate N.POLY

# IDVG curve
METHOD  DAMPED  NEWTON  CARRIERS=2

SOLVE   initial
#
solve vdrain=0 vstep=0.1 vfinal=$"Vd" name=drain
solve vgate=0 vstep=-0.1 vfinal=-1 name=gate
LOG outfile=NMOS_IDVG_Vd$"Vd".log master
solve vgate=-1 vstep=0.05 vfinal=$"Vg" name=gate
#SOLVE   Vsubstrate=0 Vsource=0 Vgate=0 Vdrain=0 ELEC=drain  VSTEP=0.05  NSTEP=20
#SOLVE   Vsubstrate=0 Vsource=0 Vdrain=1 Vgate=-1 ELEC=gate  VSTEP=0.05  NSTEP=60
log off
SAVE outfile=NMOS_IDVG_Vd$"Vd".str
OUTPUT CON.BAND VAL.BAND
SAVE outfile=NMOS_IDVG_Vd$"Vd"_Band.str
#
extract name="vt" (xintercept(maxslope(curve(abs(v."gate"),abs(i."drain")))) - abs(ave(v."drain"))/2.0)
#
extract name="subvt" 1.0/slope(maxslope(curve(abs(v."gate"),log10(abs(i."drain")))))
#

quit