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
|
//genesis
/* FILE INFORMATION
**
** Tabchannel implementation of Hodgkin-Huxley squid Na and K channels
** This version uses SI (MKS) units and is equivalent to (but faster than)
** hh_chan.g, except for the function and channel names.
*/
// CONSTANTS
float EREST_ACT = -0.060
float ENA = 0.045
float EK = -0.090
/* Square meters */
float SOMA_A = 1e-9
/*************************************************************************
The function setupalpha uses the same form for the gate variables
X and Y as the vdep_gate, namely: (A+B*V)/(C+exp((V+D)/F))
The constants above are related to the hh_channel constants by:
EXPONENTIAL :
gate variables in terms of hh-channel variables
A A
B 0
C 0
D -V0
F -B
SIGMOID :
Gate in terms of hh-ch variables
A A
B 0
C 1
D -V0
F B
LINOID :
A -A * V0
B A
C -1
D -V0
F B
*************************************************************************/
//========================================================================
// Tabchan Na channel
//========================================================================
function make_Na_hh_tchan
str chanpath = "Na_hh_tchan"
// str chanpath = "Na_squid_hh"
if ({exists {chanpath}})
return
end
create tabchannel {chanpath}
// V
// S
// A
// S
setfield ^ Ek {ENA} Gbar {1.2e3*SOMA_A} Ik 0 Gk 0 \
Xpower 3 Ypower 1 Zpower 0
setupalpha {chanpath} X {0.1e6*(0.025 + EREST_ACT)} -0.1e6 \
-1.0 {-1.0*(0.025 + EREST_ACT)} -0.01 \
4e3 0.0 0.0 {-1.0*EREST_ACT} 18e-3
setupalpha {chanpath} Y 70.0 0.0 0.0 \
{-1.0*EREST_ACT} 0.02 1.0e3 0.0 1.0 \
{-1.0*(0.030 + EREST_ACT)} -10.0e-3
end
//========================================================================
// Tabchan version of K channel
//========================================================================
function make_K_hh_tchan
str chanpath = "K_hh_tchan"
// str chanpath = "K_squid_hh"
if (({exists {chanpath}}))
return
end
create tabchannel {chanpath}
// V
// S
// A
// S
setfield ^ Ek {EK} Gbar {360.0*SOMA_A} Ik 0 Gk 0 \
Xpower 4 Ypower 0 Zpower 0
setupalpha {chanpath} X {10e3*(0.01 + EREST_ACT)} -10.0e3 \
-1.0 {-1.0*(0.01 + EREST_ACT)} -0.01 125.0 0.0 0.0 \
{-1.0*EREST_ACT} 80.0e-3
end
|
