-
Notifications
You must be signed in to change notification settings - Fork 1
/
ntIh.mod
83 lines (63 loc) · 1.25 KB
/
ntIh.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
: $Id: ntIh.mod,v 1.6 1998/08/14 20:52:38 billl Exp $
TITLE anomalous rectifier channel
:
: Anomalous Rectifier Ih - cation (Na/K) channel for geniculate interneurons
: Differential equations
:
:
: Written by Jun Zhu, Univ. Wisconsin, Jan 1996
:
INDEPENDENT {t FROM 0 TO 1 WITH 1 (ms)}
NEURON {
SUFFIX iar
USEION other WRITE iother VALENCE 1
RANGE ghbar, shift
GLOBAL h_inf, tau_h, erev, stp
}
UNITS {
(molar) = (1/liter)
(mM) = (millimolar)
(mA) = (milliamp)
(mV) = (millivolt)
(msM) = (ms mM)
}
PARAMETER {
v (mV)
erev = -44 (mV)
celsius = 36 (degC)
ghbar = .0008 (mho/cm2)
shift = -6 (mV)
stp = 7.4
}
STATE {
h
}
ASSIGNED {
i (mA/cm2)
iother (mA/cm2)
h_inf
tau_h (ms)
tadj
}
BREAKPOINT {
SOLVE state METHOD cnexp
iother = ghbar * h * (v - erev)
}
DERIVATIVE state {
evaluate_fct(v)
h' = (h_inf - h) / tau_h
}
UNITSOFF
INITIAL {
tadj = 3.0 ^ ((celsius-36)/10)
evaluate_fct(v)
h = h_inf
: Experiments of Coulter et al were at 36 deg.C
: Q10 is assumed equal to 3
:
}
PROCEDURE evaluate_fct(v (mV)) {
h_inf = 1 / ( 1 + exp((v+shift+79)/stp) )
tau_h = exp((v+shift+293.3)/29.67) / ( 1 + exp((v+shift+76.7)/7.82)) / tadj
}
UNITSON