-
Notifications
You must be signed in to change notification settings - Fork 2
/
Copy pathsimplestim.hoc
157 lines (118 loc) · 3.34 KB
/
simplestim.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
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
back=0
//////// load needed templates////////////
if(!back){ load_file("nrngui.hoc") }
if(!back){ load_file("template/ObliquePath.hoc") }
if(!back){ load_file("template/BasalPath.hoc" ) }
objref econ,f1,f2,f3,ss,cvode
f1= new File()
f2= new File()
f3= new File()
ss=new SaveState()
cvode = new CVode(1)
x=cvode.active(1)
restart=0
v_init=-70
strdef morphology_location, morpho_path, ObliqueTrunkSection, BasalTrunkSection
objref vRP, vAPEX
proc xopen_morphology(){
sprint(morpho_path,"%s/%s",morphology_location,$s1)
xopen(morpho_path)
}
// Carmen
morphology_location = "pc2b"
ObliqueTrunkSection = "trunk[17]"
BasalTrunkSection = "trunk[7]"
xopen_morphology("cell.hoc") //reads morpholofy from the file
xopen_morphology("cell-analysis-simple.hoc") //reads simplified version of CA
// --------------- Creating lists-----------------
xopen("lib/TP-lib.hoc")
Tip_sections(apical_non_trunk_list,apical_trunk_list,"Apical")
objref apical_tip_list
apical_tip_list=TP_list // Apical Tip list
print "apical_tip_list"
apical_tip_list.printnames()
print "END apical_tip_list"
objref tmp_pl[num_tips],pl[num_tips],opl[num_tips],degree_apical_tip,peri_trunk_list
objref bl[num_tips],obl[num_tips],degree_basal_tip
xopen("lib/Oblique-lib.hoc")
oblique_sections(apical_tip_list,apical_trunk_list,num_tips) // apical dendrite path lists and degree of tips
xopen("lib/vector-distance.hoc")
//-----------------------------------------------------------------------------------------------------
printf("Setting up cell\n") // load cell-setup to
xopen("cell-setup.hoc") // specify all mechanisms, membrane properties etc
///////////////////Load File//////////////////////
//if(!back) load_file("kca.ses")
if(!back) { load_file("print.ses") }
/*cvode_active(1)*/
////// Spike counter//////////
objref apc, v1
proc insert_APC() {
apc = new APCount(0.5)
apc.thresh = $1
v1 = new Vector()
apc.record(v1)
}
proc init() {
if(restart){
f1.ropen("state.old")
ss.fread(f1)
f1.close
finitialize(v_init)
ss.restore()
t=0
fcurrent()
cvode.re_init()
} else {
finitialize(v_init)
fcurrent()
}
}
////////////////////Main///////////////////
nsyn=10
objref s[nsyn], rsyn[nsyn], nc[nsyn]
objref rsynmda[nsyn], ncnmda[nsyn]
objref tvec,ampvec,iclamp
proc main(){
tstop=4000
in=0
th=-14
access soma
insert_APC(th)
f3.wopen("spike.dat")
current=0.0 // 0.18
npulse = 10
period = 10 // 10
pulsdur = 2
pulsamp = 1.0
starttime = 1000
soma {
iclamp = new IClamp( 0.5 )
tvec = new Vector( npulse * 2 )
ampvec = new Vector( npulse * 2 )
}
iclamp.del = starttime
iclamp.amp = 0
for i=0,(npulse-1){
tvec.x[i*2 ] = starttime + i*period
tvec.x[i*2+1] = starttime + i*period + pulsdur
iclamp.dur = i*period + pulsdur
ampvec.x[i*2 ] = pulsamp
ampvec.x[i*2+1] = 0
}
ampvec.play(&iclamp.amp,tvec)
access soma[0]
//////////////run////////////////////////
run()
/////////////////////////////////////
// print the spike number
if(!back) printf("\n Current: %.4f nA | # Spikes: %d\n", current, apc.n)
if(!back) printf("___________________\n")
v1.printf(f3)
f3.close
f2.wopen("state.new")
ss.save
ss.fwrite(f2)
f2.close
}
main()
/******** end file ******/