create dend1
create dend2
create SpineStalk
create SpineHead
create Soma

access dend1
 L = 630 // lambda=1261 for 1 um Dia
 diam = 1 
nseg=315

access dend2
 L =  630  
 diam = 1
nseg=315

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

access SpineHead
L=0.7
diam=1.0
nseg=1

access Soma
L=25
diam=25
nseg=3

connect Soma(1),dend1(0)
connect dend1(1), dend2(0)
connect SpineStalk(0), dend1(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
celsius = 22
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
}

access Soma
insert spike
celsius = 22
insert cad
ena = 35.0
ek = -75
taur_cad = 1.5
gcabar_spike = 0
gkbar_spike = 0 //0.012
gabar_spike = 0 //0.036
gkcbar_spike = 0 // 0.00005
gnabar_spike = 0 //0.04

access dend1
insert spike
celsius = 22
insert cad
ena = 35.0
ek = -75
taur_cad=1.5
gcabar_spike = 0
gkbar_spike = 0 //0.012
gabar_spike = 0 //0.036
gkcbar_spike = 0 // 0.00005
gnabar_spike = 0 //0.04


access dend2
insert spike
celsius = 22
insert cad
ena = 35.0
ek = -75
taur_cad=1.5
gcabar_spike = 0 // .0015
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 dend1
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")