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, 3-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, 31-1 dend2[i] {
gcabar_spike = 0.002
gkbar_spike = 0.012
gabar_spike = 0.036
gkcbar_spike = 0.000001
gnabar_spike = 0.025
}
for i = 0, 35-1 dend3[i] {
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, 1-1 dend5[i] {
gcabar_spike = 0.002
gkbar_spike = 0.012
gabar_spike = 0.036
gkcbar_spike = 0.000001
gnabar_spike = 0.025
}
for i = 0, 57-1 dend6[i] {
gcabar_spike = 0.002
gkbar_spike = 0.012
gabar_spike = 0.036
gkcbar_spike = 0.000001
gnabar_spike = 0.025
}
for i = 0, 1-1 dend7[i] {
gcabar_spike = 0.002
gkbar_spike = 0.012
gabar_spike = 0.036
gkcbar_spike = 0.000001
gnabar_spike = 0.025
}
for i = 0, 7-1 dend8[i] {
gcabar_spike = 0.002
gkbar_spike = 0.012
gabar_spike = 0.036
gkcbar_spike = 0.000001
gnabar_spike = 0.025
}

soma {
diam=13.9804
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.054
gkcbar_spike = 0.000065
gnabar_spike = 0.1
}

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

load_proc("nrnmainmenu")
access soma

// xopen("test.ses")

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