-
Notifications
You must be signed in to change notification settings - Fork 5
/
nafPR.mod
167 lines (137 loc) · 3.18 KB
/
nafPR.mod
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
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
: $Id: nafPR.mod,v 1.6 2004/02/09 21:19:55 billl Exp $
INDEPENDENT {t FROM 0 TO 1 WITH 1 (ms)}
NEURON { SUFFIX nafPR }
NEURON { USEION na WRITE ina }
ASSIGNED { ina }
PARAMETER {
erev = 60. (mV)
gmax = 0.030 (mho/cm2)
vrest = -60.0
exptemp = 37
maflag = 3
malphaA = -0.32
malphaB = -4.0
malphaV0 = 13.1
mbflag = 3
mbetaA = 0.28
mbetaB = 5.0
mbetaV0 = 40.1
mq10 = 3
mexp = -2
haflag = 1
halphaA = 0.128
halphaB = -18
halphaV0 = 17.
hbflag = 2
hbetaA = 4.
hbetaB = -5.
hbetaV0 = 40.
hq10 = 3
hexp = 1
celsius (degC)
dt (ms)
v (mV)
} : end PARAMETER
NEURON {
RANGE gmax, g, i
GLOBAL erev, Inf, Tau, vrest, qq10
} : end NEURON
CONSTANT {
FARADAY = 96489.0 : Faraday's constant
R= 8.31441 : Gas constant
} : end CONSTANT
UNITS {
(mA) = (milliamp)
(mV) = (millivolt)
(umho) = (micromho)
} : end UNITS
ASSIGNED {
i (mA/cm^2)
g (mho/cm^2)
Inf[2] : 0 = m and 1 = h
Tau[2] : 0 = m and 1 = h
qq10[2]
} : end ASSIGNED
STATE { h }
INITIAL {
mh(v)
h = Inf[1]
}
BREAKPOINT {
SOLVE states METHOD cnexp
mh(v)
g = gmax * Inf[0]*Inf[0] * h
i = g*(v-erev)
ina=i
} : end BREAKPOINT
: ASSIGNMENT PROCEDURES
: Must be given by a user routines in parameters.multi
: E.G.:
: PROCEDURE iassign () { i = g*(v-erev) ina=i }
: PROCEDURE iassign () { i = g*ghkca(v) ica=i }
:-------------------------------------------------------------------
DERIVATIVE states {
mh(v)
h' = (-h + Inf[1]) / Tau[1]
}
:-------------------------------------------------------------------
: NOTE : 0 = m and 1 = h
PROCEDURE mh (v) {
LOCAL a, b, j
qq10[0] = mq10^((celsius-exptemp)/10.)
qq10[1] = hq10^((celsius-exptemp)/10.)
: Calculater Inf and Tau values for h and m
FROM j = 0 TO 1 {
a = alpha (v, j)
b = beta (v, j)
Inf[j] = a / (a + b)
Tau[j] = 1. / (a + b) / qq10[j]
if (hexp==0) { Tau[1] = 1. Inf[1] = 1.}
}
} : end PROCEDURE mh (v)
:-------------------------------------------------------------------
FUNCTION alpha(v,j) {
LOCAL flag, A, B, V0
if (j==1 && hexp==0) {
alpha = 0
} else {
if (j == 1) {
A = halphaA B = halphaB V0 = halphaV0+vrest flag = haflag
} else {
A = malphaA B = malphaB V0 = malphaV0+vrest flag = maflag
}
if (flag == 1) { : EXPONENTIAL
alpha = A*exp((v-V0)/B)
} else if (flag == 2) { : SIGMOID
alpha = A/(exp((v-V0)/B)+1)
} else if (flag == 3) { : LINOID
if(v == V0) {
alpha = A*B
} else {
alpha = A*(v-V0)/(exp((v-V0)/B)-1) }
}
}
} : end FUNCTION alpha (v,j)
:-------------------------------------------------------------------
FUNCTION beta (v,j) {
LOCAL flag, A, B, V0
if (j==1 && hexp==0) {
beta = 1
} else {
if (j == 1) {
A = hbetaA B = hbetaB V0 = hbetaV0+vrest flag = hbflag
} else {
A = mbetaA B = mbetaB V0 = mbetaV0+vrest flag = mbflag
}
if (flag == 1) { : EXPONENTIAL
beta = A*exp((v-V0)/B)
} else if (flag == 2) { : SIGMOID
beta = A/(exp((v-V0)/B)+1)
} else if (flag == 3) { : LINOID
if(v == V0) {
beta = A*B
} else {
beta = A*(v-V0)/(exp((v-V0)/B)-1) }
}
}
} : end FUNCTION beta (v,j)