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
    
mainDT = dt
preDT = 0.001

dt = preDT
tstop = prelength
run()

if ( stoprun == 1) return

soma.switch_Kv3_AKP06_Kv3 = 1
dt = mainDT
continuerun(prelength + mainlength1)
if ( stoprun == 1) return

soma.switch_Kv3_AKP06_Kv3 = 0
continuerun(prelength + mainlength1 + mainlength2)
if ( stoprun == 1) return

soma.switch_Kv3_AKP06_Kv3 = 1
continuerun(prelength + mainlength1 + mainlength2 + mainlength3)

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()
quit()