XPP model

This model was converted from XPP ode format to SBML using sbmlutils-0.1.5a6.

# rubin_terman_dbs.ode :  For Rubin & Terman, JCNS, 2004

# identical to rubin_terman_pd.ode but with
# periodic DBS stimulation applied

# turn off stimulation by changing the initial condition
# for i1 to il=0

# note that there may be some numerical issues after 1000 msec so we
#  avoided long time simulations in the paper

# stn parameters
p vl=-60,vna=55.,vk=-80.,thetam=30,sm=15
p gl=2.25,gna=37.5,gk=45,tn=1.,th=0.05
p gahp=9.,gca=.5,vca=140.,k1=15.,eps=5e-05
p kca=22.5,thetas=39.,ss=8.,xp=1.,i=0.
p thetah=-39,sh=3.1,thetan=-32.,sn=-8.
p taun0=1,taun1=100.,thn=80.,sigman=26.
p tauh0=1,tauh1=500,thh=57.,sigmah=3.,phi=.75
p thetat=-63.,kt=-7.8,gt=.5,phir=.5
p thetar=-67,kr=2.,taur0=7.1,taur1=17.5,thr=-68,sigmar=2.2
p alpha=5,beta=1.,ab=-30.,gsyn=0.9,vsyn=-100
p rth=.25,rsig=-.07,rho1=.5,a1=.9,hets=0,inp=25
p alphai=1,betai=.05

# gp parameters
p gnag=120.,gkg=30.,gahpg=30.,gtg=.5,gcag=.1,glg=.1
p vnag=55.,vkg=-80.,vcag=120.,vlg=-55.
p thetasg=-57.,ksg=2.,thetas1g=-35.,ks1g=2.
p thetarg=-70.,krg=-2.,taurg=30.
p thetamg=-37.,sigmamg=10.
p thetang=-50.,sigmang=14.
p taun0g=.05,taun1g=.27,thng=-40,sng=-12
p thetahg=-58,sigmahg=-12
p tauh0g=.05,tauh1g=.27,thhg=-40,shg=-12
p k1g=30.,kcag=15.,epsg=0.0001
p phig=1.,deltang=.1,deltahg=.05
p iapp=2,gsyngg=.3,vsyngg=-80.
p gsyng=0.3,vsyng=0,alphag=2,betag=.04,abg=-20
p iappi=5,igl=1.,thresg=0.0,hetg=0.,betagi=.08,kcagi=15
p gsyngi=.3,gsynggi=1.5,alphaggi=1,betaggi=.1

# TC parameters
p vsyntc=-85,gsyntc=.08,asg=200,bsg=.4
p itc=6,shi=-80,shi2=-90,dur=5,dur2=10,period=25
p gnabar=3,gkbar=5,glbar=.05,ena=50,ek=-90,eleak=-70
p gtbar=5,qht=2.5,tadj=1,apr=4,apt=1

# stn functions
sinf(v)=1./(1.+exp(-(v+thetas)/ss))
minf(v)=1./(1.+exp(-(v+thetam)/sm))
hinf(v)=1./(1.+exp((v-thetah)/sh))
ninf(v)=1./(1.+exp((v-thetan)/sn))
taun(v)=taun0+taun1/(1+exp((v+thn)/sigman))
tauh(v)=tauh0+tauh1/(1+exp((v+thh)/sigmah))
rinf(v)=1/(1+exp((v-thetar)/kr))
taur(v)=taur0+taur1/(1+exp((v+thr)/sigmar))
tinf(v)=1/(1+exp((v-thetat)/kt))
rnew(r)=1/(1+exp((r-rth)/rsig))-1/(1+exp(-rth/rsig))

# gp functions
sinfg(vg)=1/(1+exp(-(vg-thetasg)/ksg))
sinf1g(vg)=1/(1+exp(-(vg-thetas1g)/ks1g))
rinfg(vg)=1/(1+exp(-(vg-thetarg)/krg))
minfg(vg)=1./(1.+exp(-(vg-thetamg)/sigmamg))
ninfg(vg)=1./(1.+exp(-(vg-thetang)/sigmang))
taung(vg)=taun0g+taun1g/(1+exp(-(vg-thng)/sng))
hinfg(vg)=1./(1.+exp(-(vg-thetahg)/sigmahg))
tauhg(vg)=tauh0g+tauh1g/(1+exp(-(vg-thhg)/shg))
hv(v)=1./(1.+exp(-v/.001))

# TC functions
inaf(v,m,h)=gnabar*m^3*h*(v-ena)
ikf(v,n)=gkbar*n^4*(v-ek)
ilf(v)=glbar*(v-eleak)
itf(v,mt,ht)=gtbar*mt^2*ht*v
minftc(v)  =  1/(1+exp(-(v+37)/7))
mtinf(v)  = 1/(1+exp(-(v+60)/6.2))
ah(v) =  0.128*exp(-(46+v)/18)
bh(v) =  apr/(1+exp(-(23+v)/5))
tauhtc(v)  =  1/(ah(v)+bh(v))
hinftc(v)  =  1/(1+exp((v+41)/4))
htinf(v)  = 1/(1+exp((v+84)/4))
tauht(v)=(28+apt*exp((v+25)/(-10.5)))

# function for cortical inputs, periodic case:
ftc(t)=itc*hv(sin(6.2831853*(t+shi)/period))*(1-hv(sin(6.2831853*(t+shi+dur)/period)))

# function used to set up windows for calculation of error index:
ftc2(t)=hv(sin(6.2831853*(t+shi2)/period))*(1-hv(sin(6.2831853*(t+shi2+dur2)/period)))

# stn currents
il(v)=gl*(v-vl)
ina(v,h)=gna*(minf(v))^3*h*(v-vna)
ik(v,n)=gk*n^4*(v-vk)
iahp(v,ca)=gahp*(v-vk)*ca/(ca+k1)
ica(v)=gca*((sinf(v))^2)*(v-vca)
it(v,r)=gt*(tinf(v)**3)*(rnew(r)^2)*(v-vca)

isyn1=gsyn*(sg5+sg2)*(v1-vsyn)
isyn2=gsyn*(sg6+sg1)*(v2-vsyn)
isyn3=gsyn*(sg8+sg4)*(v3-vsyn)
isyn4=gsyn*(sg7+sg3)*(v4-vsyn)
isyn5=gsyn*(sg2+sg6)*(v5-vsyn)
isyn6=gsyn*(sg1+sg5)*(v6-vsyn)
isyn7=gsyn*(sg3+sg8)*(v7-vsyn)
isyn8=gsyn*(sg4+sg7)*(v8-vsyn)
isyn9=gsyn*(sg13+sg10)*(v9-vsyn)
isyn10=gsyn*(sg14+sg9)*(v10-vsyn)
isyn11=gsyn*(sg16+sg12)*(v11-vsyn)
isyn12=gsyn*(sg15+sg11)*(v12-vsyn)
isyn13=gsyn*(sg10+sg14)*(v13-vsyn)
isyn14=gsyn*(sg9+sg13)*(v14-vsyn)
isyn15=gsyn*(sg11+sg16)*(v15-vsyn)
isyn16=gsyn*(sg12+sg15)*(v16-vsyn)


# gp currents
itg(vg,rg)=gtg*(sinfg(vg)^3)*rg*(vg-vcag)
inag(vg,hg)=gnag*(minfg(vg)^3)*hg*(vg-vnag)
ikg(vg,ng)=gkg*(ng^4)*(vg-vkg)
iahpg(vg,cag)=gahpg*(vg-vkg)*cag/(cag+k1g)
icag(vg)=gcag*((sinf1g(vg))^2)*(vg-vcag)
ilg(vg)=glg*(vg-vlg)

isyng1=gsyng*(s8+s4+.2*s16)*(vg1-vsyng)
isyng2=gsyng*(s7+s3+.2*s15)*(vg2-vsyng)
isyng3=gsyng*(s1+s5+.2*s9)*(vg3-vsyng)
isyng4=gsyng*(s2+s6+.2*s10)*(vg4-vsyng)
isyng5=gsyng*(s4+s8+.2*s12)*(vg5-vsyng)
isyng6=gsyng*(s3+s7+.2*s11)*(vg6-vsyng)
isyng7=gsyng*(s5+s2+.2*s13)*(vg7-vsyng)
isyng8=gsyng*(s6+s1+.2*s14)*(vg8-vsyng)

isyng9=gsyng*(s16+s12+.2*s8)*(vg9-vsyng)
isyng10=gsyng*(s15+s11+.2*s7)*(vg10-vsyng)
isyng11=gsyng*(s9+s13+.2*s1)*(vg11-vsyng)
isyng12=gsyng*(s10+s14+.2*s2)*(vg12-vsyng)
isyng13=gsyng*(s12+s16+.2*s4)*(vg13-vsyng)
isyng14=gsyng*(s11+s15+.2*s3)*(vg14-vsyng)
isyng15=gsyng*(s13+s10+.2*s5)*(vg15-vsyng)
isyng16=gsyng*(s14+s9+.2*s6)*(vg16-vsyng)


isyngg1=gsyngg*(sg2+sg3)*(vg1-vsyngg)
isyngg2=gsyngg*(sg1+sg5)*(vg2-vsyngg)
isyngg3=gsyngg*(sg4+sg8)*(vg3-vsyngg)
isyngg4=gsyngg*(sg3+sg1)*(vg4-vsyngg)
isyngg5=gsyngg*(sg6+sg7)*(vg5-vsyngg)
isyngg6=gsyngg*(sg5+sg2)*(vg6-vsyngg)
isyngg7=gsyngg*(sg8+sg3)*(vg7-vsyngg)
isyngg8=gsyngg*(sg7+sg4)*(vg8-vsyngg)
isyngg9=gsyngg*(sg10+sg11)*(vg9-vsyngg)
isyngg10=gsyngg*(sg9+sg13)*(vg10-vsyngg)
isyngg11=gsyngg*(sg12+sg16)*(vg11-vsyngg)
isyngg12=gsyngg*(sg11+sg9)*(vg12-vsyngg)
isyngg13=gsyngg*(sg14+sg15)*(vg13-vsyngg)
isyngg14=gsyngg*(sg13+sg10)*(vg14-vsyngg)
isyngg15=gsyngg*(sg16+sg11)*(vg15-vsyngg)
isyngg16=gsyngg*(sg15+sg12)*(vg16-vsyngg)


isyngi[1..16]=gsyngi*si[j]*(vgi[j]-vsyng)

isynggi[1..16]=gsynggi*sggi[j]*(vgi[j]-vsyn)

i V1=-78.,V2=-55.,V3=-78.,V4=-52.
i V5=-78.,V6=-78.,H1=0.2,H2=0.08
i H3=0.2,H4=0.08,H5=0.2,H6=0.2
i N1=0.08,N2=0.4,N3=0.08,N4=0.43
i N5=0.08,N6=0.08,R1=0.67,R2=0.38
i R3=0.67,R4=0.38,R5=0.67,R6=0.67
i CA1=0.44,CA2=0.47,CA3=0.4,CA4=0.4
i CA5=0.44,CA6=0.44,S1=0.001,S2=0.3
i S3=0.001,S4=0.35,S5=0.001,S6=0.31
i VG1=-16.,VG2=-69.,VG3=-6.8,VG4=-69.
i VG5=-8.,VG6=-69.,NG1=0.84,NG2=0.2
i NG3=0.8,NG4=0.2,NG5=0.83,NG6=0.2
i HG1=0.1,HG2=0.7,HG3=0.14,HG4=0.7
i HG5=0.15,HG6=0.7,RG1=0.6,RG2=0.47
i RG3=0.66,RG4=0.47,RG5=0.66,RG6=0.47
i CAG1=0.068,CAG2=0.06,CAG3=0.067,CAG4=0.05
i CAG5=0.067,CAG6=0.05,SG1=0.025,SG2=0.025
i SG3=0.06,SG4=0.025,SG5=0.025,SG6=0.025

i V10=-78.,V9=-512.,V11=-78.,V13=-129.
i V12=-78.,V15=-78.,H10=0.9,H9=0.08
i H11=0.9,H13=0.08,H12=0.9,H15=0.9
i N10=0.08,N9=0.13,N11=0.08,N13=0.1311
i N12=0.08,N15=0.08,R10=0.157,R9=0.118
i R11=0.157,R13=0.118,R12=0.157,R15=0.157
i CA10=0.1313,CA9=0.137,CA11=0.13,CA13=0.13
i CA12=0.1313,CA15=0.1313,S10=0.0010,S9=0.11
i S11=0.0010,S13=0.1112,S12=0.0010,S15=0.1110
i VG10=-1015.,VG9=-159.,VG11=-15.8,VG13=-159.
i VG12=-8.,VG15=-159.,NG10=0.813,NG9=0.9
i NG11=0.8,NG13=0.9,NG12=0.811,NG15=0.9
i HG10=0.10,HG9=0.7,HG11=0.1013,HG13=0.7
i HG12=0.1012,HG15=0.7,RG10=0.15,RG9=0.137
i RG11=0.1515,RG13=0.137,RG12=0.1515,RG15=0.137
i CAG10=0.0158,CAG9=0.015,CAG11=0.0157,CAG13=0.012
i CAG12=0.0157,CAG15=0.012,SG10=0.0912,SG9=0.0912
i SG11=0.015,SG13=0.0912,SG12=0.0912,SG15=0.0912

i vt1=-65.7,ht1=.008577

# stn equations
v[1..16]'=-(il(v[j])+ina(v[j],h[j])+ik(v[j],n[j])+iahp(v[j],ca[j])+ica(v[j])+it(v[j],r[j]))-isyn[j]+hets*[j]+inp+dbs(t)
h[1..16]'=phi*( hinf(v[j])-h[j] )/tauh(v[j])
n[1..16]'=phi*( ninf(v[j])-n[j] )/taun(v[j])
r[1..16]'=phir*(rinf(v[j])-r[j])/taur(v[j])
ca[1..16]'=phi*eps*(-ica(v[j])-it(v[j],r[j]) - kca*ca[j])
s[1..16]'=alpha*(1-s[j])*sinf(v[j]+ab)-beta*s[j]
si[1..16]'=alphai*(1-si[j])*sinf(v[j]+ab)-betai*si[j]

# gp equations
vg[1..16]'= -(itg(vg[j],rg[j])+inag(vg[j],hg[j])+ikg(vg[j],ng[j])+iahpg(vg[j],cag[j])+icag(vg[j])+ilg(vg[j]))+iapp-isyngg[j]-isyng[j]+hetg*[j]
ng[1..16]'= deltang*(ninfg(vg[j])-ng[j])/taung(vg[j])
hg[1..16]'= deltahg*(hinfg(vg[j])-hg[j])/tauhg(vg[j])
rg[1..16]'= phig*(rinfg(vg[j])-rg[j])/taurg
cag[1..16]'=epsg*(-icag(vg[j])-itg(vg[j],rg[j]) - kcag*cag[j])
sg[1..16]'=alphag*(1-sg[j])*sinfg(vg[j]+abg)-betag*sg[j]
sggi[1..16]'=alphaggi*(1-sggi[j])*sinfg(vg[j]+abg)-betaggi*sggi[j]

# gpi equations
vgi[1..16]'= -(itg(vgi[j],rgi[j])+inag(vgi[j],hgi[j])+ikg(vgi[j],ngi[j])+iahpg(vgi[j],cagi[j])+icag(vgi[j])+ilg(vgi[j]))+iappi-isyngi[j]-isynggi[j]
ngi[1..16]'= deltang*(ninfg(vgi[j])-ngi[j])/taung(vgi[j])
hgi[1..16]'= deltahg*(hinfg(vgi[j])-hgi[j])/tauhg(vgi[j])
rgi[1..16]'= phig*(rinfg(vgi[j])-rgi[j])/taurg
cagi[1..16]'=epsg*(-icag(vgi[j])-itg(vgi[j],rgi[j]) - kcagi*cagi[j])
sgi[1..16]'=alphag*(1-sgi[j])*sinfg(vgi[j]+abg)-betagi*sgi[j]

# TC equations
vt1'=-(ilf(vt1)+inaf(vt1,minftc(vt1),ht1)+ikf(vt1,.75*(1-ht1))+itf(vt1,mtinf(vt1),htt1))+ftc(t)-gsyntc*sga*(vt1-vsyntc)
htt1'=qht*(htinf(vt1)-htt1)/tauht(vt1)
ht1'=tadj*(hinftc(vt1)-ht1)/tauhtc(vt1)

vt2'=-(ilf(vt2)+inaf(vt2,minftc(vt2),ht2)+ikf(vt2,.75*(1-ht2))+itf(vt2,mtinf(vt2),htt2))+ftc(t)-gsyntc*sgb*(vt2-vsyntc)
htt2'=qht*(htinf(vt2)-htt2)/tauht(vt2)
ht2'=tadj*(hinftc(vt2)-ht2)/tauhtc(vt2)

# additional variables used for error index calculations and DBS
# note that i1 is amplitude of DBS:
#  with i1'=0, signal turns on or off
#  il'=f(t) can be used to specify form of signal while on
# pdbs is period of DBS
# ddbs is on duration of each DBS pulse
# both of these are constant by default but can be made time-dependent

coun1'=0
coun2'=0
pdbs'=0
ddbs'=0
i1'=0
z[1..2]'=1-hv(z[j]-50)
good[1..2]'=0

p thr0=20
i z1=50,z2=50
i good1=0,good2=0

global 1 {(vt1+thr0)*hv(t-25)*hv(ftc2(t)-0.5)} {coun1=coun1+1}
global 1 {(vt2+thr0)*hv(t-25)*hv(ftc2(t)-0.5)} {coun2=coun2+1}
global 1 {(vt1+10)*hv(t-25)*hv(ftc2(t)-0.5)*hv(z1-40)} {z1=0;good1=good1+1}
global 1 {(vt2+10)*hv(t-25)*hv(ftc2(t)-0.5)*hv(z2-40)} {z2=0;good2=good2+1}
dbs(t)=i1*hv(sin(6.2831853*t/pdbs))*(1-hv(sin(6.2831853*(t+ddbs)/pdbs)))

sga=sgi1+sgi2+sgi5+sgi6+sgi9+sgi10+sgi13+sgi14
sgb=sgi3+sgi4+sgi7+sgi8+sgi11+sgi12+sgi15+sgi16

aux sg=ftc(t)
aux sgt=ftc2(t)

aux sa=sga
aux sb=sgb


@ maxstor=100001,dt=.1,total=1000,meth=qualrk,xp=t,yp=v1,xlo=0,xhi=1000,ylo=-80,yhi=20.,bound=50000

done


This file has been produced by sbmlutils.

Terms of use

Copyright © 2017 Matthias Koenig

Redistribution and use of any part of this model, with or without modification, are permitted provided that the following conditions are met:

  1. Redistributions of this SBML file must retain the above copyright notice, this list of conditions and the following disclaimer.
  2. Redistributions in a different form must reproduce the above copyright notice, this list of conditions and the following disclaimer in the documentation and/or other materials provided with the distribution.
This model is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.


Model :

id
name
time
substance
extent
volume
area
length
Access SBML model  L3V1

FunctionDefinitions [50] name math sbo cvterm
max minimum x y x x y y
min maximum x y x x y y
heav heavyside x 0 x 0 0.5 x 0 1 x 0 0
mod modulo x y x y x y x 0 y 0 x y x y
sinf v ss thetas 1 1 v thetas ss
minf v sm thetam 1 1 v thetam sm
hinf v sh thetah 1 1 v thetah sh
ninf v sn thetan 1 1 v thetan sn
taun v sigman taun0 taun1 thn taun0 taun1 1 v thn sigman
tauh v sigmah tauh0 tauh1 thh tauh0 tauh1 1 v thh sigmah
rinf v kr thetar 1 1 v thetar kr
taur v sigmar taur0 taur1 thr taur0 taur1 1 v thr sigmar
tinf v kt thetat 1 1 v thetat kt
rnew r rsig rth 1 1 r rth rsig 1 1 rth rsig
sinfg vg ksg thetasg 1 1 vg thetasg ksg
sinf1g vg ks1g thetas1g 1 1 vg thetas1g ks1g
rinfg vg krg thetarg 1 1 vg thetarg krg
minfg vg sigmamg thetamg 1 1 vg thetamg sigmamg
ninfg vg sigmang thetang 1 1 vg thetang sigmang
taung vg sng taun0g taun1g thng taun0g taun1g 1 vg thng sng
hinfg vg sigmahg thetahg 1 1 vg thetahg sigmahg
tauhg vg shg tauh0g tauh1g thhg tauh0g tauh1g 1 vg thhg shg
hv v 1 1 v 0.001
inaf v m h ena gnabar gnabar m 3 h v ena
ikf v n ek gkbar gkbar n 4 v ek
ilf v eleak glbar glbar v eleak
itf v mt ht gtbar gtbar mt 2 ht v
minftc v 1 1 v 37 7
mtinf v 1 1 v 60 6.2
ah v 0.128 46 v 18
bh v apr apr 1 23 v 5
tauhtc v apr 1 ah v bh v apr
hinftc v 1 1 v 41 4
htinf v 1 1 v 84 4
tauht v apt 28 apt v 25 10.5
ftc t dur itc period shi itc hv 6.2831853 t shi period 1 hv 6.2831853 t shi dur period
ftc2 t dur2 period shi2 hv 6.2831853 t shi2 period 1 hv 6.2831853 t shi2 dur2 period
il v gl vl gl v vl
ina v h gna sm thetam vna gna minf v thetam sm 3 h v vna
ik v n gk vk gk n 4 v vk
iahp v ca gahp k1 vk gahp v vk ca ca k1
ica v gca ss thetas vca gca sinf v ss thetas 2 v vca
it v r gt kt rsig rth thetat vca gt tinf v kt thetat 3 rnew r rsig rth 2 v vca
itg vg rg gtg ksg thetasg vcag gtg sinfg vg ksg thetasg 3 rg vg vcag
inag vg hg gnag sigmamg thetamg vnag gnag minfg vg sigmamg thetamg 3 hg vg vnag
ikg vg ng gkg vkg gkg ng 4 vg vkg
iahpg vg cag gahpg k1g vkg gahpg vg vkg cag cag k1g
icag vg gcag ks1g thetas1g vcag gcag sinf1g vg thetas1g ks1g 2 vg vcag
ilg vg glg vlg glg vg vlg
dbs t ddbs i1 pdbs i1 hv 6.2831853 t pdbs 1 hv 6.2831853 t ddbs pdbs

Parameters [371] name constant value unit derived unit sbo cvterm
vl vl = -60 -60.0 None
vna vna = 55. 55.0 None
vk vk = -80. -80.0 None
thetam thetam = 30 30.0 None
sm sm = 15 15.0 None
gl gl = 2.25 2.25 None
gna gna = 37.5 37.5 None
gk gk = 45 45.0 None
tn tn = 1. 1.0 None
th th = 0.05 0.05 None
gahp gahp = 9. 9.0 None
gca gca = .5 0.5 None
vca vca = 140. 140.0 None
k1 k1 = 15. 15.0 None
eps eps = 5e-05 5e-05 None
kca kca = 22.5 22.5 None
thetas thetas = 39. 39.0 None
ss ss = 8. 8.0 None
xp xp = 1. 1.0 None
i i = 0. 0.0 None
thetah thetah = -39 -39.0 None
sh sh = 3.1 3.1 None
thetan thetan = -32. -32.0 None
sn sn = -8. -8.0 None
taun0 taun0 = 1 1.0 None
taun1 taun1 = 100. 100.0 None
thn thn = 80. 80.0 None
sigman sigman = 26. 26.0 None
tauh0 tauh0 = 1 1.0 None
tauh1 tauh1 = 500 500.0 None
thh thh = 57. 57.0 None
sigmah sigmah = 3. 3.0 None
phi phi = .75 0.75 None
thetat thetat = -63. -63.0 None
kt kt = -7.8 -7.8 None
gt gt = .5 0.5 None
phir phir = .5 0.5 None
thetar thetar = -67 -67.0 None
kr kr = 2. 2.0 None
taur0 taur0 = 7.1 7.1 None
taur1 taur1 = 17.5 17.5 None
thr thr = -68 -68.0 None
sigmar sigmar = 2.2 2.2 None
alpha alpha = 5 5.0 None
beta beta = 1. 1.0 None
ab ab = -30. -30.0 None
gsyn gsyn = 0.9 0.9 None
vsyn vsyn = -100 -100.0 None
rth rth = .25 0.25 None
rsig rsig = -.07 -0.07 None
rho1 rho1 = .5 0.5 None
a1 a1 = .9 0.9 None
hets hets = 0 0.0 None
inp inp = 25 25.0 None
alphai alphai = 1 1.0 None
betai betai = .05 0.05 None
gnag gnag = 120. 120.0 None
gkg gkg = 30. 30.0 None
gahpg gahpg = 30. 30.0 None
gtg gtg = .5 0.5 None
gcag gcag = .1 0.1 None
glg glg = .1 0.1 None
vnag vnag = 55. 55.0 None
vkg vkg = -80. -80.0 None
vcag vcag = 120. 120.0 None
vlg vlg = -55. -55.0 None
thetasg thetasg = -57. -57.0 None
ksg ksg = 2. 2.0 None
thetas1g thetas1g = -35. -35.0 None
ks1g ks1g = 2. 2.0 None
thetarg thetarg = -70. -70.0 None
krg krg = -2. -2.0 None
taurg taurg = 30. 30.0 None
thetamg thetamg = -37. -37.0 None
sigmamg sigmamg = 10. 10.0 None
thetang thetang = -50. -50.0 None
sigmang sigmang = 14. 14.0 None
taun0g taun0g = .05 0.05 None
taun1g taun1g = .27 0.27 None
thng thng = -40 -40.0 None
sng sng = -12 -12.0 None
thetahg thetahg = -58 -58.0 None
sigmahg sigmahg = -12 -12.0 None
tauh0g tauh0g = .05 0.05 None
tauh1g tauh1g = .27 0.27 None
thhg thhg = -40 -40.0 None
shg shg = -12 -12.0 None
k1g k1g = 30. 30.0 None
kcag kcag = 15. 15.0 None
epsg epsg = 0.0001 0.0001 None
phig phig = 1. 1.0 None
deltang deltang = .1 0.1 None
deltahg deltahg = .05 0.05 None
iapp iapp = 2 2.0 None
gsyngg gsyngg = .3 0.3 None
vsyngg vsyngg = -80. -80.0 None
gsyng gsyng = 0.3 0.3 None
vsyng vsyng = 0 0.0 None
alphag alphag = 2 2.0 None
betag betag = .04 0.04 None
abg abg = -20 -20.0 None
iappi iappi = 5 5.0 None
igl igl = 1. 1.0 None
thresg thresg = 0.0 0.0 None
hetg hetg = 0. 0.0 None
betagi betagi = .08 0.08 None
kcagi kcagi = 15 15.0 None
gsyngi gsyngi = .3 0.3 None
gsynggi gsynggi = 1.5 1.5 None
alphaggi alphaggi = 1 1.0 None
betaggi betaggi = .1 0.1 None
vsyntc vsyntc = -85 -85.0 None
gsyntc gsyntc = .08 0.08 None
asg asg = 200 200.0 None
bsg bsg = .4 0.4 None
itc itc = 6 6.0 None
shi shi = -80 -80.0 None
shi2 shi2 = -90 -90.0 None
dur dur = 5 5.0 None
dur2 dur2 = 10 10.0 None
period period = 25 25.0 None
gnabar gnabar = 3 3.0 None
gkbar gkbar = 5 5.0 None
glbar glbar = .05 0.05 None
ena ena = 50 50.0 None
ek ek = -90 -90.0 None
eleak eleak = -70 -70.0 None
gtbar gtbar = 5 5.0 None
qht qht = 2.5 2.5 None
tadj tadj = 1 1.0 None
apr apr = 4 4.0 None
apt apt = 1 1.0 None
v1 v1 = -78. -78.0 None
v2 v2 = -55. -55.0 None
v3 v3 = -78. -78.0 None
v4 v4 = -52. -52.0 None
v5 v5 = -78. -78.0 None
v6 v6 = -78. -78.0 None
h1 h1 = 0.2 0.2 None
h2 h2 = 0.08 0.08 None
h3 h3 = 0.2 0.2 None
h4 h4 = 0.08 0.08 None
h5 h5 = 0.2 0.2 None
h6 h6 = 0.2 0.2 None
n1 n1 = 0.08 0.08 None
n2 n2 = 0.4 0.4 None
n3 n3 = 0.08 0.08 None
n4 n4 = 0.43 0.43 None
n5 n5 = 0.08 0.08 None
n6 n6 = 0.08 0.08 None
r1 r1 = 0.67 0.67 None
r2 r2 = 0.38 0.38 None
r3 r3 = 0.67 0.67 None
r4 r4 = 0.38 0.38 None
r5 r5 = 0.67 0.67 None
r6 r6 = 0.67 0.67 None
ca1 ca1 = 0.44 0.44 None
ca2 ca2 = 0.47 0.47 None
ca3 ca3 = 0.4 0.4 None
ca4 ca4 = 0.4 0.4 None
ca5 ca5 = 0.44 0.44 None
ca6 ca6 = 0.44 0.44 None
s1 s1 = 0.001 0.001 None
s2 s2 = 0.3 0.3 None
s3 s3 = 0.001 0.001 None
s4 s4 = 0.35 0.35 None
s5 s5 = 0.001 0.001 None
s6 s6 = 0.31 0.31 None
vg1 vg1 = -16. -16.0 None
vg2 vg2 = -69. -69.0 None
vg3 vg3 = -6.8 -6.8 None
vg4 vg4 = -69. -69.0 None
vg5 vg5 = -8. -8.0 None
vg6 vg6 = -69. -69.0 None
ng1 ng1 = 0.84 0.84 None
ng2 ng2 = 0.2 0.2 None
ng3 ng3 = 0.8 0.8 None
ng4 ng4 = 0.2 0.2 None
ng5 ng5 = 0.83 0.83 None
ng6 ng6 = 0.2 0.2 None
hg1 hg1 = 0.1 0.1 None
hg2 hg2 = 0.7 0.7 None
hg3 hg3 = 0.14 0.14 None
hg4 hg4 = 0.7 0.7 None
hg5 hg5 = 0.15 0.15 None
hg6 hg6 = 0.7 0.7 None
rg1 rg1 = 0.6 0.6 None
rg2 rg2 = 0.47 0.47 None
rg3 rg3 = 0.66 0.66 None
rg4 rg4 = 0.47 0.47 None
rg5 rg5 = 0.66 0.66 None
rg6 rg6 = 0.47 0.47 None
cag1 cag1 = 0.068 0.068 None
cag2 cag2 = 0.06 0.06 None
cag3 cag3 = 0.067 0.067 None
cag4 cag4 = 0.05 0.05 None
cag5 cag5 = 0.067 0.067 None
cag6 cag6 = 0.05 0.05 None
sg1 sg1 = 0.025 0.025 None
sg2 sg2 = 0.025 0.025 None
sg3 sg3 = 0.06 0.06 None
sg4 sg4 = 0.025 0.025 None
sg5 sg5 = 0.025 0.025 None
sg6 sg6 = 0.025 0.025 None
v10 v10 = -78. -78.0 None
v9 v9 = -512. -512.0 None
v11 v11 = -78. -78.0 None
v13 v13 = -129. -129.0 None
v12 v12 = -78. -78.0 None
v15 v15 = -78. -78.0 None
h10 h10 = 0.9 0.9 None
h9 h9 = 0.08 0.08 None
h11 h11 = 0.9 0.9 None
h13 h13 = 0.08 0.08 None
h12 h12 = 0.9 0.9 None
h15 h15 = 0.9 0.9 None
n10 n10 = 0.08 0.08 None
n9 n9 = 0.13 0.13 None
n11 n11 = 0.08 0.08 None
n13 n13 = 0.1311 0.1311 None
n12 n12 = 0.08 0.08 None
n15 n15 = 0.08 0.08 None
r10 r10 = 0.157 0.157 None
r9 r9 = 0.118 0.118 None
r11 r11 = 0.157 0.157 None
r13 r13 = 0.118 0.118 None
r12 r12 = 0.157 0.157 None
r15 r15 = 0.157 0.157 None
ca10 ca10 = 0.1313 0.1313 None
ca9 ca9 = 0.137 0.137 None
ca11 ca11 = 0.13 0.13 None
ca13 ca13 = 0.13 0.13 None
ca12 ca12 = 0.1313 0.1313 None
ca15 ca15 = 0.1313 0.1313 None
s10 s10 = 0.0010 0.001 None
s9 s9 = 0.11 0.11 None
s11 s11 = 0.0010 0.001 None
s13 s13 = 0.1112 0.1112 None
s12 s12 = 0.0010 0.001 None
s15 s15 = 0.1110 0.111 None
vg10 vg10 = -1015. -1015.0 None
vg9 vg9 = -159. -159.0 None
vg11 vg11 = -15.8 -15.8 None
vg13 vg13 = -159. -159.0 None
vg12 vg12 = -8. -8.0 None
vg15 vg15 = -159. -159.0 None
ng10 ng10 = 0.813 0.813 None
ng9 ng9 = 0.9 0.9 None
ng11 ng11 = 0.8 0.8 None
ng13 ng13 = 0.9 0.9 None
ng12 ng12 = 0.811 0.811 None
ng15 ng15 = 0.9 0.9 None
hg10 hg10 = 0.10 0.1 None
hg9 hg9 = 0.7 0.7 None
hg11 hg11 = 0.1013 0.1013 None
hg13 hg13 = 0.7 0.7 None
hg12 hg12 = 0.1012 0.1012 None
hg15 hg15 = 0.7 0.7 None
rg10 rg10 = 0.15 0.15 None
rg9 rg9 = 0.137 0.137 None
rg11 rg11 = 0.1515 0.1515 None
rg13 rg13 = 0.137 0.137 None
rg12 rg12 = 0.1515 0.1515 None
rg15 rg15 = 0.137 0.137 None
cag10 cag10 = 0.0158 0.0158 None
cag9 cag9 = 0.015 0.015 None
cag11 cag11 = 0.0157 0.0157 None
cag13 cag13 = 0.012 0.012 None
cag12 cag12 = 0.0157 0.0157 None
cag15 cag15 = 0.012 0.012 None
sg10 sg10 = 0.0912 0.0912 None
sg9 sg9 = 0.0912 0.0912 None
sg11 sg11 = 0.015 0.015 None
sg13 sg13 = 0.0912 0.0912 None
sg12 sg12 = 0.0912 0.0912 None
sg15 sg15 = 0.0912 0.0912 None
vt1 vt1 = -65.7 -65.7 None
ht1 ht1 = .008577 0.008577 None
thr0 thr0 = 20 20.0 None
z1 z1 = 50 50.0 None
z2 z2 = 50 50.0 None
good1 good1 = 0 0.0 None
good2 good2 = 0 0.0 None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
htt1 0.0 dimensionless None
vt2 0.0 dimensionless None
htt2 0.0 dimensionless None
ht2 0.0 dimensionless None
coun1 0.0 dimensionless None
coun2 0.0 dimensionless None
pdbs 0.0 dimensionless None
ddbs 0.0 dimensionless None
i1 0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
isyn1 0.0 dimensionless None
isyn2 0.0 dimensionless None
isyn3 0.0 dimensionless None
isyn4 0.0 dimensionless None
isyn5 0.0 dimensionless None
isyn6 0.0 dimensionless None
isyn7 0.0 dimensionless None
isyn8 0.0 dimensionless None
isyn9 0.0 dimensionless None
isyn10 0.0 dimensionless None
isyn11 0.0 dimensionless None
isyn12 0.0 dimensionless None
isyn13 0.0 dimensionless None
isyn14 0.0 dimensionless None
isyn15 0.0 dimensionless None
isyn16 0.0 dimensionless None
isyng1 0.0 dimensionless None
isyng2 0.0 dimensionless None
isyng3 0.0 dimensionless None
isyng4 0.0 dimensionless None
isyng5 0.0 dimensionless None
isyng6 0.0 dimensionless None
isyng7 0.0 dimensionless None
isyng8 0.0 dimensionless None
isyng9 0.0 dimensionless None
isyng10 0.0 dimensionless None
isyng11 0.0 dimensionless None
isyng12 0.0 dimensionless None
isyng13 0.0 dimensionless None
isyng14 0.0 dimensionless None
isyng15 0.0 dimensionless None
isyng16 0.0 dimensionless None
isyngg1 0.0 dimensionless None
isyngg2 0.0 dimensionless None
isyngg3 0.0 dimensionless None
isyngg4 0.0 dimensionless None
isyngg5 0.0 dimensionless None
isyngg6 0.0 dimensionless None
isyngg7 0.0 dimensionless None
isyngg8 0.0 dimensionless None
isyngg9 0.0 dimensionless None
isyngg10 0.0 dimensionless None
isyngg11 0.0 dimensionless None
isyngg12 0.0 dimensionless None
isyngg13 0.0 dimensionless None
isyngg14 0.0 dimensionless None
isyngg15 0.0 dimensionless None
isyngg16 0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
sga 0.0 dimensionless None
sgb 0.0 dimensionless None
sg 0.0 dimensionless None
sgt 0.0 dimensionless None
sa 0.0 dimensionless None
sb 0.0 dimensionless None
t model time 0.0 dimensionless None

Rules [90]   assignment name derived units sbo cvterm
= None None
= None None
= None None
= None None
= None None
= None None
= None None
= None None
= None None
= None None
= None None
= None None
= None None
= None None
= None None
= None None
= None None
= None None
= None None
= None None
d vt1/dt = ilf vt1 eleak glbar inaf vt1 minftc vt1 ht1 ena gnabar ikf vt1 0.75 1 ht1 ek gkbar itf vt1 mtinf vt1 htt1 gtbar ftc t dur itc period shi gsyntc sga vt1 vsyntc None
d htt1/dt = qht htinf vt1 htt1 tauht vt1 apt None
d ht1/dt = tadj hinftc vt1 ht1 tauhtc vt1 apr None
d vt2/dt = ilf vt2 eleak glbar inaf vt2 minftc vt2 ht2 ena gnabar ikf vt2 0.75 1 ht2 ek gkbar itf vt2 mtinf vt2 htt2 gtbar ftc t dur itc period shi gsyntc sgb vt2 vsyntc None
d htt2/dt = qht htinf vt2 htt2 tauht vt2 apt None
d ht2/dt = tadj hinftc vt2 ht2 tauhtc vt2 apr None
d coun1/dt = 0 None
d coun2/dt = 0 None
d pdbs/dt = 0 None
d ddbs/dt = 0 None
d i1/dt = 0 None
= None None
= 0 None
isyn1 = gsyn sg5 sg2 v1 vsyn None
isyn2 = gsyn sg6 sg1 v2 vsyn None
isyn3 = gsyn sg8 sg4 v3 vsyn None
isyn4 = gsyn sg7 sg3 v4 vsyn None
isyn5 = gsyn sg2 sg6 v5 vsyn None
isyn6 = gsyn sg1 sg5 v6 vsyn None
isyn7 = gsyn sg3 sg8 v7 vsyn None
isyn8 = gsyn sg4 sg7 v8 vsyn None
isyn9 = gsyn sg13 sg10 v9 vsyn None
isyn10 = gsyn sg14 sg9 v10 vsyn None
isyn11 = gsyn sg16 sg12 v11 vsyn None
isyn12 = gsyn sg15 sg11 v12 vsyn None
isyn13 = gsyn sg10 sg14 v13 vsyn None
isyn14 = gsyn sg9 sg13 v14 vsyn None
isyn15 = gsyn sg11 sg16 v15 vsyn None
isyn16 = gsyn sg12 sg15 v16 vsyn None
isyng1 = gsyng s8 s4 0.2 s16 vg1 vsyng None
isyng2 = gsyng s7 s3 0.2 s15 vg2 vsyng None
isyng3 = gsyng s1 s5 0.2 s9 vg3 vsyng None
isyng4 = gsyng s2 s6 0.2 s10 vg4 vsyng None
isyng5 = gsyng s4 s8 0.2 s12 vg5 vsyng None
isyng6 = gsyng s3 s7 0.2 s11 vg6 vsyng None
isyng7 = gsyng s5 s2 0.2 s13 vg7 vsyng None
isyng8 = gsyng s6 s1 0.2 s14 vg8 vsyng None
isyng9 = gsyng s16 s12 0.2 s8 vg9 vsyng None
isyng10 = gsyng s15 s11 0.2 s7 vg10 vsyng None
isyng11 = gsyng s9 s13 0.2 s1 vg11 vsyng None
isyng12 = gsyng s10 s14 0.2 s2 vg12 vsyng None
isyng13 = gsyng s12 s16 0.2 s4 vg13 vsyng None
isyng14 = gsyng s11 s15 0.2 s3 vg14 vsyng None
isyng15 = gsyng s13 s10 0.2 s5 vg15 vsyng None
isyng16 = gsyng s14 s9 0.2 s6 vg16 vsyng None
isyngg1 = gsyngg sg2 sg3 vg1 vsyngg None
isyngg2 = gsyngg sg1 sg5 vg2 vsyngg None
isyngg3 = gsyngg sg4 sg8 vg3 vsyngg None
isyngg4 = gsyngg sg3 sg1 vg4 vsyngg None
isyngg5 = gsyngg sg6 sg7 vg5 vsyngg None
isyngg6 = gsyngg sg5 sg2 vg6 vsyngg None
isyngg7 = gsyngg sg8 sg3 vg7 vsyngg None
isyngg8 = gsyngg sg7 sg4 vg8 vsyngg None
isyngg9 = gsyngg sg10 sg11 vg9 vsyngg None
isyngg10 = gsyngg sg9 sg13 vg10 vsyngg None
isyngg11 = gsyngg sg12 sg16 vg11 vsyngg None
isyngg12 = gsyngg sg11 sg9 vg12 vsyngg None
isyngg13 = gsyngg sg14 sg15 vg13 vsyngg None
isyngg14 = gsyngg sg13 sg10 vg14 vsyngg None
isyngg15 = gsyngg sg16 sg11 vg15 vsyngg None
isyngg16 = gsyngg sg15 sg12 vg16 vsyngg None
= None None
= None None
sga = sgi1 sgi2 sgi5 sgi6 sgi9 sgi10 sgi13 sgi14 None
sgb = sgi3 sgi4 sgi7 sgi8 sgi11 sgi12 sgi15 sgi16 None
sg = ftc t dur itc period shi None
sgt = ftc2 t dur2 period shi2 None
sa = sga None
sb = sgb None
t = time None

Events [4] name trigger priority delay assignments sbo cvterm
e0 None
initialValue = False
persistent = True
coun1 = coun1 1
e1 None
initialValue = False
persistent = True
coun2 = coun2 1
e2 None
initialValue = False
persistent = True
z1 = 0
good1 = good1 1
e3 None
initialValue = False
persistent = True
z2 = 0
good2 = good2 1