I am pretty new to NEURON, but I am using the Impedance class tools to analyse the electrotonic structure of some neurons and I am using the following chunk of code to do that:
Code: Select all
// prepare to use Impedance class
// always a good idea to finitialize before computing impedance
v_init=-65
finitialize(v_init)
// demonstrate use of impedance class
objref zz
zz = new Impedance()
FREQ = 0 // Hz
WHERE = 0.5 // location in the soma that is the reference point
soma distance(0, WHERE) // sets origin for distance calculations
proc calcZ() {
soma zz.loc(WHERE) // sets origin for impedance calculations
zz.compute(FREQ, 1) // takes the impedance contributions of
// gating state differential equations into account
// but requires mechanisms to be compatible with CVODE
print "distance(x) ratio(x)"
forall {
print secname()
for (x) print x, x3d(x), y3d(x), z3d(x), log(1/(zz.ratio(x)))
}
}
calcZ()
Code: Select all
dend[38]
0 -15.63 3.8699999 5.23 0
0.16666667 -15.63 3.8699999 5.23 0.64765023
0.5 -15.63 3.8699999 5.23 1.3266357
0.83333333 -15.63 3.8699999 5.23 1.7573273
1 -16.469999 3.8699999 5.23 1.9388532
Is there any way to tweak the previous code to have it print only the "x=1" line?
Thank you all!