// ET_rig.ses
// studies properties of single ET cell model under voltage clanp

load_file("et.hoc")
objref et[4]
et[0] = new ET()

// add Ih 

objref ih_sects
ih_sects = new SectionList()
et[0].priden ih_sects.append()
et[0].tuftden ih_sects.append()
forsec ih_sects {
  insert Ih
  gkhbar_Ih  = 0.0024/10 // 1/10 value used in Rubin Cleland 2007 mitral cell
}

forall L = L/2

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

//Begin PointProcessGroupManager
{
load_file("pointgrp.hoc")
}
{
ocbox_ = new PointProcessGroupManager(0)
}
{object_push(ocbox_)}
{
ms = new MechanismStandard("IClamp")
ms.set("del", 500, 0)
ms.set("dur", 1000, 0)
ms.set("amp", -0.05, 0)
mnews("IClamp", ms)
select(0)
execute("ET[0].soma ocbox_.move(0.5)")
}
{object_pop() doNotify()}
{
ocbox_ = ocbox_.v1
ocbox_.map("PointProcessGroupManager", 680, 14, 556.8, 295.68)
}
objref ocbox_
//End PointProcessGroupManager

{
save_window_ = new Graph(0)
save_window_.size(0,2500,-80,40)
scene_vector_[3] = save_window_
{save_window_.view(0, -80, 2500, 120, 159, 227, 300.48, 200.32)}
graphList[0].append(save_window_)
save_window_.save_name("graphList[0].")
save_window_.addexpr("v(.5)", 1, 1, 0.8, 0.9, 2)
}
{
xpanel("RunControl", 0)
v_init = -65
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 = 2500
xvalue("t","t", 2 )
tstop = 2500
xvalue("Tstop","tstop", 1,"tstop_changed()", 0, 1 )
dt = 0.025
xvalue("dt","dt", 1,"setdt()", 0, 1 )
steps_per_ms = 40
xvalue("Points plotted/ms","steps_per_ms", 1,"setdt()", 0, 1 )
screen_update_invl = 0.05
xvalue("Scrn update invl","screen_update_invl", 1,"", 0, 1 )
realtime = 3.85
xvalue("Real Time","realtime", 0,"", 0, 1 )
xpanel(175,623)
}
objectvar scene_vector_[1]
{doNotify()}