#Date converted 11 March 2004
#*$
#.subckt fds8870  2 1 3      
# 
#*11 Mar 04
#
template FDS8870 n2,n1,n3
electrical n2,n1,n3
{
var i iscl
dp..model dbodymod =  (isl=1e-11,ikf=17,nl=1.01,rs=2.8e-3,trs1=2e-3,trs2=2e-7,cjo=1.95e-9,m=0.55,tt=9e-11,xti=2.6)
dp..model dbreakmod = (rs=8e-2,trs1=1e-3,trs2=-8.9e-6)
dp..model dplcapmod =  (cjo=1.42e-9,isl=10e-30,nl=10,m=0.38)
m..model mmedmod = (type=_n,vto=1.85,kp=15,is=1e-30, tox=1)
m..model mstrongmod = (type=_n,vto=2.2,kp=650,is=1e-30, tox=1)
m..model mweakmod = (type=_n,vto=1.48,kp=0.05,is=1e-30, tox=1,rs=0.1) 
sw_vcsp..model s1amod =  (ron=1e-5,roff=0.1,von=-5,voff=-3)
sw_vcsp..model s1bmod =  (ron=1e-5,roff=0.1,von=-3,voff=-5)
sw_vcsp..model s2amod =  (ron=1e-5,roff=0.1,von=-2,voff=-0.5)
sw_vcsp..model s2bmod =  (ron=1e-5,roff=0.1,von=-0.5,voff=-2)
c.ca n12 n8 = 2.8e-9
c.cb n15 n14 = 2.8e-9
c.cin n6 n8 = 4.3e-9

dp.dbody n7 n5 = model=dbodymod
dp.dbreak n5 n11 = model=dbreakmod
dp.dplcap n10 n5 = model=dplcapmod

spe.ebreak n11 n7 n17 n18 = 33.62
spe.eds n14 n8 n5 n8 = 1
spe.egs n13 n8 n6 n8 = 1
spe.esg n6 n10 n6 n8 = 1
spe.evthres n6 n21 n19 n8 = 1
spe.evtemp n20 n6 n18 n22 = 1

i.it n8 n17 = 1

l.lgate n1 n9 = 1e-9
l.ldrain n2 n5 = 1.0e-9
l.lsource n3 n7 = 7e-11

res.rlgate n1 n9 = 10
res.rldrain n2 n5 = 10
res.rlsource n3 n7 = 0.7

m.mmed n16 n6 n8 n8 = model=mmedmod, l=1u, w=1u
m.mstrong n16 n6 n8 n8 = model=mstrongmod, l=1u, w=1u 
m.mweak n16 n21 n8 n8 = model=mweakmod, l=1u, w=1u    

res.rbreak n17 n18  = 1, tc1=8.3e-4,tc2=-9e-7
res.rdrain n50 n16  = 3.05e-3, tc1=1.8e-3,tc2=5e-6
res.rgate n9 n20 = 2
res.rslc1 n5 n51  = 1e-6, tc1=1e-4,tc2=1e-6
res.rslc2 n5 n50 = 1e3
res.rsource n8 n7  = 9e-4, tc1=8e-3,tc2=1e-6
res.rvthres n22 n8  = 1, tc1=-1.8e-3,tc2=-9e-6
res.rvtemp n18 n19  = 1, tc1=-2.5e-3,tc2=2e-7
sw_vcsp.s1a n6 n12 n13 n8 = model=s1amod
sw_vcsp.s1b n13 n12 n13 n8 = model=s1bmod
sw_vcsp.s2a n6 n15 n14 n13 = model=s2amod
sw_vcsp.s2b n13 n15 n14 n13 = model=s2bmod

v.vbat n22 n19 = dc=1
equations {
i (n51->n50) +=iscl
iscl: v(n51,n50) = ((v(n5,n51)/(1e-9+abs(v(n5,n51))))*((abs(v(n5,n51)*1e6/500))** 10))
}
}


#.subckt fds8870t th tl
#*thermal model subcircuit
#*11 Mar 04 
# Copper area = 1 sqr inch
#
template thermal_model th tl
thermal_c th, tl
{
ctherm.ctherm1 th 8 =2.0e-3
ctherm.ctherm2 8 7 =5.0e-3
ctherm.ctherm3 7 6 =1.0e-2
ctherm.ctherm4 6 5 =4.0e-2
ctherm.ctherm5 5 4 =9.0e-2
ctherm.ctherm6 4 3 =2e-1
ctherm.ctherm7 3 2 =1e0
ctherm.ctherm8 2 tl =3e0

rtherm.rtherm1 th 8 =1e-1
rtherm.rtherm2 8 7 =5e-1
rtherm.rtherm3 7 6 =1
rtherm.rtherm4 6 5 =5
rtherm.rtherm5 5 4 =8
rtherm.rtherm6 4 3 =12
rtherm.rtherm7 3 2 =18
rtherm.rtherm8 2 tl =25
}