-
Notifications
You must be signed in to change notification settings - Fork 1
/
ctt1209a_denstim.hoc
139 lines (123 loc) · 2.36 KB
/
ctt1209a_denstim.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
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
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
create initseg,narrowr,axon
initseg {L=40}
narrowr {L=90}
axon {L=5340}
initseg nseg=int(L/15)
narrowr nseg=int(L/15)
axon nseg=int(L/15)
forall insert pas
forall insert spike
forall ena = 35.0
forall ek = -75
forall insert cad
forall g_pas = .000008
forall e_pas = -62.5
forall Ra=110
global_ra=110
celsius=22
for i = 0, 7-1 dend1[i] {
gcabar_spike = 0.002
gkbar_spike = 0.012
gabar_spike = 0.036
gkcbar_spike = 0.000001
gnabar_spike = 0.025
}
for i = 0, 11-1 dend2[i] {
gcabar_spike = 0.002
gkbar_spike = 0.012
gabar_spike = 0.036
gkcbar_spike = 0.000001
gnabar_spike = 0.025
}
dend3 {
gcabar_spike = 0.002
gkbar_spike = 0.012
gabar_spike = 0.036
gkcbar_spike = 0.000001
gnabar_spike = 0.025
}
for i = 0, 23-1 dend4[i] {
gcabar_spike = 0.002
gkbar_spike = 0.012
gabar_spike = 0.036
gkcbar_spike = 0.000001
gnabar_spike = 0.025
}
for i = 0, 47-1 dend5[i] {
gcabar_spike = 0.002
gkbar_spike = 0.012
gabar_spike = 0.036
gkcbar_spike = 0.000001
gnabar_spike = 0.025
}
soma {
diam=10.0289
gcabar_spike = 0.0015
gkbar_spike = 0.018
gabar_spike = 0.054
gkcbar_spike = 0.000065
gnabar_spike = 0.080
}
initseg {
diam=1
gcabar_spike = 0.0015
gkbar_spike = 0.018
gabar_spike = 0.054
gkcbar_spike = 0.000065
gnabar_spike = 0.150
}
narrowr {
diam=0.4
gcabar_spike = 0.0
gkbar_spike = 0.018
gabar_spike = 0.0
gkcbar_spike = 0.000065
gnabar_spike = 0.2
}
axon {
diam=1
gcabar_spike = 0.0
gkbar_spike = 0.018
gabar_spike = 0.0
gkcbar_spike = 0.000065
gnabar_spike = 0.070
}
connect soma(1),initseg(0)
connect initseg(1),narrowr(0)
connect narrowr(1),axon(0)
// objectvar clamp
// soma clamp = new IClamp(0.5)
// clamp.del = 0
// clamp.dur = 400
// clamp.amp = 0.01
objectvar clamp1
dend1[1] clamp1 = new IClamp(0.5)
clamp1.del = 0
clamp1.dur = 600
clamp1.amp = 0.014
objectvar clamp2
dend2[6] clamp2 = new IClamp(0.5)
clamp2.del = 0
clamp2.dur = 600
clamp2.amp = 0.014
objectvar clamp3
dend3[0] clamp3 = new IClamp(0.625)
clamp3.del = 0
clamp3.dur = 600
clamp3.amp = 0.014
objectvar clamp4
dend4[13] clamp4 = new IClamp(0.5)
clamp4.del = 0
clamp4.dur = 600
clamp4.amp = 0.014
objectvar clamp5
dend5[2] clamp5 = new IClamp(0.5)
clamp5.del = 0
clamp5.dur = 600
clamp5.amp = 0.014
load_proc("nrnmainmenu")
access soma
xopen("dstima.ses")
init()
forall depth_cad = diam / 2
init()