objref sr,k
objref matrix_coord,vector_L
objref TP_list
//---------------------- Tip_sections() ----------------------------------
//Outputs: matrix_coord is the a matrix with TP coord.
// TP_list is a list with the TP sections
// vector_L is a vector with the lenght of the Terminal section
// num_tips is the number of Tips
//Inputs: $o1 is object type SectionList() where are dendrites (or apical dendrites) sections
// $o2 is object type SectionList() where are soma (or trunk) sections
// $s3 is a string where has "Dendrites" or "Apical dendrites"
objref jo
proc Tip_sections(){local i, nrow
i=0
nrow=1
matrix_coord = new Matrix(nrow,3)
TP_list = new SectionList()
vector_L=new Vector(nrow)
strdef tmp_str, tmp_str2
forsec $o1{
sr=new SectionRef()
if (sr.nchild==0) {
secname() TP_list.append()
matrix_coord.resize(nrow,3)
matrix_coord.x[nrow-1][0]=x3d(n3d()-1)
matrix_coord.x[nrow-1][1]=y3d(n3d()-1)
matrix_coord.x[nrow-1][2]=z3d(n3d()-1)
vector_L.resize(nrow)
vector_L.x[nrow-1]=L
nrow+=1
}
}
// TP_list.printnames()
nrow-=1
print "\nNumber of Terminal Point in ",$s3,": ", nrow, "\n"
num_tips=nrow
}