-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathsimpleleak.mod
More file actions
62 lines (51 loc) · 1.24 KB
/
simpleleak.mod
File metadata and controls
62 lines (51 loc) · 1.24 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
TITLE simple leak current that uses a non-specific current to allow
: the easy setting of both conductance and reversal potential.
UNITS {
(pA) = (picoamp)
(molar) = (1/liter)
(mV) = (millivolt)
(S) = (siemens)
(mA) = (milliamp)
(mM) = (millimolar)
F = (faraday) (coulomb)
R = (mole k) (mV-coulomb/degC)
}
INDEPENDENT {v FROM -100 TO 50 WITH 50 (mV)}
NEURON {
SUFFIX leek
NONSPECIFIC_CURRENT ileak
RANGE ileak,gbar,eleak
}
PARAMETER {
dt (ms)
ena (mV)
nai (mM)
celsius = 35 (degC)
gkbar = 5.5e-6 (S/cm2)
gnabar = 2.375e-6 (S/cm2)
gcabar = 0.136e-6 (S/cm2)
ggabaa = 0.0e-6 (S/cm2)
gampa = 0.0e-6 (S/cm2)
eampa = 0 (mV)
eca = 120 (mV)
ecl = -70 (mV)
ek = -90.0 (mV)
nao = 145 (mM)
}
ASSIGNED {
ica (mA/cm2)
ina (mA/cm2)
ik (mA/cm2)
icl (mA/cm2)
iampa (mA/cm2)
ileak (mA/cm2)
}
BREAKPOINT {
ena = 50
ica = gcabar*(v-eca)
ina = gnabar*(v-ena)
icl = ggabaa*(v-ecl)
ik = gkbar*(v-ek)
iampa = gampa*(v-eampa)
ileak= ica + ina + ik + icl + iampa
}