create dend1
create dend2
create SpineStalk
create SpineHead

access dend1
 L = 4000
 diam = 1 
nseg=401

access dend2
 L =  4000  
 diam = 1
nseg=401

access SpineStalk
L=1
diam=.1
nseg=5

access SpineHead
L=0.5
diam=0.5
nseg=1

connect dend1(0), dend2(1)
connect SpineStalk(0), dend2(1)
connect SpineHead(0), SpineStalk(1)

forall{
// nseg=int(L/10)+1
insert pas
g_pas=1/70000
Ra=110
global_ra=110
e_pas=-70
}

CABAR=2.2
SpineHead{
insert spike
insert cad
ena = 35.0
ek = -75
taur_cad=0.5
gcabar_spike = CABAR
gkbar_spike = 0 //0.012
gabar_spike = 0 //0.036
gkcbar_spike = 0 // 0.00005
gnabar_spike = 0 //0.04
}

proc set_cabar(){
SpineHead{gcabar_spike = CABAR}
}

/*------------------------------------*/
// SPINE STMULATION
access SpineHead
objectvar salpha,IChead

// IChead= new IClamp(.5)

salpha = new AlphaSynapse(0.50)
salpha.onset = 10
salpha.tau = 10
salpha.gmax = 0.01

// ------------------------------------
// DENDRITE STIMULATION
access dend2
objectvar  dalpha,ddalpha,ICdend

// ICdend = new IClamp(1)

dalpha = new AlphaSynapse(1)
dalpha.onset = 10
dalpha.tau = 10
dalpha.gmax = 0.01

ddalpha =  new AlphaSynapse(0.95)
ddalpha.onset = 10
ddalpha.tau = 10
ddalpha.gmax = 0.01

// ------------------------------------
objectvar vbox1, tempobj1
proc show_alphas(){
        vbox1 = new VBox()
        vbox1.intercept(1)

        xpanel("")
        xvalue("Spine Head  gmax","salpha.gmax", 1, "")
        xvalue("Spine Base (a) gmax", "dalpha.gmax", 1, "")
        xvalue("Dendrite (b) gmax", "ddalpha.gmax", 1, "")
        xvalue("Ca++ channel density","CABAR", 1,"set_cabar()") 

        xpanel()
        vbox1.intercept(0)
        vbox1.map
}//

 load_proc ("nrnmainmenu")
 // Load_file("spine.ses")