A liitle update, I have got the template for my simulation, but I had problem with extracting th Vth from the simulations, anyone can help me?

go atlas

mesh space.mult=0.5

# Adjust x.mesh and y.mesh spacings to reduce node count

x.mesh loc=0.000 spac=0.001

x.mesh loc=0.010 spac=0.001

x.mesh loc=0.015 spac=0.001

x.mesh loc=0.050 spac=0.001

y.mesh loc=0.00 spac=0.001

y.mesh loc=0.002 spac=0.001

y.mesh loc=0.005 spac=0.001

y.mesh loc=0.015 spac=0.001

y.mesh loc=0.018 spac=0.001

y.mesh loc=0.020 spac=0.001

# Define regions

region num=1 material=Air x.min=0 y.min=0

region num=2 material=SiO2 x.min=0.01 x.max=0.02 y.min=0.004 y.max=0.005

region num=3 material=SiO2 x.min=0.01 x.max=0.02 y.min=0.010 y.max=0.011

region num=4 material=silicon x.min=0 x.max=0.010 y.min=0.005 y.max=0.010

region num=6 material=silicon x.min=0.01 x.max=0.02 y.min=0.005 y.max=0.010

region num=5 material=silicon x.min=0.02 x.max=0.030 y.min=0.005 y.max=0.010

# Define electrodes

electrode name=gate x.min=0.01 x.max=0.02 y.min=0.002 y.max=0.004

electrode name=gate2 x.min=0.01 x.max=0.02 y.min=0.011 y.max=0.013

electrode name=source x.min=0 x.max=0.005 y.min=0.003 y.max=0.005

electrode name=drain x.min=0.025 x.max=0.030 y.min=0.003 y.max=0.005

# Define doping

doping uniform concentration=1E18 p.type region=5

doping uniform concentration=1E20 n.type region=4

doping uniform concentration=1E20 n.type region=6

# Define interface

interf qf=3e10 y.max=0.005

# Define contacts

contact name=gate workfun=4.6

contact name=gate2 workfun=4.6 common=gate

contact name=source

contact name=drain

# Activate models

models conmob srh auger bgn fldmob print

# Define outputs

output val.band con.band qfn qfp e.field j.electron j.hole j.conduction j.total ex.field ey.field flowline e.mobility h.mobility qss e.temp h.temp j.disp

# Define solution methods

method newton trap

# Solve initial conditions and perform sweeps

solve init

solve prev

solve vdrain=0

solve vdrain=0.01

solve vdrain=0.02

solve vdrain=0.03

solve vdrain=0.04

solve vdrain=0.05

log outf=nw_gc_igdl.log

solve name=gate vgate=0 vstep=0.05 vfinal=2

Tonyplot nw_gc_igdl.log

save outf=nw_gc_igd_eq.str

Tonyplot nw_gc_igd_eq.str

extract init infile="nw_gc_igdl.log"

extract name="vt1" x.val from curve ((v."gate"), log10(abs(i."drain"))) where y.val=-22

solve init

solve prev

solve vdrain=0

solve vdrain=0.05

solve vdrain=0.1

solve vdrain=0.2

solve vdrain=0.4

solve vdrain=0.6

solve vdrain=0.8

solve vdrain=1

log outf=nw_gc_igd.log

solve ac freq=1e16 name=gate vgate=0 vstep=0.05 vfinal=1

extract init infile="nw_gc_igd.log"

extract name="vt2" x.val from curve ((v."gate"),log10(abs(i."drain"))) where y.val=-22

save outf=nw_gc_igd_eq.str

Tonyplot nw_gc_igd.log

tonyplot -overlay nw_gc_igdl.log nw_gc_igd.log -set nw_gc_igd.log

extract name="dibl" ($"vt1"-$"vt2")/(1-0.05)

extract init infile="nw_gc_igd.log"

extract name="subvt" 1.0/slope (maxslope (curve(v."gate", log10(abs(i."drain")))))

extract name="Ioff_n" y.val from curve(v."gate", i."drain") where x.val=0

extract name="Ion_n" y.val from curve(v."gate", i."drain") where x.val=1

extract name="Ion_Ioff" $Ion_n/$Ioff_n

extract name="gm" curve (v."gate",dydx (v."gate",i."drain")) outf="gm.dat"

quit