//objref stim // for removing pervious AP IClamp
// g_vgraph.view(110, -60, 40, 120, 496, 152, 300.24, 200.1)

g_vgraph.size(110,150,-60,60)
objref stim2
soma stim2 = new IClamp(0.5)
stim2.del  = 110
stim2.dur  = 15
stim2.amp  = 0.18

{
xpanel("IClamp Paramters", 0)
xvalue("Delay","stim2.del", 1,"stim2.del", 0, 0 )
xvalue("Duration","stim2.dur", 1,"stim2.dur", 0, 0 )
xvalue("Amplitude","stim2.amp", 1,"stim2.amp", 0, 0 )
xpanel(80,721)
}

tstop = 150

objref gr3, gr4, gr5 //, gr1, gr2

// {
// gr1 = new Graph(0)
// gr1.size(0,150,-90,70)
// {gr1.view(0, -90, 150, 160, 448, 61, 249.12, 179.22)}
// gr1.addvar("soma.v(0.5)", 1, 2, 0, 1, 2)
// }

// {
// gr2 = new Graph(0)
// gr2.size(0,200,136,166)
// {gr2.view(0, 136, 200, 30, 441, 441, 249.84, 189.3)}
// gr2.addexpr("soma.cai(0.5)*1e6", 2, 2, 0, 1, 2)
// }

{
gr3 = new Graph(0)
gr3.size(100,150,0.2,0.55)
{gr3.view(110, 0.2, 40, 0.35, 903, 430, 294.48, 186.42)}
graphList[2].append(gr3)
gr3.addexpr("soma.m_skca3( 0.5 )", 1, 1, 0.601956, 0.930898, 2)
}

{
gr4 = new Graph(0)
gr4.size(100,150,0.0015,0.006)
{gr4.view(110, 0, 40, 0.006, 897, 2, 297.36, 198.66)}
graphList[2].append(gr4)
gr4.addexpr("soma.oinf_skca3( 0.5 )", 1, 1, 0.57724, 0.943491, 2)
}

{
gr5 = new Graph(0)
gr5.size(100,150,0,2.5e-05)
{gr5.view(110, 0, 40, 2.5e-05, 904, 831, 295.2, 174.9)}
graphList[2].append(gr5)
gr5.addvar("soma.g_skca3( 0.5 )", 1, 1, 0.621951, 0.9, 2)
}

gr3.exec_menu("Keep Lines")
gr4.exec_menu("Keep Lines")
gr5.exec_menu("Keep Lines")

{
xpanel("SKCa Conductance", 0)
xlabel("For Fig 16A,C & D,")
xlabel("change gbar_skca3 to 0.0027 and 0.0045 mho/cm2")
xlabel("and Press Init & Run in RunControl Panel")
soma.gbar_skca3(0.5) = 0.0009
xvalue("gbar_skca3","soma.gbar_skca3(0.5)", 1,"soma.gbar_skca3(0.5)", 0, 0 )
xpanel(25,948)
}