// This experiment is to test the sensitivity of V_S
load_file("nrngui.hoc")
load_file("simulation_base.hoc")
dt = 0.001
objref cutoff_file
strdef cutoff_file_name,run_file_name
// want parameters to be set specifically by arguments,
// otherwise fail
// g_gj, vs, layer
// vary soma leak with g_pas=2e-3
// vary distance of gap junction from soma
kinetics = 1
load_gj_axons(2,layer,1,0,g_gj,kinetics)
print "loaded axons"
//setup_gj_session()
if(test_axon.stim_loc==0){ quit() }
cutoff_file = new File()
sprint(cutoff_file_name,"vs-expcellgj_J_gj%g-vs%d.out",g_gj,vs)
cutoff_file.wopen(cutoff_file_name)
set_soma_leak(axon1,vs,2e-3)
set_soma_leak(axon2,vs,2e-3)
stim1.amp=2
stim2.amp=2
d_dist = 0.02
n_iter = (dist1_end-dist1_beg)/d_dist + 1
for(j=0;j<n_iter;j=j+1){
dist = dist1_beg + j*d_dist
print "distance=", dist
cutoff_file.printf("%g\t",dist)
sprint(run_file_name,"vs-expcellgj-run_J_gj%g-vs%d-d%d.out",g_gj,vs,dist)
gj.reset_loc(dist/axon1.main_length)
run_search(cutoff_file,run_file_name,3)
//refine_search(cutoff_file,run_file_name,3)
}
//quit() //uncomment when running as a batch