XPP model

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

# Excitatory cortical neurons, 0 [Ca].
#
number Cm=1.0
number pms=3,pns=4
number VNa=55.0,t_tauh=-40.5,t_taun=-27.0
number thetaa=-50.0,sigmaa=20.0,thetab=-80.0,sigmab=-6.0,tauBs=15.0
number sigmam=9.5,sigmah=-7.0,sigman=10.0,sigmaz=5.0,sigmak=7
# 
p gNa=35.0,gKdr=6.0,gL=0.05,Iapp=0.661914
p gA=1.4,gNaP=0.3,gZ=1.0
p thetaz=-39.0,tauZs=75.0
p phi=10.0, thetah=-45.0
p thetam=-30.0,thetan=-35.0,thetap=-47.0,sigmap=3.0
p VK=-90.0,VL=-70.0
#
GAMMAF(VV,theta,sigma)=1.0/(1.0+exp(-(VV-theta)/sigma))
#ZFUNC(AA,CA,zz)=1/(1+(AA^zz/CA^zz))
#
VVs'=(-gL*(VVs-VL)-INa-INaP-IKdr-IA-Iz+Iappx)/Cm
hhs'=phi*(GAMMAF(VVs,thetah,sigmah)-hhs)/(1.0+7.5*GAMMAF(VVs,t_tauh,-6.0))
nns'=phi*(GAMMAF(VVs,thetan,sigman)-nns)/(1.0+5.0*GAMMAF(VVs,t_taun,-15.0))
bbs'=(GAMMAF(VVs,thetab,sigmab)-bbs)/tauBs
zzs'=(GAMMAF(VVs,thetaz,sigmaz)-zzs)/tauZs
#
Iappx=Iapp
#Iappx=if(t<=3.0)then(Iapp)else(0.0)
Minfs=GAMMAF(VVs,thetam,sigmam)
Pinfs=GAMMAF(VVs,thetap,sigmap)
Ainfs=GAMMAF(VVs,thetaa,sigmaa)
#
INa=gNa*(Minfs^pms)*hhs*(VVs-VNa)
INaP=gNaP*Pinfs*(VVs-VNa)
IKdr=gKdr*(nns^pns)*(VVs-VK)
IA=gA*Ainfs^3*bbs*(VVs-VK)
Iz=gZ*zzs*(VVs-VK)
#
VVs(0)=-71.81327
hhs(0)=0.98786
nns(0)=0.02457
bbs(0)=0.203517
zzs(0)=0.00141
#
@ MAXSTOR=800000
@ BACK=Black
@ XP=T
@ YP=VVs
@ AXES=2
@ TOTAL=500.0
@ DT=0.05
@ NJMP=1
@ T0=0.0
@ TRANS=0.0
@ NMESH=40
@ METH=rungekutta
@ DTMIN=0.001
@ DTMAX=1.0
@ TOLER=0.00001
@ BOUND=10000.0
@ DELAY=0
@ XLO=0.0, XHI=500.0, YLO=-90.0, YHI=30.0
@ NTST=50,NMAX=2000,NPR=50
@ DS=0.02,DSMIN=0.001,DSMAX=0.5
@ PARMIN=-10,PARMAX=50,NORMMIN=0.0,NORMMAX=10000.0
@ AUTOVAR=VVs1,AUTOXMIN=-10.0,AUTOXMAX=50.0,AUTOYMIN=-90.0,AUTOYMAX=50.0
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 [5] 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
gammaf vv theta sigma 1 1 vv theta sigma

Parameters [48] name constant value unit derived unit sbo cvterm
cm cm = 1.0 1.0 None
pms pms = 3 3.0 None
pns pns = 4 4.0 None
vna vna = 55.0 55.0 None
t_tauh t_tauh = -40.5 -40.5 None
t_taun t_taun = -27.0 -27.0 None
thetaa thetaa = -50.0 -50.0 None
sigmaa sigmaa = 20.0 20.0 None
thetab thetab = -80.0 -80.0 None
sigmab sigmab = -6.0 -6.0 None
taubs taubs = 15.0 15.0 None
sigmam sigmam = 9.5 9.5 None
sigmah sigmah = -7.0 -7.0 None
sigman sigman = 10.0 10.0 None
sigmaz sigmaz = 5.0 5.0 None
sigmak sigmak = 7 7.0 None
gna gna = 35.0 35.0 None
gkdr gkdr = 6.0 6.0 None
gl gl = 0.05 0.05 None
iapp iapp = 0.661914 0.661914 None
ga ga = 1.4 1.4 None
gnap gnap = 0.3 0.3 None
gz gz = 1.0 1.0 None
thetaz thetaz = -39.0 -39.0 None
tauzs tauzs = 75.0 75.0 None
phi phi = 10.0 10.0 None
thetah thetah = -45.0 -45.0 None
thetam thetam = -30.0 -30.0 None
thetan thetan = -35.0 -35.0 None
thetap thetap = -47.0 -47.0 None
sigmap sigmap = 3.0 3.0 None
vk vk = -90.0 -90.0 None
vl vl = -70.0 -70.0 None
vvs vvs = -71.81327 -71.81327 None
hhs hhs = 0.98786 0.98786 None
nns nns = 0.02457 0.02457 None
bbs bbs = 0.203517 0.203517 None
zzs zzs = 0.00141 0.00141 None
iappx 0.0 dimensionless None
minfs 0.0 dimensionless None
pinfs 0.0 dimensionless None
ainfs 0.0 dimensionless None
ina 0.0 dimensionless None
inap 0.0 dimensionless None
ikdr 0.0 dimensionless None
ia 0.0 dimensionless None
iz 0.0 dimensionless None
t model time 0.0 dimensionless None

Rules [15]   assignment name derived units sbo cvterm
d vvs/dt = gl vvs vl ina inap ikdr ia iz iappx cm None
d hhs/dt = phi gammaf vvs thetah sigmah hhs 1 7.5 gammaf vvs t_tauh 6 None
d nns/dt = phi gammaf vvs thetan sigman nns 1 5 gammaf vvs t_taun 15 None
d bbs/dt = gammaf vvs thetab sigmab bbs taubs None
d zzs/dt = gammaf vvs thetaz sigmaz zzs tauzs None
iappx = iapp None
minfs = gammaf vvs thetam sigmam None
pinfs = gammaf vvs thetap sigmap None
ainfs = gammaf vvs thetaa sigmaa None
ina = gna minfs pms hhs vvs vna None
inap = gnap pinfs vvs vna None
ikdr = gkdr nns pns vvs vk None
ia = ga ainfs 3 bbs vvs vk None
iz = gz zzs vvs vk None
t = time None