124063/Cone_ih.mod at master · ModelDBRepository/124063 · GitHub

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
: Cone Photoreceptor h channel by using kinetics in Barnes' paper

NEURON
{
	SUFFIX h_cone

	NONSPECIFIC_CURRENT ih

	RANGE ghbar, gh, eh


}

UNITS
{
	(mA) = (milliamp)
	(mV) = (millivolt)
	(mS) = (millimho)
}

PARAMETER
{


       : h channel
        ghbar = 3.5 (mS/cm2) <0,1e9>
        eh  = -32.5 (mV)


}

STATE
{
	nh

}

ASSIGNED
{
	v (mV)

	ih (mA/cm2)

	infh
	tauh   (ms)

	gh (mho/cm2)

}

INITIAL
{
	rate(v)
	nh  = infh
}


BREAKPOINT
{
	SOLVE states METHOD cnexp
	gh  = (0.001)*ghbar*(1-(1+3*nh)*(1-nh)^3)
	ih  = gh*(v - eh)

	: the current is in the unit of mA/cm2


}

DERIVATIVE states
{
	rate(v)
	nh'  = (infh  - nh )/tauh

}


FUNCTION alphah(v(mV)) (/ms)
{
	alphah = 0.001*18/( exp  (  ( v+88)/12 ) + 1 )
}


FUNCTION betah(v(mV)) (/ms)
{
	betah = 0.001*18/( exp  ( - ( v+18)/19 ) + 1 )
}


PROCEDURE rate(v (mV))
{
        LOCAL a, b


	a = alphah(v)
	b = betah(v)
	tauh = 1/(a + b)
	infh = a/(a + b)

}