-
Notifications
You must be signed in to change notification settings - Fork 5
/
kcRT03.mod
78 lines (61 loc) · 1.19 KB
/
kcRT03.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
: $Id: kcRT03.mod,v 1.5 2003/09/30 21:11:43 billl Exp $
TITLE Potasium C type current for RD Traub, J Neurophysiol 89:909-921, 2003
COMMENT
Implemented by Maciej Lazarewicz 2003 (mlazarew@seas.upenn.edu)
ENDCOMMENT
INDEPENDENT { t FROM 0 TO 1 WITH 1 (ms) }
UNITS {
(mV) = (millivolt)
(mA) = (milliamp)
}
NEURON {
SUFFIX kcRT03
USEION k READ ek WRITE ik
USEION ca READ cai
RANGE gmax, ik, g, i
GLOBAL alpha,beta
}
PARAMETER {
gmax = 0.0 (mho/cm2)
v ek (mV)
cai (1)
}
ASSIGNED {
g i
ik (mA/cm2)
alpha beta (/ms)
}
STATE {
m
}
BREAKPOINT {
SOLVE states METHOD cnexp
if( 0.004 * cai < 1 ) {
g = gmax * m * 0.004 * cai
} else {
g = gmax * m
}
i = g * (v-ek)
ik=i
}
INITIAL {
settables(v)
m = alpha / ( alpha + beta )
m = 0
}
DERIVATIVE states {
settables(v)
m' = alpha * ( 1 - m ) - beta * m
}
UNITSOFF
PROCEDURE settables(v) {
if( v < -10.0 ) {
alpha = 2 / 37.95 * ( exp( ( v + 50 ) / 11 - ( v + 53.5 ) / 27 ) )
: Note that there is typo in the paper - missing minus sign in the front of 'v'
beta = 2 * exp( ( - v - 53.5 ) / 27 ) - alpha
}else{
alpha = 2 * exp( ( - v - 53.5 ) / 27 )
beta = 0
}
}
UNITSON