load_file("nrngui.hoc")
load_file("morphology_mechanisms.hoc")
load_file("OFF_ON_OFF.ses")

prelength = 1000
mainlength1 = 1500
mainlength2 = 1500
mainlength3 = 1500

access soma

//**********************************************************************
proc simulate() { local preDT, mainDT, Kbin
    
Kbin = gbar_Kv3_AKP06
gbar_Kv3_AKP06 = 0

mainDT = dt
preDT = 0.025

dt = preDT
tstop = prelength
run()

gbar_Kv3_AKP06 = Kbin
if ( stoprun == 1) return

gbar_Kv3_AKP06 = 0
dt = mainDT
continuerun(prelength + mainlength1)
gbar_Kv3_AKP06 = Kbin
if ( stoprun == 1) return

gbar_Kv3_AKP06 = Kbin
continuerun(prelength + mainlength1 + mainlength2)
if ( stoprun == 1) return

gbar_Kv3_AKP06 = 0
continuerun(prelength + mainlength1 + mainlength2 + mainlength3)
gbar_Kv3_AKP06 = Kbin

if ( stoprun == 1) return

}


//*************User-Interface*******************************************

dt = 0.001
nrnsecmenu(0.5, 1)

xpanel("Spontaneous firing")
xlabel("Kbin Off...on...off")
xvalue("Time for Initialization", "prelength")
xvalue("Duration w/o Kbin", "mainlength1")
xvalue("Duration with Kbin", "mainlength2")
xvalue("Duration w/o Kbin", "mainlength3")

xvalue("dt", "dt")
xvalue("t", "t")
xlabel("")
xbutton("Start", "simulate()")
xbutton("Stop", "stoprun = 1")
xpanel()

simulate()