Nd = 10
objref f1
f1 = new File()
strdef filepath, filename
filepath = "celldata/"
objref time, Vol
time = new Vector()
Vol = new Vector()
time.record(&t, Nd*dt)
Vol.record(&Gran.soma.v(0.5), Nd*dt)
proc save_data() {
sprint(filename, "%st",filepath)
f1.wopen(filename)
time.printf(f1)
f1.close()
sprint(filename, "%sVgc",filepath)
f1.wopen(filename)
Vol.printf(f1)
f1.close()
sprint(filename, "%sGs",filepath)
f1.wopen(filename)
Gran.spiketimes.printf(f1)
f1.close()
sprint(filename, "%sGd",filepath)
f1.wopen(filename)
Gran.dendspike.printf(f1)
f1.close()
}