create initseg,narrowr,axon

initseg {L=40}
narrowr {L=90}
axon {L=5340}

initseg nseg=int(L/15)
narrowr nseg=int(L/15)
axon nseg=int(L/15)

forall insert pas
forall insert spike
forall ena = 35.0
forall ek = -75
forall insert cad
forall g_pas = .000008
forall e_pas = -62.5
forall Ra=110
global_ra=110

celsius=22

for i = 0, 7-1 dend1[i] {
gcabar_spike = 0.002
gkbar_spike = 0.012
gabar_spike = 0.036
gkcbar_spike = 0.000001
gnabar_spike = 0.025
}
for i = 0, 11-1 dend2[i] {
gcabar_spike = 0.002
gkbar_spike = 0.012
gabar_spike = 0.036
gkcbar_spike = 0.000001
gnabar_spike = 0.025
}
dend3 {
gcabar_spike = 0.002
gkbar_spike = 0.012
gabar_spike = 0.036
gkcbar_spike = 0.000001
gnabar_spike = 0.025
}
for i = 0, 23-1 dend4[i] {
gcabar_spike = 0.002
gkbar_spike = 0.012
gabar_spike = 0.036
gkcbar_spike = 0.000001
gnabar_spike = 0.025
}
for i = 0, 47-1 dend5[i] {
gcabar_spike = 0.002
gkbar_spike = 0.012
gabar_spike = 0.036
gkcbar_spike = 0.000001
gnabar_spike = 0.025
}

soma {
diam=10.0289 
gcabar_spike = 0.0015
gkbar_spike = 0.018
gabar_spike = 0.054
gkcbar_spike = 0.000065
gnabar_spike = 0.080
}

initseg {
diam=1
gcabar_spike = 0.0015
gkbar_spike = 0.018
gabar_spike = 0.054
gkcbar_spike = 0.000065
gnabar_spike = 0.150
}

narrowr {
diam=0.4
gcabar_spike = 0.0
gkbar_spike = 0.018
gabar_spike = 0.0
gkcbar_spike = 0.000065
gnabar_spike = 0.2
}

axon {
diam=1
gcabar_spike = 0.0
gkbar_spike = 0.018
gabar_spike = 0.0
gkcbar_spike = 0.000065
gnabar_spike = 0.070
}


connect soma(1),initseg(0)
connect initseg(1),narrowr(0)
connect narrowr(1),axon(0)

// objectvar clamp
// soma clamp = new IClamp(0.5)
// clamp.del = 0
// clamp.dur = 400
// clamp.amp = 0.01
objectvar clamp1
dend1[1] clamp1 = new IClamp(0.5)
clamp1.del = 0
clamp1.dur = 600
clamp1.amp = 0.014
objectvar clamp2
dend2[6] clamp2 = new IClamp(0.5)
clamp2.del = 0
clamp2.dur = 600
clamp2.amp = 0.014
objectvar clamp3
dend3[0] clamp3 = new IClamp(0.625)
clamp3.del = 0
clamp3.dur = 600
clamp3.amp = 0.014
objectvar clamp4
dend4[13] clamp4 = new IClamp(0.5)
clamp4.del = 0
clamp4.dur = 600
clamp4.amp = 0.014
objectvar clamp5
dend5[2] clamp5 = new IClamp(0.5)
clamp5.del = 0
clamp5.dur = 600
clamp5.amp = 0.014

load_proc("nrnmainmenu")
access soma

xopen("dstima.ses")

init()
forall depth_cad = diam / 2
init()