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)

//====================================================
//                Save DATA
//====================================================	
	
proc save_data() {

// Simulation time
    sprint(filename, "%st",filepath)
    f1.wopen(filename)
    time.printf(f1)
    f1.close()

// Voltage 
    sprint(filename, "%sVgc",filepath)
    f1.wopen(filename)
    Vol.printf(f1)
    f1.close()

// Somatic spike time
    sprint(filename, "%sGs",filepath)
    f1.wopen(filename)
    Gran.spiketimes.printf(f1)
    f1.close()
	
// Dendritic spike time
    sprint(filename, "%sGd",filepath)
    f1.wopen(filename)
    Gran.dendspike.printf(f1)
    f1.close()	
	
}