neuron 7.5-1 / usr / share / nrn / demo / dynclamp.ses

{load_file("nrngui.hoc")}
objectvar save_window_, rvp_
objectvar scene_vector_[10]
objectvar ocbox_, ocbox_list_, scene_, scene_list_
{ocbox_list_ = new List()  scene_list_ = new List()}
{pwman_place(0,0,0)}

//Begin CellBuild[0]
{
load_file("celbild.hoc", "CellBuild")
}
{ocbox_ = new CellBuild(1)}
{object_push(ocbox_)}
{
version(5.7)
continuous = 1
}
{object_push(topol)}
{
first = 0
slist.remove_all()
sname = "soma"
objref tobj
}
{
tobj = new CellBuildSection("soma",0, 0, tobj, 0) slist.append(tobj)
  tobj.position(0,0,15,0) tobj.lx=8.75202 tobj.ly=-19.3733 tobj.i3d=0
all_init()
}
for i=0, slist.count-1 {slist.object(i).rdses()}
{object_pop()}
{
}
{object_push(subsets)}
{first = 0}
{consist()}
{object_pop()}
{
}
{object_push(geom)}
{
first = 0
tobj = new GeoSpec(4)
tobj.value = 100
bild.subsets.snlist.object(0).geo.append(tobj)
set_default()
}
{object_pop()}
{
}
{object_push(memb)}
{first=0}
{
tobj = new FakeMechStan(0)
tobj.value = 80
tobj.set_default()
tobj = new MStanWrap(tobj, 0)
bild.subsets.snlist.object(0).ml.append(tobj)
}
{
tobj = new FakeMechStan(1)
tobj.value = 1
tobj.set_default()
tobj = new MStanWrap(tobj, 0)
bild.subsets.snlist.object(0).ml.append(tobj)
}
{
tobj = new MechanismStandard("khh")
tobj.set("gmax_khh", 0.036, 0)
tobj = new MStanWrap(tobj, 1)
bild.subsets.snlist.object(0).ml.append(tobj)
}
{
tobj = new MechanismStandard("leak")
tobj.set("gmax_leak", 0.0003, 0)
tobj.set("e_leak", -54.3, 0)
tobj = new MStanWrap(tobj, 1)
bild.subsets.snlist.object(0).ml.append(tobj)
}
{object_pop()}
{
}
{object_push(manage)}
{
first = 0
classname = "Cell"
etop=1 esub=1 egeom=1 emem=1
itop=1 isub=0 igeom=0 imem=0
bild.topol.names_off = 0
bild.topol.circles_off = 0
output_index = 0  output_x = 1
}
{object_pop()}
{
cexport()
}
{object_pop()}
{
save_window_=ocbox_.gtopol
save_window_.size(-200,200,-150,150)
scene_vector_[2] = save_window_
ocbox_.gtopol = save_window_
save_window_.save_name("ocbox_.gtopol")
}
{
ocbox_.map("CellBuild[0]", 16, 540, 609.6, 287.04)
}
objref ocbox_
//End CellBuild[0]


//Begin CellBuild[1]
{
load_file("celbild.hoc", "CellBuild")
}
{ocbox_ = new CellBuild(1)}
{object_push(ocbox_)}
{
version(5.7)
continuous = 1
}
{object_push(topol)}
{
first = 0
slist.remove_all()
sname = "synth"
objref tobj
}
{
tobj = new CellBuildSection("synth",0, 0, tobj, 0) slist.append(tobj)
  tobj.position(0,20,15,0) tobj.lx=4.86222 tobj.ly=-16.7801 tobj.i3d=0
all_init()
}
for i=0, slist.count-1 {slist.object(i).rdses()}
{object_pop()}
{
}
{object_push(subsets)}
{first = 0}
{consist()}
{object_pop()}
{
}
{object_push(geom)}
{
first = 0
tobj = new GeoSpec(4)
tobj.value = 100
bild.subsets.snlist.object(0).geo.append(tobj)
set_default()
}
{object_pop()}
{
}
{object_push(memb)}
{first=0}
{
tobj = new FakeMechStan(1)
tobj.value = 1e-06
tobj.set_default()
tobj = new MStanWrap(tobj, 0)
bild.subsets.snlist.object(0).ml.append(tobj)
}
{
tobj = new MechanismStandard("nahh")
tobj.set("gmax_nahh", 0.12, 0)
tobj = new MStanWrap(tobj, 1)
bild.subsets.snlist.object(0).ml.append(tobj)
}
{
tobj = new MechanismStandard("extracellular")
tobj.set("xraxial", 1e+09, 0)
tobj.set("xraxial", 1e+09, 1)
tobj.set("xg", 1e-09, 0)
tobj.set("xg", 1e+09, 1)
tobj.set("xc", 0, 0)
tobj.set("xc", 0, 1)
tobj.set("e_extracellular", 0, 0)
tobj = new MStanWrap(tobj, 1)
bild.topol.slist.object(0).ml.append(tobj)
}
{object_pop()}
{
}
{object_push(manage)}
{
first = 0
classname = "Cell"
etop=1 esub=1 egeom=1 emem=1
itop=1 isub=0 igeom=0 imem=0
bild.topol.names_off = 0
bild.topol.circles_off = 0
output_index = 0  output_x = 1
}
{object_pop()}
{
cexport()
}
{object_pop()}
{
save_window_=ocbox_.gtopol
save_window_.size(-200,200,-150,150)
scene_vector_[3] = save_window_
ocbox_.gtopol = save_window_
save_window_.save_name("ocbox_.gtopol")
}
{
ocbox_.map("CellBuild[1]", 11, 626, 611.52, 288.96)
}
objref ocbox_
//End CellBuild[1]


//Begin PointProcessManager
{
load_file("pointman.hoc")
}
{
soma ocbox_ = new PointProcessManager(0)
}
{object_push(ocbox_)}
{
mt.select("IClamp") i = mt.selected()
ms[i] = new MechanismStandard("IClamp")
ms[i].set("del", 0.5, 0)
ms[i].set("dur", 0.1, 0)
ms[i].set("amp", 0.3, 0)
mt.select("IClamp") i = mt.selected() maction(i)
hoc_ac_ = 0.5
sec.sec move() d1.flip_to(0)
}
{object_pop() doNotify()}
{
ocbox_ = ocbox_.v1
ocbox_.map("PointProcessManager", 137, 466, 208.32, 326.4)
}
objref ocbox_
//End PointProcessManager

//Begin LinearCircuit[0]
{
load_file("lincir.hoc", "LinearCircuit")
}
{
ocbox_ = new LinearCircuit(1)
}
{object_push(ocbox_)}
{version(2)}
{mkelm(8, 190, 20, 3, 0)}
{mklabel(0, "soma(0.5)", 0.115, -19.654)}
{mklabel(1, "Vm", -14.669, -2.8385)}
{sel.extra_info.set("soma", 0.5) sel.extra_info.name(sel)}
{mkelm(1, 220, 70, 3.60555, 0.588003)}
5
{mklabel(0, "Re2", 12.002, -9.6973)}
{mklabel(2, "Vx", -1.571, 10.548)}
{restore_ic(2, -65)}
{mkelm(7, 110, 80, 1, 3.14159)}
1
0
{mklabel(0, "Ac", 3.732, 0.8358)}
{mklabel(3, "Vc", 4.0199, 13.1556)}
{restore_ic(3, -65)}
{mkelm(1, 270, 190, 1, 0)}
1
{mklabel(0, "Ri", -0.073, -9.409)}
{mkelm(7, 270, 230, 1, 1.5708)}
1
0
{mklabel(0, "Ai", -2.301, 2.1792)}
{mkelm(0, 245, 250, 2.5, 0)}
{mkelm(7, 240, 190, 1, 0)}
100000
0
{mklabel(0, "Acmd", 1.251, 2.6388)}
{mkelm(0, 220, 225, 2.5, -1.5708)}
{mkelm(0, 290, 190, 1, 3.14159)}
{mkelm(9, 50, 135, 2.5, 1.5708)}
{mklabel(0, "synth(0.5)", 44.0199, 1.686)}
{sel.extra_info.set("synth", 0.5) sel.extra_info.name(sel)}
{mkelm(7, 110, 180, 1, 0)}
1e+06
0
{mklabel(0, "AI", -1.5878, 1.3549)}
{mklabel(3, "VIc", 2.622, -7.997)}
{mkelm(1, 110, 210, 2, 0)}
1
{mklabel(0, "RI", 0.3138, 10.3188)}
{mkelm(0, 175, 180, 4.5, 3.14159)}
{mkelm(1, 160, 70, 3.60555, -0.588003)}
100
{mklabel(0, "Re1", 16.974, 14.1211)}
{mkelm(0, 300, 140, 5, 1.5708)}
{mkelm(0, 270, 200, 1.41421, 0.785398)}
{mkelm(0, 270, 200, 1.41421, -0.785398)}
{mkelm(4, 275, 90, 2.5, 3.14159)}
{mklabel(0, "I2", -0.289, -19.8414)}
{sel.extra_info.restore()}
3
1e+09 0
1e+09 0
1e+09 0
{mkelm(4, 70, 80, 2, 3.14159)}
{mklabel(0, "I1", -1.4411, -15.5907)}
{sel.extra_info.restore()}
3
1e+09 0
1e+09 0
1e+09 0
{mkelm(6, 130, 60, 1, 0)}
{mkelm(0, 130, 195, 1.5, -1.5708)}
{mkelm(0, 90, 200, 1, 1.5708)}
{mkelm(6, 90, 160, 1, 0)}
{mkelm(0, 70, 190, 2, 0)}
{mkelm(0, 50, 175, 1.5, -1.5708)}
{mkelm(0, 50, 95, 1.5, 1.5708)}
{parasitic_ = 0  noconsist_ = 0}
{graphlist.append(new LincirGraph(this, 1))}
1
Vm (mV)
2 1 0.8 0.9 2
0 5 -80 40 // graph size
687 25 338.88 245.76 // box size
// end info
{graphlist.append(new LincirGraph(this, 1))}
2
I1 I (nA)
1 1 0.8 0.9 2
I2 I (nA)
2 1 0.8 0.9 2
0 5 -1 1 // graph size
684 293 342.72 240 // box size
// end info
{g.exec_menu("Simulate")  tool(2)}
{sel = nil}
{object_pop()}
{
{
save_window_=ocbox_.g
save_window_.size(0,300,0,200)
scene_vector_[5] = save_window_
ocbox_.g = save_window_
save_window_.save_name("ocbox_.g")
save_window_.label(190.115, 0.346, "soma(0.5)", 1, 1, 0.5, 0.5, 1)
save_window_.label(175.331, 47.1615, "Vm", 1, 1, 0.5, 0.5, 1)
save_window_.label(232.002, 60.3027, "Re2", 1, 1, 0.5, 0.5, 1)
save_window_.label(248.429, 100.548, "Vx", 1, 1, 0.5, 0.5, 1)
save_window_.label(113.732, 80.8358, "Ac", 1, 1, 0.5, 0.5, 1)
save_window_.label(94.0199, 93.1557, "Vc", 1, 1, 0.5, 0.5, 1)
save_window_.label(269.927, 180.591, "Ri", 1, 1, 0.5, 0.5, 1)
save_window_.label(267.699, 232.179, "Ai", 1, 1, 0.5, 0.5, 1)
save_window_.label(241.251, 192.639, "Acmd", 1, 1, 0.5, 0.5, 1)
save_window_.label(94.0199, 136.686, "synth(0.5)", 1, 1, 0.5, 0.5, 1)
save_window_.label(108.412, 181.355, "AI", 1, 1, 0.5, 0.5, 1)
save_window_.label(132.622, 172.003, "VIc", 1, 1, 0.5, 0.5, 1)
save_window_.label(110.314, 220.319, "RI", 1, 1, 0.5, 0.5, 1)
save_window_.label(176.974, 84.1211, "Re1", 1, 1, 0.5, 0.5, 1)
save_window_.label(274.711, 70.1586, "I2", 1, 1, 0.5, 0.5, 1)
save_window_.label(68.5589, 64.4093, "I1", 1, 1, 0.5, 0.5, 1)
}
ocbox_.map("LinearCircuit[0]", 681, 572, 477.12, 299.52)
}
objref ocbox_
//End LinearCircuit[0]

objectvar scene_vector_[1]
{doNotify()}