-
Notifications
You must be signed in to change notification settings - Fork 1
/
purkinje.hoc
113 lines (82 loc) · 1.93 KB
/
purkinje.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
load_file("nrngui.hoc")
create dend
access dend
nseg = 1
L = 20
diam = 4
cm=1.5
insert climbing
del_climbing=2
g_climbing = 0.002
icin_climbing=-0.016
e_climbing=-74
tauC_climbing=10
tauO_climbing=0.4
delf_climbing=2
ef_climbing=-74
tauCf_climbing=10
tauOf_climbing=10
fastfact_climbing=0
insert leak
gbar_leak =9e-5
insert CaT3_1
pcabar_CaT3_1=0.000
t1_CaT3_1=0.6
t2_CaT3_1=55
insert KA1
gkbar_KA1=0.03
vh1_KA1=50
vm1_KA1=27
vm2_KA1=30
insert newCaP1
pcabar_newCaP1=0.000
pp_newCaP1=0.75
insert Kdr1
gkbar_Kdr1=0.092
pp_Kdr1=8
pdr_Kdr1=0.4
vsh1_Kdr1=0
vsh2_Kdr1=0
hdr_Kdr1=1
insert KC3
gkbar_KC3=4.800
tin_KC3=4
insert SK2
gkbar_SK2 =0.005
insert cdp5
Nannuli_cdp5 = 0.326 + (1.94 * diam ) + (0.289*diam *diam ) - ((3.33e-2)*diam *diam *diam ) + ((1.55e-3)*diam *diam *diam *diam ) - (2.55e-5*diam *diam *diam *diam *diam )
Buffnull2_cdp5 = 64.2 - 57.3*exp(-diam /1.4)
rf3_cdp5 = 0.162 - 0.106*exp(-diam /2.29)
if (diam >=2) {
rf4_cdp5 = 0.000267 + 0.0167*exp(-diam /0.722) + 0.0028*exp(-diam /4)
} else {
rf4_cdp5 = 0.003
}
//Pump
kpmp1_cdp5 = 3e-3
kpmp2_cdp5 = 1.75e-5
kpmp3_cdp5 = 7.255e-5
TotalPump_cdp5 = 1e-9
//Dye Fura-ff
BTCnull_cdp5 =1
b1_cdp5 =570
b2_cdp5 = 3.5
//Buffer
DMNPEnull_cdp5=0.75
c1_cdp5=570
c2_cdp5=4
//Calbindin
CBnull_cdp5=0.16
nf1_cdp5=43.5
nf2_cdp5=3.58e-2
ns1_cdp5=5.5
ns2_cdp5=0.26e-2
//Parvalbumin
PVnull_cdp5=0.075
m1_cdp5=107
m2_cdp5=9.5e-4
p1_cdp5=0.8
p2_cdp5=2.5e-2
global_ra = 250
set_ra()
celsius=37