-
Notifications
You must be signed in to change notification settings - Fork 0
/
Sim1.hoc
67 lines (50 loc) · 1.09 KB
/
Sim1.hoc
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
load_file("nrngui.hoc")
/* axon morphology */
create s, a[10], b[10]
forall {
Ra = 110
cm = 1
}
access s
nseg = 1
diam = 10
L = 10
for i = 0, 9 a[i] {
nseg = 100
diam = 0.2
L = 100
}
for i = 0, 9 b[i] {
nseg = 4
diam = 4
L = 4
}
connect s(1), a[0](0)
for i = 0, 9 {connect b[i](0), a[i](1)
}
for i = 0, 8 {connect a[i+1](0), b[i](1)
}
/* channel */
forall {insert hhmfb
insert KIn
gnabar_hhmfb = 0.05 /* S/cm^2 */
gkbar_hhmfb = 0 /* S/cm^2 */
gl_hhmfb = 0.0001 /* S/cm^2 */
gkbar_KIn = 0.036 /* S/cm^2 */
ena = 50 /* mV */
ek = -85 /* mV */
celsius = 25
}
forsec "s" gnabar_hhmfb = 0.01 /* S/cm^2 */
// forsec "b[7]" gkbar_KIn = 0 /* S/cm^2*/
// forsec "b[8]" gkbar_KIn = 0 /* S/cm^2*/
// forsec "b[9]" gkbar_KIn = 0 /* S/cm^2*/
// forsec "a[8]" gkbar_KIn = 0 /* S/cm^2*/
// forsec "a[9]" gkbar_KIn = 0 /* S/cm^2*/
/* single shock */
objectvar stim
s stim = new IClamp(0.5)
stim.del = 50
stim.dur = 0.2
stim.amp = 1
load_file("rec1.ses")