/*
5/17/2017  Revised by N.T. Carnevale for the sake of conceptual clarity
and to facilitate attributed reuse.
In this version, the reference temperature is 23 deg C
and the value assigned to celsius is the actual operating temperature
in degrees celsius.
*/

objectvar save_window_, rvp_
objectvar scene_vector_[4]
objectvar ocbox_, ocbox_list_, scene_, scene_list_
{ocbox_list_ = new List()  scene_list_ = new List()}

//Begin I/V Clamp Electrode
{
load_file("electrod.hoc")
}
{
ocbox_=new Electrode(0)
execute("can_locate=1 sec=\"soma\" xloc=0.5 locate(0)", ocbox_)
execute("vc.dur[0]=0 vc.amp[0]=-80", ocbox_)
execute("vc.dur[1]=40 vc.amp[1]=-10", ocbox_)
execute("vc.dur[2]=0 vc.amp[2]=-65", ocbox_)
execute("stim.del=10 stim.dur=400 stim.amp=0.02", ocbox_)
execute("vcsteps=5", ocbox_)
execute("samp=stim.amp  store_vclamp() glyph()", ocbox_)
ocbox_ = ocbox_.v1
ocbox_.map("I/V Clamp Electrode", 946, 377, 237.75, 444.75)
}
objref ocbox_
//End I/V Clamp Electrode

{
xpanel("NEURON Main Panel", 0)
xcheckbox("Quiet",&stdrun_quiet,"")
realtime = 0
xvalue("Real Time","realtime", 0,"", 0, 1 )
// celsius = 6.3
celsius = 23 // actual operating temperature
xvalue("celsius","celsius", 1,"", 0, 1 )
global_ra = 200
xvalue("Axial Resistivity","global_ra", 1,"set_ra()", 1, 1 )
xbutton("RunControl","nrncontrolmenu()")
xmenu("New Graph", 0)
xbutton("Voltage axis","newPlotV()")
xbutton("Current axis","newPlotI()")
xbutton("State axis","newPlotS()")
xbutton("Shape plot","newshapeplot()")
xbutton("Vector movie","newvectorplot()")
xbutton("Phase Plane","newphaseplane()")
xbutton("Grapher","load_proc(\"makegrapher\") makegrapher()")
xmenu()
xmenu("Point Processes", 0)
xmenu("Managers", 0)
xbutton("Point Manager","load_template(\"PointProcessManager\") makeppm()")
xbutton("Point Group","load_template(\"PointProcessGroupManager\") makeppgm()")
xbutton("Electrode","load_proc(\"makeelectrode\") makeelectrode()")
xmenu()
xmenu("Viewers", 0)
xmenu("PointProcesses", 0)
xbutton("IClamp","makePointBrowser(\"IClamp\")")
xbutton("AlphaSynapse","makePointBrowser(\"AlphaSynapse\")")
xbutton("SEClamp","makePointBrowser(\"SEClamp\")")
xbutton("VClamp","makePointBrowser(\"VClamp\")")
xbutton("APCount","makePointBrowser(\"APCount\")")
xmenu()
xmenu()
xmenu()
xmenu("Distributed Mechanisms", 0)
xmenu("Managers", 0)
xbutton("Inserter","load_template(\"Inserter\") makeinserter()")
xbutton("Homogeneous Spec","load_proc(\"makeshowmechanism\") makeshowmechanism()")
xmenu()
xmenu("Viewers", 0)
xbutton("Shape Name","load_template(\"MenuExplore\") makeMenuExplore()")
xbutton("Name Values","nrnallsectionmenu()")
xmenu("Mechanisms (Globals)", 0)
xbutton("na_ion","nrnglobalmechmenu(\"na_ion\")")
xbutton("k_ion","nrnglobalmechmenu(\"k_ion\")")
xbutton("hh","nrnglobalmechmenu(\"hh\")")
xmenu()
xmenu()
xmenu()
xmenu("Miscellaneous", 0)
xmenu("Clipboard", 0)
xbutton("Save to File","load_proc(\"clipboard_save\") clipboard_save()")
xbutton("Retrieve from File","load_proc(\"clipboard_retrieve\") clipboard_retrieve()")
xmenu()
xbutton("Family","load_template(\"Family\") makeFamily()")
xbutton("Parameterized Function","load_template(\"FunctionFitter\") makefitter()")
xbutton("Run Fitter","load_template(\"RunFitter\") makerunfitter()")
xmenu("Impedance", 0)
xbutton("Frequency","load_template(\"ImpedanceRatio\") makeImpRatio()")
xbutton("Path","load_template(\"Impx\") makeImpx()")
xbutton("log(A) vs x","load_template(\"LogAvsX\") makelogax()")
xbutton("Shape","load_template(\"ImpShape\") makeImpShape()")
xmenu()
xmenu()
xpanel(-3,378)
}
{
xpanel("RunControl", 0)
v_init = -120
xvalue("Init","v_init", 1,"stdinit()", 1, 1 )
xbutton("Init & Run","run()")
xbutton("Stop","stoprun=1")
runStopAt = 5
xvalue("Continue til","runStopAt", 1,"{continuerun(runStopAt) stoprun=1}", 1, 1 )
runStopIn = 1
xvalue("Continue for","runStopIn", 1,"{continuerun(t + runStopIn) stoprun=1}", 1, 1 )
xbutton("Single Step","steprun()")
t = 10
xvalue("t","t", 2 )
tstop = 10
xvalue("Tstop","tstop", 1,"tstop_changed()", 0, 1 )
dt = 0.02
xvalue("dt","dt", 1,"setdt()", 0, 1 )
steps_per_ms = 50
xvalue("Points plotted/ms","steps_per_ms", 1,"setdt()", 0, 1 )
xpanel(275,536)
}
{
save_window_ = new Graph(0)
save_window_.size(0,10,-1,1)
scene_vector_[3] = save_window_
{save_window_.view(0, -1, 10, 2, 325, 5, 433.5, 298)}
graphList[1].append(save_window_)
save_window_.save_name("graphList[1].")
save_window_.xaxis(0,250,0,5,0,0,1)
save_window_.yaxis(-2,2,0,4,0,0,1)
save_window_.addexpr("soma.ina_B_Na( 0.5 )", 1, 1, 0.30974, 0.919914, 2)
save_window_.addexpr("soma.m_B_Na( 0.5 )", 3, 1, 0.8, 0.9, 2)
save_window_.addvar("soma.h_B_Na( 0.5 )", 2, 1, 0.8, 0.9, 2)
}
{
xpanel("soma(0 - 1) (Parameters)", 0)
xlabel("soma(0 - 1) (Parameters)")
xlabel("nseg = 1")
soma.L = 100
xvalue("L","soma.L", 1,"define_shape()", 0, 0 )
soma.diam = 10
xvalue("diam","soma.diam", 1,"", 0, 0 )
soma.cm = 1
xvalue("cm","soma.cm", 1,"", 0, 0 )
soma.gnabar_B_Na = 0.07
xvalue("gnabar_B_Na","soma.gnabar_B_Na", 1,"", 0, 0 )
soma.gkbar_KDRI = 0.0042
xvalue("gkbar_KDRI","soma.gkbar_KDRI", 1,"", 0, 0 )
soma.ek = -84
xvalue("ek","soma.ek", 1,"", 0, 0 )
soma.g_pas = 5e-05
xvalue("g_pas","soma.g_pas", 1,"", 0, 0 )
soma.e_pas = -80
xvalue("e_pas","soma.e_pas", 1,"", 0, 0 )
soma.gnabar_SS = 0.0006
xvalue("gnabar_SS","soma.gnabar_SS", 1,"", 0, 0 )
soma.ena = 53
xvalue("ena","soma.ena", 1,"", 0, 0 )
xpanel(599,519)
}
objectvar scene_vector_[1]
{doNotify()}