diff --git a/vdmos/ro_11_vdmos.cir b/vdmos/ro_11_vdmos.cir
index 3795e9f00..c5dcd6179 100644
--- a/vdmos/ro_11_vdmos.cir
+++ b/vdmos/ro_11_vdmos.cir
@@ -20,8 +20,8 @@ xinv8 10 9 1 0 inv
 xinv9 2 10 1 0 inv
 
 
-.model  N1  vdmos cgdmin=0.05p cgdmax=0.2p a=1.2 cgs=0.15p rg=10 kp=2e-5 rb=1e7 cjo=1n
-.model  P1  vdmos cgdmin=0.05p cgdmax=0.2p a=1.2 cgs=0.15p rg=10 kp=2e-5 rb=1e7 cjo=1n pchan
+.model  N1  vdmos cgdmin=0.05p cgdmax=0.2p a=1.2 cgs=0.15p rg=10 kp=2e-5 rb=1e7 cjo=1n subthres=0.2
+.model  P1  vdmos cgdmin=0.05p cgdmax=0.2p a=1.2 cgs=0.15p rg=10 kp=2e-5 rb=1e7 cjo=1n pchan subthres=0.2
 
 .tran 0.1n 1u
 
diff --git a/vdmos/vdmos-out.cir b/vdmos/vdmos-out.cir
index 140f77e19..2c4c599f9 100644
--- a/vdmos/vdmos-out.cir
+++ b/vdmos/vdmos-out.cir
@@ -1,7 +1,8 @@
 VDMOS output
 
 m1 d g s s n1
-.model n1 vdmos rb=0.05 is=10n kp=2 bv=12 rd=0.1
+*.model n1 vdmos rb=0.05 is=10n kp=2 bv=12 rd=0.1
+.model  N1  vdmos vto=1 cgdmin=0.05p cgdmax=0.2p a=1.2 cgs=0.15p rg=10 kp=2e-4 rb=1e4 is=1e-9 bv=12 cjo=1p subthres=1
 *d1 s d dmod
 *.model dmod d  is=10n rs=0.05
 
@@ -10,11 +11,13 @@ vg g 0 1
 vs s 0 0
 vb b 0 0
 
-.dc vd -2 15 0.05 vg 0 5 1
+*.dc vd -2 15 0.05 vg 0 5 1
 
 .control
-run
+dc vd -2 15 0.05 vg 0 5 1
 plot vs#branch
+dc vg 0 5 0.05 vd 0.5 2.5 0.5
+plot log(vs#branch)
 .endc
 
 .end
diff --git a/vdmos/vdmosp-out-mtr.cir b/vdmos/vdmosp-out-mtr.cir
index 6931249d2..576913937 100644
--- a/vdmos/vdmosp-out-mtr.cir
+++ b/vdmos/vdmosp-out-mtr.cir
@@ -9,18 +9,24 @@ m1 d g s s IRF7233
 m2 d g s2 s2 IRF7233_2
 .model IRF7233_2 VDMOS(pchan mtriode=2 Rg=3 Rd=8m Rs=6m Vto=-1 Kp=70 Cgdmax=2n Cgdmin=.25n Cgs=3.3n Cjo=.98n Is=98p Rb=10m mfg=International_Rectifier Vds=-12 Ron=20m Qg=49n)
 
-
+m3 d g s3 s3 IRF7233_3
+.model IRF7233_3 VDMOS(pchan mtriode=2 Rg=3 Rd=8m Rs=6m Vto=-1 Kp=70 Cgdmax=2n Cgdmin=.25n Cgs=3.3n Cjo=.98n Is=98p Rb=10m mfg=International_Rectifier Vds=-12 Ron=20m Qg=49n subthres=0.5)
 
 vd d 0 -5
 vg g 0 -5
 vs s 0 0
 vs2 s2 0 0
+vs3 s3 0 0
+
 
 .dc vd -12 1 0.05 vg 0 -5 -1
 
 .control
-run
-plot vs#branch vs2#branch
+dc vd -12 1 0.05 vg 0 -5 -1
+plot vs#branch vs2#branch vs3#branch
+dc vg 0 -4 -0.05 vd -1 -12 -2
+plot vs#branch vs2#branch vs3#branch
+plot log(-vs#branch) log(-vs2#branch) log(-vs3#branch)
 .endc
 
 .end