#nwell
region num=6 silicon x.min=0.05 x.max=0.95 y.min=0 y.max=2.5
#n+
region num=7 silicon x.min=0.15 x.max=0.25 y.min=0 y.max=0.8
region num=8 silicon x.min=0.6 x.max=0.8 y.min=0 y.max=0.8
#gate
region num=9 Al x.min=0.35 x.max=0.55 y.min=-0.5 y.max=-0.35
#D/S
region num=10 Al x.min=0.15 x.max=0.25 y.min=-0.3 y.max=-0.2
region num=11 Al x.min=0.6 x.max=0.8 y.min=-0.3 y.max=-0.2
#sub
region num=12 Al x.min=1.5 x.max=3.5 y.min=4.8 y.max=5
#psub
doping uniform conc=4e12 p.type region=1
#nwell erfc[(junc-peak)/char]=0.1 so (junc-peak)/char=1.15 and you can count the char
doping gauss conc=5e13 n.type peak=0.5 char=1.3 junc=2 region=6
#n+
doping gauss conc=5e13 n.type peak=0.25 char=0.22 junc=0.5 region=7
doping gauss conc=5e13 n.type peak=0.25 char=0.22 junc=0.5 region=8
material material=Si taun0=1.0e-7 taup0=1.0e-7
material material=Poly taun0=1.0e-7 taup0=1.0e-7
material region=1 tc.a=2.27 tc.b=0 tc.c=0
material region=3 tc.a=7.9 tc.b=0 tc.c=0