-
Notifications
You must be signed in to change notification settings - Fork 0
/
customcode_121060_chan_inKIR.hoc
87 lines (72 loc) · 3.43 KB
/
customcode_121060_chan_inKIR.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
/* Execute main.hoc to run the program. The program implements the model described in Steephen & Manchanda, 2009.
*
* Steephen, J. E., & Manchanda, R. (2009). Differences in biophysical properties of nucleus accumbens medium spiny neurons emerging from inactivation of inward rectifying potassium currents. J Comput Neurosci,
* doi:10.1007/s10827-009-0161-7
*/
//KIR
objref vecmtau_KIR, vecv_KIR
vecmtau_KIR = new Vector()
vecv_KIR = new Vector()
vecv_KIR.indgen(-120, 0, 10)
vecmtau_KIR.append(7.465, 7.465, 7.465, 8, 9.435, 10.755, 12.12, 13.795, 15.385, 14.285, 11.765, 8.89, 8) // At 35 deg C
//table_tabmtau_KIR(&vecmtau_KIR.x[0], vecv_KIR.size, &vecv_KIR.x[0])
//inKIR
objref vecv_inKIR, vechinf_inKIR, vechtau_inKIR, vecv_tau_inKIR
vecv_tau_inKIR = new Vector()
vechtau_inKIR= new Vector()
vecv_tau_inKIR.append(-120,-90, -50)
vechtau_inKIR.append(7.767, 15, 25.333) // At 35 deg C
table_tabhtau_inKIR(&vechtau_inKIR.x[0],vecv_tau_inKIR.size, &vecv_tau_inKIR.x[0])
vecv_inKIR = new Vector()
vechinf_inKIR = new Vector()
vecv_inKIR.append(-120,-90, -50)
vechinf_inKIR.append(0, 0.13, 1)
table_tabhinf_inKIR(&vechinf_inKIR.x[0], vecv_inKIR.size, &vecv_inKIR.x[0])
table_tabmtau_inKIR(&vecmtau_KIR.x[0], vecv_KIR.size, &vecv_KIR.x[0])
/*
//KRP
objref vecv_KRP, vecmtau_KRP, vechtau_KRP
vecv_KRP = new Vector()
vecv_KRP.indgen(-100, 50, 5)
vecmtau_KRP = new Vector()
vecmtau_KRP.append(13.33, 15, 16.27, 18.33, 21.47, 25, 27.97, 30, 31.17, 31.67, 31.8, 32.33, 33.07, 31.67, 26.57, 20, 14.83, 11.67, 9.77, 8.33, 6.67, 5, 3.87, 3.33, 3.2, 3.33, 3.5, 3.33, 2.67, 1.67, 1.67)
table_tabmtau_KRP(&vecmtau_KRP.x[0], vecv_KRP.size, &vecv_KRP.x[0])
vechtau_KRP = new Vector()
vechtau_KRP.append(2333.33,2333.33,2333.33,2333.33,2333.33,2333.33,2333.33,2333.33,2333.33,2333.33,2333.33,2333.33,2333.33,2333.33,2333.33,2333.33,2333.33,2333.33,2333.33,2333.33,2247.5,2000,1580.07,1166.67,927.77,833.33,778.77,733.33,694.5,666.67,666.67)
table_tabhtau_KRP(&vechtau_KRP.x[0], vecv_KRP.size, &vecv_KRP.x[0])
objref vecv_NaF, vecmtau_NaF, vechtau_NaF
vecv_NaF = new Vector()
vecv_NaF.indgen(-90, 30, 10)
vecmtau_NaF = new Vector()
vecmtau_NaF.append(0.02,0.023,0.03,0.037,0.043,0.067,0.107,0.053,0.05,0.04,0.027,0.02,0.02)
table_tabmtau_NaF(&vecmtau_NaF.x[0], vecv_NaF.size, &vecv_NaF.x[0])
vechtau_NaF = new Vector()
vechtau_NaF.append(0.433,0.433,0.433,0.433,0.433,0.433,0.433,0.283,0.167,0.15,0.107,0.1,0.093)
table_tabhtau_NaF(&vechtau_NaF.x[0], vecv_NaF.size, &vecv_NaF.x[0])
*/
/*
//CaT
objref vecv_CaT, vecmtau_CaT, vechtau_CaT
vecv_CaT = new Vector()
vecv_CaT.indgen(-65, 10, 5)
vecmtau_CaT = new Vector()
vecmtau_CaT.append(6.73,6.73,4.37,2.9,2.27,1.87,1.47,1.27,1.2,1.1,1.2,1.2,1.1,1.1,1.1,1.1)
table_tabmtau_CaT(&vecmtau_CaT.x[0], vecv_CaT.size, &vecv_CaT.x[0])
vechtau_CaT = new Vector()
vechtau_CaT.append(127.33,69.33,54,43,39.67,35.67,35.67,35.67,36,36.33,36.33,36.67,36.67,36.67,36.67,36.67)
table_tabhtau_CaT(&vechtau_CaT.x[0], vecv_CaT.size, &vecv_CaT.x[0])
//CaR
objref vecv_CaR, vechtau_CaR
vecv_CaR = new Vector()
vecv_CaR.indgen(-30, 10, 10)
vechtau_CaR = new Vector()
vechtau_CaR.append(33.33,21.67,11.67,10.00,6.67)
table_tabhtau_CaR(&vechtau_CaR.x[0], vecv_CaR.size, &vecv_CaR.x[0])
//NaP
objref vecv_NaP, vechtau_NaP
vecv_NaP = new Vector()
vecv_NaP.indgen(-100, 10, 10)
vechtau_NaP = new Vector()
vechtau_NaP.append(1500, 1583.33, 1733.33, 2033.33, 2100, 1666.67, 1416.67, 1166.67, 1000, 900, 833.33, 700)
table_tabhtau_NaP(&vechtau_NaP.x[0], vecv_NaP.size, &vecv_NaP.x[0])
*/