// This function is used to deduce the maximum AMPA coductance value
// at any location along the cell, such that a single current pulse will give
// rise to approximately 5mV local depolarization.
// tuneepsps.hoc is used is advance to calculate maximum AMPA conductance values
// for a few locations along each section (reference list). This function is used to
// deduce the AMPA conductance for ANY specified location in ANY section by finding
// the estimated AMPA value of the reference point closest to the specified point.
// written by Terrence Brannon, last modified by Yiota Poirazi, July 2001, poirazi@LNC.usc.edu
// last modified by Jose Gomez, 1st May 2006, jfcgomez@ull.es
objref find_gmax_vec, find_gmax_dest, find_gmax_index, find_gmax_tmp
func find_gmax() { local range_ref, pre_diff, pre_i, post_gmax, diff, i,R
range_ref = $1
find_gmax_vec = new Vector()
find_gmax_index = new Vector()
find_gmax_tmp = new Vector()
post_i=1
for i=0,tune_epsp_list.count()1 { // for all sections in tune_epsp_list (reference list)
if (issection(tune_epsp_list.object(i).section_name)) {
// print "\t\t tuning", secname() //
// printf("R=tune_epsp_list.object(%d).range_ref\n", i) //
R = tune_epsp_list.object(i).range_ref
// printf("R=tune_epsp_list.object(%d).range_ref == %g\n", i,R) //
if (R <= range_ref) {
pre_i=i
pre_R=R
}
if (R > range_ref) {
if (post_i < 0) {
post_i=i
post_R=R
}
}
}
}
//print "===================================== calculating ret_gmax, ret_R"
//print "jooooo, post_i", post_i
if ( (post_i <0)  (abs(pre_Rrange_ref) < abs(post_Rrange_ref)) ) {
if(pre_R<=0){pre_i=post_i //I put this if, Jose Gomez,I had wrong
//value if the first element if was ret_R=0
pre_R=post_R
}
// print "pre_i", pre_i //
ret_gmax=tune_epsp_list.object(pre_i).gbar_ampa
ret_R=pre_R
} else {
// print "post_i", post_i, "post_R ", post_R //
ret_gmax=tune_epsp_list.object(post_i).gbar_ampa
ret_R=post_R
}
// printf("%g was the closest range_ref to %g for %s\n", ret_R, range_ref, secname()) //
return(ret_gmax)
}
